The contribution of DNA metabarcoding to fungal conservation: diversity assessment, habitat partitioning and mapping red-listed fungi in protected coastal Salix repens communities in the Netherlands

Multiple graphical views for automatically generating SQL for the MycoDiversity DB; making fungal biodiversity studies accessible

Fungi is a highly diverse group of eukaryotic organisms that live under an extremely wide range of environmental conditions. Nowadays, there is a fundamental focus on observing how biodiversity varies on different spatial scales, in addition to understanding the environmental factors which drive fungal biodiversity. Metabarcoding is a high-throughput DNA sequencing technology that has positively contributed to observing fungal communities in environments. While the DNA sequencing data generated from metabarcoding studies are available in public archives, this valuable data resource is not directly usable for fungal biodiversity investigation. Additionally, due to its fragmented storage and distributed nature, it is not immediately accessible through a single user interface. We developed the MycoDiversity DataBase User Interface (https://mycodiversity.liacs.nl) to provide direct access and retrieval of fungal data that was previously inaccessible in the public domain. The user interface provides multiple graphical views of the data components used to reveal fungal biodiversity. These components include reliable geo-location terms, the reference taxonomic scientific names associated with fungal species and the standard features describing the environment where they occur. Direct observation of the public DNA sequencing data in association with fungi is accessible through SQL search queries created by interactively manipulating topological maps and dynamic hierarchical tree views. The search results are presented in configurable data table views that can be downloaded for further use. With the MycoDiversity DataBase User Interface, we make fungal biodiversity data accessible, assisting researchers and other stakeholders in using metabarcoding studies for assessing fungal biodiversity.

A case study on the application of spore sampling for the monitoring of macrofungi

Fungi play a vital role in ecosystem functioning, yet significant knowledge gaps persist in understanding their diversity and distribution leading to uncertainties about their threat status and extinction risk. This is partly owed to the difficulty of monitoring fungi using traditional fruiting body surveys. The present study evaluates airborne environmental DNA (eDNA) sampling as a monitoring tool with a focus on grassland macrofungi. We applied active and passive air sampling methods, complemented by extensive field surveys of waxcap and clavarioid fungi-species groups of high relevance for conservation. Twenty-nine species were recorded during the field surveys, 19 of which were also detectable by ITS2 metabarcoding of the collected samples. An additional 12 species from the studied genera were identified exclusively in air eDNA. We found that the patterns of species detection and read abundance in air samples reflected the abundance and occurrence of fruiting bodies on the field. Dispersal kernels fitted for the three dominant species predicted rapidly decreasing spore concentrations with increasing distance from fruitbodies. Airborne assemblages were dominated by a high diversity of common species, while rare and threatened red-listed species were under-represented, which underscores the difficulty in detecting rare species, not only in conventional surveys. Considering the benefits and drawbacks of air sampling and fruitbody surveys, we conclude that air sampling serves as a cost- and time-efficient tool to characterize local macrofungal communities, providing the potential to facilitate and improve future fungal monitoring efforts.

Comparison of survey methods for fungi using metabarcoding and fruit body inventories in an altitudinal gradient

Metabarcoding of environmental samples is nowadays an established method in biodiversity research. When it comes to studying fungal populations in various ecotypes, fruit body inventories are the traditional method to assess the diversity of fungal communities. In this study, both methods-metabarcoding of soil samples and a traditional fruit body inventory-were conducted on 144 sample plots in an altitudinal gradient in the Bavarian Forest (Germany) and the results were compared. Metabarcoding detected significantly more species than the traditional fruit body inventory. The majority of taxa recorded in the fruit body inventory belonged to the Basidiomycota, whereas in the metabarcoding data, the distribution of species between Basidiomycota and Ascomycota was approximately balanced. Species of several orders forming inconspicuous or hypogeous fruit bodies were detected only by metabarcoding, while several wood decomposers were recorded only in the fruit body inventory. The proportion of detected wood-colonising species with melanized spores was considerably higher with metabarcoding than with the fruit body inventory, where more than 70% of recorded wood-colonisers had hyaline spores. Based on the metabarcoding data, a decline of species richness with increasing altitude was evident, but this was not visible in the fruit body inventory data. Detrended correspondence analyses yielded similar results for relative species community similarities with both survey methods.

Changing Rhizosphere Microbial Community and Metabolites with Developmental Stages of Coleus barbatus

Coleus barbatus is a medicinal herb belonging to Lamiaceae. It is the only living organism known to produce forskolin, which is a labdane diterpene and is reported to activate adenylate cyclase. Microbes associated with plants play an important role in maintaining plant health. Recently, the targeted application of beneficial plant-associated microbes and their combinations in abiotic and biotic stress tolerance has gained momentum. In this work, we carried out the rhizosphere metagenome sequencing of C. barbatus at different developmental stages to understand how rhizosphere microflora are affected by and affect the metabolite content in plants. We found that the Kaistobacter genus was abundantly present in the rhizosphere of C. barbatus and its accumulation pattern appears to correlate with the quantities of forskolin in the roots at different developmental stages. Members of the Phoma genus, known for several pathogenic species, were in lower numbers in the C. barbatus rhizosphere in comparison with C. blumei. To our knowledge, this is the first metagenomic study of the rhizospheric microbiome of C. barbatus, which may help to explore and exploit the culturable and non-culturable microbial diversity present in the rhizosphere.

Abiotic environmental factors drive the diversity, compositional dynamics and habitat preference of ectomycorrhizal fungi in Pannonian forest types

Ectomycorrhizal (ECM) fungi are among the most diverse and dominant fungal groups in temperate forests and are crucial for ecosystem functioning of forests and their resilience toward disturbance. We carried out DNA metabarcoding of ECM fungi from soil samples taken at 62 sites in the Bükk Mountains in northern Hungary. The selected sampling sites represent the characteristic Pannonian forest types distributed along elevation (i.e., temperature), pH and slope aspect gradients. We compared richness and community composition of ECM fungi among forest types and explored relationships among environmental variables and ECM fungal alpha and beta diversity. The DNA sequence data generated in this study indicated strong correlations between fungal community composition and environmental variables, particularly with pH and soil moisture, with many ECM fungi showing preference for specific zonal, topographic or edaphic forest types. Several ECM fungal genera showed significant differences in richness among forest types and exhibited strong compositional differences mostly driven by differences in environmental factors. Despite the relatively high proportions of compositional variance explained by the tested environmental variables, a large proportion of the compositional variance remained unexplained, indicating that both niche (environmental filtering) and neutral (stochastic) processes shape ECM fungal community composition at landscape level. Our work provides unprecedented insights into the diversity, landscape-level distribution, and habitat preferences of ECM fungi in the Pannonian forests of Northern Hungary.

Habitat Protection Approaches Facilitate Conservation of Overlooked Fungal Diversity - A Case Study From the Norwegian Coastal Heathland System

European coastal heathlands are distinct ecosystems shaped by land use tradition and they have experienced an 80% area reduction from their historical maximum. These mosaics of mires and wind exposed patches have ericaceous shrub dominated vegetation, and soils within coastal heathlands are characterized by low pH and high levels of recalcitrant debris. Using a culture-based approach with molecular identification of isolates, we characterized root-associated fungal communities of six ericaceous species in eight heathland localities along Norway's western coast. Site-level alpha diversity ranged from 21-38 OTUs, while the total estimated gamma diversity for culturable heathland root fungi was 190-231 OTUs. Most species recovered are previously reported at low abundance in Norway, suggesting the biodiversity in this community is underreported, rather than novel for science. The fungi recovered were primarily Ascomycota, specifically endophytic Phialocephala, and Pezicula, and no host specificity was observed in the communities. The fungal communities exhibited high turnover and low nestedness, both between ericaceous hosts and across heathland sites. We observed no spatial patterns in fungal betadiversity, and this heterogeneity may be a product of the unique historic land use practices at each locality creating a distinct mycofloral "fingerprint". Robust diversity estimates will be key for managing fungal biodiversity in coastal heathlands. Our results indicate that sampling schemes that maximize the number of host plants sampled per site, rather than the number of cultures per plant yield improved alpha diversity estimates. Similarly, gamma diversity estimates are improved by maximizing the total number of localities sampled, rather than increasing the number of plants sampled per locality. We argue that while the current protected status of coastal heathland habitats and restoration efforts have knock-on effects for the conservation of fungal biodiversity, fungi have a vital functional role in the ecosystem and holistic conservation plans that consider fungal biodiversity would be beneficial.

Molecular-Based Diversity Studies and Field Surveys Are Not Mutually Exclusive: On the Importance of Integrated Methodologies in Mycological Research

Understanding and describing the diversity of living organisms is a great challenge. Fungi have for a long time been, and unfortunately still are, underestimated when it comes to taxonomic research. The foundations were laid by the first mycologists through field observations. These important fundamental works have been and remain vital reference works. Nevertheless, a non-negligible part of the studied funga escaped their attention. Thanks to modern developments in molecular techniques, the study of fungal diversity has been revolutionized in terms of tools and knowledge. Despite a number of disadvantages inherent to these techniques, traditional field-based inventory work has been increasingly superseded and neglected. This perspective aims to demonstrate the central importance of field-based research in fungal diversity studies, and encourages researchers not to be blinded by the sole use of molecular methods.

Current Insight into Traditional and Modern Methods in Fungal Diversity Estimates

Fungi are an important and diverse component in various ecosystems. The methods to identify different fungi are an important step in any mycological study. Classical methods of fungal identification, which rely mainly on morphological characteristics and modern use of DNA based molecular techniques, have proven to be very helpful to explore their taxonomic identity. In the present compilation, we provide detailed information on estimates of fungi provided by different mycologistsover time. Along with this, a comprehensive analysis of the importance of classical and molecular methods is also presented. In orderto understand the utility of genus and species specific markers in fungal identification, a polyphasic approach to investigate various fungi is also presented in this paper. An account of the study of various fungi based on culture-based and cultureindependent methods is also provided here to understand the development and significance of both approaches. The available information on classical and modern methods compiled in this study revealed that the DNA based molecular studies are still scant, and more studies are required to achieve the accurate estimation of fungi present on earth.

A long-read amplicon approach to scaling up the metabarcoding of lichen herbarium specimens

Reference sequence databases are critical to the accurate detection and identification of fungi in the environment. As repositories of large numbers of well-curated specimens, herbaria and fungal culture collections have the material resources to generate sequence data for large number of taxa, and could therefore allow filling taxonomic gaps often present in reference sequence databases. Financial resources to do that are however often lacking, so that recent efforts have focused on decreasing sequencing cost by increasing the number of multiplexed samples per sequencing run while maintaining high sequence quality. Following a previous study that aimed at decreasing sequencing cost for lichen specimens by generating fungal ITS barcodes for 96 specimens using PacBio amplicon sequencing, we present a method that further decreases lichen specimen metabarcoding costs. A total of 384 mixed DNA extracts obtained from lichen herbarium specimens, mostly from the four genera Buellia, Catillaria, Endocarpon and Parmotrema, were used to generate new fungal ITS sequences using a Sequel I sequencing platform and the PacBio M13 barcoded primers. The average success rate across all taxa was high (86.5%), with particularly high rates for the crustose saxicolous taxa (Buellia, Catillaria and others; 93.3%) and the terricolous squamulose taxa (Endocarpon and others; 96.5%). On the other hand, the success rate for the foliose genus Parmotrema was lower (60.4%). With this taxon sampling, greater specimen age did not appear to impact sequencing success. In fact, the 1966–1980 collection date category showed the highest success rate (97.3%). Compared to the previous study, the abundance-based sequence denoising method showed some limitations, but the cost of generating ITS barcodes was further decreased thanks to the higher multiplexing level. In addition to contributing new ITS barcodes for specimens of four interesting lichen genera, this study further highlights the potential and challenges of using new sequencing technologies on collection specimens to generate DNA sequences for reference databases.

Soil Fungal Community Composition and Diversity of Culturable Endophytic Fungi from Plant Roots in the Reclaimed Area of the Eastern Coast of China

As an important resource for screening microbial strains capable of conferring stress tolerance in plants, the fungal community associated with the plants grown in stressful environments has received great attention. In this study, high-throughput sequencing was employed to study the rhizosphere fungal community in the reclaimed area (i.e., sites F, H, and T) of the eastern coast of China. Moreover, endophytic fungi from the root of six plant species colonizing the investigated sites were isolated and identified. The differences in soil physicochemical parameters, fungal diversity, and community structure were detected among the sampling sites and between the seasons. Ectomycorrhizal (ECM) fungi (e.g., genera Tuber and Geopora) were dominant at site F, which was characterized by high soil total carbon (SC) and total nitrogen (SN) contents and low soil electrical conductivity (EC) value. Arbuscular mycorrhizal (AM) fungi, including genera Glomus, Rhizophagus, and Entrophospora were dominant at sites H (winter), H (summer), and T (summer), respectively. The positive relationship between the EC value and the abundance of genus Glomus indicated the ability of this AM fungus to protect plants against the salt stress. Endophytic fungi at sites F (Aspergillus and Tetracladium), H (Nigrospora), and T (Nigrospora, Coniochaeta and Zopfiella) were recognized as the biomarkers or keystone taxa, among which only genus Aspergillus was isolated from the plant roots. The aforementioned AM fungi and endophytic fungi could contribute to the promotion of plant growth in the newly reclaimed land.

Similarities and Differences among Soil Fungal Assemblages in Managed Forests and Formerly Managed Forest Reserves

Unlike the numerous works concerning the effect of management on the forest mycobiome, only a few studies have addressed how fungi from different trophic groups recover from natural and anthropogenic disturbances and develop structural features typical of unmanaged old-growth forests. Our objective is to compare the soil fungal assemblages represented by different functional/trophic groups in protected and managed stands located in European mixed forests dominated by Scots pine. Fungal communities were analyzed using high-throughput Illumina MiSeq sequencing of fungal internal transcribed spacer 1 (ITS1) amplicons. Formerly managed forest reserves (established around 50 years ago) and forests under standard forest management appeared to be similar in terms of total and mean species richness of all fungal operational taxonomic units (OTUs), as well as OTUs assigned to different functional trophic groups. Among the 599 recorded OTUs, 497 (83%) were shared between both management types, whereas 9.5% of taxa were unique to forest reserves and 7.5% were unique to managed stands. Ascomycota and Basidiomycota were the predominant phyla, comprising 88% of all identified fungi. The main functional components of soil fungal assemblages consisted of saprotrophic (42% fungal OTUs; 27% reads) and ectomycorrhizal fungi (16%; 47%). Two-way analysis of similarities (ANOSIM) revealed that both site and management strategy influenced the species composition of soil fungal communities, with site being a primary effect for saprotrophic and ectomycorrhizal fungi. Volume of coarse and very fine woody debris and soil pH significantly influenced the ectomycorrhizal fungal community, whereas saprotrophic fungi were influenced primarily by volume of coarse woody debris and soil nitrate concentration. Among the identified fungal OTUs, 18 red-listed fungal species were identified from both forest reserves and managed forests, comprising two ECM fungi and four saprotrophs from the category of endangered species. Our results suggest that the transformation of fungal diversity after cessation of forest management is rather slow, and that both forest reserves and managed forests help uphold fungal diversity.

Sareomycetes: more diverse than meets the eye

Since its resurrection, the resinicolous discomycete genus Sarea has been accepted as containing two species, one with black apothecia and pycnidia, and one with orange. We investigate this hypothesis using three ribosomal (nuITS, nuLSU, mtSSU) regions from and morphological examination of 70 specimens collected primarily in Europe and North America. The results of our analyses support separation of the traditional Sarea difformis s.lat. and Sarea resinae s.lat. into two distinct genera, Sarea and Zythia. Sarea as circumscribed is shown to conservatively comprise three phylospecies, with one corresponding to Sarea difformis s.str. and two, morphologically indistinguishable, corresponding to the newly combined Sarea coeloplata. Zythia is provisionally maintained as monotypic, containing only a genetically and morphologically variable Z. resinae. The new genus Atrozythia is erected for the new species A. klamathica. Arthrographis lignicola is placed in this genus on molecular grounds, expanding the concept of Sareomycetes by inclusion of a previously unknown type of asexual morph. Dating analyses using additional marker regions indicate the emergence of the Sareomycetes was roughly concurrent with the diversification of the genus Pinus, suggesting that this group of fungi emerged to exploit the newly-available resinous ecological niche supplied by Pinus or another, extinct group of conifers. Our phylogeographic studies also permitted us to study the introductions of these fungi to areas where they are not native, including Antarctica, Cape Verde, and New Zealand and are consistent with historical hypotheses of introduction.

A global meta-analysis of ITS rDNA sequences from material belonging to the genus Ganoderma (Basidiomycota, Polyporales) including new data from selected taxa

Ganoderma P. Karst. is a cosmopolitan genus of white-rot fungi which comprises species with highly-prized pharmaceutical properties, valuable biotechnological applications and of significant phytopathological interest. However, the status of the taxonomy within the genus is still highly controversial and ambiguous despite the progress made through molecular approaches. A metadata analysis of 3908 nuclear ribosomal internal transcribed spacer (ITS) rDNA sequences obtained from GenBank/ENA/DDBJ and UNITE was performed by targeting sequences annotated as Ganoderma, but also sequences from environmental samples and from material examined for the first time. Ganoderma taxa segregated into five main lineages (Clades A to E). Clade A corresponds to the core of laccate species and includes G. shanxiense and three major well-supported clusters: Cluster A.1 ('G. lucidum sensu lato') consists of taxa from Eurasia and North America, Cluster A.2 of material with worldwide occurrence including G. resinaceum and Cluster A.3 is composed of species originating from all continents except Europe and comprises G. lingzhi. Clade B includes G. applanatum and allied species with a Holarctic distribution. Clade C comprises taxa from Asia and Africa only. Clade D consists of laccate taxa with tropical/subtropical occurrence, while clade E harbours the highest number of non-laccate species with a cosmopolitan distribution. The 92 Ganoderma-associated names, initially used for sequences labelling, correspond to at least 80 taxa. Amongst them, 21 constitute putatively new phylospecies after our application of criteria relevant to the robustness/support of the terminal clades, intra- and interspecific genetic divergence and available biogeographic data. Moreover, several other groups or individual sequences seem to represent distinct taxonomic entities and merit further investigation. A particularly large number of the public sequences was revealed to be insufficiently and/or incorrectly identified, for example, 87% and 78% of entries labelled as G. australe and G. lucidum, respectively. In general, ITS demonstrated high efficacy in resolving relationships amongst most of the Ganoderma taxa; however, it was not equally useful at elucidating species barriers across the entire genus and such cases are outlined. Furthermore, we draw conclusions on biogeography by evaluating species occurrence on a global scale in conjunction with phylogenetic structure/patterns. The sequence variability assessed in ITS spacers could be further exploited for diagnostic purposes.

Soil stabilisation for DNA metabarcoding of plants and fungi. Implications for sampling at remote locations or via third-parties

Storage of soil samples prior to metagenomic analysis presents a problem. If field sites are remote or if samples are collected by third parties, transport to analytical laboratories may take several days or even weeks. The bulk of such samples and requirement for later homogenisation precludes the convenient use of a stabilisation buffer, so samples are usually cooled or frozen during transit. There has been limited testing of the most appropriate storage methods for later study of soil organisms by eDNA approaches. Here we tested a range of storage methods on two contrasting soils, comparing these methods to the control of freezing at -80 °C, followed by freeze-drying. To our knowledge, this is the first study to examine the effect of storage conditions on eukaryote DNA in soil, including both viable organisms (fungi) and DNA contained within dying/dead tissues (plants). For fungi, the best storage regimes (closest to the control) were storage at 4 °C (for up to 14 d) or active air-drying at room temperature. The worst treatments involved initial freezing, followed by thawing which led to significant later spoilage. The key spoilage organisms were identified as Metarhizium carneum and Mortierella spp., with a general increase in saprotrophic fungi and reduced abundances of mycorrhizal/biotrophic fungi. Plant data showed a similar pattern, but with greater variability in community structure, especially in the freeze-thaw treatments, probably due to stochastic variation in substrates for fungal decomposition, algal proliferation and some seed germination. In the absence of freeze drying facilities, samples should be shipped refrigerated, but not frozen if there is any risk of thawing.

Assessment of fungal diversity in soil rhizosphere associated with Rhazya stricta and some desert plants using metagenomics

This study aimed to compare the fungal rhizosphere communities of Rhazya stricta, Enneapogon desvauxii, Citrullus colocynthis, Senna italica, and Zygophyllum simplex, and the gut mycobiota of Poekilocerus bufonius (Orthoptera, Pyrgomorphidae, "Usherhopper"). A total of 164,485 fungal reads were observed from the five plant rhizospheres and Usherhopper gut. The highest reads were in S. italica rhizosphere (29,883 reads). Species richness in the P. bufonius gut was the highest among the six samples. Ascomycota was dominant in all samples, with the highest reads in E. desvauxii (26,734 reads) rhizosphere. Sordariomycetes and Dothideomycetes were the dominant classes detected with the highest abundance in C. colocynthis and E. desvauxii rhizospheres. Aspergillus and Ceratobasidium were the most abundant genera in the R. stricta rhizosphere, Fusarium and Penicillium in the E. desvauxii rhizosphere and P. bufonius gut, Ceratobasidium and Myrothecium in the C. colocynthis rhizosphere, Aspergillus and Fusarium in the S. italica rhizosphere, and Cochliobolus in the Z. simplex rhizosphere. Aspergillus terreus was the most abundant species in the R. stricta and S. italica rhizospheres, Fusarium sp. in E. desvauxii rhizosphere, Ceratobasidium sp. in C. colocynthis rhizosphere, Cochliobolus sp. in Z. simplex rhizosphere, and Penicillium sp. in P. bufonius gut. The phylogenetic results revealed the unclassified species were related closely to Ascomycota and the species in E. desvauxii, S. italica and Z. simplex rhizospheres were closely related, where the species in the P. bufonius gut, were closely related to the species in the R. stricta, and C. colocynthis rhizospheres.

Naming the untouchable - environmental sequences and niche partitioning as taxonomical evidence in fungi

Due to their submerged and cryptic lifestyle, the vast majority of fungal species are difficult to observe and describe morphologically, and many remain known to science only from sequences detected in environmental samples. The lack of practices to delimit and name most fungal species is a staggering limitation to communication and interpretation of ecology and evolution in kingdom Fungi. Here, we use environmental sequence data as taxonomical evidence and combine phylogenetic and ecological data to generate and test species hypotheses in the class Archaeorhizomycetes (Taphrinomycotina, Ascomycota). Based on environmental amplicon sequencing from a well-studied Swedish pine forest podzol soil, we generate 68 distinct species hypotheses of Archaeorhizomycetes, of which two correspond to the only described species in the class. Nine of the species hypotheses represent 78% of the sequenced Archaeorhizomycetes community, and are supported by long read data that form the backbone for delimiting species hypothesis based on phylogenetic branch lengths.

Soil fungal communities are shaped by environmental filtering and competitive exclusion so that closely related species are less likely to co-occur in a niche if adaptive traits are evolutionarily conserved. In soil profiles, distinct vertical horizons represent a testable niche dimension, and we found significantly differential distribution across samples for a well-supported pair of sister species hypotheses. Based on the combination of phylogenetic and ecological evidence, we identify two novel species for which we provide molecular diagnostics and propose names. While environmental sequences cannot be automatically translated to species, they can be used to generate phylogenetically distinct species hypotheses that can be further tested using sequences as ecological evidence. We conclude that in the case of abundantly and frequently observed species, environmental sequences can support species recognition in the absences of physical specimens, while rare taxa remain uncaptured at our sampling and sequencing intensity.

GlobalFungi, a global database of fungal occurrences from high-throughput-sequencing metabarcoding studies

Fungi are key players in vital ecosystem services, spanning carbon cycling, decomposition, symbiotic associations with cultivated and wild plants and pathogenicity. The high importance of fungi in ecosystem processes contrasts with the incompleteness of our understanding of the patterns of fungal biogeography and the environmental factors that drive those patterns. To reduce this gap of knowledge, we collected and validated data published on the composition of soil fungal communities in terrestrial environments including soil and plant-associated habitats and made them publicly accessible through a user interface at https://globalfungi.com . The GlobalFungi database contains over 600 million observations of fungal sequences across > 17 000 samples with geographical locations and additional metadata contained in 178 original studies with millions of unique nucleotide sequences (sequence variants) of the fungal internal transcribed spacers (ITS) 1 and 2 representing fungal species and genera. The study represents the most comprehensive atlas of global fungal distribution, and it is framed in such a way that third-party data addition is possible.

Fungal metabarcoding data integration framework for the MycoDiversity DataBase (MDDB)

Fungi have crucial roles in ecosystems, and are important associates for many organisms. They are adapted to a wide variety of habitats, however their global distribution and diversity remains poorly documented. The exponential growth of DNA barcode information retrieved from the environment is assisting considerably the traditional ways for unraveling fungal diversity and detection. The raw DNA data in association to environmental descriptors of metabarcoding studies are made available in public sequence read archives. While this is potentially a valuable source of information for the investigation of Fungi across diverse environmental conditions, the annotation used to describe environment is heterogenous. Moreover, a uniform processing pipeline still needs to be applied to the available raw DNA data. Hence, a comprehensive framework to analyses these data in a large context is still lacking. We introduce the MycoDiversity DataBase, a database which includes public fungal metabarcoding data of environmental samples for the study of biodiversity patterns of Fungi. The framework we propose will contribute to our understanding of fungal biodiversity and aims to become a valuable source for large-scale analyses of patterns in space and time, in addition to assisting evolutionary and ecological research on Fungi.

Notes, outline and divergence times of Basidiomycota

The Basidiomycota constitutes a major phylum of the kingdom Fungi and is second in species numbers to the Ascomycota. The present work provides an overview of all validly published, currently used basidiomycete genera to date in a single document. An outline of all genera of Basidiomycota is provided, which includes 1928 currently used genera names, with 1263 synonyms, which are distributed in 241 families, 68 orders, 18 classes and four subphyla. We provide brief notes for each accepted genus including information on classification, number of accepted species, type species, life mode, habitat, distribution, and sequence information. Furthermore, three phylogenetic analyses with combined LSU, SSU, 5.8s, rpb1, rpb2, and ef1 datasets for the subphyla Agaricomycotina, Pucciniomycotina and Ustilaginomycotina are conducted, respectively. Divergence time estimates are provided to the family level with 632 species from 62 orders, 168 families and 605 genera. Our study indicates that the divergence times of the subphyla in Basidiomycota are 406–430 Mya, classes are 211–383 Mya, and orders are 99–323 Mya, which are largely consistent with previous studies. In this study, all phylogenetically supported families were dated, with the families of Agaricomycotina diverging from 27–178 Mya, Pucciniomycotina from 85–222 Mya, and Ustilaginomycotina from 79–177 Mya. Divergence times as additional criterion in ranking provide additional evidence to resolve taxonomic problems in the Basidiomycota taxonomic system, and also provide a better understanding of their phylogeny and evolution.

Large-scale generation and analysis of filamentous fungal DNA barcodes boosts coverage for kingdom fungi and reveals thresholds for fungal species and higher taxon delimitation

Species identification lies at the heart of biodiversity studies that has in recent years favoured DNA-based approaches. Microbial Biological Resource Centres are a rich source for diverse and high-quality reference materials in microbiology, and yet the strains preserved in these biobanks have been exploited only on a limited scale to generate DNA barcodes. As part of a project funded in the Netherlands to barcode specimens of major national biobanks, sequences of two nuclear ribosomal genetic markers, the Internal Transcribed Spaces and 5.8S gene (ITS) and the D1/D2 domain of the 26S Large Subunit (LSU), were generated as DNA barcode data for ca. 100 000 fungal strains originally assigned to ca. 17 000 species in the CBS fungal biobank maintained at the Westerdijk Fungal Biodiversity Institute, Utrecht. Using more than 24 000 DNA barcode sequences of 12 000 ex-type and manually validated filamentous fungal strains of 7 300 accepted species, the optimal identity thresholds to discriminate filamentous fungal species were predicted as 99.6 % for ITS and 99.8 % for LSU. We showed that 17 % and 18 % of the species could not be discriminated by the ITS and LSU genetic markers, respectively. Among them, ∼8 % were indistinguishable using both genetic markers. ITS has been shown to outperform LSU in filamentous fungal species discrimination with a probability of correct identification of 82 % vs. 77.6 %, and a clustering quality value of 84 % vs. 77.7 %. At higher taxonomic classifications, LSU has been shown to have a better discriminatory power than ITS. With a clustering quality value of 80 %, LSU outperformed ITS in identifying filamentous fungi at the ordinal level. At the generic level, the clustering quality values produced by both genetic markers were low, indicating the necessity for taxonomic revisions at genus level and, likely, for applying more conserved genetic markers or even whole genomes. The taxonomic thresholds predicted for filamentous fungal identification at the genus, family, order and class levels were 94.3 %, 88.5 %, 81.2 % and 80.9 % based on ITS barcodes, and 98.2 %, 96.2 %, 94.7 % and 92.7 % based on LSU barcodes. The DNA barcodes used in this study have been deposited to GenBank and will also be publicly available at the Westerdijk Institute's website as reference sequences for fungal identification, marking an unprecedented data release event in global fungal barcoding efforts to date.

Isotopic evidence of biotrophy and unusual nitrogen nutrition in soil-dwelling Hygrophoraceae

Several lines of evidence suggest that the agaricoid, non-ectomycorrhizal members of the family Hygrophoraceae (waxcaps) are biotrophic with unusual nitrogen nutrition. However, methods for the axenic culture and lab-based study of these organisms remain to be developed, so our current knowledge is limited to field-based investigations. Addition of nitrogen, lime or organophosphate pesticide at an experimental field site (Sourhope) suppressed fruiting of waxcap basidiocarps. Furthermore, stable isotope natural abundance in basidiocarps were unusually high in 15 N and low in 13 C, the latter consistent with mycorrhizal nutritional status. Similar patterns were found in waxcap basidiocarps from diverse habitats across four continents. Additional data from 14 C analysis of basidiocarps and 13 C pulse label experiments suggest that these fungi are not saprotrophs but rather biotrophic endophytes and possibly mycorrhizal. The consistently high but variable δ15 N values (10-20‰) of basidiocarps further indicate that N acquisition or processing differ from other fungi; we suggest that N may be derived from acquisition of N via soil fauna high in the food chain.

Tackling critical parameters in metazoan meta-barcoding experiments: a preliminary study based on coxI DNA barcode

Nowadays DNA meta-barcoding is a powerful instrument capable of quickly discovering the biodiversity of an environmental sample by integrating the DNA barcoding approach with High Throughput Sequencing technologies. It mainly consists of the parallel reading of informative genomic fragment/s able to discriminate living entities. Although this approach has been widely studied, it still needs optimization in some necessary steps requested in its advanced accomplishment. A fundamental element concerns the standardization of bioinformatic analyses pipelines. The aim of the present study was to underline a number of critical parameters of laboratory material preparation and taxonomic assignment pipelines in DNA meta-barcoding experiments using the cytochrome oxidase subunit-I (coxI) barcode region, known as a suitable molecular marker for animal species identification. We compared nine taxonomic assignment pipelines, including a custom in-house method, based on Hidden Markov Models. Moreover, we evaluated the potential influence of universal primers amplification bias in qPCR, as well as the correlation between GC content with taxonomic assignment results. The pipelines were tested on a community of known terrestrial invertebrates collected by pitfall traps from a chestnut forest in Italy. Although the present analysis was not exhaustive and needs additional investigation, our results suggest some potential improvements in laboratory material preparation and the introduction of additional parameters in taxonomic assignment pipelines. These include the correct setup of OTU clustering threshold, the calibration of GC content affecting sequencing quality and taxonomic classification, as well as the evaluation of PCR primers amplification bias on the final biodiversity pattern. Thus, careful attention and further validation/optimization of the above-mentioned variables would be required in a DNA meta-barcoding experimental routine.

Comparative metagenomics approaches to characterize the soil fungal communities of western coastal region, Saudi Arabia

A total of 145007 reads were obtained from pyrosequencing for all the 4 samples. The total count ranged from 11,301,014 (Mecca old road) to 23,503,512 bp (Thuwal). A total of 460 fungal species belonging to 133 genera, 58 families, 33 orders, 13 classes and 4 phyla was identified across the four sites. The most abundant phylum at all four sites was Ascomycota followed by Basidiomycota. Four phyla (Ascomycota—99.31%, Basidiomycota—0.59%, Chytridiomycota—0.04%, Glomeromycota—0.03%) were detected in Khulais. Except for Glomeromycota, all phyla were detected at Mecca old road (Ascomycota—74.26%, Basidiomycota—25.71%, Chytridiomycota—0.01%) and Thuwal (Ascomycota—99.59%, Basidiomycota—0.40%, Chytridiomycota—0.002%); while only Ascomycota—90.98% and Basidiomycota—9.01% were detected in Asfan road. At the class level, Sordariomycetes was predominantly observed at Asfan road—59.88%, Khulais—68.26% and Thuwal—94.84%; while Pezizomycetes was dominant at Mecca old road—56.01%, was absent at Asfan road. Agaricomycetes was present only at Mecca old road—25.73%; while Tremellomycetes—5.77%, Malasseizomycetes—2.13% and Microbotryomycetes—1.10% were found only at Asfan road. The phylogenetic trees revealed that clear genus level differences are visible across all the four sites, with an overall predominance of Thielavia followed by Madurella, Aspergillus, and Gelasinospora. Chaetomium sp., Aspergillus caespitosus and Aspergillus sp. were found in moderate (Mecca old road and Thuwal) to abundant (Asfan road and Khulais) quantities. Thielavia sp., Thielavia hyalocarpa and Madurella sp. are found in moderate quantities at Khulais and Mecca old road, while in abundant levels at Asfan road and Thuwal. Fusarium equisati and F. oxysporum were detected at Thuwal and Khulais. Sordaria araneosa was present at Khulais, while Malasseiza globosa species was detected in moderate quantities across all sites except Khulais.

Mycorrhizal Associations and Trophic Modes in Coexisting Orchids: An Ecological Continuum between Auto- and Mixotrophy

Two distinct nutritional syndromes have been described in temperate green orchids. Most orchids form mycorrhizas with rhizoctonia fungi and are considered autotrophic. Some orchids, however, associate with fungi that simultaneously form ectomycorrhizas with surrounding trees and derive their carbon from these fungi. This evolutionarily derived condition has been called mixotrophy or partial mycoheterotrophy and is characterized by 13C enrichment and high N content. Although it has been suggested that the two major nutritional syndromes are clearly distinct and tightly linked to the composition of mycorrhizal communities, recent studies have challenged this assumption. Here, we investigated whether mycorrhizal communities and nutritional syndromes differed between seven green orchid species that co-occur under similar ecological conditions (coastal dune slacks). Our results showed that mycorrhizal communities differed significantly between orchid species. Rhizoctonia fungi dominated in Dactylorhiza sp., Herminium monorchis, and Epipactis palustris, which were autotrophic based on 13C and N content. Conversely, Liparis loeselii and Epipactis neerlandica associated primarily with ectomycorrhizal fungi but surprisingly, 13C and N content supported mixotrophy only in E. neerlandica. This, together with the finding of some ectomycorrhizal fungi in rhizoctonia-associated orchids, suggests that there exists an ecological continuum between the two syndromes. The presence of a large number of indicator species associating with individual orchid species further confirms previous findings that mycorrhizal fungi may be important factors driving niche-partitioning in terrestrial orchids and therefore contribute to orchid coexistence.

Nestedness in Arbuscular Mycorrhizal Fungal Communities along Soil pH Gradients in Early Primary Succession: Acid-Tolerant Fungi Are pH Generalists

Soil acidity is a major constraint on plant productivity. Arbuscular mycorrhizal (AM) fungi support plant colonization in acidic soil, but soil acidity also constrains fungal growth and diversity. Fungi in extreme environments generally evolve towards specialists, suggesting that AM fungi in acidic soil are acidic-soil specialists. In our previous surveys, however, some AM fungi detected in strongly acidic soils could also be detected in a soil with moderate pH, which raised a hypothesis that the fungi in acidic soils are pH generalists. To test the hypothesis, we conducted a pH-manipulation experiment and also analyzed AM fungal distribution along a pH gradient in the field using a synthesized dataset of the previous and recent surveys. Rhizosphere soils of the generalist plant Miscanthus sinensis were collected both from a neutral soil and an acidic soil, and M. sinensis seedlings were grown at three different pH. For the analysis of field communities, rhizosphere soils of M. sinensis were collected from six field sites across Japan, which covered a soil pH range of 3.0-7.4, and subjected to soil trap culture. AM fungal community compositions were determined based on LSU rDNA sequences. In the pH-manipulation experiment the acidification of medium had a significant impact on the compositions of the community from the neutral soil, but the neutralization of the medium had no effect on those of the community from the acidic soil. Furthermore, the communities in lower -pH soils were subsets of (nested in) those in higher-pH soils. In the field communities a significant nestedness pattern was observed along the pH gradient. These observations suggest that the fungi in strongly acidic soils are pH generalists that occur not only in acidic soil but also in wide ranges of soil pH. Nestedness in AM fungal community along pH gradients may have important implications for plant community resilience and early primary succession after disturbance in acidic soils.

Long-term increase in snow depth leads to compositional changes in arctic ectomycorrhizal fungal communities

Many arctic ecological processes are regulated by soil temperature that is tightly interconnected with snow cover distribution and persistence. Recently, various climate-induced changes have been observed in arctic tundra ecosystems, e.g. shrub expansion, resulting in reduction in albedo and greater C fixation in aboveground vegetation as well as increased rates of soil C mobilization by microbes. Importantly, the net effects of these shifts are unknown, in part because our understanding of belowground processes is limited. Here, we focus on the effects of increased snow depth, and as a consequence, increased winter soil temperature on ectomycorrhizal (ECM) fungal communities in dry and moist tundra. We analyzed deep DNA sequence data from soil samples taken at a long-term snow fence experiment in Northern Alaska. Our results indicate that, in contrast with previously observed responses of plants to increased snow depth at the same experimental site, the ECM fungal community of the dry tundra was more affected by deeper snow than the moist tundra community. In the dry tundra, both community richness and composition were significantly altered while in the moist tundra, only community composition changed significantly while richness did not. We observed a decrease in richness of Tomentella, Inocybe and other taxa adapted to scavenge the soil for labile N forms. On the other hand, richness of Cortinarius, and species with the ability to scavenge the soil for recalcitrant N forms, did not change. We further link ECM fungal traits with C soil pools. If future warmer atmospheric conditions lead to greater winter snow fall, changes in the ECM fungal community will likely influence C emissions and C fixation through altering N plant availability, fungal biomass and soil-plant C-N dynamics, ultimately determining important future interactions between the tundra biosphere and atmosphere.

Differences in mycorrhizal communities between Epipactis palustris, E. helleborine and its presumed sister species E. neerlandica

BACKGROUND AND AIMS

In orchid species that have populations occurring in strongly contrasting habitats, mycorrhizal divergence and other habitat-specific adaptations may lead to the formation of reproductively isolated taxa and ultimately to species formation. However, little is known about the mycorrhizal communities associated with recently diverged sister taxa that occupy different habitats.

METHODS

In this study, 454 amplicon pyrosequencing was used to investigate mycorrhizal communities associating with Epipactis helleborine in its typical forest habitat and with its presumed sister species E. neerlandica that almost exclusively occurs in coastal dune habitats. Samples of the phylogenetically more distant E. palustris, which co-occurred with E. neerlandica, were also included to investigate the role of habitat-specific conditions on mycorrhizal communities.

RESULTS

A total of 105 operational taxonomic units (OTUs) of putative orchid mycorrhizal fungi were observed in the three studied species. The majority of these fungi were endophytic fungi of Helotiales and ectomycorrhizal fungi belonging to Thelephoraceae, Sebacinaceae and Inocybaceae. In addition, a large number of other ectomycorrhizal taxa were detected, including Cortinarius, Cenococcum, Tuber, Geopora, Wilcoxina, Meliniomyces, Hebeloma, Tricholoma, Russula and Peziza Mycorrhizal communities differed significantly between the three species, but differences were most pronounced between the forest species (E. helleborine) and the two dune slack species (E. neerlandica and E. palustris).

CONCLUSION

The results clearly showed that recently diverged orchid species that occupy different habitats were characterized by significantly different mycorrhizal communities and call for more detailed experiments that aim at elucidating the contribution of habitat-specific adaptations in general and mycorrhizal divergence in particular to the process of speciation in orchids.

The utility of DNA metabarcoding for studying the response of arthropod diversity and composition to land-use change in the tropics

Metabarcoding potentially offers a rapid and cheap method of monitoring biodiversity, but real-world applications are few. We investigated its utility in studying patterns of litter arthropod diversity and composition in the tropics. We collected litter arthropods from 35 matched forest-plantation sites across Xishuangbanna, southwestern China. A new primer combination and the MiSeq platform were used to amplify and sequence a wide variety of litter arthropods using simulated and real-world communities. Quality filtered reads were clustered into 3,624 MOTUs at ≥97% similarity and the taxonomy of each MOTU was predicted. We compared diversity and compositional differences between forests and plantations (rubber and tea) for all MOTUs and for eight arthropod groups. We obtained ~100% detection rate after in silico sequencing six mock communities with known arthropod composition. Ordination showed that rubber, tea and forest communities formed distinct clusters. α-diversity declined significantly between forests and adjacent plantations for more arthropod groups in rubber than tea, and diversity of order Orthoptera increased significantly in tea. Turnover was higher in forests than plantations, but patterns differed among groups. Metabarcoding is useful for quantifying diversity patterns of arthropods under different land-uses and the MiSeq platform is effective for arthropod metabarcoding in the tropics.

Characterization of rhizosphere and endophytic fungal communities from roots of Stipa purpurea in alpine steppe around Qinghai Lake

Stipa purpurea is among constructive endemic species in the alpine steppe on the Qinghai-Xizang Plateau. To reveal the fungal community structure and diversity in the rhizosphere and roots of this important grass and to analyze the potential influence of different habitats on the structure of fungal communities, we explored the root endophyte and the directly associated rhizosphere communities of S. purpurea by using internal transcribed spacer rRNA cloning and sequencing methods. We found that the roots of S. purpurea are associated with a diverse consortium of Basidiomycota (59.8%) and Ascomycota (38.5%). Most fungi obtained from rhizosphere soil in S. purpurea have been identified as Ascomycetes, while the high proportion detected in roots were basidiomycetous endophytes. The species richness, diversity, and evenness of fungal assemblages were higher in roots than in the rhizosphere soil. Fungi inhabiting the rhizosphere and roots of S. purpurea are significantly different, and the rhizosphere and endophyte communities are largely independent with little overlap in the dominant phyla or operational taxonomic units. Taken together, these results suggested that a wide variety of fungal communities are associated with the roots and rhizosphere soil of S. purpurea and that the fungal assemblages are strongly influenced by different habitats.

Long-term warming alters richness and composition of taxonomic and functional groups of arctic fungi

Fungi, including symbionts, pathogens and decomposers, play crucial roles in community dynamics and nutrient cycling in terrestrial ecosystems. Despite their ecological importance, the response of most arctic fungi to climate warming is unknown, so are their potential roles in driving the observed and predicted changes in tundra communities. We carried out deep DNA sequencing of soil samples to study the long-term effects of experimental warming on fungal communities in dry heath and moist tussock tundra in Arctic Alaska. The data presented here indicate that fungal community composition responds strongly to warming in the moist tundra, but not in the dry tundra. While total fungal richness was not significantly affected by warming, there were clear correlations among operational taxonomic unit richness of various ecological and taxonomic groups and long-term warming. Richness of ectomycorrhizal, ericoid mycorrhizal and lichenized fungi generally decreased with warming, while richness of saprotrophic, plant and animal pathogenic, and root endophytic fungi tended to increase in the warmed plots. More importantly, various taxa within these functional guilds followed opposing trends that highlight the importance of species-specific responses to warming. We recommend that species-level ecological differences be taken into account in climate change and nutrient cycling studies that involve arctic fungi.

Mycorrhizal fungal communities in coastal sand dunes and heaths investigated by pyrosequencing analyses

Maritime sand dunes and coastal ericaceous heaths are unstable and dynamic habitats for mycorrhizal fungi. Creeping willow (Salix repens) is an important host plant in these habitats in parts of Europe. In this study, we wanted to assess which mycorrhizal fungi are associated with S. repens in four different coastal vegetation types in Southern Norway, three types from sand dunes and one from heaths. Moreover, we investigated which ecological factors are important for the fungal community structure in these vegetation types. Mycorrhizal fungi on S. repens root samples were identified by 454 pyrosequencing of tag-encoded internal transcribed spacer 1 (ITS1) amplicons. Significantly higher fungal richness was observed in hummock dunes and dune slacks compared to eroded dune vegetation. The compositional variation was mainly accounted for by location (plot) and vegetation type and was significantly correlated to content of carbon, nitrogen and phosphorus in soil. The investigated maritime sand dunes and coastal ericaceous heaths hosted mycorrhizal taxa mainly associated with Helotiales, Sebacinales, Thelephorales and Agaricales.

Phylogenic diversity and tissue specificity of fungal endophytes associated with the pharmaceutical plant, Stellera chamaejasme L. revealed by a cultivation-independent approach

The fungal endophytes associated with medicinal plants have been demonstrated as a reservoir with novel natural products useful in medicine and agriculture. It is desirable to explore the species composition, diversity and tissue specificity of endophytic fungi that inhabit in different tissues of medicinal plants. In this study, a culture-independent survey of fungal diversity in the rhizosphere, leaves, stems and roots of a toxic medicinal plant, Stellera chamaejasme L., was conducted by sequence analysis of clone libraries of the partial internal transcribed spacer region. Altogether, 145 fungal OTUs (operational taxonomic units), represented by 464 sequences, were found in four samples, of these 109 OTUs (75.2 %) belonging to Ascomycota, 20 (13.8 %) to Basidiomycota, 14 (9.7 %) to Zygomycota, 1 (0.7 %) to Chytridiomycota, and 1 (0.7 %) to Glomeromycota. The richness and diversity of fungal communities were strongly influenced by plant tissue environments, and the roots are associated with a surprisingly rich endophyte community. The endophyte assemblages associated with S. chamaejasme were strongly shaped by plant tissue environments, and exhibited a certain degree of tissue specificity. Our results suggested that a wide variety of fungal assemblages inhabit in S. chamaejasme, and plant tissue environments conspicuously influence endophyte community structure.

Summer temperature increase has distinct effects on the ectomycorrhizal fungal communities of moist tussock and dry tundra in Arctic Alaska

Arctic regions are experiencing the greatest rates of climate warming on the planet and marked changes have already been observed in terrestrial arctic ecosystems. While most studies have focused on the effects of warming on arctic vegetation and nutrient cycling, little is known about how belowground communities, such as fungi root-associated, respond to warming. Here, we investigate how long-term summer warming affects ectomycorrhizal (ECM) fungal communities. We used Ion Torrent sequencing of the rDNA internal transcribed spacer 2 (ITS2) region to compare ECM fungal communities in plots with and without long-term experimental warming in both dry and moist tussock tundra. Cortinarius was the most OTU-rich genus in the moist tundra, while the most diverse genus in the dry tundra was Tomentella. On the diversity level, in the moist tundra we found significant differences in community composition, and a sharp decrease in the richness of ECM fungi due to warming. On the functional level, our results indicate that warming induces shifts in the extramatrical properties of the communities, where the species with medium-distance exploration type seem to be favored with potential implications for the mobilization of different nutrient pools in the soil. In the dry tundra, neither community richness nor community composition was significantly altered by warming, similar to what had been observed in ECM host plants. There was, however, a marginally significant increase in OTUs identified as ECM fungi with the medium-distance exploration type in the warmed plots. Linking our findings of decreasing richness with previous results of increasing ECM fungal biomass suggests that certain ECM species are favored by warming and may become more abundant, while many other species may go locally extinct due to direct or indirect effects of warming. Such compositional shifts in the community might affect nutrient cycling and soil organic C storage.