Habitat is more important than climate for structuring soil fungal communities associated in truffle sites.

The ectomycorrhizal fungi Tuber melanosporum Vittad. and Tuber aestivum Vittad. produce highly valuable truffles, but little is known about the soil fungal communities associated with these truffle species in places where they co-occur. Here, we compared soil fungal communities present in wild and planted truffle sites, in which T. melanosporum and T. aestivum coexist, in Mediterranean and temperate regions over three sampling seasons spanning from 2018 to 2019. We showed that soil fungal community composition and ectomycorrhizal species composition are driven by habitat type rather than climate regions. Also, we observed the influence of soil pH, organic matter content and C:N ratio structuring total and ectomycorrhizal fungal assemblages. Soil fungal communities in wild sites revealed more compositional variability than those of plantations. Greater soil fungal diversity was found in temperate compared to Mediterranean sites when considering all fungal guilds. Ectomycorrhizal diversity was significantly higher in wild sites compared to plantations. Greater mould abundance at wild sites than those on plantation was observed while tree species and seasonal effects were not significant predictors in fungal community structure. Our results suggested a strong influence of both ecosystem age and management on the fungal taxa composition in truffle habitats.

Contrasting plant-soil-microbial feedbacks stabilize vegetation types and uncouple topsoil C and N stocks across a subarctic-alpine landscape.

Global vegetation regimes vary in belowground carbon (C) and nitrogen (N) dynamics. However, disentangling large-scale climatic controls from the effects of intrinsic plant-soil-microbial feedbacks on belowground processes is challenging. In local gradients with similar pedo-climatic conditions, effects of plant-microbial feedbacks may be isolated from large-scale drivers. Across a subarctic-alpine mosaic of historic grazing fields and surrounding heath and birch forest, we evaluated whether vegetation-specific plant-microbial feedbacks involved contrasting N cycling characteristics and C and N stocks in the organic topsoil. We sequenced soil fungi, quantified functional genes within the inorganic N cycle, and measured 15 N natural abundance. In grassland soils, large N stocks and low C : N ratios associated with fungal saprotrophs, archaeal ammonia oxidizers, and bacteria capable of respiratory ammonification, indicating maintained inorganic N cycling a century after abandoned reindeer grazing. Toward forest and heath, increasing abundance of mycorrhizal fungi co-occurred with transition to organic N cycling. However, ectomycorrhizal fungal decomposers correlated with small soil N and C stocks in forest, while root-associated ascomycetes associated with small N but large C stocks in heath, uncoupling C and N storage across vegetation types. We propose that contrasting, positive plant-microbial feedbacks stabilize vegetation trajectories, resulting in diverging soil C : N ratios at the landscape scale.

Silvicultural management and altitude prevail on soil properties and fungal community in shaping understorey plant communities in a Mediterranean pine forest.

Understorey vegetation plays a key role in Mediterranean forest ecosystem functioning. However, we still lack a thorough understanding of the patterns and drivers of understorey composition and diversity. As a result, understoreys are often ignored during assessments of forest functioning under climate change. Here we studied the effect of silvicultural management, topography, soil fungal community composition and soil physical and chemical properties on understorey community composition and diversity. The plant cover and number of individuals of understorey perennial plants, shrubs and non-dominant trees was recorded on 24 plots (paired: control-thinned) in a Mediterranean pine-dominated mountainous area in Northeast Spain. The study area represented a broad thinning intensity gradient (from 0 to 70 % in removed stand basal area) along a 400-m altitudinal range (from 609 m to 1013 m). Our results showed that thinning intensity and topography explained the greatest proportion of the total variance in the understorey species composition, i.e., 18 % and 16 %, respectively. Interestingly, the effects of the silvicultural treatments were significant only when considering the altitudinal effect, so that, the main impacts of thinning on the understorey community composition occurred at low altitudes (between 609 m and 870 m). Moreover, we found a significant decrease in both richness and abundance of understorey species in both the control and thinned plots with increasing altitude, with thinned plots being significantly richer in species compared to the control plots. The difference in the understorey community sensitivity to forest thinning along the altitudinal gradient suggests changes in factors that limit plant growth. Low elevation plots were restrained by light availability while high altitudes plots limited by winter freezing temperature.

