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

Abstract

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.

Intrahepatic T‐Cell Receptor β Immune Repertoire Is Essential for Liver Regeneration

“…V, J, and CDR3 regions of TCRβ consensus sequences were identified using BLAST Plus in the international ImMunoGeneTics (IMGT) information system (IMGT: http://www.imgt.org/) by a standard algorithm . The diversity of TCRβ was measured by Gini coefficient, normalized Shannon diversity entropy, and Rank‐abundance, which have been widely used for assessing the richness and diversity of TCR as described …”

Ecotoxicological assessment of antibiotics: A call for improved consideration of microorganisms

“…The recently standardized DNA extraction protocol (ISO 11063) is an important step forward, but this protocol will not be the most appropriate protocol in all cases (Plassart et al, 2012;Terrat et al, 2014). In principle, it should also be possible to standardize downstream molecular analysis of the extracted DNA involving steps such as PCR, sequencing, and bioinformatic analysis, but currently the proper selection of specific methods and protocols to standardize is not obvious, as the respective techniques and instrument platforms are rapidly developing (Birtel et al, 2015). Molecular community analysis requires sophisticated statistical analysis to interpret the data.…”

“…Molecular community analysis requires sophisticated statistical analysis to interpret the data. A wide range of statistical approaches are currently available (Haegeman et al, 2013;Ramette, 2007), but standardized protocols for data analysis have yet to be developed (Birtel et al, 2015).…”

Maternal microbiomes in preterm birth: Recent progress and analytical pipelines

“…Each method has its own formula and calculation that must be considered when deciding which one to use for analysis. For example, Bray Curtis accounts for both the abundance of microbiota 54 . In contrast, unweighted UniFrac uses only the phylogenetic distance between the different taxa 55 .…”