Body size determines soil community assembly in a tropical forest

Rainfall seasonality shapes microbial assembly and niche characteristics in Yunnan Plateau lakes, China

Microorganisms are crucial components of freshwater ecosystems. Understanding the microbial community assembly processes and niche characteristics in freshwater ecosystems, which are poorly understood, is crucial for evaluating microbial ecological roles. The Yunnan Plateau lakes in China represent a freshwater ecosystem that is experiencing eutrophication due to anthropogenic activities. Here, variation in the assembly and niche characteristics of both prokaryotic and microeukaryotic communities was explored in Yunnan Plateau lakes across two seasons (dry season and rainy season) to determine the impacts of rainfall and environmental conditions on the microbial community and niche. The results showed that the environmental heterogeneity of the lakes decreased in the rainy season compared to the dry season. The microbial (bacterial and microeukaryotic) α-diversity significantly decreased during the rainy season. Deterministic processes were found to dominate microbial community assembly in both seasons. β-Diversity decomposition analysis revealed that microbial community compositional dissimilarities were dominated by species replacement processes. The co-occurrence networks indicated reduced species complexity for microbes and a destabilized network for prokaryotes prior to rainfall, while the opposite was found for microeukaryotes following rainfall. Microbial niche breadth decreased significantly in the rainy season. In addition, lower prokaryotic niche overlap, but greater microeukaryotic niche overlap, was observed after rainfall. Rainfall and environmental conditions significantly affected the microbial community assembly and niche characteristics. It can be concluded that rainfall and external pollutant input during the seasonal transition alter the lake environment, thereby regulating the microbial community and niche in these lakes. Our findings offer new insight into microbiota assembly and niche patterns in plateau lakes, further deepening the understanding of freshwater ecosystem functioning.

Are plant traits drivers of endophytic communities in seasonally flooded tropical forests?

PREMISE

In the Amazon basin, seasonally flooded (SF) forests offer varying water constraints, providing an excellent way to investigate the role of habitat selection on microbial communities within plants. However, variations in the microbial community among host plants cannot solely be attributed to environmental factors, and how plant traits contribute to microbial assemblages remains an open question.

METHODS

We described leaf- and root-associated microbial communities using ITS2 and 16 S high-throughput sequencing and investigated the stochastic-deterministic balance shaping these community assemblies using two null models. Plant ecophysiological functioning was evaluated by focusing on 10 leaf and root traits in 72 seedlings, belonging to seven tropical SF tree species in French Guiana. We then analyzed how root and leaf traits drove the assembly of endophytic communities.

RESULTS

While both stochastic and deterministic processes governed the endophyte assembly in the leaves and roots, stochasticity prevailed. Discrepancies were found between fungi and bacteria, highlighting that these microorganisms have distinct ecological strategies within plants. Traits, especially leaf traits, host species and spatial predictors better explained diversity than composition, but they were modest predictors overall.

CONCLUSIONS

This study widens our knowledge about tree species in SF forests, a habitat sensitive to climate change, through the combined analyses of their associated microbial communities with functional traits. We emphasize the need to investigate other plant traits to better disentangle the drivers of the relationship between seedlings and their associated microbiomes, ultimately enhancing their adaptive capacities to climate change.

The rhizosphere microbiome of 51 potato cultivars with diverse plant growth characteristics

Rhizosphere microbial communities play a substantial role in plant productivity. We studied the rhizosphere bacteria and fungi of 51 distinct potato cultivars grown under similar greenhouse conditions using a metabarcoding approach. As expected, individual cultivars were the most important determining factor of the rhizosphere microbial composition; however, differences were also obtained when grouping cultivars according to their growth characteristics. We showed that plant growth characteristics were related to deterministic and stochastic assembly processes of bacterial and fungal communities, respectively. The bacterial genera Arthrobacter and Massilia (known to produce indole acetic acid and siderophores) exhibited greater relative abundance in high- and medium-performing cultivars. Bacterial co-occurrence networks were larger in the rhizosphere of these cultivars and were characterized by a distinctive combination of plant beneficial Proteobacteria and Actinobacteria along with a module of diazotrophs namely Azospira, Azoarcus, and Azohydromonas. Conversely, the network within low-performing cultivars revealed the lowest nodes, hub taxa, edges density, robustness, and the highest average path length resulting in reduced microbial associations, which may potentially limit their effectiveness in promoting plant growth. Our findings established a clear pattern between plant productivity and the rhizosphere microbiome composition and structure for the investigated potato cultivars, offering insights for future management practices.

Ecological restoration enhanced the stability of epiphytic microbial food webs of submerged macrophytes: Insights from predation characteristics of epiphytic predators

The application of various submerged macrophytes for ecological restoration has gained increasing attention in urban lake ecosystems. The multitrophic microbial communities that colonized in various submerged macrophytes constitute microbial food webs through trophic cascade effects, which affect the biogeochemical cycles of the lake ecosystem and directly determine the effects of ecological restoration. Therefore, it is essential to reveal the diversity, composition, assembly processes, and stability of the microbial communities within epiphytic food webs of diverse submerged macrophytes under eutrophication and ecological restoration scenarios. In this study, we explored the epiphytic microbial food webs of Vallisneria natans and Hydrilla verticillata in both eutrophic and ecological restoration regions. The obtained results indicated that the two regions with different nutrient levels remarkably affected the diversity and composition of epiphytic multitrophic microbial communities of submerged macrophytes, among them, the community composition of epiphytic predators were more prone to change. Secondly, environmental filtering effects played a more important role in driving the community assembly of epiphytic predators than that of prey. Furthermore, the generality and intraguild predation of epiphytic predators were significantly improved within ecological restoration regions, which increased the stability of epiphytic microbial food webs. Additionally, compared with Hydrilla verticillata, the epiphytic microbial food webs of Vallisneria natans exhibited higher multitrophic diversity and higher network stability regardless of regions. Overall, this study focused on the role of the epiphytic microbial food webs of submerged macrophytes in ecological restoration and uncovered the potential of epiphytic predators to enhance the stability of microbial food webs, which may provide new insights into the development of ecological restoration strategies.

Stable body sizes in soil nematodes across altitudes: The role of intrageneric variation in community assembly

Animal body size exhibits rapid responses to environmental variations and displays considerable variability across ecological scales, significantly influencing ecological community assembly. However, our understanding of the extent of body size variation and its responses to environmental differences within soil fauna remains limited, impeding a comprehensive grasp of soil fauna's functional ecology. Here, we aim to investigate the magnitude of intrageneric body size variation and its implications for soil nematode community assembly along an altitudinal gradient. We examined soil nematode body size responses along an altitudinal gradient spanning from 3136 to 4128 m in an alpine mountain region of the eastern Tibetan Plateau. We assessed the contributions of intra- and intergeneric variations in body size, both within and among communities, using individual body size values. The implications of these variations for community assembly processes were determined through phenotypic variance ratios employing permutation tests. Our analyses did not reveal statistically significant correlations between altitude and the community-weighted mean body mass, regardless of considering intrageneric trait variation (IGTV). Approximately 15% of the variation in body size among communities and a substantial 72% of the variation in body size within communities can be attributed to IGTV. Altitude did not significantly affect IGTV within or among communities. Furthermore, our results underscored the dominant role of internal filtering within the community in governing nematode community assembly, with external filtering outside the community playing a limited role within our altitudinal range. Our findings emphasize the dominant role of body size variation within communities rather than among communities, attributable to strong internal filtering processes. These findings advance our understanding of body size variation in soil nematodes across ecological scales and highlight the pivotal role of intrageneric variation in shaping the functional ecology of soil fauna.

The joint effects of local, climatic, and spatial variables determine soil oribatid mite community assembly along a temperate forest elevational gradient