Lack of Phylogenetic Differences in Ectomycorrhizal Fungi among Distinct Mediterranean Pine Forest Habitats.

Understanding whether the occurrences of ectomycorrhizal species in a given tree host are phylogenetically determined can help in assessing different conservational needs for each fungal species. In this study, we characterized ectomycorrhizal phylogenetic composition and phylogenetic structure in 42 plots with five different Mediterranean pine forests: i.e., pure forests dominated by P. nigra, P. halepensis, and P. sylvestris, and mixed forests of P. nigra-P. halepensis and P. nigra-P. sylvestris, and tested whether the phylogenetic structure of ectomycorrhizal communities differs among these. We found that ectomycorrhizal communities were not different among pine tree hosts neither in phylogenetic composition nor in structure and phylogenetic diversity. Moreover, we detected a weak abiotic filtering effect (4%), with pH being the only significant variable influencing the phylogenetic ectomycorrhizal community, while the phylogenetic structure was slightly influenced by the shared effect of stand structure, soil, and geographic distance. However, the phylogenetic community similarity increased at lower pH values, supporting that fewer, closely related species were found at lower pH values. Also, no phylogenetic signal was detected among exploration types, although short and contact were the most abundant types in these forest ecosystems. Our results demonstrate that pH but not tree host, acts as a strong abiotic filter on ectomycorrhizal phylogenetic communities in Mediterranean pine forests at a local scale. Finally, our study shed light on dominant ectomycorrhizal foraging strategies in drought-prone ecosystems such as Mediterranean forests.

Sampling forest soils to describe fungal diversity and composition. Which is the optimal sampling size in mediterranean pure and mixed pine oak forests?

Soil sampling is a critical step affecting perceived fungal diversity, however sampling optimization for high-throughput-DNA sequencing studies have never been tested in Mediterranean forest ecosystems. We identified the minimum number of pooled samples needed to obtain a reliable description of fungal communities in terms of diversity and composition in three different Mediterranean forests (pine, oak, and mixed-pine-oak). Twenty soil samples were randomly selected in each of the three plots per type. Samples obtained in 100 m2 plots were pooled to obtain mixtures of 3, 6, 10, 15, 20 samples, and sequenced using Illumina MiSeq of fungal ITS2 amplicons. Pooling three soil samples in Pinus and Quercus stands provided consistent richness estimations, while at least six samples were needed in mixed-stands. ß-diversity decreased with increasing sample pools in monospecific-stands, while there was no effect of sample pool size on mixed-stands. Soil sample pooling had no effect over species composition. We estimate that three samples would be already optimal to describe fungal richness and composition in Mediterranean pure stands, while at least six samples would be needed in mixed stands.

Reindeer control over subarctic treeline alters soil fungal communities with potential consequences for soil carbon storage.

The climate-driven encroachment of shrubs into the Arctic is accompanied by shifts in soil fungal communities that could contribute to a net release of carbon from tundra soils. At the same time, arctic grazers are known to prevent the establishment of deciduous shrubs and, under certain conditions, promote the dominance of evergreen shrubs. As these different vegetation types associate with contrasting fungal communities, the belowground consequences of climate change could vary among grazing regimes. Yet, at present, the impact of grazing on soil fungal communities and their links to soil carbon have remained speculative. Here we tested how soil fungal community composition, diversity and function depend on tree vicinity and long-term reindeer grazing regime and assessed how the fungal communities relate to organic soil carbon stocks in an alpine treeline ecotone in Northern Scandinavia. We determined soil carbon stocks and characterized soil fungal communities directly underneath and >3 m away from mountain birches (Betula pubescens ssp. czerepanovii) in two adjacent 55-year-old grazing regimes with or without summer grazing by reindeer (Rangifer tarandus). We show that the area exposed to year-round grazing dominated by evergreen dwarf shrubs had higher soil C:N ratio, higher fungal abundance and lower fungal diversity compared with the area with only winter grazing and higher abundance of mountain birch. Although soil carbon stocks did not differ between the grazing regimes, stocks were positively associated with root-associated ascomycetes, typical to the year-round grazing regime, and negatively associated with free-living saprotrophs, typical to the winter grazing regime. These findings suggest that when grazers promote dominance of evergreen dwarf shrubs, they induce shifts in soil fungal communities that increase soil carbon sequestration in the long term. Thus, to predict climate-driven changes in soil carbon, grazer-induced shifts in vegetation and soil fungal communities need to be accounted for.

