Fungi play important roles in forest ecosystems and understanding fungal diversity is crucial to address essential questions about species conservation and ecosystems management. Changes in fungal diversity can have severe impacts on ecosystem functionality. Unfortunately, little is known about fungal diversity in northern temperate and boreal forests, and we have yet to understand how abiotic variables shape fungal richness and composition. Our objectives were to make an overview of the fungal richness and the community composition in the region and identify their major abiotic drivers. We sampled 262 stands across the northern temperate and boreal Quebec forest located in the region of Abitibi-Témiscamingue, Mauricie, and Haute-Mauricie. At each site, we characterized fungal composition using Illumina sequencing, as well as several potential abiotic drivers (e.g., humus thickness, soil pH, vegetation cover, etc.). We tested effects of abiotic drivers on species richness using generalized linear models, while difference in fungal composition between stands was analyzed with permutational multivariate analysis of variance and beta-diversity partitioning analyses. Fungi from the order Agaricales, Helotiales, and Russulales were the most frequent and sites from the north of Abitibi-Témiscamingue showed the highest OTUs (Operational Taxonomic Unit) richness. Stand age and moss cover were the best predictors of fungal richness. On the other hand, the strongest drivers of fungal community structure were soil pH, average cumulative precipitation, and stand age, although much of community variance was left unexplained in our models. Overall, our regional metacommunity was characterized by high turnover rate, even when rare OTUs were removed. This may indicate strong environmental filtering by several unmeasured abiotic filters, or stronger than expected dispersal limitations in soil fungal communities. Our results show how difficult it can be to predict fungal community assembly even with high replication and efforts to include several biologically relevant explanatory variables.

Boreal forests provide important ecosystem services, most notably being the mitigation of increasing atmospheric CO2 emissions. Microbial biodiversity, particularly the local diversity of fungi, has been shown to promote multiple functions of the boreal forests of Northeastern China. However, this microbial biodiversity-multifunctionality relationship has yet to be explored in Northeastern Canada, where historical environment have shaped a different regional pool of microbial diversity. This study focuses on the relationship between the soil microbiome and ecosystem multifunctionality, as well as the influence of pH and redox potential (Eh) on the regulation of such relationship. Structural equation modelling (SEM) was used to explore the different causal relationships existing in the studied ecosystems. In a managed part of the Canadian boreal forest, 156 forest polygons were sampled to (1) estimate the α- and β-diversity of fungal and bacterial communities and (2) measure 12 ecosystem functions mainly related to soil nutrient storage and cycling. Both bacteria and fungi influenced ecosystem multifunctionality, but on their own respective functions. Bacterial β-diversity was the most important factor increasing primary productivity and soil microbial biomass, while reducing soil emitted atmospheric CO2. Environmental characteristics, particularly low levels of organic matter in soil, were shown to have the strongest positive impact on boreal ecosystem multifunctionality. Overall, our results were consistent with those obtained in Northeastern China; however, some differences need to be further explored especially considering the history of forest management in Northeastern Canada.