Global hotspots for soil nature conservation.

Soils are the foundation of all terrestrial ecosystems1. However, unlike for plants and animals, a global assessment of hotspots for soil nature conservation is still lacking2. This hampers our ability to establish nature conservation priorities for the multiple dimensions that support the soil system: from soil biodiversity to ecosystem services. Here, to identify global hotspots for soil nature conservation, we performed a global field survey that includes observations of biodiversity (archaea, bacteria, fungi, protists and invertebrates) and functions (critical for six ecosystem services) in 615 composite samples of topsoil from a standardized survey in all continents. We found that each of the different ecological dimensions of soils-that is, species richness (alpha diversity, measured as amplicon sequence variants), community dissimilarity and ecosystem services-peaked in contrasting regions of the planet, and were associated with different environmental factors. Temperate ecosystems showed the highest species richness, whereas community dissimilarity peaked in the tropics, and colder high-latitudinal ecosystems were identified as hotspots of ecosystem services. These findings highlight the complexities that are involved in simultaneously protecting multiple ecological dimensions of soil. We further show that most of these hotspots are not adequately covered by protected areas (more than 70%), and are vulnerable in the context of several scenarios of global change. Our global estimation of priorities for soil nature conservation highlights the importance of accounting for the multidimensionality of soil biodiversity and ecosystem services to conserve soils for future generations.

The UNITE database for molecular identification and taxonomic communication of fungi and other eukaryotes: sequences, taxa and classifications reconsidered.