Numerous factors influence mountain biodiversity variation across elevational gradients and recognizing the relative importance is vital for understanding species distribution mechanisms. We examined oribatid mites at nine elevations (from 600 to 2200 m a.s.l) and four vegetation types from mixed coniferous and broad-leaved forests to alpine tundra on Changbai Mountain. We assessed the contribution of environmental factors (climatic and local factors) and spatial processes (geographic or elevation distances) to oribatid mite community assembly and identified 59 oribatid mite species from 38 families and 51 genera. With increasing elevation, species richness and the Shannon index declined significantly, whereas abundance followed a hump-shaped trend. Soil TP, NH4 +-N, MAT, MAP, and elevation were the critical variables shaping oribatid mite communities based on random forest analysis. Moreover, environmental and spatial factors, and oribatid mite communities were significantly correlated based on Mantel and partial Mantel tests. Local characteristics (3.9%), climatic factors (1.9%), and spatial filtering (8.8%) played crucial roles in determining oribatid mite communities across nine elevational bands (based on variation partitioning analyses of abundance data). Within the same vegetation types, spatial processes had relatively little effects, with local characteristics the dominant drivers of oribatid mite community variation. Environmental and spatial filters together shape oribatid mite community assembly and their relative roles varied with elevation and vegetation type. These findings are crucial for the conservation, restoration, and management of Changbai mountain ecosystems in the context of climate change, along with the prediction of future vertical biotic gradient pattern evolution.

Historical and dispersal processes drive community assembly of multiple aquatic taxa in glacierized catchments in the Qinghai-Tibet plateau

Understanding the relative role of dispersal dynamics and niche constraints is not only a core task in community ecology, but also becomes an important prerequisite for bioassessment. Despite the recent progress in our knowledge of community assembly in space and time, patterns and processes underlying biotic communities in alpine glacierized catchments remain mostly ignored. To fill this knowledge gap, we combined the recently proposed dispersal-niche continuum index (DNCI) with traditional constrained ordinations and idealized patterns of species distributions to unravel community assembly mechanisms of different key groups of primary producers and consumers (i.e., phytoplankton, epiphytic algae, zooplankton, macroinvertebrates, and fishes) in rivers in the Qinghai-Tibet Plateau, the World's Third Pole. We tested whether organismal groups with contrasting body sizes differed in their assembly processes, and discussed their applicability in bioassessment in alpine zones. We found that community structure of alpine river biotas was always predominantly explained in terms of dispersal dynamics and historical biogeography. These patterns are most likely the result of differences in species-specific functional attributes, the stochastic colonization-extinction dynamics driven by multi-year glacier disturbances and the repeated hydrodynamic separation among alpine catchments after the rising of the Qilian mountains. Additionally, we found that the strength of dispersal dynamics and niche constraints was partially mediated by organismal body sizes, with dispersal processes being more influential for microscopic primary producers. Finding that zooplankton and macroinvertebrate communities followed clumped species replacement structures (i.e., Clementsian gradients) supports the notion that environmental filtering also contributes to the structure of high-altitude animal communities in glacierized catchments. In terms of the applied fields, we argue that freshwater bioassessment in glacierized catchments can benefit from incorporating the metacommunity perspective and applying novel approaches to (i) detect the optimal spatial scale for species sorting and (ii) identify and eliminate the species that are sensitive to dispersal-related processes.

eDNA Metabarcoding Reveals the Species-Area Relationship of Amphibians on the Zhoushan Archipelago

The species-area relationship is important for understanding species diversity patterns at spatial scales, but few studies have examined the relationship using environmental DNA (eDNA) techniques. We investigated amphibian diversity on 21 islands of the Zhoushan Archipelago and nearby mainland areas in China using the combination of eDNA metabarcoding and the traditional line transect method (TLTM) and identified the species-area relationship for amphibians on the islands. The mean detection probability of eDNA is 0.54, while the mean detection probability of TLTM is 0.24. The eDNA metabarcoding detected eight amphibian species on the islands and nine species in the mainland areas, compared with seven species on the islands and nine species in the mainland areas that were identified by TLTM. Amphibian richness on the islands increased with island area and habitat diversity. The species-area relationship for amphibians in the archipelago was formulated as the power function (S = 0.47A0.21) or exponential function (S = 2.59 + 2.41 (logA)). Our results suggested that eDNA metabarcoding is more sensitive for the detection of amphibian species. The combined use of eDNA metabarcoding and the traditional line transect method may optimize the survey results for amphibians.

Environmental DNA methods for biosecurity and invasion biology in terrestrial ecosystems: Progress, pitfalls, and prospects

Analysis of environmental DNA (eDNA) enables indirect detection of species without the need to directly observe and sample them. For biosecurity and invasion biology, eDNA-based methods are useful to address biological invasions at all phases, from detecting arrivals to confirming eradication of past invasions. We conducted a systematic review of the literature and found that in biosecurity and invasion biology, eDNA has primarily been used to detect new incursions and monitor spread in marine and freshwater ecosystems, with much slower uptake in terrestrial ecosystems, reflecting a broader trend common to the usage of eDNA tools. In terrestrial ecosystems, eDNA research has mostly focussed on the use of eDNA metabarcoding to characterise biodiversity, rather than targeting biosecurity threats or non-native populations. We discuss how eDNA-based methods are being applied to terrestrial ecosystems for biosecurity and managing non-native populations at each phase of the invasion continuum: transport, introduction, establishment, and spread; across different management options: containment, control, and eradication; and for detecting the impact of non-native organisms. Finally, we address some of the current technical issues and caveats of eDNA-based methods, particularly for terrestrial ecosystems, and how these might be solved. As eDNA-based methods improve, they will play an increasingly important role in the early detection and adaptive management of biological invasions, and the implementation of effective biosecurity controls.

Body size mediates the functional potential of soil organisms by diversity and community assembly across soil aggregates

Body size is an important life-history trait that affects organism niche occupancy and ecological interactions. However, it is still unclear to what extent the assembly process of organisms with different body sizes affects soil biogeochemical cycling processes at the aggregate level. Here, we examined the diversity and community assembly of soil microorganisms (bacteria, fungi, and protists) and microfauna (nematodes) with varying body sizes. The microbial functional potential associated with carbon, nitrogen, phosphorus, and sulfur metabolism within three soil aggregate sizes (large macroaggregates, > 2 mm; small macroaggregates, 0.25-2 mm; and microaggregates, < 0.25 mm) were determined by metagenomics. We found that the smallest microbes (bacteria) had higher α-diversity and lower β-diversity and were mostly structured by stochastic processes, while all larger organisms (fungi, protists, and nematodes) had lower α-diversity and were relatively more influenced by deterministic processes. Structural equation modeling indicated that the microbial functional potential associated with carbon, nitrogen, phosphorus, and sulfur metabolism was mainly influenced by the bacterial and protist diversity in microaggregates. In contrast, the microbial functional potential was primarily mediated by the assembly processes of four organism groups, especially the nematode community in macroaggregates. This study reveals the important roles of soil organisms with different body sizes in the functional potential related to nutrient cycling, and provides new insights into the ecological processes structuring the diversity and community assembly of organisms of different body sizes at the soil aggregate level, with implications for soil nutrient cycling dynamics.

Improving the efficiency of adaptive management methods in multiple fishways using environmental DNA

Dams and weirs impede the continuity of rivers and transit of migratory fish. To overcome this obstacle, fishways are installed worldwide; however, management after installation is important. The Miyanaka Intake Dam has three fish ladders with different flow velocities and discharges and has been under adaptive management since 2012. Fish catch surveys, conducted as an adaptive management strategy, place a heavy burden on fish. Furthermore, a large number of investigators must be mobilized during the 30-day investigation period. Thus, a monitoring method using environmental DNA that exerts no burden on fish and requires only a few surveyors (to obtain water samples) and an in-house analyst was devised; however, its implementation in a fishway away from the point of analysis and with limited flow space and its effective water sampling frequency have not been reported. Therefore, in 2019, we started a trial aiming to evaluate the methods and application conditions of environmental DNA surveys for the continuous and long-term monitoring of various fish fauna upstream and downstream of the Miyanaka Intake Dam. To evaluate the fish fauna, the results of an environmental DNA survey (metabarcoding method) for 2019 to 2022 were compared to those of a catch survey in the fishway from 2012 to 2022. The results confirmed the use of environmental DNA surveys in evaluating the contribution of fishways to biodiversity under certain conditions and introduced a novel method for sample collection.

Coupling remote sensing and eDNA to monitor environmental impact: A pilot to quantify the environmental benefits of sustainable agriculture in the Brazilian Amazon

