Microbiome Analysis of Stool Samples from African Americans with Colon Polyps

Background

Colonic polyps are common tumors occurring in ~50% of Western populations with ~10% risk of malignant progression. Dietary agents have been considered the primary environmental exposure to promote colorectal cancer (CRC) development. However, the colonic mucosa is permanently in contact with the microbiota and its metabolic products including toxins that also have the potential to trigger oncogenic transformation.

Aim

To analyze fecal DNA for microbiota composition and functional potential in African Americans with pre-neoplastic lesions.

Materials & Methods

We analyzed the bacterial composition of stool samples from 6 healthy individuals and 6 patients with colon polyps using 16S ribosomal RNA-based phylogenetic microarray; the Human intestinal Tract Chip (HITChip) and 16S rRNA gene barcoded 454 pyrosequencing. The functional potential was determined by sequence-based metagenomics using 454 pyrosequencing.

Results

Fecal microbiota profiling of samples from the healthy and polyp patients using both a phylogenetic microarraying (HITChip) and barcoded 454 pyrosequencing generated similar results. A distinction between both sets of samples was only obtained when the analysis was performed at the sub-genus level. Most of the species leading to the dissociation were from the Bacteroides group. The metagenomic analysis did not reveal major differences in bacterial gene prevalence/abundances between the two groups even when the analysis and comparisons were restricted to available Bacteroides genomes.

Conclusion

This study reveals that at the pre-neoplastic stages, there is a trend showing microbiota changes between healthy and colon polyp patients at the sub-genus level. These differences were not reflected at the genome/functions levels. Bacteria and associated functions within the Bacteroides group need to be further analyzed and dissected to pinpoint potential actors in the early colon oncogenic transformation in a large sample size.     

16S rRNA基因在微生物生态学中的应用

16S rRNA(Small subunit ribosomal RNA)基因是对原核微生物进行系统进化分类研究时最常用的分子标志物(Biomarker),广泛应用于微生物生态学研究中。近些年来随着高通量测序技术及数据分析方法等的不断进步,大量基于16S rRNA基因的研究使得微生物生态学得到了快速发展,然而使用16S rRNA基因作为分子标志物时也存在诸多问题,比如水平基因转移、多拷贝的异质性、基因扩增效率的差异、数据分析方法的选择等,这些问题影响了微生物群落组成和多样性分析时的准确性。对当前使用16S rRNA基因分析微生物群落组成和多样性的进展情况做一总结,重点讨论当前存在的主要问题以及各种分析方法的发展,尤其是与高通量测序技术有关的实验和数据处理问题。     

Semi-automated extraction of microbial DNA from feces for qPCR and phylogenetic microarray analysis

The human gastrointestinal tract (GI-tract) harbors a complex microbial ecosystem, largely composed of so far uncultured species, which can be detected only by using techniques such as PCR and by different hybridization techniques including phylogenetic microarrays. Manual DNA extraction from feces is laborious and is one of the bottlenecks holding up the application of microarray and other DNA-based techniques in large cohort studies. In order to enhance the DNA extraction step we combined mechanical disruption of microbial cells by repeated bead-beating (RBB) with two automated DNA extraction methods, KingFisher with InviMag Stool DNA kit (KF) and NucliSENS easyMAG (NeM). The semi-automated DNA extraction methods, RBB combined with either KF or NeM, were compared to the manual extraction method currently considered the most suited method for fecal DNA extraction by assessing the yield of 16S rRNA gene copies by qPCR and total microbiota composition by the HITChip, a phylogenetic microarray. Parallel DNA extractions from infant fecal samples by using the three methods showed that the KF and manual methods gave comparable yields of 16S rRNA gene copies as assessed by qPCR, whereas NeM showed a significantly lower yield. All three methods showed highly similar microbiota profiles in HITChip. Both KF and NeM were found to be suitable methods for DNA extraction from fecal samples after the mechanical disruption of microbial cells by bead-beating. The semi-automated methods could be performed in half of the time required for the manual protocol, while being comparable to the manual method in terms of reagent costs.     

Comparative Analysis of Pyrosequencing and a Phylogenetic Microarray for Exploring Microbial Community Structures in the Human Distal Intestine

Background

