Diversified pattern of the human colorectal cancer microbiome
© Geng et al; licensee BioMed Central Ltd. 2013
Received: 28 January 2013
Accepted: 5 March 2013
Published: 7 March 2013
The aim of this study is to expand existing knowledge about the CRC-associated microbiome among Han Chinese, and to further discover the variation pattern of the human CRC microbiome across all population.
Using pyrosequencing-based molecular monitoring of bacterial 16S rRNA gene from eight tumor/normal tissue pairs of eight Chinese CRC patients, we analyzed and characterized the basic features of the CRC-associated microbiome. Firstly, we discovered an increasing diversity among tumor-associated bacterial communities. Secondly, in 50% of Chinese CRC patients, we found a significant increase of Roseburia (P = 0.017), and a concurrent decrease of both Microbacterium (P = 0.009) and Anoxybacillus (P = 0.009) in tumor tissue.
We discovered a novel CRC microbiome pattern in Chinese. Both the over-represented Roseburia bacteria at tumor sites and the over-represented Microbacterium and Anoxybacillus bacteria away from tumor sites were both closely related in Chinese CRC patients. Across several populations reported in this study and previously, we observed both common and distinctive patterns of human CRC microbiome’s association with a high-risk of CRC.
KeywordsChinese CRC Microbiome Pyrosequencing
Bacterial infections play a potentially significant role in the pathogenesis of colorectal cancer (CRC) [1–5], though there are a variety of dietary, genetic, and environmental factors that add to CRC risk . Previous studies however have reported that various infectious agents (e.g., Fusobacterium spp., Bacteroides fragilis, and Escherichia coli) are related to a high-risk of CRC across several different populations. These findings prompted us ask whether such associations between CRC and microbiome patterns are common or population-specific. Such a finding to this question is essential in developing personalized medicine strategies and treatment options for CRC patients.
The aim of this study was to quantitatively evaluate the differences of bacterial communities and compositions between eight tumor/normal pairs from eight Chinese CRC patients, as well as to characterize both the common and various patterns of the human CRC microbiome among different populations, ultimately to contribute towards a preliminary understanding of the bacterial driving forces at play in CRC.
Patients, sample collection, and DNA extraction
Eight Chinese CRC patients from Kunming, China were used in this study with (56.9±14.4) (SD) average age, (22.97±1.56) body mass index (BMI), and (1:1) male/female ratio. Four patients were diagnosed as having rectal cancers while the other half had colon cancers (ascending colon, transverse colon, descending colon, and sigmoid). From the eight patients, we obtained a total 16 tissue samples, including eight cancerous tissues and their matched adjacent normal tissues obtained via colonoscopy following the Standard Operating Procedures at the First People’s Hospital of Yunnan Province of China, Kunming. Participants were enrolled in the study prior to colonoscopy, and written informed consent was obtained from the patient for publication of this report and any accompanying images. Samples were removed endoscopically and immediately transferred from the colonoscopy room to the pathology suite and subsequently evaluated by the pathologist. Where possible one fragment of healthy tissue and one of tumor were chosen and placed in a cryotube, then frozen immediately in liquid nitrogen. All specimens were stored in their original tubes at -80°C prior to DNA extraction. This study conformed to the ethical guidelines outlined in the 1975 Declaration of Helsinki as reflected by a priori approval from the Medical Ethics Board of the First People's Hospital of Yunnan Province of China.
While frozen, an aliquot (~25 mg) of each specimen was suspended in a solution containing 200 ul buffer ATL (QIAGEN Kit Buffer for tissue lysis) and 200 ul of a slurry of 0.1-mm-diameter zirconia/silica beads (BioSpec Products, Bartlesville, OK). The mixed sample was then lysed by mechanical disruption with a bead beater (BioSpec Products), set on high for 2 min (20°C), followed by extraction with the QIAamp® DNA Mini Kit (Qiagen, Inc., Valencia CA). DNA from tissue was eluted in a final volume of 200 ul elution buffer and stored at -20°C. Tubes containing only QIAamp® DNA Mini Kit extraction controls were included throughout lysis and PCR to serve as negative controls.
PCR amplification of V1-V2 hypervariable regions of 16S rRNA gene and pyrosequencing
The forward primer included the 454 Life Sciences primer B sequence (5′ - CTATGCGCCTTGCCAGCCCGCTCAG -3′) and the broadly conserved bacterial primer 27 F (5′-AGAGTTTGATCCTGGCTCAG-3′). The reverse primer included the 454 Life Sciences primer A sequence (5′-CGTATCGCCTCCCTCGCGCCATCAG- 3′), a unique 10-nt barcode used to tag each PCR product, and the broad-range bacterial primer 338R (5′- TGCTGCCTCCCGTAGGAGT-3′). Using the primer pair described above, triplicate PCR reactions were performed on each sample. Each resulting 25 ul reaction contained 0.2 μM forward and reverse primers, 3 μl template DNA, and 2.5 ul 10X PCR buffer plus Mg2+ (TaKaRa), 2.0 ul dNTP (2.5 mM each) (TaKaRa), 0.75 ul DMSO (100%), 0.25 ul TaKaRa Taq™ (5 U/μl). Thermal cycling was conducted at 95°C for 10 minutes, followed by 30 cycles of 95°C for 30 seconds, 52°C for 30 seconds, and 72°C for 90 seconds, with a final extension of 10 minutes at 72°C. Replicate amplicons were pooled and visualized on 1.5% agarose gels using EB stain in 0.5X TE. Amplicons were cleaned using MinElute® Gel Extraction Kit (Qiagen) according to the manufacturer’s instructions. Amplicon DNA concentrations were determined using the Quant-iT PicoGreen dsDNA reagent and kit (Invitrogen). Using the amplicon pool, we carried out pyrosequencing using primer A and Titanium chemistry on a 454 Life Sciences Genome Sequencer FLX instrument (Roche) at the DNA Sequencing Facility of the Kunming Institute of Zoology, Chinese Academy of Sciences.
Sequences were processed and analyzed using Qiime 1.4 . Sequences were assigned to each sample by examining the 10-nt barcode based on the following criteria: a minimum (min) sequence length of 200 nt; maximum (max) sequence length of 400 nt; min qual score of 25; maximum number of errors in barcode of 0; maximum length of homopolymer run of 6; number of mismatches in primer of 0; excluding ambiguous and unassigned characters; and removed 454-adaptor B and 27 F bacterial primer from all assigned sequence data. Sequencing errors were removed from filtered sequences using denoiser 0.91 . Using Chimera Slayer , chimera sequences arising from the PCR amplification were detected and excluded from the denoised sequences. The chimera-free sequences were then clustered into operational taxonomic units (OTUs) using CD-hit  with a criterion of a minimum identity of 97%. Representative sequences per OTU were classified using BLAST with default parameters in Qiime 1.4. All OTUs found in at least two samples were retained for performing the following further analyses.
General characteristics were expressed as median and mean or percentages. Comparisons were performed between tumor and normal samples using the Mann-Whitney rank sum test or t-test in SigmaPlot 12.0 (Systat Software, Inc.). Statistical significance was set at P < 0.05.
Our observations over the course of the study suggest that we have discovered a novel pattern to the human CRC microbiome among Han Chinese. According to the bacterial driver-passenger model for CRC put forward by Tjalsma et al (2012) , among 50% of Chinese CRC patients in this study, the over-represented Roseburia bacteria at tumor sites should in fact be considered as ‘passenger bacteria’ for CRC (previously known in other populations, but not previously observed in Chinese) and the over-represented Microbacterium and Anoxybacillus bacteria away from tumor sites—i.e., adjacent non-malignant tissue—as ‘driver bacteria’ for CRC (novel among all populations).
Consistent with several previous reports [3, 4, 11–13], we also found an increasing trend of Fusobacterium spp. in tumors among 87.5% of patients, though we did not discover a significant increase. Across several of the different, previously studied populations including Europeans, Americans, and Asians (in particular, Chinese and Vietnamese), a similar pattern was observed; such a consistent overrepresentation of Fusobacterium spp. in tumor tissue suggests that there exists a common mechanism of gut microbial disorder connected with CRC. This finding implies that it is crucial to begin defining an underlying association of the gut passenger bacteria Fusobacterium with an increasing risk of CRC for most populations, despite the lack of relative reports for many different populations to date.
We also found another gut bacterium, Roseburia, is potentially associated with an increased risk of CRC, due to the overrepresentation of Roseburia in tumor tissue among Chinese, as we noted in this study, and as previously reported among Dutch . This finding is contrary to the study done by Wang et al study on another Chinese population  that compared the difference of stool flora between CRC patients and healthy subjects. Among the Chinese, the distinct discovery of Roseburia and its potential association with CRC may be due to different sample types, as another study done by Chen et al indicated a different microbial structure between the intestinal lumen and cancerous tissue in Chinese CRC patients. The microbial structure difference between the intestinal lumen and mucosa tissue in healthy subjects was further confirmed by Eckburg et al. Additionally, two factors —diet and genetics—may have minor effects on the differentiation of Roseburia for CRC among Chinese, as we observed consistent overrepresentation of Roseburia in tumor tissues between the Chinese and Dutch, two populations with higher divergence of diet and genetics than within the Chinese population. Accordingly, the potential role of candidate gut passenger bacteria Roseburia should be emphasized more heavily during the occurrence of CRC, regardless of the population (e.g. Chinese or Dutch) and further investigated.
Similar to the relationship between CRC and gut bacteria Fusobacterium and Roseburia described above, the ‘passenger bacteria’ role of Bacteroides for CRC was supported by studies on the Chinese , Dutch , and French . And the ‘driver bacteria’ role of Bacteroides was likewise confirmed by two other studies that examined a similar pattern among the Spanish, American and Vietnamese [3, 4]. Furthermore, Wu et al confirmed that gut bacteria Bacteroides fragilis enables the promotion of colon tumorigenesis . Similarly, based on in vivo experiments of mice model, E. coli-induced colitis is a driving factor of colorectal cancer , and the ‘passenger bacteria’ role of Escherichia species for CRC among the Chinese was also implied by the study done by Wang et al.. Moreover, besides the common pattern of the human CRC microbiome represented by gut passenger bacteria Fusobacterium within and between populations, there exists a diversified pattern in the human CRC microbiome due to three possible factors. One may be due largely to high variation of the normal human gut microbiome , potentially associated with diet [17, 18], age [18, 19], sample type (mucosa or stool) , host genetic factors , or other factors, such as antibiotic abuse . Another factor may result from the different stage of tumor progression that is randomly selected by different studies, as the CRC microbiome variance may be temporally associated with developing tumors . The last possible factor is that most findings to date have only been derived from bacterial 16S rRNA-based analyses, though some evidence from a metagenomic approach with a markedly more powerful ability to decipher the landscape of human CRC microbiome are intriguing [3, 4]. Further studies of this kind will be helpful in confirming and elucidating the potential associations we have outlined in the present study.
In summary, in the present study we presented some initial findings the lead towards a deeper and more comprehensive view of the human CRC microbiome. The existing findings are suggestive of further research, and underscore the necessity of borrowing from both high-throughput meta-genomic or transcriptomic data and (animal) model experiments that will better define and validate the association of high-risk microbial populations with occurrence of CRCs across different populations.
JG, HF, XT, HZ: Department of Gastroenterology, The First People’s Hospital of Yunnan Province, Kunming, China. ZZ: State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, China.
Operational taxonomic units
The authors would like to thank Z.D. and S.X. (Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, China) for their superb technical assistance and Andrew Willden (Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, China) for improving the language of the manuscript. This research was supported by grants from Natural Science Foundation of Yunnan Province of China (Grant No. 2010CD191 and Grant No. 2011FA035) and from National Natural Science Foundation of China (NSFC) (Grant No. 31100916).
The 454 sequence date set was deposited and available in the NCBI Sequence Read Archive (SRA) (accession number SRA060808).
- Wu S, Rhee K-J, Albesiano E, Rabizadeh S, Wu X, Yen H-R, Huso DL, Brancati FL, Wick E, McAllister F: A human colonic commensal promotes colon tumorigenesis via activation of T helper type 17 T cell responses. Nat Med. 2009, 15 (9): 1016-1022. 10.1038/nm.2015.PubMed CentralView ArticlePubMedGoogle Scholar
- Arthur JC, Perez-Chanona E, Mühlbauer M, Tomkovich S, Uronis JM, Fan T-J, Campbell BJ, Abujamel T, Dogan B, Rogers AB: Intestinal inflammation targets cancer-inducing activity of the microbiota. Science. 2012, 338 (6103): 120-123. 10.1126/science.1224820.PubMed CentralView ArticlePubMedGoogle Scholar
- Castellarin M, Warren RL, Freeman JD, Dreolini L, Krzywinski M, Strauss J, Barnes R, Watson P, Allen-Vercoe E, Moore RA: Fusobacterium nucleatum infection is prevalent in human colorectal carcinoma. Genome Res. 2012, 22 (2): 299-306. 10.1101/gr.126516.111.PubMed CentralView ArticlePubMedGoogle Scholar
- Kostic AD, Gevers D, Pedamallu CS, Michaud M, Duke F, Earl AM, Ojesina AI, Jung J, Bass AJ, Tabernero J: Genomic analysis identifies association of Fusobacterium with colorectal carcinoma. Genome Res. 2012, 22 (2): 292-298. 10.1101/gr.126573.111.PubMed CentralView ArticlePubMedGoogle Scholar
- Tjalsma H, Boleij A, Marchesi JR, Dutilh BE: A bacterial driver–passenger model for colorectal cancer: beyond the usual suspects. Nat Rev Micro. 2012, 10 (8): 575-582. 10.1038/nrmicro2819.View ArticleGoogle Scholar
- Ahmed FE: Gene-gene, gene-environment & multiple interactions in colorectal cancer. J. Environ. Sci. Health. 2006, 24 (1): 1-101.View ArticleGoogle Scholar
- Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, Fierer N, Pena AG, Goodrich JK, Gordon JI: QIIME allows analysis of high-throughput community sequencing data. Nat Meth. 2010, 7 (5): 335-336. 10.1038/nmeth.f.303.View ArticleGoogle Scholar
- Reeder J, Knight R: Rapidly denoising pyrosequencing amplicon reads by exploiting rank-abundance distributions. Nat Meth. 2010, 7 (9): 668-669.View ArticleGoogle Scholar
- Haas BJ, Gevers D, Earl A, Feldgarden M, Ward DV, Giannokous G, Ciulla D, Tabbaa D, Highlander SK, Sodergren E: Chimeric 16S rRNA sequence formation and detection in Sanger and 454-pyrosequenced PCR amplicons. Genome Res. 2011, 21: 494-504. 10.1101/gr.112730.110.PubMed CentralView ArticlePubMedGoogle Scholar
- Li W, Godzik A: Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics. 2006, 22 (13): 1658-1659. 10.1093/bioinformatics/btl158.View ArticlePubMedGoogle Scholar
- Chen W, Liu F, Ling Z, Tong X, Xiang C: Human intestinal lumen and mucosa-associated microbiota in patients with colorectal cancer. PLoS One. 2012, 7 (6): e39743-10.1371/journal.pone.0039743.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang T, Cai G, Qiu Y, Fei N, Zhang M, Pang X, Jia W, Cai S, Zhao L: Structural segregation of gut microbiota between colorectal cancer patients and healthy volunteers. ISME J. 2012, 6 (2): 320-329. 10.1038/ismej.2011.109.PubMed CentralView ArticlePubMedGoogle Scholar
- Marchesi JR, Dutilh BE, Hall N, Peters WHM, Roelofs R, Boleij A, Tjalsma H: Towards the human colorectal cancer microbiome. PLoS One. 2011, 6 (5): e20447-10.1371/journal.pone.0020447.PubMed CentralView ArticlePubMedGoogle Scholar
- Eckburg PB, Bik EM, Bernstein CN, Purdom E, Dethlefsen L, Sargent M, Gill SR, Nelson KE, Relman DA: Diversity of the human intestinal microbial flora. Science. 2005, 308 (5728): 1635-1638. 10.1126/science.1110591.PubMed CentralView ArticlePubMedGoogle Scholar
- Sobhani I, Tap J, Roudot-Thoraval F, Roperch JP, Letulle S, Langella P, Corthier G, Van Nhieu JT, Furet JP: Microbial dysbiosis in colorectal cancer (CRC) patients. PLoS One. 2011, 6 (1): e16393-10.1371/journal.pone.0016393.PubMed CentralView ArticlePubMedGoogle Scholar
- Schloissnig S, Arumugam M, Sunagawa S, Mitreva M, Tap J, Zhu A, Waller A, Mende DR, Kultima JR, Martin J: Genomic variation landscape of the human gut microbiome. Nature. 2013, 493 (7430): 45-50.PubMed CentralView ArticlePubMedGoogle Scholar
- Muegge BD, Kuczynski J, Knights D, Clemente JC, González A, Fontana L, Henrissat B, Knight R, Gordon JI: Diet drives convergence in gut microbiome functions across mammalian phylogeny and within humans. Science. 2011, 332 (6032): 970-974. 10.1126/science.1198719.PubMed CentralView ArticlePubMedGoogle Scholar
- Claesson MJ, Jeffery IB, Conde S, Power SE, O/'Connor EM, Cusack S, Harris HMB, Coakley M, Lakshminarayanan B, O/'Sullivan O: Gut microbiota composition correlates with diet and health in the elderly. Nature. 2012, 488 (7410): 178-184. 10.1038/nature11319.View ArticlePubMedGoogle Scholar
- Claesson MJ, Cusack S, O'Sullivan O, Greene-Diniz R, De Weerd H, Flannery E, Marchesi JR, Falush D, Dinan T, Fitzgerald G: Composition, variability, and temporal stability of the intestinal microbiota of the elderly. Proc Natl Acad USA. 2011, 108 (Supplement 1): 4586-4591.View ArticleGoogle Scholar
- Benson AK, Kelly SA, Legge R, Ma F, Low SJ, Kim J, Zhang M, Oh PL, Nehrenberg D, Hua K: Individuality in gut microbiota composition is a complex polygenic trait shaped by multiple environmental and host genetic factors. Proc Natl Acad USA. 2010, 107 (44): 18933-18938. 10.1073/pnas.1007028107.View ArticleGoogle Scholar
- Dethlefsen L, Relman DA: Incomplete recovery and individualized responses of the human distal gut microbiota to repeated antibiotic perturbation. Proc Natl Acad USA. 2011, 108 (Supplement 1): 4554-4561.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.