Utility of large subunit for environmental sequencing of arbuscular mycorrhizal fungi: a new reference database and pipeline

Tracking arbuscular mycorrhizal fungi to their source: active inoculation and passive dispersal differentially affect community assembly in urban soils

Communities of arbuscular mycorrhizal (AM) fungi assemble passively over time via biotic and abiotic mechanisms. In degraded soils, AM fungal communities can assemble actively when humans manage mycorrhizas for ecosystem restoration. We investigated mechanisms of urban AM fungal community assembly in a 2-yr green roof experiment. We compared AM fungal communities in inoculated and uninoculated trays to samples from two potential sources: the inoculum and air. Active inoculation stimulated more distinct and diverse AM fungal communities, an effect that intensified over time. In the treatment trays, 45% of AM fungal taxa were detected in the inoculum, 2% were detected in aerial samples, 23% were detected in both inoculum and air, and 30% were not detected in either source. Passive dispersal of AM fungi likely resulted in the successful establishment of a small number of species, but active inoculation with native AM fungal species resulted in an immediate shift to a diverse and unique fungal community. When urban soils are constructed or modified by human activity, this is an opportunity for intervention with AM fungi that will persist and add diversity to that system.

Three new species of arbuscular mycorrhizal fungi (Glomeromycota) and Acaulospora gedanensis revised

Studies of the morphology and the 45S nuc rDNA phylogeny of three potentially undescribed arbuscular mycorrhizal fungi (phylum Glomeromycota) grown in cultures showed that one of these fungi is a new species of the genus Diversispora in the family Diversisporaceae; the other two fungi are new Scutellospora species in Scutellosporaceae. Diversispora vistulana sp. nov. came from maritime sand dunes of the Vistula Spit in northern Poland, and S. graeca sp. nov. and S. intraundulata sp. nov. originally inhabited the Mediterranean dunes of the Peloponnese Peninsula, Greece. In addition, the morphological description of spores of Acaulospora gedanensis, originally described in 1988, was emended based on newly found specimens, and the so far unknown phylogeny of this species was determined. The phylogenetic analyses of 45S sequences placed this species among Acaulospora species with atypical phenotypic and histochemical features of components of the two inner germinal walls.

Dilution of specialist pathogens drives productivity benefits from diversity in plant mixtures

Productivity benefits from diversity can arise when compatible pathogen hosts are buffered by unrelated neighbors, diluting pathogen impacts. However, the generality of pathogen dilution has been controversial and rarely tested within biodiversity manipulations. Here, we test whether soil pathogen dilution generates diversity- productivity relationships using a field biodiversity-manipulation experiment, greenhouse assays, and feedback modeling. We find that the accumulation of specialist pathogens in monocultures decreases host plant yields and that pathogen dilution predicts plant productivity gains derived from diversity. Pathogen specialization predicts the strength of the negative feedback between plant species in greenhouse assays. These feedbacks significantly predict the overyielding measured in the field the following year. This relationship strengthens when accounting for the expected dilution of pathogens in mixtures. Using a feedback model, we corroborate that pathogen dilution drives overyielding. Combined empirical and theoretical evidence indicate that specialist pathogen dilution generates overyielding and suggests that the risk of losing productivity benefits from diversity may be highest where environmental change decouples plant-microbe interactions.

Interactions with fungi vary among Tripsacum dactyloides genotypes from across a precipitation gradient