Monitoring is essential to ensure that environmental goals are being achieved, including those of sustainable agriculture. Growing interest in environmental monitoring provides an opportunity to improve monitoring practices. Approaches that directly monitor land cover change and biodiversity annually by coupling the wall-to-wall coverage from remote sensing and the site-specific community composition from environmental DNA (eDNA) can provide timely, relevant results for parties interested in the success of sustainable agricultural practices. To ensure that the measured impacts are due to the environmental projects and not exogenous factors, sites where projects have been implemented should be benchmarked against counterfactuals (no project) and control (natural habitat) sites. Results can then be used to calculate diverse sets of indicators customized to monitor different projects. Here, we report on our experience developing and applying one such approach to assess the impact of shaded cocoa projects implemented by the Instituto de Manejo e Certificação Florestal e Agrícola (IMAFLORA) near São Félix do Xingu, in Pará, Brazil. We used the Continuous Degradation Detection (CODED) and LandTrendr algorithms to create a remote sensing-based assessment of forest disturbance and regeneration, estimate carbon sequestration, and changes in essential habitats. We coupled these remote sensing methods with eDNA analyses using arthropod-targeted primers by collecting soil samples from intervention and counterfactual pasture field sites and a control secondary forest. We used a custom set of indicators from the pilot application of a coupled monitoring framework called TerraBio. Our results suggest that, due to IMAFLORA's shaded cocoa projects, over 400 acres were restored in the intervention area and the community composition of arthropods in shaded cocoa is closer to second-growth forests than that of pastures. In reviewing the coupled approach, we found multiple aspects worked well, and we conclude by presenting multiple lessons learned.

Community metabarcoding reveals the relative role of environmental filtering and spatial processes in metacommunity dynamics of soil microarthropods across a mosaic of montane forests

Disentangling the relative role of environmental filtering and spatial processes in driving metacommunity structure across mountainous regions remains challenging, as the way we quantify spatial connectivity in topographically and environmentally heterogeneous landscapes can influence our perception of which process predominates. More empirical data sets are required to account for taxon- and context-dependency, but relevant research in understudied areas is often compromised by the taxonomic impediment. Here we used haplotype-level community DNA metabarcoding, enabled by stringent filtering of amplicon sequence variants (ASVs), to characterize metacommunity structure of soil microarthropod assemblages across a mosaic of five forest habitats on the Troodos mountain range in Cyprus. We found similar β diversity patterns at ASV and species (OTU, operational taxonomic unit) levels, which pointed to a primary role of habitat filtering resulting in the existence of largely distinct metacommunities linked to different forest types. Within-habitat turnover was correlated to topoclimatic heterogeneity, again emphasizing the role of environmental filtering. However, when integrating landscape matrix information for the highly fragmented Quercus alnifolia habitat, we also detected a major role of spatial isolation determined by patch connectivity, indicating that stochastic and niche-based processes synergistically govern community assembly. Alpha diversity patterns varied between ASV and OTU levels, with OTU richness decreasing with elevation and ASV richness following a longitudinal gradient, potentially reflecting a decline of genetic diversity eastwards due to historical pressures. Our study demonstrates the utility of haplotype-level community metabarcoding for characterizing metacommunity structure of complex assemblages and improving our understanding of biodiversity dynamics across mountainous landscapes worldwide.

Body size determines multitrophic soil microbiota community assembly associated with soil and plant attributes in a tropical seasonal rainforest

To understand soil biodiversity we need to know how soil communities are assembled. However, the relationship between soil community assembly and environmental factors, and the linkages between soil microbiota taxonomic groups and their body sizes, remain unexplored in tropical seasonal rainforests. Systematic and stratified random sampling was used to collect 243 soil and organism samples across a 20-ha plot in a tropical seasonal rainforest in southwestern China. High-throughput sequencing, variation analysis and principal coordinates of neighbourhood matrices were performed. Soil community composition, spatial distribution and assembly processes based on propagule size (including archaea, bacteria, fungi and nematodes) were investigated. The results showed that: (i) the community assembly of small soil microorganisms (bacteria, fungi) was mostly influenced by stochastic processes while that of larger soil organisms (nematodes) was more deterministic; (ii) the independent effects of habitat (including soil and topographic variables) and its interaction with plant attributes for community structure significantly decreased with increasing body size; and (iii) plant leaf phosphorus directly influenced the spatial distribution of soil-available phosphorus, which indicates their indirect impact on the assembly of the soil communities. Our data suggest that the assembly of multitrophic soil communities can be explained to some extent by changes in above-ground plant attributes. This highlights the importance of above- and below-ground linkages in influencing multitrophic soil microbiota community assembly.

Metabarcoding data reveal vertical multitaxa variation in topsoil communities during the colonization of deglaciated forelands

Ice-free areas are expanding worldwide due to dramatic glacier shrinkage and are undergoing rapid colonization by multiple lifeforms, thus representing key environments to study ecosystem development. It has been proposed that the colonization dynamics of deglaciated terrains is different between surface and deep soils but that the heterogeneity between communities inhabiting surface and deep soils decreases through time. Nevertheless, tests of this hypothesis remain scarce, and it is unclear whether patterns are consistent among different taxonomic groups. Here, we used environmental DNA metabarcoding to test whether community diversity and composition of six groups (Eukaryota, Bacteria, Mycota, Collembola, Insecta, and Oligochaeta) differ between the surface (0-5 cm) and deeper (7.5-20 cm) soil at different stages of development and across five Alpine glaciers. Taxonomic diversity increased with time since glacier retreat and with soil evolution. The pattern was consistent across groups and soil depths. For Eukaryota and Mycota, alpha-diversity was highest at the surface. Time since glacier retreat explained more variation of community composition than depth. Beta-diversity between surface and deep layers decreased with time since glacier retreat, supporting the hypothesis that the first 20 cm of soil tends to homogenize through time. Several molecular operational taxonomic units of bacteria and fungi were significant indicators of specific depths and/or soil development stages, confirming the strong functional variation of microbial communities through time and depth. The complexity of community patterns highlights the importance of integrating information from multiple taxonomic groups to unravel community variation in response to ongoing global changes.

Polar lake microbiomes have distinct evolutionary histories

Toward the poles, life on land is increasingly dominated by microorganisms, yet the evolutionary origin of polar microbiomes remains poorly understood. Here, we use metabarcoding of Arctic, sub-Antarctic, and Antarctic lacustrine benthic microbial communities to test the hypothesis that high-latitude microbiomes are recruited from a globally dispersing species pool through environmental selection. We demonstrate that taxonomic overlap between the regions is limited within most phyla, even at higher-order taxonomic levels, with unique deep-branching phylogenetic clades being present in each region. We show that local and regional taxon richness and net diversification rate of regionally restricted taxa differ substantially between polar regions in both microeukaryotic and bacterial biota. This suggests that long-term evolutionary divergence resulting from low interhemispheric dispersal and diversification in isolation has been a prominent process shaping present-day polar lake microbiomes. Our findings illuminate the distinctive biogeography of polar lake ecosystems and underscore that conservation efforts should include their unique microbiota.

Biological Interactions and Environmental Influences Shift Microeukaryotes in Permafrost Active Layer Soil Across the Qinghai-Tibet Plateau

Permafrost active layer soils are harsh environments with thaw/freeze cycles and sub-zero temperatures, harboring diverse microorganisms. However, the distribution patterns, assembly mechanism, and driving forces of soil microeukaryotes in permafrost remain largely unknown. In this study, we investigated microeukaryotes in permafrost active layer across the Qinghai-Tibet Plateau (QTP) using 18S rRNA gene sequencing. The results showed that the microbial eukaryotic communities were dominated by Nematozoa, Ciliophora, Ascomycota, Cercozoa, Arthropoda, and Basidiomycota in terms of relative abundance and operational taxonomic unit (OTU) richness. Nematozoa had the highest relative abundance, while Ciliophora had the highest OTU richness. These phyla had strong interactions between each other. Their alpha diversity and community structure were differently influenced by the factors associated to location, climate, and soil properties, particularly the soil properties. Significant but weak distance-decay relationships with different slopes were established for the communities of these dominant phyla, except for Basidiomycota. According to the null model, community assemblies of Nematozoa and Cercozoa were dominated by heterogeneous selection, Ciliophora and Ascomycota were dominated by dispersal limitation, while Arthropoda and Basidiomycota were highly dominated by non-dominant processes. The assembly mechanisms can be jointly explained by biotic interactions, organism treats, and environmental influences. Modules in the co-occurrence network of the microeukaryotes were composed by members from different taxonomic groups. These modules also had interactions and responded to different environmental factors, within which, soil properties had strong influences on these modules. The results suggested the importance of biological interactions and soil properties in structuring microbial eukaryotic communities in permafrost active layer soil across the QTP.

