High-throughput identification technologies provide efficient tools for understanding the ecology and functioning of microorganisms. Yet, these methods have been only rarely used for monitoring and testing ecological hypotheses in plant pathogens and pests in spite of their immense importance in agriculture, forestry and plant community dynamics. The main objectives of this manuscript are the following: (a) to provide a comprehensive overview about the state-of-the-art high-throughput quantification and molecular identification methods used to address population dynamics, community ecology and host associations of microorganisms, with a specific focus on antagonists such as pathogens, viruses and pests; (b) to compile available information and provide recommendations about specific protocols and workable primers for bacteria, fungi, oomycetes and insect pests; and (c) to provide examples of novel methods used in other microbiological disciplines that are of great potential use for testing specific biological hypotheses related to pathology. Finally, we evaluate the overall perspectives of the state-of-the-art and still evolving methods for diagnostics and population- and community-level ecological research of pathogens and pests.

Arbuscular mycorrhizal fungi (AMF) strongly affect ecosystem functioning. To understand and quantify the mechanisms of this control, knowledge about the relationship between the actual abundance and community composition of AMF in the soil and in plant roots is needed. We collected soil and root samples in a natural dune grassland to test whether, across a plant community, the abundance of AMF in host roots (measured as the total length of roots colonized) is related to soil AMF abundance (using the neutral lipid fatty acids (NLFA) 16:1ω5 as proxy). Next-generation sequencing was used to explore the role of community composition in abundance patterns. We found a strong positive relationship between the total length of roots colonized by AMF and the amount of NLFA 16:1ω5 in the soil. We provide the first field-based evidence of proportional biomass allocation between intra-and extraradical AMF mycelium, at ecosystem level. We suggest that this phenomenon is made possible by compensatory colonization strategies of individual fungal species. Finally, our findings open the possibility of using AMF total root colonization as a proxy for soil AMF abundances, aiding further exploration of the AMF impacts on ecosystems functioning.

Soil microbiome has a pivotal role in ecosystem functioning, yet little is known about its build-up from local to regional scales. In a multi-year regional-scale survey involving 1251 plots and long-read third-generation sequencing, we found that soil pH has the strongest effect on the diversity of fungi and its multiple taxonomic and functional groups. The pH effects were typically unimodal, usually both direct and indirect through tree species, soil nutrients or mold abundance. Individual tree species, particularly Pinus sylvestris, Picea abies, and Populus x wettsteinii, and overall ectomycorrhizal plant proportion had relatively stronger effects on the diversity of biotrophic fungi than saprotrophic fungi. We found strong temporal sampling and investigator biases for the abundance of molds, but generally all spatial, temporal and microclimatic effects were weak. Richness of fungi and several functional groups was highest in woodlands and around ruins of buildings but lowest in bogs, with marked group-specific trends. In contrast to our expectations, diversity of soil fungi tended to be higher in forest island habitats potentially due to the edge effect, but fungal richness declined with island distance and in response to forest fragmentation. Virgin forests supported somewhat higher fungal diversity than old non-pristine forests, but there were no differences in richness between natural and anthropogenic habitats such as parks and coppiced gardens. Diversity of most fungal groups suffered from management of seminatural woodlands and parks and thinning of forests, but especially for forests the results depended on fungal group and time since partial harvesting. We conclude that the positive effects of tree diversity on overall fungal richness represent a combined niche effect of soil properties and intimate associations.

