Large-scale fungal strain sequencing unravels the molecular diversity in mating loci maintained by long-term balancing selection.

Balancing selection, an evolutionary force that retains genetic diversity, has been detected in multiple genes and organisms, such as the sexual mating loci in fungi. However, to quantify the strength of balancing selection and define the mating-related genes require a large number of strains. In tetrapolar basidiomycete fungi, sexual type is determined by two unlinked loci, MATA and MATB. Genes in both loci define mating type identity, control successful mating and completion of the life cycle. These loci are usually highly diverse. Previous studies have speculated, based on culture crosses, that species of the non-model genus Trichaptum (Hymenochaetales, Basidiomycota) possess a tetrapolar mating system, with multiple alleles. Here, we sequenced a hundred and eighty strains of three Trichaptum species. We characterized the chromosomal location of MATA and MATB, the molecular structure of MAT regions and their allelic richness. The sequencing effort was sufficient to molecularly characterize multiple MAT alleles segregating before the speciation event of Trichaptum species. Analyses suggested that long-term balancing selection has generated trans-species polymorphisms. Mating sequences were classified in different allelic classes based on an amino acid identity (AAI) threshold supported by phylogenetics. 17,550 mating types were predicted based on the allelic classes. In vitro crosses allowed us to support the degree of allelic divergence needed for successful mating. Even with the high amount of divergence, key amino acids in functional domains are conserved. We conclude that the genetic diversity of mating loci in Trichaptum is due to long-term balancing selection, with limited recombination and duplication activity. The large number of sequenced strains highlighted the importance of sequencing multiple individuals from different species to detect the mating-related genes, the mechanisms generating diversity and the evolutionary forces maintaining them.

Mast seeding in European beech (Fagus sylvatica L.) is associated with reduced fungal sporocarp production and community diversity.

Mast seeding is a well-documented phenomenon across diverse forest ecosystems. While its effect on aboveground food webs has been thoroughly studied, how it impacts the soil fungi that drive soil carbon and nutrient cycling has not yet been explored. To evaluate the relationship between mast seeding and fungal resource availability, we paired a Swiss 29-year fungal sporocarp census with contemporaneous seed production for European beech (Fagus sylvatica L.). On average, mast seeding was associated with a 55% reduction in sporocarp production and a compositional community shift towards drought-tolerant taxa across both ectomycorrhizal and saprotrophic guilds. Among ectomycorrhizal fungi, traits associated with carbon cost did not explain species' sensitivity to seed production. Together, our results support a novel hypothesis that mast seeding limits annual resource availability and reproductive investment in soil fungi, creating an ecosystem 'rhythm' to forest processes that is synchronized above- and belowground.

Global field collection data confirm an affinity of brown rot fungi for coniferous habitats and substrates.

Unlike 'white rot' (WR) wood-decomposing fungi that remove lignin to access cellulosic sugars, 'brown rot' (BR) fungi selectively extract sugars and leave lignin behind. The relative frequency and distribution of these fungal types (decay modes) have not been thoroughly assessed at a global scale; thus, the fate of one-third of Earth's aboveground carbon, wood lignin, remains unclear. Using c. 1.5 million fungal sporocarp and c. 30 million tree records from publicly accessible databases, we mapped and compared decay mode and tree type (conifer vs angiosperm) distributions. Additionally, we mined fungal record metadata to assess substrate specificity per decay mode. The global average for BR fungi proportion (BR/(BR + WR records)) was 13% and geographic variation was positively correlated (R2 = 0.45) with conifer trees proportion (conifer/(conifer + angiosperm records)). Most BR species (61%) were conifer, rather than angiosperm (22%), specialists. The reverse was true for WR (conifer: 19%; angiosperm: 62%). Global BR proportion patterns were predicted with greater accuracy using the relative distributions of individual tree species (R2 = 0.82), rather than tree type. Fungal decay mode distributions can be explained by tree type and, more importantly, tree species distributions, which our data suggest is due to strong substrate specificities.

Symbiotic status alters fungal eco-evolutionary offspring trajectories.

