Optimal sequence similarity thresholds for clustering of molecular operational taxonomic units in DNA metabarcoding studies

Dynamics and drivers of mycorrhizal fungi after glacier retreat

The development of terrestrial ecosystems depends greatly on plant mutualists such as mycorrhizal fungi. The global retreat of glaciers exposes nutrient-poor substrates in extreme environments and provides a unique opportunity to study early successions of mycorrhizal fungi by assessing their dynamics and drivers. We combined environmental DNA metabarcoding and measurements of local conditions to assess the succession of mycorrhizal communities during soil development in 46 glacier forelands around the globe, testing whether dynamics and drivers differ between mycorrhizal types. Mycorrhizal fungi colonized deglaciated areas very quickly (< 10 yr), with arbuscular mycorrhizal fungi tending to become more diverse through time compared to ectomycorrhizal fungi. Both alpha- and beta-diversity of arbuscular mycorrhizal fungi were significantly related to time since glacier retreat and plant communities, while microclimate and primary productivity were more important for ectomycorrhizal fungi. The richness and composition of mycorrhizal communities were also significantly explained by soil chemistry, highlighting the importance of microhabitat for community dynamics. The acceleration of ice melt and the modifications of microclimate forecasted by climate change scenarios are expected to impact the diversity of mycorrhizal partners. These changes could alter the interactions underlying biotic colonization and belowground-aboveground linkages, with multifaceted impacts on soil development and associated ecological processes.

The importance of species addition 'versus' replacement varies over succession in plant communities after glacier retreat

The mechanisms underlying plant succession remain highly debated. Due to the local scope of most studies, we lack a global quantification of the relative importance of species addition 'versus' replacement. We assessed the role of these processes in the variation (β-diversity) of plant communities colonizing the forelands of 46 retreating glaciers worldwide, using both environmental DNA and traditional surveys. Our findings indicate that addition and replacement concur in determining community changes in deglaciated sites, but their relative importance varied over time. Taxa addition dominated immediately after glacier retreat, as expected in harsh environments, while replacement became more important for late-successional communities. These changes were aligned with total β-diversity changes, which were more pronounced between early-successional communities than between late-successional communities (>50 yr since glacier retreat). Despite the complexity of community assembly during plant succession, the observed global pattern suggests a generalized shift from the dominance of facilitation and/or stochastic processes in early-successional communities to a predominance of competition later on.

Environmental DNA and visual encounter surveys for amphibian biomonitoring in aquatic environments of the Ecuadorian Amazon

Background

The development of anthropogenic activities has generated a decline in aquatic fauna populations, and amphibians have been the most affected. The decline of batrachofauna is concerning, as 41% of all species worldwide are endangered. For this reason, rapid, efficient, and non-invasive biodiversity monitoring techniques are needed, and environmental DNA (eDNA) is one such tool that has been sparsely applied in Ecuador. This technique has allowed scientists generates information on species diversity and amphibian community composition from a water sample. This study applied eDNA-based biomonitoring analyses and visual encounter surveys (VES) as inventory techniques to identify the diversity of aquatic amphibians in the Tena River micro-basin (TRMB).

Methods

The experimental design was divided into three components: (1) fieldwork: all amphibians were recorded by the VES technique and water samples were collected; (2) laboratory work: DNA isolation from amphibian tissue samples and eDNA-containing filters, amplification, electrophoresis, and sequencing were performed; (3) Data analysis: a local DNA reference database was constructed, and eDNA sequence data were processed for classification, taxonomic assignment, and ecological interpretation.

Results

Using both eDNA and VES, we detected 33 amphibian species (13 with eDNA only, five with VES only, and 15 with both methods). These species belonged to six amphibian families: Hylidae being the richest with 14 species (three eDNA, one VES, and 10 with both methods), followed by Strabomantidae with nine species (six eDNA, one VES, and two with both methods). All families were detected with both methods, except for the Aromobatidae, having one single record (Allobates aff. insperatus) by VES. Individually, eDNA detected 28 species and had a detection probability (DP) of 0.42 CI [0.40-0.45], while VES recorded 20 species with a DP of 0.17 CI [0.14-0.20]. Similarly, using VES, Cochranella resplendens was detected for the first time in TRMB, while with eDNA, four mountain frogs Pristimantis acerus, Pristimantis eriphus, Pristimantis mallii, and Pristimantis sp. (INABIO 15591) previously recorded at 1,518 m.a.s.l. at altitudes below 600 m.a.s.l. were detected.

Conclusions

Results obtained in this study showed that eDNA-based detection had a greater capacity to detect amphibians in aquatic environments compared to VES. The combination of VES and eDNA improves the sensitivity of species detection and provides more reliable, robust, and detailed information. The latter is essential for developing conservation strategies in the Ecuadorian Amazon.

Microbial Diversity Biased Estimation Caused by Intragenomic Heterogeneity and Interspecific Conservation of 16S rRNA Genes

The 16S rRNA gene has been extensively used as a molecular marker to explore evolutionary relationships and profile microbial composition throughout various environments. Despite its convenience and prevalence, limitations are inevitable. Variable copy numbers, intragenomic heterogeneity, and low taxonomic resolution have caused biases in estimating microbial diversity. Here, analysis of 24,248 complete prokaryotic genomes indicated that the 16S rRNA gene copy number ranged from 1 to 37 in bacteria and 1 to 5 in archaea, and intragenomic heterogeneity was observed in 60% of prokaryotic genomes, most of which were below 1%. The overestimation of microbial diversity caused by intragenomic variation and the underestimation introduced by interspecific conservation were calculated when using full-length or partial 16S rRNA genes. Results showed that, at the 100% threshold, microbial diversity could be overestimated by as much as 156.5% when using the full-length gene. The V4 to V5 region-based analyses introduced the lowest overestimation rate (4.4%) but exhibited slightly lower species resolution than other variable regions under the 97% threshold. For different variable regions, appropriate thresholds rather than the canonical value 97% were proposed for minimizing the risk of splitting a single genome into multiple clusters and lumping together different species into the same cluster. This study has not only updated the 16S rRNA gene copy number and intragenomic variation information for the currently available prokaryotic genomes, but also elucidated the biases in estimating prokaryotic diversity with quantitative data, providing references for choosing amplified regions and clustering thresholds in microbial community surveys. IMPORTANCE Microbial diversity is typically analyzed using marker gene-based methods, of which 16S rRNA gene sequencing is the most widely used approach. However, obtaining an accurate estimation of microbial diversity remains a challenge, due to the intragenomic variation and low taxonomic resolution of 16S rRNA genes. Comprehensive examination of the bias in estimating such prokaryotic diversity using 16S rRNA genes within ever-increasing prokaryotic genomes highlights the importance of the choice of sequencing regions and clustering thresholds based on the specific research objectives.