UNITE (https://unite.ut.ee) is a web-based database and sequence management environment for molecular identification of eukaryotes. It targets the nuclear ribosomal internal transcribed spacer (ITS) region and offers nearly 10 million such sequences for reference. These are clustered into ∼2.4M species hypotheses (SHs), each assigned a unique digital object identifier (DOI) to promote unambiguous referencing across studies. UNITE users have contributed over 600 000 third-party sequence annotations, which are shared with a range of databases and other community resources. Recent improvements facilitate the detection of cross-kingdom biological associations and the integration of undescribed groups of organisms into everyday biological pursuits. Serving as a digital twin for eukaryotic biodiversity and communities worldwide, the latest release of UNITE offers improved avenues for biodiversity discovery, precise taxonomic communication and integration of biological knowledge across platforms.

Interaction effects of pH and land cover on soil microbial diversity are climate-dependent.

Factors regulating the diversity and composition of soil microbial communities include soil properties, land cover and climate. How these factors interact at large scale remains poorly investigated. Here, we used an extensive dataset including 715 locations from 24 European countries to investigate the interactive effects of climatic region, land cover and pH on soil bacteria and fungi. We found that differences in microbial diversity and community composition between land cover types depended on the climatic region. In Atlantic, Boreal and Continental regions, microbial richness was higher in croplands and grasslands than woodlands while richness in Mediterranean areas did not vary significantly among land cover types. These differences were further related to soil pH, as a driver of bacterial and fungal richness in most climatic regions, but the interaction of pH with land cover depended on the region. Microbial community composition differed the most between croplands and woodlands in all regions, mainly due to differences in pH. In the Mediterranean region, bacterial communities in woodlands and grasslands were the most similar, whereas in other regions, grassland and cropland-associated bacteria showed more similarity. Overall, we showed that key factors interact in shaping soil microbial communities in a climate-dependent way at large scale.

Unearthing the soil-borne microbiome of land plants.

Plant-soil biodiversity interactions are fundamental for the functioning of terrestrial ecosystems. Yet, the existence of a set of globally distributed topsoil microbial and small invertebrate organisms consistently associated with land plants (i.e., their consistent soil-borne microbiome), together with the environmental preferences and functional capabilities of these organisms, remains unknown. We conducted a standardized field survey under 150 species of land plants, including 58 species of bryophytes and 92 of vascular plants, across 124 locations from all continents. We found that, despite the immense biodiversity of soil organisms, the land plants evaluated only shared a small fraction (less than 1%) of all microbial and invertebrate taxa that were present across contrasting climatic and soil conditions and vegetation types. These consistent taxa were dominated by generalist decomposers and phagotrophs and their presence was positively correlated with the abundance of functional genes linked to mineralization. Finally, we showed that crossing environmental thresholds in aridity (aridity index of 0.65, i.e., the transition from mesic to dry ecosystems), soil pH (5.5; i.e., the transition from acidic to strongly acidic soils), and carbon (less than 2%, the lower limit of fertile soils) can result in drastic disruptions in the associations between land plants and soil organisms, with potential implications for the delivery of soil ecosystem processes under ongoing global environmental change.

Structure and function of the global topsoil microbiome.

Soils harbour some of the most diverse microbiomes on Earth and are essential for both nutrient cycling and carbon storage. To understand soil functioning, it is necessary to model the global distribution patterns and functional gene repertoires of soil microorganisms, as well as the biotic and environmental associations between the diversity and structure of both bacterial and fungal soil communities1-4. Here we show, by leveraging metagenomics and metabarcoding of global topsoil samples (189 sites, 7,560 subsamples), that bacterial, but not fungal, genetic diversity is highest in temperate habitats and that microbial gene composition varies more strongly with environmental variables than with geographic distance. We demonstrate that fungi and bacteria show global niche differentiation that is associated with contrasting diversity responses to precipitation and soil pH. Furthermore, we provide evidence for strong bacterial-fungal antagonism, inferred from antibiotic-resistance genes, in topsoil and ocean habitats, indicating the substantial role of biotic interactions in shaping microbial communities. Our results suggest that both competition and environmental filtering affect the abundance, composition and encoded gene functions of bacterial and fungal communities, indicating that the relative contributions of these microorganisms to global nutrient cycling varies spatially.

Effects of nitrogen deposition on carbon and nutrient cycling along a natural soil acidity gradient as revealed by metagenomics.

Nitrogen (N) deposition and soil acidification are environmental challenges affecting ecosystem functioning, health, and biodiversity, but their effects on functional genes are poorly understood. Here, we utilized metabarcoding and metagenomics to investigate the responses of soil functional genes to N deposition along a natural soil pH gradient. Soil N content was uncorrelated with pH, enabling us to investigate their effects separately. Soil acidity strongly and negatively affected the relative abundances of most cluster of orthologous gene categories of the metabolism supercategory. Similarly, soil acidity negatively affected the diversity of functional genes related to carbon and N but not phosphorus cycling. Multivariate analyses showed that soil pH was the most important factor affecting microbial and functional gene composition, while the effects of N deposition were less important. Relative abundance of KEGG functional modules related to different parts of the studied cycles showed variable responses to soil acidity and N deposition. Furthermore, our results suggested that the diversity-function relationship reported for other organisms also applies to soil microbiomes. Since N deposition and soil pH affected microbial taxonomic and functional composition to a different extent, we conclude that N deposition effects might be primarily mediated through soil acidification in forest ecosystems.

Traces of Hymenoscyphus fraxineus in Northeastern Europe Extend Further Back in History than Expected.

Herbaria are a promising but still poorly applied information source for retrospective microbiological studies. In order to find any evidence of the virulent European origin of ash dieback agent Hymenoscyphus fraxineus and other fungal pathogens, we analyzed 109 leaf samples from three different Estonian botanical herbaria, sampled during 171 years from 20 ash species and cultivars, using a PacBio third-generation sequencing of the fungal internal transcribed spacer ITS1-5.8S-ITS2 ribosomal DNA region. We identified a large amount of saprotrophic fungi naturally colonizing ash leaves. Hymenoscyphus fraxineus colonized a Fraxinus chinensis subsp. rhynchophylla specimen and a F. chinensis specimen collected from Tallinn Botanic Garden in July 1978 and July 1992, respectively. The samples originated from trees grown in this garden from seeds collected from Shamora, Far-East Russia, in 1961 and from a Beijing botanical garden in eastern China in 1985, respectively. Repeated subsequent DNA extraction, real-time quantitative PCR, and Sanger and Illumina sequencing confirmed our findings of these apparently oldest cases of the ash dieback agent in Europe. These results show that H. fraxineus evidently was present in Estonia 19 years earlier than our previous data from fungal herbaria documented and 14 years before the first visible damage of ash trees was registered in Poland. Because we found no evidence of the saprotrophic H. albidus from earlier mycological and botanical herbarium specimens, the presence of H. albidus in Estonia remains questionable.

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.

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.

Community dynamics of soil-borne fungal communities along elevation gradients in neotropical and palaeotropical forests.

Because of their steep gradients in abiotic and biotic factors, mountains offer an ideal setting to illuminate the mechanisms that underlie patterns of species distributions and community assembly. We compared the composition of taxonomically and functionally diverse fungal communities in soils along five elevational gradients in mountains of the Neo- and Palaeotropics (northern Argentina, southern Brazil, Panama, Malaysian Borneo and Papua New Guinea). Both the richness and composition of soil fungal communities reflect environmental factors, particularly temperature and soil pH, with some shared patterns among neotropical and palaeotropical regions. Community dynamics are characterized by replacement of species along elevation gradients, implying a relatively narrow elevation range for most fungi, which appears to be driven by contrasting environmental preferences among both functional and taxonomic groups. For functional groups dependent on symbioses with plants (especially ectomycorrhizal fungi), the distribution of host plants drives richness and community composition, resulting in important differences in elevational patterns between neotropical and palaeotropical montane communities. The pronounced compositional and functional turnover along elevation gradients implies that tropical montane forest fungi will be sensitive to climate change, resulting in shifts in composition and functionality over time.

Global patterns in endemicity and vulnerability of soil fungi.

Fungi are highly diverse organisms, which provide multiple ecosystem services. However, compared with charismatic animals and plants, the distribution patterns and conservation needs of fungi have been little explored. Here, we examined endemicity patterns, global change vulnerability and conservation priority areas for functional groups of soil fungi based on six global surveys using a high-resolution, long-read metabarcoding approach. We found that the endemicity of all fungi and most functional groups peaks in tropical habitats, including Amazonia, Yucatan, West-Central Africa, Sri Lanka, and New Caledonia, with a negligible island effect compared with plants and animals. We also found that fungi are predominantly vulnerable to drought, heat and land-cover change, particularly in dry tropical regions with high human population density. Fungal conservation areas of highest priority include herbaceous wetlands, tropical forests, and woodlands. We stress that more attention should be focused on the conservation of fungi, especially root symbiotic arbuscular mycorrhizal and ectomycorrhizal fungi in tropical regions as well as unicellular early-diverging groups and macrofungi in general. Given the low overlap between the endemicity of fungi and macroorganisms, but high conservation needs in both groups, detailed analyses on distribution and conservation requirements are warranted for other microorganisms and soil organisms.

Global taxonomic and phylogenetic assembly of AM fungi.

Arbuscular mycorrhizal (AM) fungi are a ubiquitous group of plant symbionts, yet processes underlying their global assembly - in particular the roles of dispersal limitation and historical drivers - remain poorly understood. Because earlier studies have reported niche conservatism in AM fungi, we hypothesized that variation in taxonomic community composition (i.e., unweighted by taxon relatedness) should resemble variation in phylogenetic community composition (i.e., weighted by taxon relatedness) which reflects ancestral adaptations to historical habitat gradients. Because of the presumed strong dispersal ability of AM fungi, we also anticipated that the large-scale structure of AM fungal communities would track environmental conditions without regional discontinuity. We used recently published AM fungal sequence data (small-subunit ribosomal RNA gene) from soil samples collected worldwide to reconstruct global patterns in taxonomic and phylogenetic community variation. The taxonomic structure of AM fungal communities was primarily driven by habitat conditions, with limited regional differentiation, and there were two well-supported clusters of communities - occurring in cold and warm conditions. Phylogenetic structure was driven by the same factors, though all relationships were markedly weaker. This suggests that niche conservatism with respect to habitat associations is weakly expressed in AM fungal communities. We conclude that the composition of AM fungal communities tracks major climatic and edaphic gradients, with the effects of dispersal limitation and historic factors considerably less apparent than those of climate and soil.

Perspectives and Benefits of High-Throughput Long-Read Sequencing in Microbial Ecology.

Short-read, high-throughput sequencing (HTS) methods have yielded numerous important insights into microbial ecology and function. Yet, in many instances short-read HTS techniques are suboptimal, for example, by providing insufficient phylogenetic resolution or low integrity of assembled genomes. Single-molecule and synthetic long-read (SLR) HTS methods have successfully ameliorated these limitations. In addition, nanopore sequencing has generated a number of unique analysis opportunities, such as rapid molecular diagnostics and direct RNA sequencing, and both Pacific Biosciences (PacBio) and nanopore sequencing support detection of epigenetic modifications. Although initially suffering from relatively low sequence quality, recent advances have greatly improved the accuracy of long-read sequencing technologies. In spite of great technological progress in recent years, the long-read HTS methods (PacBio and nanopore sequencing) are still relatively costly, require large amounts of high-quality starting material, and commonly need specific solutions in various analysis steps. Despite these challenges, long-read sequencing technologies offer high-quality, cutting-edge alternatives for testing hypotheses about microbiome structure and functioning as well as assembly of eukaryote genomes from complex environmental DNA samples.

Seasonal dynamics of mycoplankton in the Yellow Sea reflect the combined effect of riverine inputs and hydrographic conditions.

Little is known about how multiple factors including land-based inputs and ocean currents affect the spatiotemporal distribution of the mycoplankton in coastal regions. To explore the seasonal changes of mycoplanktonic communities and potential environmental drivers, we collected water samples from the Yellow Sea, used here as a model for subtropical sea habitats, in different seasons over two years. Compared with winter and spring, summer exhibited higher levels of fungal richness and community heterogeneity in the water column. The seasonal shifts in mycoplankton diversity and community composition were mainly ascribed to freshwater inputs, the Cold Water Mass and invasion of the Yellow Sea Warm Current. Among the physicochemical variables tested, temperature was the primary determinant of fungal diversity and showed contrasting influences on fungal richness in the surface and bottom waters during summer. In addition, we provide evidence for the community similarity and dissolved nutrients of different water bodies to highlight the potential origin of the Cold Water Mass. Our findings bring new understanding on the factors determining the dynamics of mycoplankton communities by modelling the influence of physicochemical variables and tracking the geographical distribution of certain fungal taxa.

The UNITE database for molecular identification of fungi: handling dark taxa and parallel taxonomic classifications.

UNITE (https://unite.ut.ee/) is a web-based database and sequence management environment for the molecular identification of fungi. It targets the formal fungal barcode-the nuclear ribosomal internal transcribed spacer (ITS) region-and offers all ∼1 000 000 public fungal ITS sequences for reference. These are clustered into ∼459 000 species hypotheses and assigned digital object identifiers (DOIs) to promote unambiguous reference across studies. In-house and web-based third-party sequence curation and annotation have resulted in more than 275 000 improvements to the data over the past 15 years. UNITE serves as a data provider for a range of metabarcoding software pipelines and regularly exchanges data with all major fungal sequence databases and other community resources. Recent improvements include redesigned handling of unclassifiable species hypotheses, integration with the taxonomic backbone of the Global Biodiversity Information Facility, and support for an unlimited number of parallel taxonomic classification systems.

Surveying for Batrachochytrium dendrobatidis and B. salamandrivorans in wild and captive amphibian populations in Estonia and Latvia.

The pathogenic chytrid fungi Batrachochytrium dendrobatidis (Bd) and B. salamandrivorans (Bsal) cause infections that have become primary drivers of amphibian biodiversity loss. While globally widespread, the distribution margins of Bd and Bsal have not been determined, and the presence of these pathogens has probably gone unnoticed in many areas, especially in northern Eurasia. To better understand the presence and distribution of both pathogens in the northern temperate and boreal forest biomes, 243 individuals were sampled from 8 native amphibian species across Estonia. Additionally, 68 amphibians were sampled from captive collections in Estonia and Latvia. Pathogen infection was assessed using metabarcoding of the ITS2 marker. No positive matches for Bsal infection were found. Bd was detected in 13 specimens, 3 of which were sampled at the Riga Zoo (with a prevalence of 5.2%) and 10 in natural environments in Estonia (3.3%). The infected wild individuals belonged to 6 amphibian species and were detected throughout the mainland of Estonia, but not on islands. Prevalence of infection with Bd ranged between 3.1 and 12.5% among native species. In addition, we found molecular evidence for a potentially new sister species to Bd in nature. Although outbreaks of chytridiomycosis have never been observed in Estonia, it cannot be excluded that the dynamics of local amphibian populations are affected by Bd infections. Therefore, further work, including capture-mark-recapture studies and long-term monitoring, are required to clarify the impact of Bd on amphibians in Northern Europe.

Diversity of arbuscular mycorrhizal fungi and its chemical drivers across dryland habitats.

Qatar is largely characterized by a hyper-arid climate and low soil fertility which create a stressful soil environment for arbuscular mycorrhizal (AM) fungi. In a study of AM fungal communities and their relationship with soil chemical characteristics, we used a high-throughput sequencing technique to explore AM fungal diversity and community composition in different habitats across Qatar. We identified a total of 79 AM fungal taxa, over 77% of which were species from the Glomeraceae family. The lowest AM fungal diversity was observed in saltmarsh and in one rawdha site, while the highest richness, effective number of species, and diversity were observed in rawdha and sabkha communities. NMDS and multiple regression analyses showed that AM fungi were negatively correlated with soil pH and TC, but positively correlated with K and NO3-. AM fungi also were positively correlated with Cd, with the latter suggesting that very low levels of heavy metals may not always be harmful to AM fungi. These findings provide baseline information on AM fungal assemblages and the chemical drivers of diversity across communities in Qatar. This work partly compensates for the current lack of broad-scale studies in the Arabian Peninsula by providing understanding of overall patterns of AM fungi and their drivers in the region.

Plant nutrient-acquisition strategies drive topsoil microbiome structure and function.

Plant nutrient-acquisition strategies drive soil processes and vegetation performance, but their effect on the soil microbiome remains poorly understood. This knowledge is important to predict the shifts in microbial diversity and functions due to increasing changes in vegetation traits under global change. Here we documented the topsoil microbiomes of 145 boreal and temperate terrestrial sites in the Baltic region that broadly differed in vegetation type and nutritional traits, such as mycorrhizal types and symbiotic nitrogen-fixation. We found that sites dominated by arbuscular mycorrhizal (AM) vegetation harbor relatively more AM fungi, bacteria, fungal saprotrophs, and pathogens in the topsoil compared with sites dominated by ectomycorrhizal (EM) plants. These differences in microbiome composition reflect the rapid nutrient cycling and negative plant-soil feedback in AM soils. Lower fungal diversity and bacteria : fungi ratios in EM-dominated habitats are driven by monodominance of woody vegetation as well as soil acidification by EM fungi, which are associated with greater diversity and relative abundance of carbohydrate-active enzymes. Our study suggests that shifts in vegetation related to global change and land use may strongly alter the topsoil microbiome structure and function.

FungalRoot: global online database of plant mycorrhizal associations.

Testing of ecological, biogeographical and phylogenetic hypotheses of mycorrhizal traits requires a comprehensive reference dataset about plant mycorrhizal associations. Here we present a database, FungalRoot, which summarizes publicly available data about vascular plant mycorrhizal type and intensity of root colonization by mycorrhizal fungi, accompanied with rich metadata. We compiled and digitized data about plant mycorrhizal colonization in nine widespread languages. The present version of the FungalRoot database contains 36 303 species-by-site observations for 14 870 plant species, tripling the previously available compiled information about plant mycorrhizal associations. Based on these data, we provide a recommended list of genus-level plant mycorrhizal associations, based on the majority of data for species and careful analysis of conflicting data. The majority of ectomycorrhizal and ericoid mycorrhizal plants are trees (92%) and shrubs (85%), respectively. The majority of arbuscular and nonmycorrhizal plant species are herbaceous (50% and 70%, respectively). Our publicly available database is a powerful resource for mycorrhizal scientists and ecologists. It features possibilities for dynamic updating and addition of data about plant mycorrhizal associations. The new database will promote research on plant and fungal biogeography and evolution, and on links between above- and belowground biodiversity and ecosystem functioning.

Distinct fungal successional trajectories following wildfire between soil horizons in a cold-temperate forest.

Soil fungi represent a major component of below-ground biodiversity that determines the succession and recovery of forests after disturbance. However, their successional trajectories and driving mechanisms following wildfire remain unclear. We examined fungal biomass, richness, composition and enzymes across three soil horizons (Oe, A1 and A2) along a near-complete fire chronosequence (1, 2, 8, 14, 30, 49 and c. 260 yr) in cold-temperate forests of the Great Khingan Mountains, China. The importance of soil properties, spatial distance and tree composition were also tested. Ectomycorrhizal fungal richness and ß-glucosidase activity were strongly reduced by burning and significantly increased with 'time since fire' in the Oe horizon but not in the mineral horizons. Time since fire and soil C : N ratio were the primary drivers of fungal composition in the Oe and A1/A2 horizons, respectively. Ectomycorrhizal fungal composition was remarkably sensitive to fire history in the Oe horizon, while saprotroph community was strongly affected by time since fire in the deeper soil horizon and this effect emerged 18 years after fire in the A2 horizon. Our study demonstrates pronounced horizon-dependent successional trajectories following wildfire and indicates interactive effects of time since fire, soil stoichiometry and spatial distance in the reassembly of below-ground fungal communities in a cold and fire-prone region.