Despite host-fungal symbiotic interactions being ubiquitous in all ecosystems, understanding how symbiosis has shaped the ecology and evolution of fungal spores that are involved in dispersal and colonization of their hosts has been ignored in life-history studies. We assembled a spore morphology database covering over 26,000 species of free-living to symbiotic fungi of plants, insects and humans and found more than eight orders of variation in spore size. Evolutionary transitions in symbiotic status correlated with shifts in spore size, but the strength of this effect varied widely among phyla. Symbiotic status explained more variation than climatic variables in the current distribution of spore sizes of plant-associated fungi at a global scale while the dispersal potential of their spores is more restricted compared to free-living fungi. Our work advances life-history theory by highlighting how the interaction between symbiosis and offspring morphology shapes the reproductive and dispersal strategies among living forms.

Mycena species can be opportunist-generalist plant root invaders.

Traditional strict separation of fungi into ecological niches as mutualist, parasite or saprotroph is increasingly called into question. Sequences of assumed saprotrophs have been amplified from plant root interiors, and several saprotrophic genera can invade and interact with host plants in laboratory growth experiments. However, it is uncertain if root invasion by saprotrophic fungi is a widespread phenomenon and if laboratory interactions mirror field conditions. Here, we focused on the widespread and speciose saprotrophic genus Mycena and performed (1) a systematic survey of their occurrences (in ITS1/ITS2 datasets) in mycorrhizal roots of 10 plant species, and (2) an analysis of natural abundances of 13 C/15 N stable isotope signatures of Mycena basidiocarps from five field locations to examine their trophic status. We found that Mycena was the only saprotrophic genus consistently found in 9 out of 10 plant host roots, with no indication that the host roots were senescent or otherwise vulnerable. Furthermore, Mycena basidiocarps displayed isotopic signatures consistent with published 13 C/15 N profiles of both saprotrophic and mutualistic lifestyles, supporting earlier laboratory-based studies. We argue that Mycena are widespread latent invaders of healthy plant roots and that Mycena species may form a spectrum of interactions besides saprotrophy also in the field.

DNA metabarcoding reveals host-specific communities of arthropods residing in fungal fruit bodies.

Biological communities within living organisms are structured by their host's traits. How host traits affect biodiversity and community composition is poorly explored for some associations, such as arthropods within fungal fruit bodies. Using DNA metabarcoding, we characterized the arthropod communities in living fruit bodies of 11 wood-decay fungi from boreal forests and investigated how they were affected by different fungal traits. Arthropod diversity was higher in fruit bodies with a larger surface area-to-volume ratio, suggesting that colonization is crucial to maintain arthropod populations. Diversity was not higher in long-lived fruit bodies, most likely because these fungi invest in physical or chemical defences against arthropods. Arthropod community composition was structured by all measured host traits, namely fruit body size, thickness, surface area, morphology and toughness. Notably, we identified a community gradient where soft and short-lived fruit bodies harboured more true flies, while tougher and long-lived fruit bodies had more oribatid mites and beetles, which might reflect different development times of the arthropods. Ultimately, close to 75% of the arthropods were specific to one or two fungal hosts. Besides revealing surprisingly diverse and host-specific arthropod communities within fungal fruit bodies, our study provided insight into how host traits structure communities.

The Indoor Mycobiomes of Daycare Centers Are Affected by Occupancy and Climate.

Many children spend considerable time in daycare centers and may be influenced by the indoor microorganisms there, including fungi. In this study, we investigate the indoor mycobiomes of 125 daycare centers distributed along strong environmental gradients throughout Norway. Dust samples were collected from doorframes outside and inside buildings using a community science sampling approach. Fungal communities in the dust samples were analyzed using DNA metabarcoding of the internal transcribed spacer 2 (ITS2) region. We observed a marked difference between the outdoor and indoor mycobiomes. The indoor mycobiomes included considerably more yeasts and molds than the outdoor samples, with Saccharomyces, Mucor, Malassezia, and Penicillium being among the most dominant fungal genera. Changes in the indoor fungal richness and composition correlated with numerous variables related to both outdoor and indoor conditions; there was a clear geographic structure in the indoor mycobiome composition that mirrored the outdoor climate, ranging from humid areas in western Norway to drier and colder areas in eastern Norway. Moreover, the number of children in the daycare centers, as well as various building features, influenced the indoor mycobiome composition. We conclude that the indoor mycobiomes in Norwegian daycare centers are structured by multiple factors and are dominated by yeasts and molds. This study exemplifies how community science sampling enables DNA-based analyses of a high number of samples covering wide geographic areas. IMPORTANCE With an alarming increase in chronic diseases like childhood asthma and allergies, there is an increased focus on the exposure of young children to indoor biological and chemical air pollutants. Our study of 125 daycares throughout Norway demonstrates that the indoor mycobiome not only reflects cooccurring outdoor fungi but also includes a high abundance of yeast and mold fungi with an affinity for indoor environments. A multitude of factors influence the indoor mycobiomes in daycares, including the building type, inhabitants, as well as the outdoor environment. Many of the detected yeasts and molds are likely associated with the human body, where some have been coupled with allergies and respiratory problems. Our results call for further studies investigating the potential impact of the identified daycare-associated mycobiomes on children's health.