Variations in the composition of the human intestinal microbiota are linked to diverse health conditions. High-throughput molecular technologies have recently elucidated microbial community structure at much higher resolution than was previously possible. Here we compare two such methods, pyrosequencing and a phylogenetic array, and evaluate classifications based on two variable 16S rRNA gene regions.

Methods and Findings

Over 1.75 million amplicon sequences were generated from the V4 and V6 regions of 16S rRNA genes in bacterial DNA extracted from four fecal samples of elderly individuals. The phylotype richness, for individual samples, was 1,400–1,800 for V4 reads and 12,500 for V6 reads, and 5,200 unique phylotypes when combining V4 reads from all samples. The RDP-classifier was more efficient for the V4 than for the far less conserved and shorter V6 region, but differences in community structure also affected efficiency. Even when analyzing only 20% of the reads, the majority of the microbial diversity was captured in two samples tested. DNA from the four samples was hybridized against the Human Intestinal Tract (HIT) Chip, a phylogenetic microarray for community profiling. Comparison of clustering of genus counts from pyrosequencing and HITChip data revealed highly similar profiles. Furthermore, correlations of sequence abundance and hybridization signal intensities were very high for lower-order ranks, but lower at family-level, which was probably due to ambiguous taxonomic groupings.

Conclusions

The RDP-classifier consistently assigned most V4 sequences from human intestinal samples down to genus-level with good accuracy and speed. This is the deepest sequencing of single gastrointestinal samples reported to date, but microbial richness levels have still not leveled out. A majority of these diversities can also be captured with five times lower sampling-depth. HITChip hybridizations and resulting community profiles correlate well with pyrosequencing-based compositions, especially for lower-order ranks, indicating high robustness of both approaches. However, incompatible grouping schemes make exact comparison difficult.     

Comparative analysis of human gut microbiota by barcoded pyrosequencing

Humans host complex microbial communities believed to contribute to health maintenance and, when in imbalance, to the development of diseases. Determining the microbial composition in patients and healthy controls may thus provide novel therapeutic targets. For this purpose, high-throughput, cost-effective methods for microbiota characterization are needed. We have employed 454-pyrosequencing of a hyper-variable region of the 16S rRNA gene in combination with sample-specific barcode sequences which enables parallel in-depth analysis of hundreds of samples with limited sample processing. In silico modeling demonstrated that the method correctly describes microbial communities down to phylotypes below the genus level. Here we applied the technique to analyze microbial communities in throat, stomach and fecal samples. Our results demonstrate the applicability of barcoded pyrosequencing as a high-throughput method for comparative microbial ecology.     

Optimizing the analysis of human intestinal microbiota with phylogenetic microarray

Phylogenetic microarrays present an attractive strategy to high-throughput interrogation of complex microbial communities. In this work, we present several approaches to optimize the analysis of intestinal microbiota with the recently developed Microbiota Array. First, we determined how 16S rDNA-specific PCR amplification influenced bacterial detection and the consistency of measured abundance values. Bacterial detection improved with an increase in the number of PCR amplification cycles, but 25 cycles were sufficient to achieve the maximum possible detection. A PCR-caused deviation in the measured abundance values was also observed. We also developed two mathematical algorithms that aimed to account for a predicted cross-hybridization of 16S rDNA fragments among different species, and to adjust the measured hybridization signal based on the number of 16S rRNA gene copies per species genome. The 16S rRNA gene copy adjustment indicated that the presence of members of the class Clostridia might be overestimated in some 16S rDNA-based studies. Finally, we show that the examination of total community RNA with phylogenetic microarray can provide estimates of the relative metabolic activity of individual community members. Complementary profiling of genomic DNA and total RNA isolated from the same sample presents an opportunity to assess population structure and activity in the same microbial community.     

基于不同通用引物比较分析肠道微生物群落变化

