Variations in 16S rRNA-based microbiome profiling between pyrosequencing runs and between pyrosequencing facilities

Pyrosequencing of 16S rRNA gene amplicons on the 454 FLX Titanium platform has been widely used to analyze microbiomes in various environments. However, different results may stem from variations among sequencing runs or among sequencing facilities. This study aimed to evaluate these variations between different pyrosequencing runs by sequencing 16S rRNA gene amplicon libraries generated from three sets of rumen samples twice each on the 454 FLX Titanium system at two independent sequencing facilities. Similar relative abundances were found for predominant taxa represented by large numbers of sequence reads but not for minor taxa represented by small numbers of sequence reads. The two sequencing facilities revealed different bacterial profiles with respect to both predominant taxa and minor taxa, including the most predominant genus Prevotella, the family Lachnospiraceae, and the phylum Proteobacteria. Differences in primers used to generate amplicon libraries may be a major source of variations in microbiome profiling. Because different primers and regions of 16S rRNA genes are often used by different researchers, significant variations likely exist among studies. Quantitative interpretation for relative abundance of taxa, especially minor taxa, from prevalence of sequence reads and comparisons of results from different studies should be done with caution.     

Genetic and phenotypic aspects of foot lesion scores in sheep of different breeds and ages

Footrot is a costly endemic disease of sheep. This study investigates the potential to decrease its prevalence through selective breeding for decreased lesion score. Pedigreed mule and Scottish Blackface (SBF) ewes were scored for lesions on each hoof on a 0 to 4 scale for up to 2 (SBF ewes) or 4 (mules) times over 2 years. One score was obtained for SBF lambs. An animal was deemed to have lesions (severe lesions) if at least one hoof had a score of at least 1 (2). The prevalence of lesions was 34% in lambs, 17% in SBF ewes and 51% in mules. The heritability of lesions (severe lesions) analysed as repeated measurements of the same trait in a threshold model was 0.19 (0.26) in SBF ewes and 0.12 (0.19) in mules. Estimates for the sum and maximum of scores as well as the number of feet affected were much lower, as were estimates for permanent animal effects (i.e. non-genetic effects associated with an animal). When successive scores on the same animal were analysed as correlated traits, heritability estimates for most traits tended to be higher, except for severe footrot in mules where estimates varied greatly over time. The phenotypic correlations between successive scores in SBF ewes were close to 0, genetic correlations were moderately positive (0.18 to 0.55). Correlations in mules were generally of a similar size, but some genetic correlations were higher (up to 0.92). There was a clear trend for heritabilities for lesions and severe lesions to increase with higher prevalence of lesions, even when analysed in a threshold model. Heritability estimates for traits that combine scores over several events in mules, identifying the more persistently affected animals, ranged from 0.12 to 0.23 with the highest estimates for the average number of feet that were (severely) affected in animals scored for a minimum at two events. The heritability of all lesion traits in lambs was estimated as 0. It is concluded that selection for lower lesions is possible in ewes but not lambs, and that a simple binary score at an animal level is at least as effective as a comprehensive score at hoof level. Given the low repeatability of lesion scores, repeated measures over time will improve effectiveness of selection. Selection across environments (flocks, seasons) with different prevalences of lesions scores will need to take account of variation in the heritability.     

16S rRNA-based analysis of microbiota from the cecum of broiler chickens.