Shift in tree species changes the belowground biota of boreal forests.

The replacement of native birch with Norway spruce has been initiated in Norway to increase long-term carbon storage in forests. However, there is limited knowledge on the impacts that aboveground changes will have on the belowground microbiota. We examined which effects a tree species shift from birch to spruce stands has on belowground microbial communities, soil fungal biomass and relationships with vegetation biomass and soil organic carbon (SOC). Replacement of birch with spruce negatively influenced soil bacterial and fungal richness and strongly altered microbial community composition in the forest floor layer, most strikingly for fungi. Tree species-mediated variation in soil properties was a major factor explaining variation in bacterial communities. For fungi, both soil chemistry and understorey vegetation were important community structuring factors, particularly for ectomycorrhizal fungi. The relative abundance of ectomycorrhizal fungi and the ectomycorrhizal : saprotrophic fungal ratio were higher in spruce compared to birch stands, particularly in the deeper mineral soil layers, and vice versa for saprotrophs. The positive relationship between ergosterol (fungal biomass) and SOC stock in the forest floor layer suggests higher carbon sequestration potential in spruce forest soil, alternatively, that the larger carbon stock leads to an increase in soil fungal biomass.

Warming drives a 'hummockification' of microbial communities associated with decomposing mycorrhizal fungal necromass in peatlands.

Dead fungal mycelium (necromass) represents a critical component of soil carbon (C) and nutrient cycles. Assessing how the microbial communities associated with decomposing fungal necromass change as global temperatures rise will help in determining how these belowground organic matter inputs contribute to ecosystem responses. In this study, we characterized the structure of bacterial and fungal communities associated with multiple types of decaying mycorrhizal fungal necromass incubated within mesh bags across a 9°C whole ecosystem temperature enhancement in a boreal peatland. We found major taxonomic and functional shifts in the microbial communities present on decaying mycorrhizal fungal necromass in response to warming. These changes were most pronounced in hollow microsites, which showed convergence towards the necromass-associated microbial communities present in unwarmed hummocks. We also observed a high colonization of ericoid mycorrhizal fungal necromass by fungi from the same genera as the necromass. These results indicate that microbial communities associated with mycorrhizal fungal necromass decomposition are likely to change significantly with future climate warming, which may have strong impacts on soil biogeochemical cycles in peatlands. Additionally, the high enrichment of congeneric fungal decomposers on ericoid mycorrhizal necromass may help to explain the increase in ericoid shrub dominance in warming peatlands.

Legacies of invertebrate exclusion and tree secondary metabolites control fungal communities in dead wood.

During decomposition of organic matter, microbial communities may follow different successional trajectories depending on the initial environment and colonizers. The timing and order of the species arrival (assembly history) can lead to divergent communities through priority effects. We explored how assembly history and resource quality affected fungal communities and decay rate of decomposing wood, 1.5 and 4.5 years after tree felling. Additionally, we investigated the effect of invertebrate exclusion during the first two summers. We measured initial resource quality of bark and wood of aspen (Populus tremula) logs and surveyed the fungal communities by DNA metabarcoding at different times during succession. We found that gradients in fungal community composition were related to resource quality and we discuss how this may reflect different fungal life history strategies. As with previous studies, the initial amount of bark tannins was negatively correlated with wood decomposition rate over 4.5 years. The initial fungal community explained variation in community composition after 1.5, but not 4.5, years of succession. Although the assembly history of initial colonizers may cause alternative trajectories in successional communities, our results indicate that the communities may converge with the arrival of secondary colonizers. We also identified a strong legacy of invertebrate exclusion on fungal communities, even after 4.5 years of succession, thereby adding crucial knowledge on the importance of invertebrates in affecting fungal community development. By measuring and manipulating aspects of assembly history and resource quality that have rarely been studied, we expand our understanding of the complexity of fungal community dynamics.