White mulch and irrigation increase black truffle soil mycelium when competing with summer truffle in young truffle orchards.

The black truffle (Tuber melanosporum Vittad.) and the summer truffle (Tuber aestivum Vittad.) are two of the most appreciated edible fungi worldwide. The natural distributions of both species partially overlap. However, the interspecific interactions between these truffles and how irrigation and mulching techniques impact the dynamics between them are still unknown. Here, an experimental truffle plantation with Quercus ilex was established in Maials (Catalonia, Spain), combining three soil mulch treatments (white mulch, black mulch and bare soil as a control) and two irrigation regimes (irrigated and non-irrigated as a control) to investigate truffle mycelial dynamics in soil when both truffle species co-occur. The development of truffle mycelium in two different seasons (spring and autumn) in two consecutive years (2017 and 2018) was quantified using qPCR. Truffle mycelia of both species showed greatest development under white mulch. When mycelia of both truffle species co-occurred in soil, irrigation combined with white mulch resulted in greater quantities of T. melanosporum mycelial biomass, whereas the control irrigation treatment favoured the development of T. aestivum. Mulch treatments were also advantageous for seedling growth, which was expressed as root collar diameter and its increment during the study period. Significant relationships between root collar diameter and root growth and the amount of mycelial biomass in the soil were observed for both truffle species. Our results indicate the potential advantages of using white mulch to support irrigation in truffle plantations located in areas with dry Mediterranean climatic conditions to promote the development of Tuber mycelium.

Production and turnover of mycorrhizal soil mycelium relate to variation in drought conditions in Mediterranean Pinus pinaster, Pinus sylvestris and Quercus ilex forests.

In forests, ectomycorrhizal mycelium is pivotal for driving soil carbon and nutrient cycles, but how ectomycorrhizal mycelial dynamics vary in ecosystems with drought periods is unknown. We quantified the production and turnover of mycorrhizal mycelium in Mediterranean Pinus pinaster, Pinus sylvestris and Quercus ilex forests and related the estimates to standardised precipitation index (SPI), to study how mycelial dynamics relates to tree species and drought-moisture conditions. Production and turnover of mycelium was estimated between July and February, by quantifying the fungal biomass (ergosterol) in ingrowth mesh bags and using statistical modelling. SPI for time scales of 1-3 months was calculated from precipitation records and precipitation data over the study period. Forests dominated by Pinus trees displayed higher biomass but were seasonally more variable, as opposed to Q. ilex forests where the mycelial biomass remained lower and stable over the season. Production and turnover, respectively, varied between 1.4-5.9 kg ha-1  d-1 and 7.2-9.9 times yr-1 over the different forest types and were positively correlated with 2-month and 3-month SPI over the study period. Our results demonstrated that mycorrhizal mycelial biomass varied with season and tree species and we speculate that production and turnover are related to physiology and plant host performance during drought.

Insect defoliation is linked to a decrease in soil ectomycorrhizal biomass and shifts in needle endophytic communities.

Insect outbreaks of increasing frequency and severity in forests are predicted due to climate change. Insect herbivory is known to promote physiological changes in forest trees. However, little is known about whether these plant phenotypic adjustments have cascading effects on tree microbial symbionts such as fungi in roots and foliage. We studied the impact of defoliation by the pine processionary moth in two infested Pinus nigra forests through a multilevel sampling of defoliated and non-defoliated trees. We measured tree growth, nutritional status and carbon allocation to chemical defenses. Simultaneously, we analysed the putative impact of defoliation on the needle endophytes and on the soil fungal communities. Higher concentrations of chemical defenses were found in defoliated trees, likely as a response to defoliation; however, no differences in non-structural carbohydrate reserves were found. In parallel to the reductions in tree growth and changes in chemical defenses, we observed shifts in the composition of needle endophytic and soil fungal communities in defoliated trees. Defoliated trees consistently corresponded with a lower biomass of ectomycorrhizal fungi in both sites, and a higher alpha diversity and greater relative abundance of belowground saprotrophs and pathogens. However, ectomycorrhizal alpha diversity was similar between non-defoliated and defoliated trees. Specific needle endophytes in old needles were strongly associated with non-defoliated trees. The potential role of these endophytic fungi in pine resistance should be further investigated. Our study suggests that lower biomass of ectomycorrhizal fungi in defoliated trees might slow down tree recovery since fungal shifts might affect tree-mycorrhizal feedbacks and can potentially influence carbon and nitrogen cycling in forest soils.