肠道微生物对于人体健康的重要作用已经得到广泛证实,目前,对肠道微生物的研究大多采用基于扩增细菌16S rRNA基因V3-V4区的高通量测序分析,对古菌的关注较少。本研究选择了一对可以同时扩增细菌和古菌16S rRNA基因的引物,通过比较人为干扰肠道微生物前后的群落变化,说明这对引物适宜分析人类肠道细菌和古菌群落变化并具有一定优越性。采集志愿者粪便样品,同时用仅能扩增细菌引物 (B引物) 和细菌古菌通用引物 (AB引物) 进行扩增和高通量测序;使用几个常用的rRNA数据库判断引物对细菌的覆盖度和对古菌的扩增能力。结果表明,AB引物在可以展示B引物扩增出的细菌群落的基础上,可以得到肠道中常见的产甲烷古菌的序列,同时也展示出人为干扰肠道微生物前后的群落结构变化。AB引物可以仅通过一次扩增和测序同时分析肠道中的细菌和古菌群落,更加全面展示肠道微生物群落结构,适用于肠道微生物相关研究。     

Assessing the Fecal Microbiota: An Optimized Ion Torrent 16S rRNA Gene-Based Analysis Protocol

Assessing the distribution of 16S rRNA gene sequences within a biological sample represents the current state-of-the-art for determination of human gut microbiota composition. Advances in dissecting the microbial biodiversity of this ecosystem have very much been dependent on the development of novel high-throughput DNA sequencing technologies, like the Ion Torrent. However, the precise representation of this bacterial community may be affected by the protocols used for DNA extraction as well as by the PCR primers employed in the amplification reaction. Here, we describe an optimized protocol for 16S rRNA gene-based profiling of the fecal microbiota.     

Novel high-rank phylogenetic lineages within a sulfur spring (Zodletone Spring, Oklahoma), revealed using a combined pyrosequencing-sanger approach

The utilization of high-throughput sequencing technologies in 16S rRNA gene-based diversity surveys has indicated that within most ecosystems, a significant fraction of the community could not be assigned to known microbial phyla. Accurate determination of the phylogenetic affiliation of such sequences is difficult due to the short-read-length output of currently available high-throughput technologies. This fraction could harbor multiple novel phylogenetic lineages that have so far escaped detection. Here we describe our efforts in accurate assessment of the novelty and phylogenetic affiliation of selected unclassified lineages within a pyrosequencing data set generated from source sediments of Zodletone Spring, a sulfide- and sulfur-rich spring in southwestern Oklahoma. Lineage-specific forward primers were designed for 78 putatively novel lineages identified within the pyrosequencing data set, and representative nearly full-length small-subunit (SSU) rRNA gene sequences were obtained by pairing those primers with reverse universal bacterial primers. Of the 78 lineages tested, amplifiable products were obtained for 52, 32 of which had at least one nearly full-length sequence that was representative of the lineage targeted. Analysis of phylogenetic affiliation of the obtained Sanger sequences identified 5 novel candidate phyla and 10 novel candidate classes (within Fibrobacteres, Planctomycetes, and candidate phyla BRC1, GN12, TM6, TM7, LD1, WS2, and GN06) in the data set, in addition to multiple novel orders and families. The discovery of multiple novel phyla within a pilot study of a single ecosystem clearly shows the potential of the approach in identifying novel diversities within the rare biosphere.     

Mathematical modeling of 16S ribosomal DNA amplification reveals optimal conditions for the interrogation of complex microbial communities with phylogenetic microarrays

MOTIVATION: Many current studies of complex microbial communities rely on the isolation of community genomic DNA, amplification of 16S ribosomal RNA genes (rDNA) and subsequent examination of community structure through interrogation of the amplified 16S rDNA pool by high-throughput sequencing, phylogenetic microarrays or quantitative PCR. RESULTS: Here we describe the development of a mathematical model aimed to simulate multitemplate amplification of 16S ribosomal DNA sample and subsequent detection of these amplified 16S rDNA species by phylogenetic microarray. Using parameters estimated from the experimental results obtained in the analysis of intestinal microbial communities with Microbiota Array, we show that both species detection and the accuracy of species abundance estimates depended heavily on the number of PCR cycles used to amplify 16S rDNA. Both parameters initially improved with each additional PCR cycle and reached optimum between 15 and 20 cycles of amplification. The use of more than 20 cycles of PCR amplification and/or more than 50 ng of starting genomic DNA template was, however, detrimental to both the fraction of detected community members and the accuracy of abundance estimates. Overall, the outcomes of the model simulations matched well available experimental data. Our simulations also showed that species detection and the accuracy of abundance measurements correlated positively with the higher sample-wide PCR amplification rate, lower template-to-template PCR bias and lower number of species in the interrogated community. The developed model can be easily modified to simulate other multitemplate DNA mixtures as well as other microarray designs and PCR amplification protocols.     