Population genomics of a forest fungus reveals high gene flow and climate adaptation signatures.

Genome sequencing of spatially distributed individuals sheds light on how evolution structures genetic variation. Populations of Phellopilus nigrolimitatus, a red-listed wood-inhabiting fungus associated with old-growth coniferous forests, have decreased in size over the last century due to a loss of suitable habitats. We assessed the population genetic structure and investigated local adaptation in P. nigrolimitatus, by establishing a reference genome and genotyping 327 individuals sampled from 24 locations in Northern Europe by RAD sequencing. We revealed a shallow population genetic structure, indicating large historical population sizes and high levels of gene flow. Despite this weak substructuring, two genetic groups were recognized; a western group distributed mostly in Norway and an eastern group covering most of Finland, Poland and Russia. This substructuring may reflect coimmigration with the main host, Norway spruce (Picea abies), into Northern Europe after the last ice age. We found evidence of low levels of genetic diversity in southwestern Finland, which has a long history of intensive forestry and urbanization. Numerous loci were significantly associated with one or more environmental factors, indicating adaptation to specific environments. These loci clustered into two groups with different associations with temperature and precipitation. Overall, our findings indicate that the current population genetic structure of P. nigrolimitatus results from a combination of gene flow, genetic drift and selection. The acquisition of similar knowledge especially over broad geographic scales, linking signatures of adaptive genetic variation to evolutionary processes and environmental variation, for other fungal species will undoubtedly be useful for assessment of the combined effects of habitat fragmentation and climate change on fungi strongly bound to old-growth forests.

Fungal community dynamics across a forest-alpine ecotone.

Climate change is causing upward shift of forest lines worldwide, with consequences for soil biota and carbon (C) sequestration. We here analyse compositional changes in the soil biota across the forest line ecotone, an important transition zone between different ecosystems. We collected soil samples along transects stretching from subalpine mountain birch forests to alpine heath. Soil fungi and micro-eukaryotes were surveyed using DNA metabarcoding of the ITS2 and 18S markers, while ergosterol was used to quantify fungal biomass. We observed a strong shift in the soil biota across the forest line ecotone: Below the forest line, there were higher proportions of basidiomycetes and mucoromycetes, including ectomycorrhizal and saprotrophic fungi. Above it, we observed relatively more root-associated ascomycetes, including Archaeorhizomycetes, ericoid mycorrhizal fungi and dark septate endophytes. Ergosterol and percentage C content in soil correlated strongly and positively with the abundance of root-associated ascomycetes. The predominance of ectomycorrhizal and saprotrophic fungi below the forest line probably promote high C turnover, while root-associated ascomycetes above the forest line may enhance C sequestration. With further rise in forest lines, there will be a corresponding shift in the below-ground biota, probably leading to enhanced release of soil C.

Analysing indoor mycobiomes through a large-scale citizen science study in Norway.

In the built environment, fungi can cause important deterioration of building materials and have adverse health effects on occupants. Increased knowledge about indoor mycobiomes from different regions of the world, and their main environmental determinants, will enable improved indoor air quality management and identification of health risks. This is the first citizen science study of indoor mycobiomes at a large geographical scale in Europe, including 271 houses from Norway and 807 dust samples from three house compartments: outside of the building, living room and bathroom. The fungal community composition determined by DNA metabarcoding was clearly different between indoor and outdoor samples, but there were no significant differences between the two indoor compartments. The 32 selected variables, related to the outdoor environment, building features and occupant characteristics, accounted for 15% of the overall variation in community composition, with the house compartment as the key factor (7.6%). Next, climate was the main driver of the dust mycobiomes (4.2%), while building and occupant variables had significant but minor influences (1.4% and 1.1%, respectively). The house-dust mycobiomes were dominated by ascomycetes (⁓70%) with Capnodiales and Eurotiales as the most abundant orders. Compared to the outdoor samples, the indoor mycobiomes showed higher species richness, which is probably due to the mixture of fungi from outdoor and indoor sources. The main indoor indicator fungi belonged to two ecological groups with allergenic potential: xerophilic moulds and skin-associated yeasts. Our results suggest that citizen science is a successful approach for unravelling the built microbiome at large geographical scales.

Contrasting demographic histories revealed in two invasive populations of the dry rot fungus Serpula lacrymans.