The microbiota of the intestinal tract of chickens plays an important role in inhibiting the establishment of intestinal pathogens. Earlier culturing and microscopic examinations indicated that only a fraction of the bacteria in the cecum of chickens could be grown in the laboratory. Therefore, a survey of cecal bacteria was done by retrieval of 16S rRNA gene sequences from DNA isolated from the cecal content and the cecal mucosa. The ribosomal gene sequences were amplified with universal primers and cloned or subjected to temporal temperature gradient gel electrophoresis (TTGE). Partial 16S rRNA gene sequences were determined from the clones and from the major bands in TTGE gels. A total of 1,656 partial 16S rRNA gene sequences were obtained and compared to sequences in the GenBank. The comparison indicated that 243 different sequences were present in the samples. Overall, sequences representing 50 phylogenetic groups or subgroups of bacteria were found, but approximately 89% of the sequences represented just four phylogenetic groups (Clostridium leptum, Sporomusa sp., Clostridium coccoides, and enterics). Sequences of members of the Bacteroides group, the Bifidobacterium infantis subgroup, and of Pseudomonas sp. each accounted for less than 2% of the total. Sequences related to those from the Escherichia sp. subgroup and from Lactobacillus, Pseudomonas, and Bifidobacterium spp. were generally between 98 and 100% identical to sequences already deposited in the GenBank. Sequences most closely related to those of the other bacteria were generally 97% or less identical to those in the databases and therefore might be from currently unknown species. TTGE and random cloning indicated that certain phylogenetic subgroups were common to all birds analyzed, but sequence data from random cloning also provided evidence for qualitative and quantitative differences among the cecal microbiota of individual birds reared under very similar conditions.     

基于16S rDNA高通量测序的短尾猴口腔微生物多样性

口腔微生物群落结构是维持机体健康的重要因素,了解动物口腔微生物多样性有助于认识和理解动物的生态学适应。本研究以栖息于安徽黄山的短尾猴为研究对象,采用非损伤性取样法收集了鱼鳞坑YA1群体中19个短尾猴口腔样本,采用改进高盐提取法提取微生物DNA,利用Illumina Miseq测序平台对微生物16S rDNA V3-V4区扩增产物进行双端测序,分析微生物群落结构多样性。研究共获得206 533条优质序列,发现4 685个OTU,归属20个门、310个属。结果表明：短尾猴口腔微生物物种丰富,以变形菌门(Proteobacteria,占总条带44.58%)、厚壁菌门(Firmicutes, 30.28%)、拟杆菌门(Bacteroidetes, 12.27%)、梭杆菌门(Fusobacteria, 7.72%) 和放线菌门(Actinobacteria,3.70%)为主;24个微生物属在所有样本中均有分布,为其核心微生物属;短尾猴口腔中存在大量与口腔疾病相关的微生物和多种低丰富度的潜在病原菌。本研究为进一步研究短尾猴口腔微生物群落结构形成与适应性提供了基础资料,也提示在保护和管理野生猴群中需要应对人畜共患病传播的潜在风险。     

Polyphyly in 16S rRNA-based LVTree Versus Monophyly in Whole-genome-based CVTree

We report an important but long-overlooked manifestation of low-resolution power of 16S rRNA sequence analysis at the species level, namely, in 16S rRNA-based phylogenetic trees polyphyletic placements of closely-related species are abundant compared to those in genome-based phylogeny. This phenomenon makes the demarcation of genera within many families ambiguous in the 16S rRNA-based taxonomy. In this study, we reconstructed phylogenetic relationship for more than ten thousand prokaryote genomes using the CVTree method, which is based on whole-genome information. And many such genera, which are polyphyletic in 16S rRNA-based trees, are well resolved as monophyletic clusters by CVTree. We believe that with genome sequencing of prokaryotes becoming a commonplace, genome-based phylogeny is doomed to play a definitive role in the construction of a natural and objective taxonomy.     

Designation of Streptomycete 16S and 23S rRNA-based target regions for oligonucleotide probes

The 16S and 23S rRNA of various Streptomyces species were partially sequenced and screened for the presence of stretches that could define all members of the genus, groups of species, or individual species. Nucleotide 929 (Streptomyces ambofaciens nomenclature [J.L. Pernodet, M.T. Alegre, F. Boccard, and M. Guerineau, Gene 79:33-46, 1989]) is a nucleotide highly unique to Streptomyces species which, in combination with flanking regions, allowed the designation of a genus-specific probe. Regions 158 through 203 of the 16S rRNA and 1518 through 1645 of the 23S rRNA (helix 54 [Pernodet et al., Gene 79:33-46, 1989]) have a high potential to define species, whereas the degree of variation in regions 982 through 998 and 1102 through 1122 of the 16S rRNA is less pronounced but characteristic for at least certain species. Alone or in combination with each other, these regions may serve as target sites for synthetic oligonucleotide probes and primers to be used in the determination of pure cultures and in the characterization of community structures. The specificity of several probes is demonstrated by dot blot hybridization.     