Optimized metabarcoding with Pacific biosciences enables semi-quantitative analysis of fungal communities.

Recent studies have questioned the use of high-throughput sequencing of the nuclear ribosomal internal transcribed spacer (ITS) region to derive a semi-quantitative representation of fungal community composition. However, comprehensive studies that quantify biases occurring during PCR and sequencing of ITS amplicons are still lacking. We used artificially assembled communities consisting of 10 ITS-like fragments of varying lengths and guanine-cytosine (GC) contents to evaluate and quantify biases during PCR and sequencing with Illumina MiSeq, PacBio RS II and PacBio Sequel I technologies. Fragment length variation was the main source of bias in observed community composition relative to the template, with longer fragments generally being under-represented for all sequencing platforms. This bias was three times higher for Illumina MiSeq than for PacBio RS II and Sequel I. All 10 fragments in the artificial community were recovered when sequenced with PacBio technologies, whereas the three longest fragments (> 447 bases) were lost when sequenced with Illumina MiSeq. Fragment length bias also increased linearly with increasing number of PCR cycles but could be mitigated by optimization of the PCR setup. No significant biases related to GC content were observed. Despite lower sequencing output, PacBio sequencing was better able to reflect the community composition of the template than Illumina MiSeq sequencing.

Soil microclimate changes affect soil fungal communities in a Mediterranean pine forest.

Soil microclimate is a potentially important regulator of the composition of plant-associated fungal communities in climates with significant drought periods. Here, we investigated the spatio-temporal dynamics of soil fungal communities in a Mediterranean Pinus pinaster forest in relation to soil moisture and temperature. Fungal communities in 336 soil samples collected monthly over 1 year from 28 long-term experimental plots were assessed by PacBio sequencing of ITS2 amplicons. Total fungal biomass was estimated by analysing ergosterol. Community changes were analysed in the context of functional traits. Soil fungal biomass was lowest during summer and late winter and highest during autumn, concurrent with a greater relative abundance of mycorrhizal species. Intra-annual spatio-temporal changes in community composition correlated significantly with soil moisture and temperature. Mycorrhizal fungi were less affected by summer drought than free-living fungi. In particular, mycorrhizal species of the short-distance exploration type increased in relative abundance under dry conditions, whereas species of the long-distance exploration type were more abundant under wetter conditions. Our observations demonstrate a potential for compositional and functional shifts in fungal communities in response to changing climatic conditions. Free-living fungi and mycorrhizal species with extensive mycelia may be negatively affected by increasing drought periods in Mediterranean forest ecosystems.

Diplodia Tip Blight on Its Way to the North: Drivers of Disease Emergence in Northern Europe.

Disease emergence in northern and boreal forests has been mostly due to tree-pathogen encounters lacking a co-evolutionary past. However, outbreaks involving novel interactions of the host or the pathogen with the environment have been less well documented. Following an increase of records in Northern Europe, the first large outbreak of Diplodia sapinea on Pinus sylvestris was discovered in Sweden in 2016. By reconstructing the development of the epidemic, we found that the attacks started approx. 10 years back from several isolated trees in the stand and ended up affecting almost 90% of the trees in 2016. Limited damage was observed in other plantations in the surroundings of the affected stand, pointing to a new introduced pathogen as the cause of the outbreak. Nevertheless, no genetic differences based on SSR markers were found between isolates of the outbreak area and other Swedish isolates predating the outbreak or from other populations in Europe and Asia Minor. On a temporal scale, we saw that warm May and June temperatures were associated with higher damage and low tree growth, while cold and rainy conditions seemed to favor growth and deter disease. At a spatial scale, we saw that spread occurred predominantly in the SW aspect-area of the stand. Within that area and based on tree-ring and isotope (δ13C) analyses, we saw that disease occurred on trees that over the years had shown a lower water-use efficiency (WUE). Spore traps showed that highly infected trees were those producing the largest amount of inoculum. D. sapinea impaired latewood growth and reduced C reserves in needles and branches. D. sapinea attacks can cause serious economic damage by killing new shoots, disrupting the crown, and affecting the quality of stems. Our results show that D. sapinea has no limitations in becoming a serious pathogen in Northern Europe. Management should focus on reducing inoculum, especially since climate change may bring more favorable conditions for this pathogen. Seedlings for planting should be carefully inspected as D. sapinea may be present in a latent stage in asymptomatic tissues.