European ash (Fraxinus excelsior) is threatened by the invasive ascomycete Hymenoscyphus fraxineus originating from Asia. Ash leaf tissues serve as a route for shoot infection but also as a sporulation substrate for this pathogen. Knowledge of the leaf niche partitioning by indigenous fungi and H. fraxineus is needed to understand the fungal community receptiveness to the invasion. We subjected DNA extracted from unwashed and washed leaflets of healthy and diseased European ash to PacBio sequencing of the fungal ITS1-5.8S-ITS2 rDNA region. Leaflets from co-inhabiting rowan trees (Sorbus aucuparia) served as a reference. The overlap in leaflet mycobiomes between ash and rowan was remarkably high, but unlike in rowan, in ash leaflets the sequence read proportion, and the qPCR-based DNA amount estimates of H. fraxineus increased vigorously towards autumn, concomitant with a significant decline in overall fungal richness. The niche of ash and rowan leaves was dominated by epiphytic propagules (Vishniacozyma yeasts, the dimorphic fungus Aureobasidion pullulans and the dematiaceous hyphomycete Cladosporium ramotenellum and H. fraxineus), and endophytic thalli of biotrophs (Phyllactinia and Taphrina species), the indigenous necrotroph Venturia fraxini and H. fraxineus. Mycobiome comparison between healthy and symptomatic European ash leaflets revealed no significant differences in relative abundance of H. fraxineus, but A. pullulans was more prevalent in symptomatic trees. The impacts of host specificity, spatiotemporal niche partitioning, species carbon utilization profiles and life cycle traits are discussed to understand the ecological success of H. fraxineus in Europe. Further, the inherent limitations of different experimental approaches in the profiling of foliicolous fungi are addressed.

Tree species identity is one of the key factors driving ectomycorrhizal (EcM) fungal richness and community composition in boreal and temperate forest ecosystems, but little is known about the influence of tree species combinations and their neighborhood effects on EcM communities. To advance our understanding of host plant effects on EcM fungi, the roots of silver birch, Scots pine, and Norway spruce were analyzed using high-throughput sequencing across mature boreal forest exploratory plots of monocultures and two- and three-species mixtures in Finland. Our analyses revealed that tree species identity was an important determinant of EcM fungal community composition, but tree species richness had no significant influence on EcM fungal richness and community composition. We found that EcM fungal community composition associated with spruce depends on neighboring tree species. Our study suggests that at a regional-scale tree species identity is the primary factor determining community composition of root-associated EcM fungi alongside with tree species composition effects on EcM fungal community of spruce in mixed stands.

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.

The relationship between the ecological success of needle pathogens of forest trees and species richness of co-inhabiting endophytic fungi is poorly understood. One of the most dangerous foliar pathogens of pine is Dothistroma septosporum, which is a widely spread threat to northern European forests. We sampled two Pinus sylvestris sites in Estonia and two in Norway in order to analyse the relations between the abundance of D. septosporum and overall fungal richness, specific fungal species composition, time of season, needle age and position in the canopy. In both countries, the overall species richness of fungi was highest in autumn, showing a trend of increase with needle age. The overall species richness in the second-year needles in Estonia and third-year needles in Norway was similar, suggesting that a critical colonization threshold for needle shed in P. sylvestris is breached earlier in Estonia than in Norway. The fungal species richness in P. sylvestris needles was largely affected by Lophodermium conigenum. Especially in older needles, the relative abundance of L. conigenum was significantly higher in spring compared to summer or autumn. The timing of recruitment and colonization mechanisms of different foliage endophytes are shortly discussed.

Increasing the number of environmental stressors could decrease ecosystem functioning in soils. Yet this relationship has never been globally assessed outside laboratory experiments. Here, using two independent global standardized field surveys, and a range of natural and human factors, we test the relationship between the number of environmental stressors exceeding different critical thresholds and the maintenance of multiple ecosystem services across biomes. Our analysis shows that, multiple stressors, from medium levels (>50%), negatively and significantly correlates with impacts on ecosystem services, and that multiple stressors crossing a high-level critical threshold (over 75% of maximum observed levels), reduces soil biodiversity and functioning globally. The number of environmental stressors >75% threshold was consistently seen as an important predictor of multiple ecosystem services, therefore improving prediction of ecosystem functioning. Our findings highlight the need to reduce the dimensionality of the human footprint on ecosystems to conserve biodiversity and function.