Defining a core microbial necrobiome associated with decomposing fungal necromass

Despite growing interest in fungal necromass decomposition due to its importance in soil carbon retention, whether a consistent group of microorganisms is associated with decomposing necromass remains unresolved. Here, we synthesize knowledge on the composition of the bacterial and fungal communities present on decomposing fungal necromass from a variety of fungal species, geographic locations, habitats, and incubation times. We found that there is a core group of both bacterial and fungal genera (i.e. a core fungal necrobiome), although the specific size of the core depended on definition. Based on a metric that included both microbial frequency and abundance, we demonstrate that the core is taxonomically and functionally diverse, including bacterial copiotrophs and oligotrophs as well as fungal saprotrophs, ectomycorrhizal fungi, and both fungal and animal parasites. We also show that the composition of the core necrobiome is notably dynamic over time, with many core bacterial and fungal genera having specific associations with the early, middle, or late stages of necromass decomposition. While this study establishes the existence of a core fungal necrobiome, we advocate that profiling the composition of fungal necromass decomposer communities in tropical environments and other terrestrial biomes beyond forests is needed to fill key knowledge gaps regarding the global nature of the fungal necrobiome.

Validation of the shotgun metabarcoding approach for comprehensively identifying herbal products containing plant, fungal, and animal ingredients

Identifying plant, fungal, and animal ingredients in a specific mixture remains challenging during the limitation of PCR amplification and low specificity of traditional methods. Genomic DNA was extracted from mock and pharmaceutical samples. Four type of DNA barcodes were generated from shotgun sequencing dataset with the help of a local bioinformatic pipeline. Taxa of each barcode was assigned by blast to TCM-BOL, BOLD, and GenBank. Traditional methods including microscopy, thin layer chromatography (TLC), and high-performance liquid chromatography (HPLC) were carried out according to Chinese pharmacopoeia. On average, 6.8 Gb shotgun reads were sequenced from genomic DNA of each sample. Then, 97, 11, 10, 14, and one operational taxonomic unit (OTU) were generated for ITS2, psbA-trnH, rbcL, matK, and COI, respectively. All the labeled ingredients including eight plant, one fungal, and one animal species were successfully detected in both the mock and pharmaceutical samples, in which Chebulae Fructus, Poria, and Fritilariae Thunbergia Bulbus were identified via mapping reads to organelle genomes. In addition, four unlabeled plant species were detected from pharmaceutical samples, while 30 genera of fungi, such as Schwanniomyces, Diaporthe, Fusarium were detected from mock and pharmaceutical samples. Furthermore, the microscopic, TLC, and HPLC analysis were all in accordance with the standards stipulated by Chinese Pharmacopoeia. This study indicated that shotgun metabarcoding could simultaneously identified plant, fungal, and animal ingredients in herbal products, which has the ability to serve as a valuable complement to traditional methods.

Weak environmental adaptation of rare phylotypes sustaining soil multi-element cycles in response to decades-long fertilization

Deciphering the ecological role of soil communities in the maintenance of multiple ecosystem functions is pivotal for the conservation and sustainability of soil biodiversity. However, few studies have investigated niche differentiation of abundant and rare microbiota, as well as their contributions to multiple soil elemental cycles, particularly in agroecosystems that have received decades of intense fertilization. Here, we characterized the environmental thresholds and phylogenetic signals for the environmental adaptation of both abundant and rare microbial subcommunities via amplicon sequencing and metagenomic sequencing and explored their importance in sustaining soil multiple nutrient cycling in agricultural fields that were fertilized for two decades. The results showed that rare taxa exhibited narrower niche breadths and weaker phylogenetic signals than abundant species. The assembly of abundant subcommunity was shaped predominantly by dispersal limitation (explained 71.1 % of the variation in bacteria) and undominated processes (explained 75 % of the variation in fungi), whereas the assembly of rare subcommunity was dominated by homogeneous selection process (explained 100 % of the variation in bacteria and 60 % of the variation in fungi). Soil ammonia nitrogen was the leading factor mediating the balance between stochastic and deterministic processes in both abundant (R² = 0.15, P < 0.001) and rare (R² = 0.08, P < 0.001) bacterial communities. Notably, the rare biosphere largely contributed to key soil processes such as carbon (R²_bacteria = 0.03, P < 0.05; R²_fungi = 0.05, P < 0.05) and nitrogen (R²_bacteria = 0.03, P < 0.05; R²_fungi = 0.17, P < 0.001) cycling. Collectively, these findings facilitate our understanding of the maintenance of rhizosphere bacterial and fungal diversity in response to agricultural fertilization and highlight the key role of rare taxa in sustaining agricultural ecosystem functions.

Comparative analysis of fish environmental DNA reveals higher sensitivity achieved through targeted sequence-based metabarcoding

Environmental DNA (eDNA)-based methods of species detection are enabling various applications in ecology and conservation including large-scale biomonitoring efforts. qPCR is widely used as the standard approach for species-specific detection, often targeting a fish species of interest from aquatic eDNA. However, DNA metabarcoding has the potential to displace qPCR in certain eDNA applications. In this study, we compare the sensitivity of the latest Illumina NovaSeq 6000 NGS platform to qPCR TaqMan assays by measuring limits of detection and by analysing eDNA from water samples collected from Churchill River and Lake Melville, NL, Canada. Species-specific, targeted next generation sequencing (NGS) assays had significantly higher sensitivity than qPCR, with limits of detection 14- to 29-fold lower. For example, when analysing eDNA, qPCR detected Gadus ogac (Greenland cod) in 21% of samples, but targeted NGS detected this species in 29% of samples. General NGS assays were as sensitive as qPCR, while simultaneously detecting 15 fish species from eDNA samples. With over 34,000 fish species on the planet, parallel and sensitive methods such as NGS will be required to support effective biomonitoring at both regional and global scales.

DNA High-Throughput Sequencing for Arthropod Gut Content Analysis to Evaluate Effectiveness and Safety of Biological Control Agents

The search for effective biological control agents without harmful non-target effects has been constrained by the use of impractical (field direct observation) or imprecise (cage experiments) methods. While advances in the DNA sequencing methods, more specifically the development of high-throughput sequencing (HTS), have been quickly incorporated in biodiversity surveys, they have been slow to be adopted to determine arthropod prey range, predation rate and food web structure, and critical information to evaluate the effectiveness and safety of a biological control agent candidate. The lack of knowledge on how HTS methods could be applied by ecological entomologists constitutes part of the problem, although the lack of expertise and the high cost of the analysis also are important limiting factors. In this review, we describe how the latest HTS methods of metabarcoding and Lazaro, a method to identify prey by mapping unassembled shotgun reads, can serve biological control research, showing both their power and limitations. We explain how they work to determine prey range and also how their data can be used to estimate predation rates and subsequently be translated into food webs of natural enemy and prey populations helping to elucidate their role in the community. We present a brief history of prey detection through molecular gut content analysis and also the attempts to develop a more precise formula to estimate predation rates, a problem that still remains. We focused on arthropods in agricultural ecosystems, but most of what is covered here can be applied to natural systems and non-arthropod biological control candidates as well.

Beyond biogeographic patterns: Processes shaping the microbial landscape in soils and sediments along the Yangtze River