Design and Experimental Application of a Novel Non-Degenerate Universal Primer Set that Amplifies Prokaryotic 16S rRNA Genes with a Low Possibility to Amplify Eukaryotic rRNA Genes

The deep sequencing of 16S rRNA genes amplified by universal primers has revolutionized our understanding of microbial communities by allowing the characterization of the diversity of the uncultured majority. However, some universal primers also amplify eukaryotic rRNA genes, leading to a decrease in the efficiency of sequencing of prokaryotic 16S rRNA genes with possible mischaracterization of the diversity in the microbial community. In this study, we compared 16S rRNA gene sequences from genome-sequenced strains and identified candidates for non-degenerate universal primers that could be used for the amplification of prokaryotic 16S rRNA genes. The 50 identified candidates were investigated to calculate their coverage for prokaryotic and eukaryotic rRNA genes, including those from uncultured taxa and eukaryotic organelles, and a novel universal primer set, 342F-806R, covering many prokaryotic, but not eukaryotic, rRNA genes was identified. This primer set was validated by the amplification of 16S rRNA genes from a soil metagenomic sample and subsequent pyrosequencing using the Roche 454 platform. The same sample was also used for pyrosequencing of the amplicons by employing a commonly used primer set, 338F-533R, and for shotgun metagenomic sequencing using the Illumina platform. Our comparison of the taxonomic compositions inferred by the three sequencing experiments indicated that the non-degenerate 342F-806R primer set can characterize the taxonomic composition of the microbial community without substantial bias, and is highly expected to be applicable to the analysis of a wide variety of microbial communities.     

Coamplification of Eukaryotic DNA with 16S rRNA Gene-Based PCR Primers: Possible Consequences for Population Fingerprinting of Complex Microbial Communities

The main aim of this study was to evaluate the specificity of three commonly used 16S rRNA gene-based polymerase chain reaction (PCR) primer sets for bacterial community analysis of samples contaminated with eukaryotic DNA. The specificity of primer sets targeting the V₃, V₃-V₅, and V₆-V₈ hypervariable regions of the 16S rRNA gene was investigated in silico and by community fingerprinting of human and fish intestinal samples. Both in silico and PCR-based analysis revealed cross-reactivity of the V₃ and V₃-V₅ primers with the 18S rRNA gene of human and sturgeon. The consequences of this primer anomaly were illustrated by denaturing gradient gel electrophoresis (DGGE) profiling of microbial communities in human feces and mixed gut of Siberian sturgeon. DGGE profiling indicated that the cross-reactivity of 16S rRNA gene primers with nontarget eukaryotic DNA might lead to an overestimation of bacterial biodiversity. This study has confirmed previous sporadic indications in literature indicating that several commonly applied 16S rRNA gene primer sets lack specificity toward bacteria in the presence of eukaryotic DNA. The phenomenon of cross-reactivity is a potential source of systematic error in all biodiversity studies where no subsequent analysis of individual community amplicons by cloning and sequencing is performed.     

The Microbial Community in the Feces of the White Rhinoceros (Ceratotherium simum) as Determined by Barcoded Pyrosequencing Analysis

As a non-ruminant herbivore, the white rhinoceros has the ability to utilize fibrous plant matter through microbial fermentation in the hindgut. So far, there has been no report using molecular techniques to study the gut microbiota of the white rhinoceros. We used barcoded pyrosequencing to characterize 105,651 sequences of 16S rRNA genes obtained from fecal samples from five white rhinoceroses. Results showed that Firmicutes and Bacteroidetes were the predominant phyla in the samples, which were comprised largely of unclassified bacteria. The microbiota of one animal treated with drug therapy differed from those in other healthy animals, and was dominated by Aerococcus -related bacteria. The core microbiota in the healthy rhinoceros were dominated by phyla Firmicutes and Bacteroidetes, represented by the Ruminococcaceae, Lachnospiraceae, Rikenellaceae and Prevotellaceae families. The present work provides a phylogenetic framework for understanding the complex microbial community of the rhinoceros; however, further studies are required to link the distinctive microbiota with their digestive role in the hindgut of the white rhinoceros.     

Bacterial Community Mapping of the Mouse Gastrointestinal Tract