DNA phylogenetic comparisons have shown that morphology-based species recognition often underestimates fungal diversity. Therefore, the need for accurate DNA sequence data, tied to both correct taxonomic names and clearly annotated specimen data, has never been greater. Furthermore, the growing number of molecular ecology and microbiome projects using high-throughput sequencing require fast and effective methods for en masse species assignments. In this article, we focus on selecting and re-annotating a set of marker reference sequences that represent each currently accepted order of Fungi. The particular focus is on sequences from the internal transcribed spacer region in the nuclear ribosomal cistron, derived from type specimens and/or ex-type cultures. Re-annotated and verified sequences were deposited in a curated public database at the National Center for Biotechnology Information (NCBI), namely the RefSeq Targeted Loci (RTL) database, and will be visible during routine sequence similarity searches with NR_prefixed accession numbers. A set of standards and protocols is proposed to improve the data quality of new sequences, and we suggest how type and other reference sequences can be used to improve identification of Fungi. Database URL: http://www.ncbi.nlm.nih.gov/bioproject/PRJNA177353.

Truffles have evolved from epigeous (aboveground) ancestors in nearly every major lineage of fleshy fungi. Because accelerated rates of morphological evolution accompany the transition to the truffle form, closely related epigeous ancestors remain unknown for most truffle lineages. This is the case for the quintessential truffle genus Tuber, which includes species with socio-economic importance and esteemed culinary attributes. Ecologically, Tuber spp. form obligate mycorrhizal symbioses with diverse species of plant hosts including pines, oaks, poplars, orchids, and commercially important trees such as hazelnut and pecan. Unfortunately, limited geographic sampling and inconclusive phylogenetic relationships have obscured our understanding of their origin, biogeography, and diversification. To address this problem, we present a global sampling of Tuberaceae based on DNA sequence data from four loci for phylogenetic inference and molecular dating. Our well-resolved Tuberaceae phylogeny shows high levels of regional and continental endemism. We also identify a previously unknown epigeous member of the Tuberaceae--the South American cup-fungus Nothojafnea thaxteri (E.K. Cash) Gamundí. Phylogenetic resolution was further improved through the inclusion of a previously unrecognized Southern hemisphere sister group of the Tuberaceae. This morphologically diverse assemblage of species includes truffle (e.g. Gymnohydnotrya spp.) and non-truffle forms that are endemic to Australia and South America. Southern hemisphere taxa appear to have diverged more recently than the Northern hemisphere lineages. Our analysis of the Tuberaceae suggests that Tuber evolved from an epigeous ancestor. Molecular dating estimates Tuberaceae divergence in the late Jurassic (~156 million years ago), with subsequent radiations in the Cretaceous and Paleogene. Intra-continental diversification, limited long-distance dispersal, and ecological adaptations help to explain patterns of truffle evolution and biodiversity.

Vegetation impacts on ecosystem functioning are mediated by mycorrhizas, plant-fungal associations formed by most plant species. Ecosystems dominated by distinct mycorrhizal types differ strongly in their biogeochemistry. Quantitative analyses of mycorrhizal impacts on ecosystem functioning are hindered by the scarcity of information on mycorrhizal distributions. Here we present global, high-resolution maps of vegetation biomass distribution by dominant mycorrhizal associations. Arbuscular, ectomycorrhizal, and ericoid mycorrhizal vegetation store, respectively, 241 ± 15, 100 ± 17, and 7 ± 1.8 GT carbon in aboveground biomass, whereas non-mycorrhizal vegetation stores 29 ± 5.5 GT carbon. Soil carbon stocks in both topsoil and subsoil are positively related to the community-level biomass fraction of ectomycorrhizal plants, though the strength of this relationship varies across biomes. We show that human-induced transformations of Earth's ecosystems have reduced ectomycorrhizal vegetation, with potential ramifications to terrestrial carbon stocks. Our work provides a benchmark for spatially explicit and globally quantitative assessments of mycorrhizal impacts on ecosystem functioning and biogeochemical cycling.

The structure and function of the soil microbiome of urban greenspaces remain largely undetermined. We conducted a global field survey in urban greenspaces and neighboring natural ecosystems across 56 cities from six continents, and found that urban soils are important hotspots for soil bacterial, protist and functional gene diversity, but support highly homogenized microbial communities worldwide. Urban greenspaces had a greater proportion of fast-growing bacteria, algae, amoebae, and fungal pathogens, but a lower proportion of ectomycorrhizal fungi than natural ecosystems. These urban ecosystems also showed higher proportions of genes associated with human pathogens, greenhouse gas emissions, faster nutrient cycling, and more intense abiotic stress than natural environments. City affluence, management practices, and climate were fundamental drivers of urban soil communities. Our work paves the way toward a more comprehensive global-scale perspective on urban greenspaces, which is integral to managing the health of these ecosystems and the well-being of human populations.