Retrieval of nearly complete 16S rRNA gene sequences from environmental DNA following 16S rRNA-based community fingerprinting

16S rRNA-based fingerprinting techniques allow rapid analyses of overall bacterial community structure but suffer from a lack of phylogenetic information hitherto retrievable from the short 16S rRNA gene sequences obtained from excised bands. An approach is presented that allows nearly complete 16S rRNA gene sequences to be retrieved for abundant components of the bacterial community as obtained by the community fingerprint, i.e. those reflected by major fingerprint bands. This was achieved by designing a pair of highly specific primers derived from the sequence of an excised band. Combined with universal 16S rRNA primers, these specific primers were applied directly to environmental DNA serving as template. This procedure allowed the generation of a nearly complete 16S rRNA gene sequence of the target taxon by specific polymerase chain reaction (PCR) followed by cycle sequencing down to a relative abundance of at least 1.5% of the environmental DNA. The procedure was exemplified for an epsilonproteobacterium related to Thiomicrospira denitrificans occurring in the central Baltic Sea. This approach is based only on PCR without any cloning step involved. It allows focussing on specific target taxa and is thus rather efficient. This approach should be applicable in general to 16S rRNA or 16S rRNA gene-based fingerprinting techniques and their respective environmental DNA.     

Novel primers for 16S rRNA-based archaeal community analyses in environmental samples

Next generation sequencing technologies for in depth analyses of complex microbial communities rely on rational primer design based on up-to-date reference databases. Most of the 16S rRNA-gene based analyses of environmental Archaea community composition use PCR primers developed from small data sets several years ago, making an update long overdue. Here we present a new set of archaeal primers targeting the 16S rRNA gene designed from 8500 aligned archaeal sequences in the SILVA database. The primers 340F-1000R showed a high archaeal specificity (< 1% bacteria amplification) covering 93 and 97% of available sequences for Crenarchaeota and Euryarchaeota respectively. In silico tests of the primers revealed at least 38% higher coverage for Archaea compared to other commonly used primers. Empirical tests with clone libraries confirmed the high specificity of the primer pair to Archaea in three biomes: surface waters in the Arctic Ocean, the pelagic zone of a temperate lake and a methanogenic bioreactor. The clone libraries featured both Euryarchaeota and Crenarchaeota in variable proportions and revealed dramatic differences in the archaeal community composition and minimal phylogenetic overlap between samples.     

Designation of Streptomycete 16S and 23S rRNA-based target regions for oligonucleotide probes.

The 16S and 23S rRNA of various Streptomyces species were partially sequenced and screened for the presence of stretches that could define all members of the genus, groups of species, or individual species. Nucleotide 929 (Streptomyces ambofaciens nomenclature [J.L. Pernodet, M.T. Alegre, F. Boccard, and M. Guerineau, Gene 79:33-46, 1989]) is a nucleotide highly unique to Streptomyces species which, in combination with flanking regions, allowed the designation of a genus-specific probe. Regions 158 through 203 of the 16S rRNA and 1518 through 1645 of the 23S rRNA (helix 54 [Pernodet et al., Gene 79:33-46, 1989]) have a high potential to define species, whereas the degree of variation in regions 982 through 998 and 1102 through 1122 of the 16S rRNA is less pronounced but characteristic for at least certain species. Alone or in combination with each other, these regions may serve as target sites for synthetic oligonucleotide probes and primers to be used in the determination of pure cultures and in the characterization of community structures. The specificity of several probes is demonstrated by dot blot hybridization.     

Reducing the effects of PCR amplification and sequencing artifacts on 16S rRNA-based studies