Plant-associated microbes, specifically fungal endophytes, augment the ability of many grasses to adapt to extreme environmental conditions. Tripsacum dactyloides (Eastern gamagrass) is a perennial, drought-tolerant grass native to the tallgrass prairies of the central USA. The extent to which the microbiome of T. dactyloides contributes to its drought tolerance is unknown. Ninety-seven genotypes of T. dactyloides were collected from native populations across an east-west precipitation gradient in Kansas, Oklahoma and Texas, and then grown together in a common garden for over 20 years. Root and leaf samples were visually examined for fungal density. Because fungal endophytes confer drought-tolerant capabilities to their host plants, we expected to find higher densities of fungal endophytes in plants from western, drier regions, compared to plants from eastern, wetter regions. Results confirmed a negative correlation between endophyte densities in roots and precipitation at the genotype's original location (r = -0.21 P = 0.04). Our analyses reveal that the host genotype's origin along the precipitation gradient predicts the absolute abundance of symbionts in the root, but not the relative abundances of particular organisms or the overall community composition. Overall, these results demonstrate that genetic variation for plant-microbe interactions can reflect historical environment, and reinforce the importance of considering plant genotype in conservation and restoration work in tallgrass prairie ecosystems.

GlobalAMFungi: a global database of arbuscular mycorrhizal fungal occurrences from high-throughput sequencing metabarcoding studies

