Amplification by PCR Artificially Reduces the Proportion of the Rare Biosphere in Microbial Communities

The microbial world has been shown to hold an unimaginable diversity. The use of rRNA genes and PCR amplification to assess microbial community structure and diversity present biases that need to be analyzed in order to understand the risks involved in those estimates. Herein, we show that PCR amplification of specific sequence targets within a community depends on the fractions that those sequences represent to the total DNA template. Using quantitative, real-time, multiplex PCR and specific Taqman probes, the amplification of 16S rRNA genes from four bacterial species within a laboratory community were monitored. Results indicate that the relative amplification efficiency for each bacterial species is a nonlinear function of the fraction that each of those taxa represent within a community or multispecies DNA template. Consequently, the low-proportion taxa in a community are under-represented during PCR-based surveys and a large number of sequences might need to be processed to detect some of the bacterial taxa within the ‘rare biosphere’. The structure of microbial communities from PCR-based surveys is clearly biased against low abundant taxa which are required to decipher the complete extent of microbial diversity in nature

Phylogeny-Aware Analysis of Metagenome Community Ecology Based on Matched Reference Genomes while Bypassing Taxonomy

We introduce the operational genomic unit (OGU) method, a metagenome analysis strategy that directly exploits sequence alignment hits to individual reference genomes as the minimum unit for assessing the diversity of microbial communities and their relevance to environmental factors. This approach is independent of taxonomic classification, granting the possibility of maximal resolution of community composition, and organizes features into an accurate hierarchy using a phylogenomic tree. The outputs are suitable for contemporary analytical protocols for community ecology, differential abundance, and supervised learning while supporting phylogenetic methods, such as UniFrac and phylofactorization, that are seldom applied to shotgun metagenomics despite being prevalent in 16S rRNA gene amplicon studies. As demonstrated in two real-world case studies, the OGU method produces biologically meaningful patterns from microbiome data sets. Such patterns further remain detectable at very low metagenomic sequencing depths. Compared with taxonomic unit-based analyses implemented in currently adopted metagenomics tools, and the analysis of 16S rRNA gene amplicon sequence variants, this method shows superiority in informing biologically relevant insights, including stronger correlation with body environment and host sex on the Human Microbiome Project data set and more accurate prediction of human age by the gut microbiomes of Finnish individuals included in the FINRISK 2002 cohort. We provide Woltka, a bioinformatics tool to implement this method, with full integration with the QIIME 2 package and the Qiita web platform, to facilitate adoption of the OGU method in future metagenomics studies. IMPORTANCE Shotgun metagenomics is a powerful, yet computationally challenging, technique compared to 16S rRNA gene amplicon sequencing for decoding the composition and structure of microbial communities. Current analyses of metagenomic data are primarily based on taxonomic classification, which is limited in feature resolution. To solve these challenges, we introduce operational genomic units (OGUs), which are the individual reference genomes derived from sequence alignment results, without further assigning them taxonomy. The OGU method advances current read-based metagenomics in two dimensions: (i) providing maximal resolution of community composition and (ii) permitting use of phylogeny-aware tools. Our analysis of real-world data sets shows that it is advantageous over currently adopted metagenomic analysis methods and the finest-grained 16S rRNA analysis methods in predicting biological traits. We thus propose the adoption of OGUs as an effective practice in metagenomic studies.Peer reviewe

Analyses of the Microbial Diversity across the Human Microbiome