Deciphering biogeographic patterns of microorganisms is important for evaluating the maintenance of microbial diversity with respect to the ecosystem functions they drives. However, ecological processes shaping distribution patterns of microorganisms across large spatial-scale watersheds remain largely unknown. Using Illumina sequencing and multiple statistical methods, we characterized distribution patterns and maintenance diversity of microorganisms (i.e., archaea, bacteria, and fungi) in soils and sediments along the Yangtze River. Distinct microbial distribution patterns were found between soils and sediments, and microbial community similarity significantly decreased with increasing geographical distance. Physicochemical properties showed a larger effect on microbial community composition than geospatial and climatic factors. Archaea and fungi displayed stronger species replacements and weaker environmental constraints in soils than that in sediments, but opposite for bacteria. Archaea, bacteria, and fungi in soils showed broader environmental breadths and stronger phylogenetic signals compared to those in sediments, suggesting stronger environmental adaptation. Stochasticity dominated community assemblies of archaea and fungi in soils and sediments, whereas determinism dominated bacterial community assembly. Our results have therefore highlighted distinct microbial distribution patterns and diversity maintenance mechanisms between soils and sediments, and emphasized important roles of species replacement, environmental adaptability, and ecological assembly processes on microbial landscape. Our findings are helpful in predicting loss of microbial diversity in the Yangtze River Basin, and might assist the establishment of environmental policies for protecting fragile watersheds.

Rethinking the nature of intraspecific variability and its consequences on species coexistence

Intraspecific variability (IV) has been proposed to explain species coexistence in diverse communities. Assuming, sometimes implicitly, that conspecific individuals can perform differently in the same environment and that IV increases niche overlap, previous studies have found contrasting results regarding the effect of IV on species coexistence. We aim at showing that the large IV observed in data does not mean that conspecific individuals are necessarily different in their response to the environment and that the role of high-dimensional environmental variation in determining IV has largely remained unexplored in forest plant communities. We first used a simulation experiment where an individual attribute is derived from a high-dimensional model, representing "perfect knowledge" of individual response to the environment, to illustrate how large observed IV can result from "imperfect knowledge" of the environment. Second, using growth data from clonal Eucalyptus plantations in Brazil, we estimated a major contribution of the environment in determining individual growth. Third, using tree growth data from long-term tropical forest inventories in French Guiana, Panama and India, we showed that tree growth in tropical forests is structured spatially and that despite a large observed IV at the population level, conspecific individuals perform more similarly locally than compared with heterospecific individuals. As the number of environmental dimensions that are well quantified at fine scale is generally lower than the actual number of dimensions influencing individual attributes, a great part of observed IV might be represented as random variation across individuals when in fact it is environmentally driven. This mis-representation has important consequences for inference about community dynamics. We emphasize that observed IV does not necessarily impact species coexistence per se but can reveal species response to high-dimensional environment, which is consistent with niche theory and the observation of the many differences between species in nature.

Differentiation strategies for planktonic bacteria and eukaryotes in response to aggravated algal blooms in urban lakes

Aggravated algal blooms potentially decreased environmental heterogeneity. Different strategies of planktonic bacteria and eukaryotes in response to aggravated algal blooms. Environmental constraints of plankton showed different patterns over time.

Habitat type and host grazing regimen influence the soil microbial diversity and communities within potential biting midge larval habitats

BACKGROUND

Biting midges (Culicoides spp.) are important vectors of diverse microbes such as viruses, protozoa, and nematodes that cause diseases in wild and domestic animals. However, little is known about the role of microbial communities in midge larval habitat utilization in the wild. In this study, we characterized microbial communities (bacterial, protistan, fungal and metazoan) in soils from disturbed (bison and cattle grazed) and undisturbed (non-grazed) pond and spring potential midge larval habitats. We evaluated the influence of habitat and grazing disturbance and their interaction on microbial communities, diversity, presence of midges, and soil properties.

RESULTS

Bacterial, protistan, fungal and metazoan community compositions were significantly influenced by habitat and grazing type. Irrespective of habitat and grazing type, soil communities were dominated by phyla Acidobacteria, Actinobacteria, Bacteroidetes, Chloroflexi, Firmicutes, Proteobacteria (Bacteria); Apicomplexa, Cercozoa, Ciliophora, Ochrophyta (Protists); Chytridiomycota, Cryptomycota (Fungi) and Nematoda, Arthropoda (Metazoa). The relative abundance of Acidobacteria, Actinobacteria, Bacteroidetes, Chloroflexi, Firmicutes, Proteobacteria, Verrucomicrobia (Bacteria); Apicomplexa, Lobosa (Protists); Ascomycota, Blastomycotina, Cryptomycota (Fungi); and Platyhelminthes (Metazoa) were significantly affected by grazing type. Of note, midge prevalence was higher in grazed sites (67-100%) than non-grazed (25%). Presence of midges in the soil was negatively correlated with bacterial, protistan, fungal and metazoan beta diversities and metazoan species richness but positively correlated with protistan and fungal species richness. Moreover, total carbon (TC), nitrogen (TN) and organic matter (OM) were negatively correlated with the presence of midges and relative abundances of unclassified Solirubrobacterales (Bacteria) and Chlamydomonadales (Protists) but positively with Proteobacteria and unclassified Burkholderiales (Bacteria).

CONCLUSIONS

Habitat and grazing type shaped the soil bacterial, protistan, fungal and metazoan communities, their compositions and diversities, as well as presence of midges. Soil properties (TN, TC, OM) also influenced soil microbial communities, diversities and the presence of midges. Prevalence of midges mainly in grazed sites indicates that midges prefer to breed and shelter in a habitat with abundant hosts, probably due to greater accessibility of food (blood meals). These results provide a first glimpse into the microbial communities, soil properties and prevalence of midges in suspected midge larval habitats at a protected natural prairie site.

Faecal metabarcoding provides improved detection and taxonomic resolution for non-invasive monitoring of gastrointestinal nematode parasites in wild moose populations

BACKGROUND

Although wild ungulate populations are heavily monitored throughout Europe, we understand little of how parasites affect population dynamics, and there is no systematic, long-term monitoring of parasite diversity and parasite loads. Such monitoring is in part hampered by a lack of time- and cost-effective assay methodologies with high sensitivity and good taxonomic resolution. DNA metabarcoding has been successfully used to characterize the parasitic nemabiome with high taxonomic resolution in a variety of wild and domestic hosts. However, in order to implement this technique in large-scale, potentially non-invasive monitoring of gastrointestinal parasitic nematodes (GIN), protocol optimization is required to maximize biodiversity detection, whilst maintaining time- and cost-effectiveness.

METHODS

Faecal samples were collected from a wild moose population and GIN communities were characterized and quantified using both parasitological techniques (egg and larva counting) and DNA metabarcoding of the ITS2 region of rDNA. Three different isolation methods were compared that differed in the volume of starting material and cell lysis method.

RESULTS

Similar nematode faunas were recovered from all samples using both parasitological and metabarcoding methods, and the approaches were largely congruent. However, metabarcoding assays showed better taxonomic resolution and slightly higher sensitivity than egg and larvae counts. The metabarcoding was not strictly quantitative, but the proportion of target nematode sequences recovered was correlated with the parasitologically determined parasite load. Species detection rates in the metabarcoding assays were maximized using a DNA isolation method that included mechanical cell disruption and maximized the starting material volume.

CONCLUSIONS

DNA metabarcoding is a promising technique for the non-invasive, large-scale monitoring of parasitic GINs in wild ungulate populations, owing to its high taxonomic resolution, increased assay sensitivity, and time- and cost-effectiveness. Although metabarcoding is not a strictly quantitative method, it may nonetheless be possible to create a management- and conservation-relevant index for the host parasite load from this data. To optimize the detection rates and time- and cost-effectiveness of metabarcoding assays, we recommend choosing a DNA isolation method that involves mechanical cell disruption and maximizes the starting material volume.

Biogeographic Pattern and Network of Rhizosphere Fungal and Bacterial Communities in Panicum miliaceum Fields: Roles of Abundant and Rare Taxa

Unraveling how microbial interactions and assembly process regulate the rhizosphere abundant and rare taxa is crucial for determining how species diversity affects rhizosphere microbiological functions. We assessed the rare and abundant taxa of rhizosphere fungal and bacterial communities in proso millet agroecosystems to explore their biogeographic patterns and co-occurrence patterns based on a regional scale. The taxonomic composition was significantly distinct between the fungal and bacterial abundant and rare taxa. Additionally, the rare taxa of bacteria and fungi exhibited higher diversity and stronger phylogenetic clustering than those of the abundant ones. The phylogenetic turnover rate of abundant taxa of bacteria was smaller than that of rare ones, whereas that of fungi had the opposite trend. Environmental variables, particularly mean annual temperature (MAT) and soil pH, were the crucial factors of community structure in the rare and abundant taxa. Furthermore, a deterministic process was relatively more important in governing the assembly of abundant and rare taxa. Our network analysis suggested that rare taxa of fungi and bacteria were located at the core of maintaining ecosystem functions. Interestingly, MAT and pH were also the important drivers controlling the main modules of abundant and rare taxa. Altogether, these observations revealed that rare and abundant taxa of fungal and bacterial communities showed obvious differences in biogeographic distribution, which were based on the dynamic interactions between assembly processes and co-occurrence networks.