Arbuscular mycorrhizal (AM) fungi are crucial mutualistic symbionts of the majority of plant species, with essential roles in plant nutrient uptake and stress mitigation. The importance of AM fungi in ecosystems contrasts with our limited understanding of the patterns of AM fungal biogeography and the environmental factors that drive those patterns. This article presents a release of a newly developed global AM fungal dataset (GlobalAMFungi database, https://globalamfungi.com) that aims to reduce this knowledge gap. It contains almost 50 million observations of Glomeromycotinian AM fungal amplicon DNA sequences across almost 8500 samples with geographical locations and additional metadata obtained from 100 original studies. The GlobalAMFungi database is built on sequencing data originating from AM fungal taxon barcoding regions in: i) the small subunit rRNA (SSU) gene; ii) the internal transcribed spacer 2 (ITS2) region; and iii) the large subunit rRNA (LSU) gene. The GlobalAMFungi database is an open source and open access initiative that compiles the most comprehensive atlas of AM fungal distribution. It is designed as a permanent effort that will be continuously updated by its creators and through the collaboration of the scientific community. This study also documented applicability of the dataset to better understand ecology of AM fungal taxa.

Above- and belowground fungal biodiversity of Populus trees on a continental scale

Understanding drivers of terrestrial fungal communities over large scales is an important challenge for predicting the fate of ecosystems under climate change and providing critical ecological context for bioengineering plant-microbe interactions in model systems. We conducted an extensive molecular and microscopy field study across the contiguous United States measuring natural variation in the Populus fungal microbiome among tree species, plant niche compartments and key symbionts. Our results show clear biodiversity hotspots and regional endemism of Populus-associated fungal communities explained by a combination of climate, soil and geographic factors. Modelling climate change impacts showed a deterioration of Populus mycorrhizal associations and an increase in potentially pathogenic foliar endophyte diversity and prevalence. Geographic differences among these symbiont groups in their sensitivity to environmental change are likely to influence broader forest health and ecosystem function. This dataset provides an above- and belowground atlas of Populus fungal biodiversity at a continental scale.

Arbuscular mycorrhizal fungi associated with the rhizosphere of an endemic terrestrial bromeliad and a grass in the Brazilian neotropical dry forest

Arbuscular mycorrhizal fungi form symbiotic associations with 80-90% of all known plants, allowing the fungi to acquire plant-synthesized carbon, and confer an increased capacity for nutrient uptake by plants, improving tolerance to abiotic and biotic stresses. We aimed at characterizing the mycorrhizal community in the rhizosphere of Neoglaziovia variegata (so-called `caroa`) and Tripogonella spicata (so-called resurrection plant), using high-throughput sequencing of the partial 18S rRNA gene. Both plants are currently undergoing a bioprospecting program to find microbes with the potential of helping plants tolerate water stress. Sampling was carried out in the Caatinga biome, a neotropical dry forest, located in northeastern Brazil. Illumina MiSeq sequencing of 37 rhizosphere samples (19 for N. variegata and 18 for T. spicata) revealed a distinct mycorrhizal community between the studied plants. According to alpha diversity analyses, T. spicata showed the highest richness and diversity based on the Observed ASVs and the Shannon index, respectively. On the other hand, N. variegata showed higher modularity of the mycorrhizal network compared to T. spicata. The four most abundant genera found (higher than 10%) were Glomus, Gigaspora, Acaulospora, and Scutellospora, with Glomus being the most abundant in both plants. Nonetheless, Gigaspora, Diversispora, and Ambispora were found only in the rhizosphere of N. variegata, whilst Scutellospora, Paraglomus, and Archaeospora were exclusive to the rhizosphere of T. spicata. Therefore, the community of arbuscular mycorrhizal fungi of the rhizosphere of each plant encompasses a unique composition, structure and modularity, which can differentially assist them in the hostile environment.

Phenol and Polyaromatic Hydrocarbons Are Stronger Drivers Than Host Plant Species in Shaping the Arbuscular Mycorrhizal Fungal Component of the Mycorrhizosphere

Changes in soil microbial communities in response to hydrocarbon pollution are critical indicators of disturbed ecosystem conditions. A core component of these communities that is functionally adjusted to the life-history traits of the host and environmental factors consists of arbuscular mycorrhizal fungi (AMF). AMF communities associated with Poa trivialis and Phragmites australis growing at a phenol and polynuclear aromatic hydrocarbon (PAH)-contaminated site and at an uncontaminated site were compared based on LSU rDNA sequencing. Dissimilarities in species composition and community structures indicated soil pollution as the main factor negatively affecting the AMF diversity. The AMF communities at the contaminated site were dominated by fungal generalists (Rhizophagus, Funneliformis, Claroideoglomus, Paraglomus) with wide ecological tolerance. At the control site, the AMF communities were characterized by higher taxonomic and functional diversity than those exposed to the contamination. The host plant identity was the main driver distinguishing the two AMF metacommunities. The AMF communities at the uncontaminated site were represented by Polonospora, Paraglomus, Oehlia, Nanoglomus, Rhizoglomus, Dominikia, and Microdominikia. Polonosporaceae and Paraglomeraceae were particularly dominant in the Ph. australis mycorrhizosphere. The high abundance of early diverging AMF could be due to the use of primers able to detect lineages such as Paraglomeracae that have not been recognized by previously used 18S rDNA primers.

Analyzing diversification dynamics using barcoding data: the case of an obligate mycorrhizal symbiont

Analysing diversification dynamics is key to understanding the past evolutionary history of clades that led to present‐day biodiversity patterns. While such analyses are widespread in well‐characterized groups of species, they are much more challenging in groups for which diversity is mostly known through molecular techniques. Here, we use the largest global database on the small subunit (SSU) rRNA gene of Glomeromycotina, a subphylum of microscopic arbuscular mycorrhizal fungi that provide mineral nutrients to most land plants by forming one of the oldest terrestrial symbioses, to analyse the diversification dynamics of this clade in the past 500 million years. We perform a range of sensitivity analyses and simulations to control for potential biases linked to the nature of the data. We find that Glomeromycotina tend to have low speciation rates compared to other eukaryotes. After a peak of speciations between 200 and 100 million years ago, they experienced an important decline in speciation rates toward the present. Such a decline could be at least partially related to a shrinking of their mycorrhizal niches and to their limited ability to colonize new niches. Our analyses identify patterns of diversification in a group of obligate symbionts of major ecological and evolutionary importance and illustrate that short molecular markers combined with intensive sensitivity analyses can be useful for studying diversification dynamics in microbial groups.

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.

Environmental identification of arbuscular mycorrhizal fungi using the LSU rDNA gene region: an expanded database and improved pipeline

Arbuscular mycorrhizal fungi (AMF; Glomeromycota) are difficult to culture; therefore, establishing a robust amplicon-based approach to taxa identification is imperative to describe AMF diversity. Further, due to low and biased sampling of AMF taxa, molecular databases do not represent the breadth of AMF diversity, making database matching approaches suboptimal. Therefore, a full description of AMF diversity requires a tool to determine sequence-based placement in the Glomeromycota clade. Nonetheless, commonly used gene regions, including the SSU and ITS, do not enable reliable phylogenetic placement. Here, we present an improved database and pipeline for the phylogenetic determination of AMF using amplicons from the large subunit (LSU) rRNA gene. We improve our database and backbone tree by including additional outgroup sequences. We also improve an existing bioinformatics pipeline by aligning forward and reverse reads separately, using a universal alignment for all tree building, and implementing a BLAST screening prior to tree building to remove non-homologous sequences. Finally, we present a script to extract AMF belonging to 11 major families as well as an amplicon sequencing variant (ASV) version of our pipeline. We test the utility of the pipeline by testing the placement of known AMF, known non-AMF, and Acaulospora sp. spore sequences. This work represents the most comprehensive database and pipeline for phylogenetic placement of AMF LSU amplicon sequences within the Glomeromycota clade.

Site-Dependent Relationships Between Fungal Community Composition, Plant Genotypic Diversity and Environmental Drivers in a Salix Biomass System

Soil fungi are strongly affected by plant species or genotypes since plants modify their surrounding environment, but the effects of plant genotype diversity on fungal diversity and function have not been extensively studied. The interactive responses of fungal community composition to plant genotypic diversity and environmental drivers were investigated in Salix biomass systems, posing questions about: (1) How fungal diversity varies as a function of plant genotype diversity; (2) If plant genotype identity is a strong driver of fungal community composition also in plant mixtures; (3) How the fungal communities change through time (seasonally and interannually)?; and (4) Will the proportion of ECM fungi increase over the rotation? Soil samples were collected over 4 years, starting preplanting from two Salix field trials, including four genotypes with contrasting phenology and functional traits, and genotypes were grown in all possible combinations (four genotypes in Uppsala, Sweden, two in Rostock, Germany). Fungal communities were identified, using Pacific Biosciences sequencing of fungal ITS2 amplicons. We found some site-dependent relationships between fungal community composition and genotype or diversity level, and site accounted for the largest part of the variation in fungal community composition. Rostock had a more homogenous community structure, with significant effects of genotype, diversity level, and the presence of one genotype ("Loden") on fungal community composition. Soil properties and plant and litter traits contributed to explaining the variation in fungal species composition. The within-season variation in composition was of a similar magnitude to the year-to-year variation. The proportion of ECM fungi increased over time irrespective of plant genotype diversity, and, in Uppsala, the 4-mixture showed a weaker response than other combinations. Species richness was generally higher in Uppsala compared with that in Rostock and increased over time, but did not increase with plant genotype diversity. This significant site-specificity underlines the need for consideration of diverse sites to draw general conclusions of temporal variations and functioning of fungal communities. A significant increase in ECM colonization of soil under the pioneer tree Salix on agricultural soils was evident and points to changed litter decomposition and soil carbon dynamics during Salix growth.

PacBio sequencing of Glomeromycota rDNA: a novel amplicon covering all widely used ribosomal barcoding regions and its applicability in taxonomy and ecology of arbuscular mycorrhizal fungi

There is no consensus barcoding region for determination of arbuscular mycorrhizal fungal (AMF) taxa. To overcome this obstacle, we have developed an approach to sequence an AMF marker within the ribosome-encoding operon (rDNA) that covers all three widely applied variable molecular markers. Using a nested PCR approach specific to AMF, we amplified a part (c. 2.5 kb) of the rDNA spanning the majority of the small subunit rRNA (SSU) gene, the complete internal transcribed spacer (ITS) region and a part of the large subunit (LSU) rRNA gene. The PCR products were sequenced on the PacBio platform utilizing Single Molecule Real Time (SMRT) sequencing. Employing this method for selected environmental DNA samples, we were able to describe complex AMF communities consisting of various glomeromycotan lineages. We demonstrate the applicability of this new 2.5 kb approach to provide robust phylogenetic assignment of AMF lineages without known sequences from pure cultures and to consolidate information about AMF taxon distributions coming from three widely used barcoding regions into one integrative dataset.