Analysis of human body microbial diversity is fundamental to understanding community structure, biology and ecology. The National Institutes of Health Human Microbiome Project (HMP) has provided an unprecedented opportunity to examine microbial diversity within and across body habitats and individuals through pyrosequencing-based profiling of 16 S rRNA gene sequences (16 S) from habits of the oral, skin, distal gut, and vaginal body regions from over 200 healthy individuals enabling the application of statistical techniques. In this study, two approaches were applied to elucidate the nature and extent of human microbiome diversity. First, bootstrap and parametric curve fitting techniques were evaluated to estimate the maximum number of unique taxa, Smax, and taxa discovery rate for habitats across individuals. Next, our results demonstrated that the variation of diversity within low abundant taxa across habitats and individuals was not sufficiently quantified with standard ecological diversity indices. This impact from low abundant taxa motivated us to introduce a novel rank-based diversity measure, the Tail statistic, (“τ”), based on the standard deviation of the rank abundance curve if made symmetric by reflection around the most abundant taxon. Due to τ’s greater sensitivity to low abundant taxa, its application to diversity estimation of taxonomic units using taxonomic dependent and independent methods revealed a greater range of values recovered between individuals versus body habitats, and different patterns of diversity within habitats. The greatest range of τ values within and across individuals was found in stool, which also exhibited the most undiscovered taxa. Oral and skin habitats revealed variable diversity patterns, while vaginal habitats were consistently the least diverse. Collectively, these results demonstrate the importance, and motivate the introduction, of several visualization and analysis methods tuned specifically for next-generation sequence data, further revealing that low abundant taxa serve as an important reservoir of genetic diversity in the human microbiome

Gene expression of bacterial collagenolytic proteases in root caries

Objective: It is unknown whether bacteria play a role in the collagen matrix degradation that occurs during caries progression. Our aim was to characterize the expression level of genes involved in bacterial collagenolytic proteases in root biofilms with and without caries. Method: we collected samples from active cavitated root caries lesions (RC, n = 30) and from sound root surfaces (SRS, n = 10). Total microbial RNA was isolated and cDNA sequenced on the Illumina Hi-Seq2500. Reads were mapped to 162 oral bacterial reference genomes. Genes encoding putative bacterial collagenolytic proteases were identified. Normalization and differential expression analysis was performed on all metatranscriptomes (FDR8) but none in SRS were Pseudoramibacter alactolyticus [HMPREF0721_RS02020; HMPREF0721_RS04640], Scardovia inopinata [SCIP_RS02440] and Olsenella uli DSM7084 [OLSU_RS02990]. Conclusion: Our findings suggest that the U32 proteases may be related to carious dentine. The contribution of a small number of species to dentine degradation should be further investigated. These proteases may have potential in future biotechnological and medical applications, serving as targets for the development of therapeutic agents

Differential preservation of endogenous human and microbial DNA in dental calculus and dentin

Dental calculus (calcified dental plaque) is prevalent in archaeological skeletal collections and is a rich source of oral microbiome and host-derived ancient biomolecules. Recently, it has been proposed that dental calculus may provide a more robust environment for DNA preservation than other skeletal remains, but this has not been systematically tested. In this study, shotgun-sequenced data from paired dental calculus and dentin samples from 48 globally distributed individuals are compared using a metagenomic approach. Overall, we find DNA from dental calculus is consistently more abundant and less contaminated than DNA from dentin. The majority of DNA in dental calculus is microbial and originates from the oral microbiome; however, a small but consistent proportion of DNA (mean 0.08 ± 0.08%, range 0.007–0.47%) derives from the host genome. Host DNA content within dentin is variable (mean 13.70 ± 18.62%, range 0.003–70.14%), and for a subset of dentin samples (15.21%), oral bacteria contribute \u3e 20% of total DNA. Human DNA in dental calculus is highly fragmented, and is consistently shorter than both microbial DNA in dental calculus and human DNA in paired dentin samples. Finally, we find that microbial DNA fragmentation patterns are associated with guanine-cytosine (GC) content, but not aspects of cellular structure

A randomised clinical study to determine the effect of a toothpaste containing enzymes and proteins on plaque oral microbiome ecology

The numerous species that make up the oral microbiome are now understood to play a key role in establishment and maintenance of oral health. The ability to taxonomically identify community members at the species level is important to elucidating its diversity and association to health and disease. We report the overall ecological effects of using a toothpaste containing enzymes and proteins compared to a control toothpaste on the plaque microbiome. The results reported here demonstrate that a toothpaste containing enzymes and proteins can augment natural salivary defences to promote an overall community shift resulting in an increase in bacteria associated with gum health and a concomitant decrease in those associated with periodontal disease. Statistical analysis shows significant increases in 12 taxa associated with gum health including Neisseria spp. and a significant decrease in 10 taxa associated with periodontal disease including Treponema spp. The results demonstrate that a toothpaste containing enzymes and proteins can significantly shift the ecology of the oral microbiome (at species level) resulting in a community with a stronger association to health