The advent of next generation sequencing has coincided with a growth in interest in using these approaches to better understand the role of the structure and function of the microbial communities in human, animal, and environmental health. Yet, use of next generation sequencing to perform 16S rRNA gene sequence surveys has resulted in considerable controversy surrounding the effects of sequencing errors on downstream analyses. We analyzed 2.7×10(6) reads distributed among 90 identical mock community samples, which were collections of genomic DNA from 21 different species with known 16S rRNA gene sequences; we observed an average error rate of 0.0060. To improve this error rate, we evaluated numerous methods of identifying bad sequence reads, identifying regions within reads of poor quality, and correcting base calls and were able to reduce the overall error rate to 0.0002. Implementation of the PyroNoise algorithm provided the best combination of error rate, sequence length, and number of sequences. Perhaps more problematic than sequencing errors was the presence of chimeras generated during PCR. Because we knew the true sequences within the mock community and the chimeras they could form, we identified 8% of the raw sequence reads as chimeric. After quality filtering the raw sequences and using the Uchime chimera detection program, the overall chimera rate decreased to 1%. The chimeras that could not be detected were largely responsible for the identification of spurious operational taxonomic units (OTUs) and genus-level phylotypes. The number of spurious OTUs and phylotypes increased with sequencing effort indicating that comparison of communities should be made using an equal number of sequences. Finally, we applied our improved quality-filtering pipeline to several benchmarking studies and observed that even with our stringent data curation pipeline, biases in the data generation pipeline and batch effects were observed that could potentially confound the interpretation of microbial community data.     

Application and validation of DNA microarrays for the 16S rRNA-based analysis of marine bacterioplankton

An oligonucleotide probe-based DNA microarray was evaluated for its ability to detect 16S rRNA targets in marine bacterioplankton samples without prior amplification by polymerase chain reaction (PCR). The results obtained were compared with those of quantitative fluorescence in situ hybridization (FISH). For extraction and direct labelling of total RNA, a fast and efficient protocol based on commercially available kits was established. A set of redundant and hierarchically structured probes was applied, and specificity of hybridization was assessed by additional control oligonucleotides comprising single central mismatches. The protocol was initially tested by microarray analysis of bacterial pure cultures. Complete discrimination of all control oligonucleotides was achieved, indicating a high degree of hybridization specificity. In a co-culture, abundant members were detected by microarray analysis, but signal ratios of positive probes did not correlate well with quantitative data from FISH experiments. A marine picoplankton sample from the German Bight was analysed. Bacterial populations with relative abundances of at least 5% were detected by hybridizing 0.1 microg of total RNA extracted from a sample of 375 ml equivalent to 4.1 x 10(8) cells. Our results demonstrate that major populations of marine bacterioplankton can be identified by microarray analysis in a fast and reliable way, even in relatively low volumes of sea water.     

中国传统发酵食品微生物组研究进展

传统发酵食品风味独特、营养丰富,多采用自然接种方式进行生产,部分类型的生产技艺已具有数千年的历史。近年来应用新技术手段对传统发酵食品发酵过程进行的研究表明,传统发酵食品微生物组具有丰富的科学内涵和重要的应用价值。本文就我国传统发酵食品微生物组的基本特征、研究进展和发展方向进行了简要的评述。     

Development and validation of a prototype 16S rRNA-based taxonomic microarray for Alphaproteobacteria

The microarray approach has been proposed for high throughput analysis of the microbial community by providing snapshots of the microbial diversity under different environmental conditions. For this purpose, a prototype of a 16S rRNA-based taxonomic microarray was developed and evaluated for assessing bacterial community diversity. The prototype microarray is composed of 122 probes that target bacteria at various taxonomic levels from phyla to species (mostly Alphaproteobacteria). The prototype microarray was first validated using bacteria in pure culture. Differences in the sequences of probes and potential target DNAs were quantified as weighted mismatches (WMM) in order to evaluate hybridization reliability. As a general feature, probes having a WMM > 2 with target DNA displayed only 2.8% false positives. The prototype microarray was subsequently tested with an environmental sample, which consisted of an Agrobacterium-related polymerase chain reaction amplicon from a maize rhizosphere bacterial community. Microarray results were compared to results obtained by cloning-sequencing with the same DNA. Microarray analysis enabled the detection of all 16S rRNA gene sequences found by cloning-sequencing. Sequences representing only 1.7% of the clone library were detected. In conclusion, this prototype 16S rRNA-based taxonomic microarray appears to be a promising tool for the analysis of Alphaproteobacteria in complex ecosystems.     