Fungi are a diverse eukaryotic group of degraders, pathogens, and symbionts, with many lineages known only from DNA sequences in soil, sediments, air, and water.

With the developments in DNA nanoball sequencing technologies and the emergence of new platforms, there is an increasing interest in their performance in comparison with the widely used sequencing-by-synthesis methods. Here, we test the consistency of metabarcoding results from DNBSEQ-G400RS (DNA nanoball sequencing platform by MGI-Tech) and NovaSeq 6000 (sequencing-by-synthesis platform by Illumina) platforms using technical replicates of DNA libraries that consist of COI gene amplicons from 120 soil DNA samples. By subjecting raw sequencing data from both platforms to a uniform bioinformatics processing, we found that the proportion of high-quality reads passing through the filtering steps was similar in both datasets. Per-sample operational taxonomic unit (OTU) and amplicon sequence variant (ASV) richness patterns were highly correlated, but sequencing data from DNBSEQ-G400RS harbored a higher number of OTUs. This may be related to the lower dominance of most common OTUs in DNBSEQ data set (thus revealing higher richness by detecting rare taxa) and/or to a lower effective read quality leading to generation of spurious OTUs. However, there was no statistical difference in the ASV and post-clustered ASV richness between platforms, suggesting that additional denoising step in the ASV workflow had effectively removed the 'noisy' reads. Both OTU-based and ASV-based composition were strongly correlated between the sequencing platforms, with essentially interchangeable results. Therefore, we conclude that DNBSEQ-G400RS and NovaSeq 6000 are both equally efficient high-throughput sequencing platforms to be utilized in studies aiming to apply the metabarcoding approach, but the main benefit of the former is related to lower sequencing cost.

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.

Novel high-throughput sequencing methods outperform earlier approaches in terms of resolution and magnitude. They enable identification and relative quantification of community members and offer new insights into fungal community ecology. These methods are currently taking over as the primary tool to assess fungal communities of plant-associated endophytes, pathogens, and mycorrhizal symbionts, as well as free-living saprotrophs. Taking advantage of the collective experience of six research groups, we here review the different stages involved in fungal community analysis, from field sampling via laboratory procedures to bioinformatics and data interpretation. We discuss potential pitfalls, alternatives, and solutions. Highlighted topics are challenges involved in: obtaining representative DNA/RNA samples and replicates that encompass the targeted variation in community composition, selection of marker regions and primers, options for amplification and multiplexing, handling of sequencing errors, and taxonomic identification. Without awareness of methodological biases, limitations of markers, and bioinformatics challenges, large-scale sequencing projects risk yielding artificial results and misleading conclusions.

Coarse woody debris (CWD) provides food and shelter to a large proportion of forest biota and is considered vital for biodiversity during periods of harsh weather. However, its importance in long-term stressed ecosystems remains largely unknown. In this work, we explored the contribution of CWD to fungal diversity along the gradient of boreal forest degradation caused by 77 years of heavy industrial emissions. We analyzed the diversity and composition of fungi in 270 samples of well-decayed Picea abies and Abies sibirica logs, as well as forest litter both adjacent to and distant from the logs. Compared with forest litter, the wood had higher water content and possessed substantially lower concentrations of heavy metals, which suggests its potential favorability for biota in polluted areas. The pollution-induced loss of fungal diversity in forest litter reached 34% and was stronger in the microhabitats not influenced by CWD. Meanwhile, wood fungal communities lost less than 10% of their total richness and even increased in alpha diversity. These processes led to the diversity and compositional convergence of fungal communities from different microhabitats and substrates in polluted areas. Despite this, the importance of wood and CWD-influenced microhabitats for fungal diversity maintenance was low. Apart from wood-associated fungi, the taxa whose diversity increased in the wood of polluted areas were ectomycorrhizal fungi and eurytopic soil saprotrophs (Mucoromycota, Mortierellomycota, Eurotiomycetes, and Helotiales) that easily tolerate highly toxic litter. Within the majority of pollution-sensitive soil saprotrophic groups, only terricolous Tricholomataceae benefit from CWD as microrefugia. Upon considering the ecological variability within low-rank taxa, the importance of decayed logs as safe sites can be high for certain soil-inhabiting fungal groups in polluted areas.