The community ecology perspective of omics data

The measurement of uncharacterized pools of biological molecules through techniques such as metabarcoding, metagenomics, metatranscriptomics, metabolomics, and metaproteomics produces large, multivariate datasets. Analyses of these datasets have successfully been borrowed from community ecology to characterize the molecular diversity of samples (ɑ-diversity) and to assess how these profiles change in response to experimental treatments or across gradients (β-diversity). However, sample preparation and data collection methods generate biases and noise which confound molecular diversity estimates and require special attention. Here, we examine how technical biases and noise that are introduced into multivariate molecular data affect the estimation of the components of diversity (i.e., total number of different molecular species, or entities; total number of molecules; and the abundance distribution of molecular entities). We then explore under which conditions these biases affect the measurement of ɑ- and β-diversity and highlight how novel methods commonly used in community ecology can be adopted to improve the interpretation and integration of multivariate molecular data. Video Abstract.

Spatial assortment of soil organisms supports the size-plasticity hypothesis

The size-plasticity hypothesis posits that larger size organisms are less plastic in their metabolic rates and, therefore, are more strongly environmental-filtered than smaller organisms. Many studies have supported this hypothesis by evaluating the relative roles of environmental filtration and dispersal for different taxonomic groups of soil organisms. Most observations are made at large spatial scales, which are assumed to have a wide array of varying habitats. However, since urbanization causes habitat fragmentation at smaller regional scales, testing the size-plasticity hypothesis at this scale would help better understand the spatial assortment of urban soil organisms which, in turn, would help to develop improved management and conservation strategies for urban soil health. Here, we used DNA metabarcoding on five groups of soil biota (bacteria, fungi, protists, nematodes, and invertebrates) to assess the relative importance of dispersal and environmental filters to examine the size-plasticity hypothesis at this spatial scale in an urban environment. We observed strong distance-decay of community similarities associated with higher levels of stochastic changes in bacteria, nematode, and protist communities but not fungal or invertebrate communities. Bacterial communities occupied the widest niche followed by protists and nematodes, potentially because of their higher dispersal abilities compared to the larger soil organisms. Null deviation of communities varied with taxonomic groups where bacteria and nematodes were mainly driven by homogenizing dispersal, protists and fungi by drift, and soil invertebrates by environmental selection. We further identified a small percentage of locally-adapted taxa (2.1%) that could be focal taxa for conservation and restoration efforts by, for example, restoring their habitats and enhancing their regional connectivity. These results support the size-plasticity hypothesis at the relatively unexplored regional scale in an urbanization context, and provide new information for improving urban soil health and sustainable city models.

Soil microbial community assembly and stability are associated with potato (Solanum tuberosum L.) fitness under continuous cropping regime

Continuous cropping obstacles caused by the over-cultivation of a single crop trigger soil degradation, yield reduction and the occurrence of plant disease. However, the relationships among stability, complexity and the assembly process of soil microbial community with continuous cropping obstacles remains unclear. In this study, molecular ecological networks analysis (MENs) and inter-domain ecological networks analysis (IDENs), and a new index named cohesion tools were used to calculate the stability and complexity of soil microbial communities from eight potato cultivars grown under a continuous cropping regime by using the high-throughput sequencing data. The results showed that the stability (i.e., robustness index) of the bacterial and fungal communities for cultivar ZS5 was significantly higher, and that the complexity (i.e., cohesion values) was also significantly higher in the bacterial, fungal and inter-domain communities (i.e., bacterial-fungal community) of cultivar ZS5 than other cultivars. Network analysis also revealed that Actinobacteria and Ascomycota were the dominant phyla within intra-domain networks of continuous cropping potato soil communities, while the phyla Proteobacteria and Ascomycota dominated the correlation of the bacterial-fungal network. Infer community assembly mechanism by phylogenetic-bin-based null model analysis (iCAMP) tools were used to calculate the soil bacterial and fungal communities' assembly processes of the eight potato cultivars under continuous cropping regime, and the results showed that the bacterial community was mainly dominated by deterministic processes (64.19% - 81.31%) while the fungal community was mainly dominated by stochastic processes (78.28% - 98.99%), indicating that the continuous-cropping regime mainly influenced the potato soil bacterial community assembly process. Moreover, cultivar ZS5 possessed a relatively lower homogeneous selection, and a higher TP, TN, AP and yield than other cultivars. Our results indicated that the soil microbial network stability and complexity, and community assemble might be associated with yield and soil properties, which would be helpful in the study for resistance to potato continuous cropping obstacles.

Plant traits and community composition drive the assembly processes of abundant and rare fungi across deserts

The difference in community assembly mechanisms between rare and abundant fungi in deserts remains unknown. Hence, we compared the distribution patterns of abundant and rare fungi, and assessed the factors driving their assembly mechanisms across major vegetation types (shrubby desert, semi-shrubby and dwarf semi-shrubby desert, dwarf semi-arboreous desert, and shrubby steppe desert) of Chinese deserts. We assessed abundant and rare fungal subcommunities base on the sequencing data of fungal ITS data. Abundant fungal assembly was more affected by neutral processes than the rare. Null model and VPA analysis indicated that heterogeneous selection dominated rare sub-communities, whereas abundant fungal assembly was mainly determined by heterogeneous selection, dispersal limitation and other, unknown processes together. As a result, abundant sub-communities exhibited a higher species turnover rate than the rare. Hierarchical partitioning analysis indicated that soil conditions and plant attributes drove the assembly processes of abundant and rare fungi, respectively. Meanwhile, the relative strength of different assembly processes differed significantly among four vegetation types. In addition, we found that plant functional traits and composition played more critical roles in shaping the assembly processes of rare fungi than those of abundant fungi. Taken together, our findings collectively suggest that rare and abundant fungi exhibit differential ecological patterns that are driven by distinct assembly processes in deserts. We emphasize that the assembly processes of abundant and rare fungi are dependent on different abiotic and biotic factors in desert ecosystems.

Dispersal limitation dominates the community assembly of abundant and rare fungi in dryland montane forests

The assembly mechanisms and drivers of abundant and rare fungi in dryland montane forest soils remain underexplored. Therefore, in this study, we compared the distribution patterns of abundant and rare fungi and explored the factors determining their assembly processes in a dryland montane forest in China. Stronger distance-decay relationships (DDRs) were found in abundant sub-communities than in rare sub-communities. In addition, abundant fungi exhibited greater presence and wider habitat niche breadth than rare fungi. Both the null model and variation partitioning analysis indicated that dispersal limitation and environmental selection work together to govern both abundant and rare fungal assembly, while dispersal limitation plays a dominant role. Meanwhile, the relative influence of dispersal limitation and environmental selection varied between abundant and rare sub-communities, where dispersal limitation showed greater dominance in abundant fungal assembly. Mantel tests demonstrated that soil pH and phosphorus played critical roles in mediating abundant and rare fungi assembly processes, respectively. Our findings highlight that the distinct biogeographic patterns of abundant and rare fungi are driven by different assembly mechanisms, and the assembly processes of abundant and rare fungi are determined by diverse ecological drivers in dryland montane forest soils.

Environmental DNA Metabarcoding: A Novel Contrivance for Documenting Terrestrial Biodiversity

Simple Summary

The innovative concept of environmental DNA has found its foot in aquatic ecosystems but remains an unexplored area of research concerning terrestrial ecosystems. When making management choices, it is important to understand the rate of eDNA degradation, the persistence of DNA in terrestrial habitats, and the variables affecting eDNA detectability for a target species. Therefore an attempt has been made to provide comprehensive information regarding the exertion of eDNA in terrestrial ecosystems from 2012 to 2022. The information provided will assist ecologists, researchers and decision-makers in developing a holistic understanding of environmental DNA and its applicability as a biodiversity monitoring contrivance.