Power law analysis of the human microbiome

Taylor's (1961, Nature, 189:732) power law, a power function (V = am^b) describing the scaling relationship between the mean and variance of population abundances of organisms, has been found to govern the population abundance distributions of single species in both space and time in macroecology. It is regarded as one of few generalities in ecology, and its parameter b has been widely applied to characterize spatial aggregation (i.e. heterogeneity) and temporal stability of single‐species populations. Here, we test its applicability to bacterial populations in the human microbiome using extensive data sets generated by the US‐NIH Human Microbiome Project (HMP). We further propose extending Taylor's power law from the population to the community level, and accordingly introduce four types of power‐law extensions (PLEs): type I PLE for community spatial aggregation (heterogeneity), type II PLE for community temporal aggregation (stability), type III PLE for mixed‐species population spatial aggregation (heterogeneity) and type IV PLE for mixed‐species population temporal aggregation (stability). Our results show that fittings to the four PLEs with HMP data were statistically extremely significant and their parameters are ecologically sound, hence confirming the validity of the power law at both the population and community levels. These findings not only provide a powerful tool to characterize the aggregations of population and community in both time and space, offering important insights into community heterogeneity in space and/or stability in time, but also underscore the three general properties of power laws (scale invariance, no average and universality) and their specific manifestations in our four PLEs.     

CORE: a phylogenetically-curated 16S rDNA database of the core oral microbiome

Comparing bacterial 16S rDNA sequences to GenBank and other large public databases via BLAST often provides results of little use for identification and taxonomic assignment of the organisms of interest. The human microbiome, and in particular the oral microbiome, includes many taxa, and accurate identification of sequence data is essential for studies of these communities. For this purpose, a phylogenetically curated 16S rDNA database of the core oral microbiome, CORE, was developed. The goal was to include a comprehensive and minimally redundant representation of the bacteria that regularly reside in the human oral cavity with computationally robust classification at the level of species and genus. Clades of cultivated and uncultivated taxa were formed based on sequence analyses using multiple criteria, including maximum-likelihood-based topology and bootstrap support, genetic distance, and previous naming. A number of classification inconsistencies for previously named species, especially at the level of genus, were resolved. The performance of the CORE database for identifying clinical sequences was compared to that of three publicly available databases, GenBank nr/nt, RDP and HOMD, using a set of sequencing reads that had not been used in creation of the database. CORE offered improved performance compared to other public databases for identification of human oral bacterial 16S sequences by a number of criteria. In addition, the CORE database and phylogenetic tree provide a framework for measures of community divergence, and the focused size of the database offers advantages of efficiency for BLAST searching of large datasets. The CORE database is available as a searchable interface and for download at http://microbiome.osu.edu.     

Environmental microbiome engineering for the mitigation of climate change

Environmental microbiome engineering is emerging as a potential avenue for climate change mitigation. In this process, microbial inocula are introduced to natural microbial communities to tune activities that regulate the long-term stabilization of carbon in ecosystems. In this review, we outline the process of environmental engineering and synthesize key considerations about ecosystem functions to target, means of sourcing microorganisms, strategies for designing microbial inocula, methods to deliver inocula, and the factors that enable inocula to establish within a resident community and modify an ecosystem function target. Recent work, enabled by high-throughput technologies and modeling approaches, indicate that microbial inocula designed from the top-down, particularly through directed evolution, may generally have a higher chance of establishing within existing microbial communities than other historical approaches to microbiome engineering. We address outstanding questions about the determinants of inocula establishment and provide suggestions for further research about the possibilities and challenges of environmental microbiome engineering as a tool to combat climate change.     

Anatomy promotes neutral coexistence of strains in the human skin microbiome

1. Download : Download high-res image (201KB)
2. Download : Download full-size image