Globalization and international trade have impacted organisms around the world leading to a considerable number of species establishing in new geographic areas. Many organisms have taken advantage of human-made environments, including buildings. One such species is the dry rot fungus Serpula lacrymans, which is the most aggressive wood-decay fungus in indoor environments in temperate regions. Using population genomic analyses of 36 full genome sequenced isolates, we demonstrated that European and Japanese isolates are highly divergent and the populations split 3000-19,000 generations ago, probably predating human influence. Approximately 250 generations ago, the European population went through a tight bottleneck, probably corresponding to the fungus colonization of the built environment in Europe. The demographic history of these populations, probably lead to low adaptive potential. Only two loci under selection were identified using a Fst outlier approach, and selective sweep analyses identified three loci with extended haplotype homozygosity. The selective sweep analyses found signals in genes possibly related to decay of various substrates in Japan and in genes involved DNA replication and protein modification in Europe. Our results suggest that the dry rot fungus independently established in indoor environments in Europe and Japan and that invasive species can potentially establish large populations in new habitats based on a few colonizing individuals.

Altitudinal upwards shifts in fungal fruiting in the Alps.

Many plant and animal species are changing their latitudinal and/or altitudinal distributions in response to climate change, but whether fungi show similar changes is largely unknown. Here, we use historical fungal fruit body records from the European Alps to assess altitudinal changes in fungal fruiting between 1960 and 2010. We observe that many fungal species are fruiting at significantly higher elevations in 2010 compared to 1960, and especially so among soil-dwelling fungi. Wood-decay fungi, being dependent on the presence of one or a few host trees, show a slower response. Species growing at higher elevations changed their altitudinal fruiting patterns significantly more than lowland species. Environmental changes in high altitudes may lead to proportionally stronger responses, since high-altitude species live closer to their physiological limit. These aboveground changes in fruiting patterns probably mirror corresponding shifts in belowground fungal communities, suggesting parallel shifts in important ecosystem functions.

In vitro evidence of root colonization suggests ecological versatility in the genus Mycena.

The root-associated habit has evolved on numerous occasions in different fungal lineages, suggesting a strong evolutionary pressure for saprotrophic fungi to switch to symbiotic associations with plants. Species within the ubiquitous, saprotrophic genus Mycena are frequently major components in molecular studies of root-associated fungal communities, suggesting that an evaluation of their trophic status is warranted. Here, we report on interactions between a range of Mycena species and the plant Betula pendula. In all, 17 Mycena species were inoculated onto B. pendula seedlings. Physical interactions between hyphae and fine roots were examined using differential staining and fluorescence microscopy. Physiological interactions were investigated using 14 C and 32 P to show potential transfer between symbionts. All Mycena species associated closely with fine roots, showing hyphal penetration into the roots, which in some cases were intracellular. Seven species formed mantle-like structures around root tips, but none formed a Hartig net. Mycena pura and Mycena galopus both enhanced seedling growth, with M. pura showing significant transfer of 32 P to the seedlings. Our results support the view that several Mycena species can associate closely with plant roots and some may potentially occupy a transitional state between saprotrophy and biotrophy.

The Inhalable Mycobiome of Sawmill Workers: Exposure Characterization and Diversity.

Exposure to fungal spores has been associated with respiratory symptoms and allergic alveolitis among sawmill workers, but the complexity of sawmill workers' fungal exposure has been poorly studied. We characterized the fungal diversity in air samples from sawmill workers' breathing zones and identified differences in the richness, diversity, and taxonomic composition between companies, departments, wood types, and seasons. Full-shift personal inhalable dust samples (n = 86) collected from 11 industrial sawmill, sorting mill, and planer mill companies processing spruce and/or pine were subjected to DNA metabarcoding using the fungal internal transcribed spacer (ITS) region 2. The workers were exposed to a higher total number of operational taxonomic units (OTUs) in summer than in winter and when processing spruce than when processing pine. Workers in the saw department had the richest fungal exposure, followed by workers in the planing department and sorting of dry timber department. Sawmills explained 11% of the variation in the fungal community composition of the exposure, followed by season (5%) and department (3%). The fungal compositions of the exposures also differed between seasons, sawmills, wood types, and departments at the taxonomic level, ranging from the phylum to the species level. The differences in exposure diversity suggest that the potential health effects of fungal inhalation may also be different; hence, a risk assessment based on the fungal diversity differences should be performed. This study may serve as a basis for establishing a fungal profile of signature species that are specific for sawmills and that can be measured quantitatively in future risk assessments of sawmill workers.IMPORTANCE To gain more knowledge about exposure-response relationships, it is important to improve exposure characterization by comprehensively identifying the temporal and spatial fungal composition and diversity of inhalable dust at workplaces. The variation in the diverse fungal communities to which individuals are exposed in different seasons and sawmills suggests that variations in exposure-related health effects between seasons and companies can be expected. More importantly, the distinct fungal profiles between departments across companies indicate that workers in different job groups are differently exposed and that health risks can be department specific. DNA metabarcoding provides insight into a broad spectrum of airborne fungi that may serve as a basis for obtaining important knowledge about the fungi to which workers are exposed.