Abstract

The dearth of cardinal data on species presence, dispersion, abundance, and habitat prerequisites, besides the threats impeded by escalating human pressure has enormously affected biodiversity conservation. The innovative concept of eDNA, has been introduced as a way of overcoming many of the difficulties of rigorous conventional investigations, and is hence becoming a prominent and novel method for assessing biodiversity. Recently the demand for eDNA in ecology and conservation has expanded exceedingly, despite the lack of coordinated development in appreciation of its strengths and limitations. Therefore it is pertinent and indispensable to evaluate the extent and significance of eDNA-based investigations in terrestrial habitats and to classify and recognize the critical considerations that need to be accounted before using such an approach. Presented here is a brief review to summarize the prospects and constraints of utilizing eDNA in terrestrial ecosystems, which has not been explored and exploited in greater depth and detail in such ecosystems. Given these obstacles, we focused primarily on compiling the most current research findings from journals accessible in eDNA analysis that discuss terrestrial ecosystems (2012–2022). In the current evaluation, we also review advancements and limitations related to the eDNA technique.

Stronger responses of soil protistan communities to legacy mercury pollution than bacterial and fungal communities in agricultural systems

Soil pollution is an important stressor affecting biodiversity and ecosystem functioning. However, we lack a holistic understanding of how soil microbial communities respond to heavy metal pollution in agricultural ecosystems. Here, we explored the distribution patterns and inter-kingdom interactions of entire soil microbiome (including bacteria, fungi, and protists) in 47 paired paddy and upland fields along a gradient of legacy mercury (Hg) pollution. We found that the richness and composition of protistan community had stronger responses to Hg pollution than those of bacterial and fungal communities in both paddy and upland soils. Mercury polluted soils harbored less protistan phototrophs but more protistan consumers. We further revealed that long-term Hg pollution greatly increased network complexity of protistan community than that of bacterial and fungal communities, as well as intensified the interactions between protists and the other microorganisms. Moreover, our results consistently indicated that protistan communities had stronger responses to long-term Hg pollution than bacterial and fungal communities in agricultural soils based on structural equation models and random forest analyses. Our study highlights that soil protists can be used as bioindicators of Hg pollution, with important implications for the assessment of contaminated farmlands and the sustainable management of agricultural ecosystems.

Long-term watermelon continuous cropping leads to drastic shifts in soil bacterial and fungal community composition across gravel mulch fields

Despite the known influence of continuous cropping on soil microorganisms, little is known about the associated difference in the effects of continuous cropping on the community compositions of soil bacteria and fungi. Here, we assessed soil physicochemical property, as well as bacterial and fungal compositions across different years (Uncropped control, 1, 6, 11, 16, and 21 years) and in the watermelon system of a gravel mulch field in the Loess Plateau of China. Our results showed that long-term continuous cropping led to substantial shifts in soil bacterial and fungal compositions. The relative abundances of dominant bacterial and fungal genera (average relative abundance > 1.0%) significantly varied among different continuous cropping years (P < 0.05). Structural equation models demonstrated that continuous cropping alter soil bacterial and fungal compositions mainly by causing substantial variations in soil attributes. Variations in soil pH, nutrient, salinity, and moisture content jointly explained 73% and 64% of the variation in soil bacterial and fungal compositions, respectively. Variations in soil moisture content and pH caused by continuous cropping drove the shifts in soil bacterial and fungal compositions, respectively (Mantel R = 0.74 and 0.54, P < 0.01). Furthermore, the variation in soil bacterial and fungal composition showed significant correlation with watermelon yield reduction (P < 0.01). Together, long-term continuous cropping can alter soil microbial composition, and thereby influencing watermelon yield. Our findings are useful for alleviating continuous cropping obstacles and guiding agricultural production.

Shotgun metagenomics of soil invertebrate communities reflects taxonomy, biomass, and reference genome properties

Metagenomics - shotgun sequencing of all DNA fragments from a community DNA extract - is routinely used to describe the composition, structure, and function of microorganism communities. Advances in DNA sequencing and the availability of genome databases increasingly allow the use of shotgun metagenomics on eukaryotic communities. Metagenomics offers major advances in the recovery of biomass relationships in a sample, in comparison to taxonomic marker gene-based approaches (metabarcoding). However, little is known about the factors which influence metagenomics data from eukaryotic communities, such as differences among organism groups, the properties of reference genomes, and genome assemblies.We evaluated how shotgun metagenomics records composition and biomass in artificial soil invertebrate communities at different sequencing efforts. We generated mock communities of controlled biomass ratios from 28 species from all major soil mesofauna groups: mites, springtails, nematodes, tardigrades, and potworms. We shotgun sequenced these communities and taxonomically assigned them with a database of over 270 soil invertebrate genomes.We recovered over 95% of the species, and observed relatively high false-positive detection rates. We found strong differences in reads assigned to different taxa, with some groups (e.g., springtails) consistently attracting more hits than others (e.g., enchytraeids). Original biomass could be predicted from read counts after considering these taxon-specific differences. Species with larger genomes, and with more complete assemblies, consistently attracted more reads than species with smaller genomes. The GC content of the genome assemblies had no effect on the biomass-read relationships. Results were similar among different sequencing efforts.The results show considerable differences in taxon recovery and taxon specificity of biomass recovery from metagenomic sequence data. The properties of reference genomes and genome assemblies also influence biomass recovery, and they should be considered in metagenomic studies of eukaryotes. We show that low- and high-sequencing efforts yield similar results, suggesting high cost-efficiency of metagenomics for eukaryotic communities. We provide a brief roadmap for investigating factors which influence metagenomics-based eukaryotic community reconstructions. Understanding these factors is timely as accessibility of DNA sequencing and momentum for reference genomes projects show a future where the taxonomic assignment of DNA from any community sample becomes a reality.

Comparison of traditional and DNA metabarcoding samples for monitoring tropical soil arthropods (Formicidae, Collembola and Isoptera)

The soil fauna of the tropics remains one of the least known components of the biosphere. Long-term monitoring of this fauna is hampered by the lack of taxonomic expertise and funding. These obstacles may potentially be lifted with DNA metabarcoding. To validate this approach, we studied the ants, springtails and termites of 100 paired soil samples from Barro Colorado Island, Panama. The fauna was extracted with Berlese-Tullgren funnels and then either sorted with traditional taxonomy and known, individual DNA barcodes ("traditional samples") or processed with metabarcoding ("metabarcoding samples"). We detected 49 ant, 37 springtail and 34 termite species with 3.46 million reads of the COI gene, at a mean sequence length of 233 bp. Traditional identification yielded 80, 111 and 15 species of ants, springtails and termites, respectively; 98%, 37% and 100% of these species had a Barcode Index Number (BIN) allowing for direct comparison with metabarcoding. Ants were best surveyed through traditional methods, termites were better detected by metabarcoding, and springtails were equally well detected by both techniques. Species richness was underestimated, and faunal composition was different in metabarcoding samples, mostly because 37% of ant species were not detected. The prevalence of species in metabarcoding samples increased with their abundance in traditional samples, and seasonal shifts in species prevalence and faunal composition were similar between traditional and metabarcoding samples. Probable false positive and negative species records were reasonably low (13-18% of common species). We conclude that metabarcoding of samples extracted with Berlese-Tullgren funnels appear suitable for the long-term monitoring of termites and springtails in tropical rainforests. For ants, metabarcoding schemes should be complemented by additional samples of alates from Malaise or light traps.

Environmental DNA of insects and springtails from caves reveals complex processes of eDNA transfer in soils