Keeping mammalian gastrointestinal (GI) tract communities in balance is crucial for host health maintenance. However, our understanding of microbial communities in the GI tract is still very limited. In this study, samples taken from the GI tracts of C57BL/6 mice were subjected to 16S rRNA gene sequence-based analysis to examine the characteristic bacterial communities along the mouse GI tract, including those present in the stomach, duodenum, jejunum, ileum, cecum, colon and feces. Further analyses of the 283,234 valid sequences obtained from pyrosequencing revealed that the gastric, duodenal, large intestinal and fecal samples had higher phylogenetic diversity than the jejunum and ileum samples did. The microbial communities found in the small intestine and stomach were different from those seen in the large intestine and fecal samples. A greater proportion of Lactobacillaceae were found in the stomach and small intestine, while a larger proportion of anaerobes such as Bacteroidaceae, Prevotellaceae, Rikenellaceae, Lachnospiraceae, and Ruminococcaceae were found in the large intestine and feces. In addition, inter-mouse variations of microbiota were observed between the large intestinal and fecal samples, which were much smaller than those between the gastric and small intestinal samples. As far as we can ascertain, ours is the first study to systematically characterize bacterial communities from the GI tracts of C57BL/6 mice.     

Oral microbiome profiles: 16S rRNA pyrosequencing and microarray assay comparison

Objectives

The human oral microbiome is potentially related to diverse health conditions and high-throughput technology provides the possibility of surveying microbial community structure at high resolution. We compared two oral microbiome survey methods: broad-based microbiome identification by 16S rRNA gene sequencing and targeted characterization of microbes by custom DNA microarray.

Methods

Oral wash samples were collected from 20 individuals at Memorial Sloan-Kettering Cancer Center. 16S rRNA gene survey was performed by 454 pyrosequencing of the V3–V5 region (450 bp). Targeted identification by DNA microarray was carried out with the Human Oral Microbe Identification Microarray (HOMIM). Correlations and relative abundance were compared at phylum and genus level, between 16S rRNA sequence read ratio and HOMIM hybridization intensity.

Results

The major phyla, Firmicutes, Proteobacteria, Bacteroidetes, Actinobacteria, and Fusobacteria were identified with high correlation by the two methods (r = 0.70∼0.86). 16S rRNA gene pyrosequencing identified 77 genera and HOMIM identified 49, with 37 genera detected by both methods; more than 98% of classified bacteria were assigned in these 37 genera. Concordance by the two assays (presence/absence) and correlations were high for common genera (Streptococcus, Veillonella, Leptotrichia, Prevotella, and Haemophilus; Correlation = 0.70–0.84).

Conclusion

Microbiome community profiles assessed by 16S rRNA pyrosequencing and HOMIM were highly correlated at the phylum level and, when comparing the more commonly detected taxa, also at the genus level. Both methods are currently suitable for high-throughput epidemiologic investigations relating identified and more common oral microbial taxa to disease risk; yet, pyrosequencing may provide a broader spectrum of taxa identification, a distinct sequence-read record, and greater detection sensitivity.     

Evaluation of new gyrB-based real-time PCR system for the detection of B. fragilis as an indicator of human-specific fecal contamination

A rapid and specific gyrB-based real-time PCR system has been developed for detecting Bacteroides fragilis as a human-specific marker of fecal contamination. Its specificity and sensitivity was evaluated by comparison with other 16S rRNA gene-based primers using closely related Bacteroides and Prevotella. Many studies have used 16S rRNA gene-based method targeting Bacteroides because this genus is relatively abundant in human feces and is useful for microbial source tracking. However, 16S rRNA gene-based primers are evolutionarily too conserved among taxa to discriminate between human-specific species of Bacteroides and other closely related genera, such as Prevotella. Recently, one of the housekeeping genes, gyrB, has been used as an alternative target in multilocus sequence analysis (MLSA) to provide greater phylogenetic resolution. In this study, a new B. fragilis-specific primer set (Bf904F/Bf958R) was designed by alignments of 322 gyrB genes and was compared with the performance of the 16S rRNA gene-based primers in the presence of B. fragilis, Bacteroides ovatus and Prevotella melaninogenica. Amplicons were sequenced and a phylogenetic tree was constructed to confirm the specificity of the primers to B. fragilis. The gyrB-based primers successfully discriminated B. fragilis from B. ovatus and P. melaninogenica. Real-time PCR results showed that the gyrB primer set had a comparable sensitivity in the detection of B. fragilis when compared with the 16S rRNA primer set. The host-specificity of our gyrB-based primer set was validated with human, pig, cow, and dog fecal samples. The gyrB primer system had superior human-specificity. The gyrB-based system can rapidly detect human-specific fecal source and can be used for improved source tracking of human contamination.     