Mushroom Emergence Detected by Combining Spore Trapping with Molecular Techniques.

Obtaining reliable and representative mushroom production data requires time-consuming sampling schemes. In this paper, we assessed a simple methodology to detect mushroom emergence by trapping the fungal spores of the fruiting body community in plots where mushroom production was determined weekly. We compared the performance of filter paper traps with that of funnel traps and combined these spore trapping methods with species-specific quantitative real-time PCR and Illumina MiSeq to determine the spore abundance. Significantly more MiSeq proportional reads were generated for both ectomycorrhizal and saprotrophic fungal species using filter traps than were obtained using funnel traps. The spores of 37 fungal species that produced fruiting bodies in the study plots were identified. Spore community composition changed considerably over time due to the emergence of ephemeral fruiting bodies and rapid spore deposition (lasting from 1 to 2 weeks), which occurred in the absence of rainfall events. For many species, the emergence of epigeous fruiting bodies was followed by a peak in the relative abundance of their airborne spores. There were significant positive relationships between fruiting body yields and spore abundance in time for five of seven fungal species. There was no relationship between fruiting body yields and their spore abundance at plot level, indicating that some of the spores captured in each plot were arriving from the surrounding areas. Differences in fungal detection capacity by spore trapping may indicate different dispersal ability between fungal species. Further research can help to identify the spore rain patterns for most common fungal species.IMPORTANCE Mushroom monitoring represents a serious challenge in economic and logistical terms because sampling approaches demand extensive field work at both the spatial and temporal scales. In addition, the identification of fungal taxa depends on the expertise of experienced fungal taxonomists. Similarly, the study of fungal dispersal has been constrained by technological limitations, especially because the morphological identification of spores is a challenging and time-consuming task. Here, we demonstrate that spores from ectomycorrhizal and saprotrophic fungal species can be identified using simple spore traps together with either MiSeq fungus-specific amplicon sequencing or species-specific quantitative real-time PCR. In addition, the proposed methodology can be used to characterize the airborne fungal community and to detect mushroom emergence in forest ecosystems.

Soil drying procedure affects the DNA quantification of Lactarius vinosus but does not change the fungal community composition.

Drying soil samples before DNA extraction is commonly used for specific fungal DNA quantification and metabarcoding studies, but the impact of different drying procedures on both the specific fungal DNA quantity and the fungal community composition has not been analyzed. We tested three different drying procedures (freeze-drying, oven-drying, and room temperature) on 12 different soil samples to determine (a) the soil mycelium biomass of the ectomycorrhizal species Lactarius vinosus using qPCR with a specifically designed TaqMan® probe and (b) the fungal community composition and diversity using the PacBio® RS II sequencing platform. Mycelium biomass of L. vinosus was significantly greater in the freeze-dried soil samples than in samples dried at oven and room temperature. However, drying procedures had no effect on fungal community composition or on fungal diversity. In addition, there were no significant differences in the proportions of fungi according to their functional roles (moulds vs. mycorrhizal species) in response to drying procedures. Only six out of 1139 operational taxonomic units (OTUs) had increased their relative proportions after soil drying at room temperature, with five of these OTUs classified as mould or yeast species. However, the magnitude of these changes was small, with an overall increase in relative abundance of these OTUs of approximately 2 %. These results suggest that DNA degradation may occur especially after drying soil samples at room temperature, but affecting equally nearly all fungi and therefore causing no significant differences in diversity and community composition. Despite the minimal effects caused by the drying procedures at the fungal community composition, freeze-drying resulted in higher concentrations of L. vinosus DNA and prevented potential colonization from opportunistic species.