A single ectomycorrhizal plant root system includes a diverse and spatially structured fungal community.

Although only a relatively small proportion of plant species form ectomycorrhizae with fungi, it is crucial for growth and survival for a number of widespread woody plant species. Few studies have attempted to investigate the fine scale spatial structure of entire root systems of adult ectomycorrhizal (EcM) plants. Here, we use the herbaceous perennial Bistorta vivipara to map the entire root system of an adult EcM plant and investigate the spatial structure of its root-associated fungi. All EcM root tips were sampled, mapped and identified using a direct PCR approach and Sanger sequencing of the internal transcribed spacer region. A total of 32.1% of all sampled root tips (739 of 2302) were successfully sequenced and clustered into 41 operational taxonomic units (OTUs). We observed a clear spatial structuring of the root-associated fungi within the root system. Clusters of individual OTUs were observed in the younger parts of the root system, consistent with observations of priority effects in previous studies, but were absent from the older parts of the root system. This may suggest a succession and fragmentation of the root-associated fungi even at a very fine scale, where competition likely comes into play at different successional stages within the root system.

Revealing hidden insect-fungus interactions; moderately specialized, modular and anti-nested detritivore networks.

Ecological networks are composed of interacting communities that influence ecosystem structure and function. Fungi are the driving force for ecosystem processes such as decomposition and carbon sequestration in terrestrial habitats, and are strongly influenced by interactions with invertebrates. Yet, interactions in detritivore communities have rarely been considered from a network perspective. In the present study, we analyse the interaction networks between three functional guilds of fungi and insects sampled from dead wood. Using DNA metabarcoding to identify fungi, we reveal a diversity of interactions differing in specificity in the detritivore networks, involving three guilds of fungi. Plant pathogenic fungi were relatively unspecialized in their interactions with insects inhabiting dead wood, while interactions between the insects and wood-decay fungi exhibited the highest degree of specialization, which was similar to estimates for animal-mediated seed dispersal networks in previous studies. The low degree of specialization for insect symbiont fungi was unexpected. In general, the pooled insect-fungus networks were significantly more specialized, more modular and less nested than randomized networks. Thus, the detritivore networks had an unusual anti-nested structure. Future studies might corroborate whether this is a common aspect of networks based on interactions with fungi, possibly owing to their often intense competition for substrate.

Explaining European fungal fruiting phenology with climate variability.

Here we assess the impact of geographically dependent (latitude, longitude, and altitude) changes in bioclimatic (temperature, precipitation, and primary productivity) variability on fungal fruiting phenology across Europe. Two main nutritional guilds of fungi, saprotrophic and ectomycorrhizal, were further separated into spring and autumn fruiters. We used a path analysis to investigate how biogeographic patterns in fungal fruiting phenology coincided with seasonal changes in climate and primary production. Across central to northern Europe, mean fruiting varied by approximately 25 d, primarily with latitude. Altitude affected fruiting by up to 30 d, with spring delays and autumnal accelerations. Fruiting was as much explained by the effects of bioclimatic variability as by their large-scale spatial patterns. Temperature drove fruiting of autumnal ectomycorrhizal and saprotrophic groups as well as spring saprotrophic groups, while primary production and precipitation were major drivers for spring-fruiting ectomycorrhizal fungi. Species-specific phenology predictors were not stable, instead deviating from the overall mean. There is significant likelihood that further climatic change, especially in temperature, will impact fungal phenology patterns at large spatial scales. The ecological implications are diverse, potentially affecting food webs (asynchrony), nutrient cycling and the timing of nutrient availability in ecosystems.