Substitution by Inosine at the 3��-Ultimate and Penultimate Positions of 16S rRNA Gene Universal Primers

Universal 16S rRNA gene primers (8F and 518R) bearing inosine substitutions at either the 3??-ultimate or the 3??-ultimate and penultimate base positions were exploited for the first time to study the bacterial community associated with coral polymicrobial Black Band Disease (BBD). Inosine-modified universal primer pairs display some shifting in the composition of 16S rRNA gene libraries, as well as expanding the observed diversity of a BBD bacterial community at the family/class level. Possible explanations for the observed shifts are discussed. These results thus point to the need for adopting multiple approaches in designing 16S rRNA universal primers for PCR amplification and subsequent construction of 16S rRNA gene libraries or pyrosequencing in the exploration of complex microbial communities.     

Estimating Bacterial Diversity for Ecological Studies: Methods,Metrics, and Assumptions

Methods to estimate microbial diversity have developed rapidly in an effort to understand the distribution and diversity of microorganisms in natural environments. For bacterial communities, the 16S rRNA gene is the phylogenetic marker gene of choice, but most studies select only a specific region of the 16S rRNA to estimate bacterial diversity. Whereas biases derived from from DNA extraction, primer choice and PCR amplification are well documented, we here address how the choice of variable region can influence a wide range of standard ecological metrics, such as species richness, phylogenetic diversity, β-diversity and rank-abundance distributions. We have used Illumina paired-end sequencing to estimate the bacterial diversity of 20 natural lakes across Switzerland derived from three trimmed variable 16S rRNA regions (V3, V4, V5). Species richness, phylogenetic diversity, community composition, β-diversity, and rank-abundance distributions differed significantly between 16S rRNA regions. Overall, patterns of diversity quantified by the V3 and V5 regions were more similar to one another than those assessed by the V4 region. Similar results were obtained when analyzing the datasets with different sequence similarity thresholds used during sequences clustering and when the same analysis was used on a reference dataset of sequences from the Greengenes database. In addition we also measured species richness from the same lake samples using ARISA Fingerprinting, but did not find a strong relationship between species richness estimated by Illumina and ARISA. We conclude that the selection of 16S rRNA region significantly influences the estimation of bacterial diversity and species distributions and that caution is warranted when comparing data from different variable regions as well as when using different sequencing techniques.     

Biased Diversity Metrics Revealed by Bacterial 16S Pyrotags Derived from Different Primer Sets

In recent years, PCR-based pyrosequencing of 16S rRNA genes has continuously increased our understanding of complex microbial communities in various environments of the Earth. However, there is always concern on the potential biases of diversity determination using different 16S rRNA gene primer sets and covered regions. Here, we first report how bacterial 16S rRNA gene pyrotags derived from a series of different primer sets resulted in the biased diversity metrics. In total, 14 types of pyrotags were obtained from two-end pyrosequencing of 7 amplicon pools generated by 7 primer sets paired by 1 of 4 forward primers (V1F, V3F, V5F, and V7F) and 1 of 4 reverse primers (V2R, V4R, V6R, and V9R), respectively. The results revealed that: i) the activated sludge exhibited a large bacterial diversity that represented a broad range of bacterial populations and served as a good sample in this methodology research; ii) diversity metrics highly depended on the selected primer sets and covered regions; iii) paired pyrotags obtained from two-end pyrosequencing of each short amplicon displayed different diversity metrics; iv) relative abundance analysis indicated the sequencing depth affected the determination of rare bacteria but not abundant bacteria; v) the primer set of V1F and V2R significantly underestimated the diversity of activated sludge; and vi) the primer set of V3F and V4R was highly recommended for future studies due to its advantages over other primer sets. All of these findings highlight the significance of this methodology research and offer a valuable reference for peer researchers working on microbial diversity determination.