Bacterial community development in experimental gingivitis.

Current knowledge of the microbial composition of dental plaque in early gingivitis is based largely on microscopy and cultural methods, which do not provide a comprehensive description of oral microbial communities. This study used 454-pyrosequencing of the V1-V3 region of 16S rRNA genes (approximately 500 bp), and bacterial culture, to characterize the composition of plaque during the transition from periodontal health to gingivitis. A total of 20 healthy volunteers abstained from oral hygiene for two weeks, allowing plaque to accumulate and gingivitis to develop. Plaque samples were analyzed at baseline, and after one and two weeks. In addition, plaque samples from 20 chronic periodontitis patients were analyzed for cross-sectional comparison to the experimental gingivitis cohort. All of the healthy volunteers developed gingivitis after two weeks. Pyrosequencing yielded a final total of 344,267 sequences after filtering, with a mean length of 354 bases, that were clustered into an average of 299 species-level Operational Taxonomic Units (OTUs) per sample. Principal coordinates analysis (PCoA) plots revealed significant shifts in the bacterial community structure of plaque as gingivitis was induced, and community diversity increased significantly after two weeks. Changes in the relative abundance of OTUs during the transition from health to gingivitis were correlated to bleeding on probing (BoP) scores and resulted in the identification of new health- and gingivitis-associated taxa. Comparison of the healthy volunteers to the periodontitis patients also confirmed the association of a number of putative periodontal pathogens with chronic periodontitis. Taxa associated with gingivitis included Fusobacterium nucleatum subsp. polymorphum, Lachnospiraceae [G-2] sp. HOT100, Lautropia sp. HOTA94, and Prevotella oulorum, whilst Rothia dentocariosa was associated with periodontal health. Further study of these taxa is warranted and may lead to new therapeutic approaches to prevent periodontal disease.BBSRC Industrial Case Studentship ref no. BB/G01714X/1 in collaboration with GlaxoSmithKline

Phylogeny-Aware Analysis of Metagenome Community Ecology Based on Matched Reference Genomes while Bypassing Taxonomy

We introduce the operational genomic unit (OGU) method, a metagenome analysis strategy that directly exploits sequence alignment hits to individual reference genomes as the minimum unit for assessing the diversity of microbial communities and their relevance to environmental factors. This approach is independent of taxonomic classification, granting the possibility of maximal resolution of community composition, and organizes features into an accurate hierarchy using a phylogenomic tree. The outputs are suitable for contemporary analytical protocols for community ecology, differential abundance, and supervised learning while supporting phylogenetic methods, such as UniFrac and phylofactorization, that are seldom applied to shotgun metagenomics despite being prevalent in 16S rRNA gene amplicon studies. As demonstrated in two real-world case studies, the OGU method produces biologically meaningful patterns from microbiome data sets. Such patterns further remain detectable at very low metagenomic sequencing depths. Compared with taxonomic unit-based analyses implemented in currently adopted metagenomics tools, and the analysis of 16S rRNA gene amplicon sequence variants, this method shows superiority in informing biologically relevant insights, including stronger correlation with body environment and host sex on the Human Microbiome Project data set and more accurate prediction of human age by the gut microbiomes of Finnish individuals included in the FINRISK 2002 cohort. We provide Woltka, a bioinformatics tool to implement this method, with full integration with the QIIME 2 package and the Qiita web platform, to facilitate adoption of the OGU method in future metagenomics studies.IMPORTANCE Shotgun metagenomics is a powerful, yet computationally challenging, technique compared to 16S rRNA gene amplicon sequencing for decoding the composition and structure of microbial communities. Current analyses of metagenomic data are primarily based on taxonomic classification, which is limited in feature resolution. To solve these challenges, we introduce operational genomic units (OGUs), which are the individual reference genomes derived from sequence alignment results, without further assigning them taxonomy. The OGU method advances current read-based metagenomics in two dimensions: (i) providing maximal resolution of community composition and (ii) permitting use of phylogeny-aware tools. Our analysis of real-world data sets shows that it is advantageous over currently adopted metagenomic analysis methods and the finest-grained 16S rRNA analysis methods in predicting biological traits. We thus propose the adoption of OGUs as an effective practice in metagenomic studies.</p