Subterranean environments host a substantial amount of biodiversity, however assessing the distribution of species living underground is still extremely challenging. Environmental DNA (eDNA) metabarcoding is a powerful tool to estimate biodiversity in poorly known environments and has excellent performance for soil organisms. Here, we tested 1) whether eDNA metabarcoding from cave soils/sediments allows to successfully detect springtails (Hexapoda: Collembola) and insects (Hexapoda: Insecta); 2) whether eDNA mostly represents autochthonous (cave-dwelling) organisms or it also incorporates information from species living in surface environments; 3) whether eDNA detection probability changes across taxa with different ecology. Environmental DNA metabarcoding analyses detected a large number of Molecular Operational Taxonomic Units (MOTUs) for both insects and springtails. For springtails, detection probability was high, with a substantial proportion of hypogean species, suggesting that eDNA provides good information on the distribution of these organisms in caves. Conversely, for insects most of MOTUs represented taxa living outside caves, and the majority of them represented taxa/organisms living in freshwater environments (Ephemeroptera, Plecoptera and Trichoptera). The eDNA of freshwater insects was particularly abundant in deep sectors of caves, far from the entrance. Furthermore, average detection probability of insects was significantly lower than the one of springtails. This suggests that cave soils/sediments act as "conveyer belts of biodiversity information", possibly because percolating water lead to the accumulation of eDNA of organisms living in nearby areas. Cave soils hold a complex mix of autochthonous and allochthonous eDNA. eDNA provided unprecedented information on the understudied subterranean cave organisms; analyses of detection probability and occupancy can help teasing apart local eDNA from the eDNA representing spatially-integrated biodiversity for whole landscape.

Best practices in metabarcoding of fungi: From experimental design to results

The development of high-throughput sequencing (HTS) technologies has greatly improved our capacity to identify fungi and unveil their ecological roles across a variety of ecosystems. Here we provide an overview of current best practices in metabarcoding analysis of fungal communities, from experimental design through molecular and computational analyses. By reanalysing published data sets, we demonstrate that operational taxonomic units (OTUs) outperform amplified sequence variants (ASVs) in recovering fungal diversity, a finding that is particularly evident for long markers. Additionally, analysis of the full-length ITS region allows more accurate taxonomic placement of fungi and other eukaryotes compared to the ITS2 subregion. Finally, we show that specific methods for compositional data analyses provide more reliable estimates of shifts in community structure. We conclude that metabarcoding analyses of fungi are especially promising for integrating fungi into the full microbiome and broader ecosystem functioning context, recovery of novel fungal lineages and ancient organisms as well as barcoding of old specimens including type material.

Unearthing soil arthropod diversity through DNA metabarcoding

DNA metabarcoding has the potential to greatly advance understanding of soil biodiversity, but this approach has seen limited application for the most abundant and species-rich group of soil fauna-the arthropods. This study begins to address this gap by comparing information on species composition recovered from metabarcoding two types of bulk samples (specimens, soil) from a temperate zone site and from bulk soil samples collected at eight sites in the Arctic. Analysis of 22 samples (3 specimen, 19 soil) revealed 410 arthropod OTUs belonging to 112 families, 25 orders, and nine classes. Studies at the temperate zone site revealed little overlap in species composition between soil and specimen samples, but more overlap at higher taxonomic levels (families, orders) and congruent patterns of α- and β-diversity. Expansion of soil analyses to the Arctic revealed locally rich, highly dissimilar, and spatially structured assemblages compatible with dispersal limited and environmentally driven assembly. The current study demonstrates that DNA metabarcoding of bulk soil enables rapid, large-scale assessments of soil arthropod diversity. However, deep sequence coverage is required to adequately capture the species present in these samples, and expansion of the DNA barcode reference library is necessary to improve taxonomic resolution of the sequences recovered through this approach.

Monitoring fish communities through environmental DNA metabarcoding in the fish pass system of the second largest hydropower plant in the world

The Itaipu Hydroelectric Power Plant is the second largest in the world in power generation. The artificial barrier created by its dam imposes an obstacle for fish migration. Thus, in 2002, a fish pass system, named Piracema Channel, was built to allow fish to access areas upstream of the reservoir. We tested the potential of environmental DNA metabarcoding to monitor the impact of both the dam and associated fish pass system in the Paraná River fish communities and to compare it with traditional monitoring methods. Using a fragment of the 12S gene, we characterized richness and community composition based on amplicon sequence variants, operational taxonomic units, and zero-radius OTUs. We combined GenBank and in-house data for taxonomic assignment. We found that different bioinformatics approaches showed similar results. Also, we found a decrease in fish diversity from 2019 to 2020 probably due to the recent extreme drought experienced in southeastern Brazil. The highest alpha diversity was recorded in the mouth of the fish pass system, located in a protected valley with the highest environmental heterogeneity. Despite the clear indication that the reference databases need to be continuously improved, our results demonstrate the analytical efficiency of the metabarcoding to monitor fish species.

Decay of similarity across tropical forest communities: integrating spatial distance with soil nutrients

Understanding the mechanisms that drive the change of biotic assemblages over space and time is the main quest of community ecology. Assessing the relative importance of dispersal and environmental species selection in a range of organismic sizes and motilities has been a fruitful strategy. A consensus for whether spatial and environmental distances operate similarly across spatial scales and taxa, however, has yet to emerge. We used censuses of four major groups of organisms (soil bacteria, fungi, ground insects, and trees) at two observation scales (1-m² sampling point vs. 2,500-m² plots) in a topographically standardized sampling design replicated in two tropical rainforests with contrasting relationships between spatial distance and nutrient availability. We modeled the decay of assemblage similarity for each taxon set and site to assess the relative contributions of spatial distance and nutrient availability distance. Then, we evaluated the potentially structuring effect of tree composition over all other taxa. The similarity of nutrient content in the litter and topsoil had a stronger and more consistent selective effect than did dispersal limitation, particularly for bacteria, fungi, and trees at the plot level. Ground insects, the only group assessed with the capacity of active dispersal, had the highest species turnover and the flattest nonsignificant distance-decay relationship, suggesting that neither dispersal limitation nor nutrient availability were fundamental drivers of their community assembly at this scale of analysis. Only the fungal communities at one of our study sites were clearly coordinated with tree composition. The spatial distance at the smallest scale was more important than nutrient selection for the bacteria, fungi, and insects. The lower initial similarity and the moderate variation in composition identified by these distance-decay models, however, suggested that the effects of stochastic sampling were important at this smaller spatial scale. Our results highlight the importance of nutrients as one of the main environmental drivers of rainforest communities irrespective of organismic or propagule size and how the overriding effect of the analytical scale influences the interpretation, leading to the perception of greater importance of dispersal limitation and ecological drift over selection associated with environmental niches at decreasing observation scales.

Community Assembly and Co-occurrence Patterns Underlying the Core and Satellite Bacterial Sub-communities in the Tibetan Lakes

Microbial communities normally comprise a few core species and large numbers of satellite species. These two sub-communities have different ecological and functional roles in natural environments, but knowledge on the assembly processes and co-occurrence patterns of the core and satellite species in Tibetan lakes is still sparse. Here, we investigated the ecological processes and co-occurrence relationships of the core and satellite bacterial sub-communities in the Tibetan lakes via 454 sequencing of 16S rRNA gene. Our studies indicated that the core and satellite bacterial sub-communities have similar dominant phyla (Proteobacteria, Bacteroidetes, and Actinobacteria). But the core sub-communities were less diverse and exhibited a stronger distance-decay relationship than the satellite sub-communities. In addition, topological properties of nodes in the network demonstrated that the core sub-communities had more complex and stable co-occurrence associations and were primarily driven by stochastic processes (58.19%). By contrast, the satellite sub-communities were mainly governed by deterministic processes (62.17%). Overall, this study demonstrated the differences in the core and satellite sub-community assembly and network stability, suggesting the importance of considering species traits to understand the biogeographic distribution of bacterial communities in high-altitude lakes.

Towards revealing the global diversity and community assembly of soil eukaryotes

Soil fungi, protists, and animals (i.e., the eukaryome) play a critical role in key ecosystem functions in terrestrial ecosystems. Yet, we lack a holistic understanding of the processes shaping the global distribution of the eukaryome. We conducted a molecular analysis of 193 composite soil samples spanning the world's major biomes. Our analysis showed that the importance of selection processes was higher in the community assemblage of smaller-bodied and wider niche breadth organisms. Soil pH and mean annual precipitation were the primary determinants of the community structure of eukaryotic microbes and animals, respectively. We further found contrasting latitudinal diversity patterns and strengths for soil eukaryotic microbes and animals. Our results point to a potential link between body size and niche breadth of soil eukaryotes and the relative effect of ecological processes and environmental factors in driving their biogeographic patterns.