High richness of ectomycorrhizal fungi and low host specificity in a coastal sand dune ecosystem revealed by network analysis

Ectomycorrhizal (EM) fungi are ubiquitous in temperate and boreal forests, comprising over 20,000 species forming root symbiotic associations with Pinaceae and woody angiosperms. As much as 100 different EM fungal species can coexist and interact with the same tree species, forming complex multispecies networks in soils. The degree of host specificity and structural properties of these interaction networks (e.g., nestedness and modularity) may influence plant and fungal community assembly and species coexistence, yet their structure has been little studied in northern coniferous forests, where trees depend on EM fungi for nutrient acquisition. We used high‐throughput sequencing to characterize the composition and diversity of bulk soil and root‐associated fungal communities in four co‐occurring Pinaceae in a relic foredune plain located at Îles de la Madeleine, Québec, Canada. We found high EM fungal richness across the four hosts, with a total of 200 EM operational taxonomic units (OTUs), mainly belonging to the Agaricomycetes. Network analysis revealed an antinested pattern in both bulk soil and roots EM fungal communities. However, there was no detectable modularity (i.e., subgroups of interacting species) in the interaction networks, indicating a low level of specificity in these EM associations. In addition, there were no differences in EM fungal OTU richness or community structure among the four tree species. Limited shared resources and competitive exclusion typically restrict the number of taxa coexisting within the same niche. As such, our finding of high EM fungal richness and low host specificity highlights the need for further studies to determine the mechanisms enabling such a large number of EM fungal species to coexist locally on the same hosts.