Fungi are ecologically important in several ecosystem processes, yet their community composition, ecophysiological roles, and responses to changing environmental factors in historical sediments are rarely studied. Here we explored ancient fungal DNA from lake Lielais Svētiņu sediment throughout the Holocene (10.5 kyr) using the ITS metabarcoding approach. Our data revealed diverse fungal taxa and smooth community changes during most of the Holocene with rapid changes occurring in the last few millennia. More precisely, plankton parasitic fungi became more diverse from the Late Holocene (2-4 kyr) which could be related to a shift towards a cooler climate. The Latest Holocene (~2 kyr) showed a distinct increase in the richness of plankton parasites, mycorrhizal, and plant pathogenic fungi which can be associated with an increased transfer rate of plant material into the lake and blooms of planktonic organisms influenced by increased, yet moderate, human impact. Thus, major community shifts in plankton parasites and mycorrhizal fungi could be utilized as potential paleo-variables that accompany host-substrate dynamics. Our work demonstrates that fungal aDNA with predicted ecophysiology and host specificity can be employed to reconstruct both aquatic and surrounding terrestrial ecosystems and to estimate the influence of environmental change.

Mycorrhizal fungi are mutualists that play crucial roles in nutrient acquisition in terrestrial ecosystems. Mycorrhizal symbioses arose repeatedly across multiple lineages of Mucoromycotina, Ascomycota, and Basidiomycota. Considerable variation exists in the capacity of mycorrhizal fungi to acquire carbon from soil organic matter. Here, we present a combined analysis of 135 fungal genomes from 73 saprotrophic, endophytic and pathogenic species, and 62 mycorrhizal species, including 29 new mycorrhizal genomes. This study samples ecologically dominant fungal guilds for which there were previously no symbiotic genomes available, including ectomycorrhizal Russulales, Thelephorales and Cantharellales. Our analyses show that transitions from saprotrophy to symbiosis involve (1) widespread losses of degrading enzymes acting on lignin and cellulose, (2) co-option of genes present in saprotrophic ancestors to fulfill new symbiotic functions, (3) diversification of novel, lineage-specific symbiosis-induced genes, (4) proliferation of transposable elements and (5) divergent genetic innovations underlying the convergent origins of the ectomycorrhizal guild.

Here, we describe the taxon hypothesis (TH) paradigm, which covers the construction, identification, and communication of taxa as datasets. Defining taxa as datasets of individuals and their traits will make taxon identification and most importantly communication of taxa precise and reproducible. This will allow datasets with standardized and atomized traits to be used digitally in identification pipelines and communicated through persistent identifiers. Such datasets are particularly useful in the context of formally undescribed or even physically undiscovered species if data such as sequences from samples of environmental DNA (eDNA) are available. Implementing the TH paradigm will to some extent remove the impediment to hastily discover and formally describe all extant species in that the TH paradigm allows discovery and communication of new species and other taxa also in the absence of formal descriptions. The TH datasets can be connected to a taxonomic backbone providing access to the vast information associated with the tree of life. In parallel to the description of the TH paradigm, we demonstrate how it is implemented in the UNITE digital taxon communication system. UNITE TH datasets include rich data on individuals and their rDNA ITS sequences. These datasets are equipped with digital object identifiers (DOI) that serve to fix their identity in our communication. All datasets are also connected to a GBIF taxonomic backbone. Researchers processing their eDNA samples using UNITE datasets will, thus, be able to publish their findings as taxon occurrences in the GBIF data portal. UNITE species hypothesis (species level THs) datasets are increasingly utilized in taxon identification pipelines and even formally undescribed species can be identified and communicated by using UNITE. The TH paradigm seeks to achieve unambiguous, unique, and traceable communication of taxa and their properties at any level of the tree of life. It offers a rapid way to discover and communicate undescribed species in identification pipelines and data portals before they are lost to the sixth mass extinction.