The biological controls of soil carbon accumulation following wildfire and harvest in boreal forests: A review.

Boreal forests are frequently subjected to disturbances, including wildfire and clear-cutting. While these disturbances can cause soil carbon (C) losses, the long-term accumulation dynamics of soil C stocks during subsequent stand development is controlled by biological processes related to the balance of net primary production (NPP) and outputs via heterotrophic respiration and leaching, many of which remain poorly understood. We review the biological processes suggested to influence soil C accumulation in boreal forests. Our review indicates that median C accumulation rates following wildfire and clear-cutting are similar (0.15 and 0.20 Mg ha-1 year-1, respectively), however, variation between studies is extremely high. Further, while many individual studies show linear increases in soil C stocks through time after disturbance, there are indications that C stock recovery is fastest early to mid-succession (e.g. 15-80 years) and then slows as forests mature (e.g. >100 years). We indicate that the rapid build-up of soil C in younger stands appears not only driven by higher plant production, but also by a high rate of mycorrhizal hyphal production, and mycorrhizal suppression of saprotrophs. As stands mature, the balance between reductions in plant and mycorrhizal production, increasing plant litter recalcitrance, and ectomycorrhizal decomposers and saprotrophs have been highlighted as key controls on soil C accumulation rates. While some of these controls appear well understood (e.g. temporal patterns in NPP, changes in aboveground litter quality), many others remain research frontiers. Notably, very little data exists describing and comparing successional patterns of root production, mycorrhizal functional traits, mycorrhizal-saprotroph interactions, or C outputs via heterotrophic respiration and dissolved organic C following different disturbances. We argue that these less frequently described controls require attention, as they will be key not only for understanding ecosystem C balances, but also for representing these dynamics more accurately in soil organic C and Earth system models.

Modelling optimal ligninolytic activity during plant litter decomposition.

A large fraction of plant litter comprises recalcitrant aromatic compounds (lignin and other phenolics). Quantifying the fate of aromatic compounds is difficult, because oxidative degradation of aromatic carbon (C) is a costly but necessary endeavor for microorganisms, and we do not know when gains from the decomposition of aromatic C outweigh energetic costs. To evaluate these tradeoffs, we developed a litter decomposition model in which the aromatic C decomposition rate is optimized dynamically to maximize microbial growth for the given costs of maintaining ligninolytic activity. We tested model performance against > 200 litter decomposition datasets collected from published literature and assessed the effects of climate and litter chemistry on litter decomposition. The model predicted a time-varying ligninolytic oxidation rate, which was used to calculate the lag time before the decomposition of aromatic C is initiated. Warmer conditions increased decomposition rates, shortened the lag time of aromatic C oxidation, and improved microbial C-use efficiency by decreasing the costs of oxidation. Moreover, a higher initial content of aromatic C promoted an earlier start of aromatic C decomposition under any climate. With this contribution, we highlight the application of eco-evolutionary approaches based on optimized microbial life strategies as an alternative parametrization scheme for litter decomposition models.

Ectomycorrhizal fungi are more sensitive to high soil nitrogen levels in forests exposed to nitrogen deposition.

Ectomycorrhizal fungi are essential for nitrogen (N) cycling in many temperate forests and responsive to anthropogenic N addition, which generally decreases host carbon (C) allocation to the fungi. In the boreal region, however, ectomycorrhizal fungal biomass has been found to correlate positively with soil N availability. Still, responses to anthropogenic N input, for instance through atmospheric deposition, are commonly negative. To elucidate whether variation in N supply affects ectomycorrhizal fungi differently depending on geographical context, we investigated ectomycorrhizal fungal communities along fertility gradients located in two nemo-boreal forest regions with similar ranges in soil N : C ratios and inorganic N availability but contrasting rates of N deposition. Ectomycorrhizal biomass and community composition remained relatively stable across the N gradient with low atmospheric N deposition, but biomass decreased and the community changed more drastically with increasing N availability in the gradient subjected to higher rates of N deposition. Moreover, potential activities of enzymes involved in ectomycorrhizal mobilisation of organic N decreased as N availability increased. In forests with low external input, we propose that stabilising feedbacks in tree-fungal interactions maintain ectomycorrhizal fungal biomass and communities even in N-rich soils. By contrast, anthropogenic N input seems to impair ectomycorrhizal functions.

Fire severity as a key determinant of aboveground and belowground biological community recovery in managed even-aged boreal forests.

Changes in fire regime of boreal forests in response to climate warming are expected to impact postfire recovery. However, quantitative data on how managed forests sustain and recover from recent fire disturbance are limited.Two years after a large wildfire in managed even-aged boreal forests in Sweden, we investigated how recovery of aboveground and belowground communities, that is, understory vegetation and soil microbial and faunal communities, responded to variation in the severity of soil (i.e., consumption of soil organic matter) and canopy fires (i.e., tree mortality).While fire overall enhanced diversity of understory vegetation through colonization of fire adapted plant species, it reduced the abundance and diversity of soil biota. We observed contrasting effects of tree- and soil-related fire severity on survival and recovery of understory vegetation and soil biological communities. Severe fires that killed overstory Pinus sylvestris promoted a successional stage dominated by the mosses Ceratodon purpureus and Polytrichum juniperinum, but reduced regeneration of tree seedlings and disfavored the ericaceous dwarf-shrub Vaccinium vitis-idaea and the grass Deschampsia flexuosa. Moreover, high tree mortality from fire reduced fungal biomass and changed fungal community composition, in particular that of ectomycorrhizal fungi, and reduced the fungivorous soil Oribatida. In contrast, soil-related fire severity had little impact on vegetation composition, fungal communities, and soil animals. Bacterial communities responded to both tree- and soil-related fire severity. Synthesis: Our results 2 years postfire suggest that a change in fire regime from a historically low-severity ground fire regime, with fires that mainly burns into the soil organic layer, to a stand-replacing fire regime with a high degree of tree mortality, as may be expected with climate change, is likely to impact the short-term recovery of stand structure and above- and belowground species composition of even-aged P. sylvestris boreal forests.

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.

Sample Preparation for Fungal Community Analysis by High-Throughput Sequencing of Barcode Amplicons.

Fungal species participate in vast numbers of processes in the landscape around us. However, their cryptic mycelial growth, inside various substrates and in highly diverse species assemblages, has been a major obstacle to thorough analysis of fungal communities, hampering exhaustive description of the fungal kingdom. Technological developments allowing rapid, high-throughput sequencing of mixed communities from many samples at once are currently having a tremendous impact in fungal community ecology. Universal DNA extraction followed by amplification and sequencing of fungal species-level barcodes such as the nuclear internal transcribed spacer (ITS) region now enables identification and relative quantification of fungal community members across well-replicated experimental settings.Here, we present the sample preparation procedure presently used in our laboratory for fungal community analysis by high-throughput sequencing of amplified ITS2 markers. We focus on the procedure optimized for studies of total fungal communities in humus-rich soils, wood, and litter. However, this procedure can be applied to other sample types and markers. We focus on the laboratory-based part of sample preparation, i.e., the procedure from the point where samples enter the laboratory until amplicons are submitted for sequencing. Our procedure comprises four main parts: (1) universal DNA extraction, (2) optimization of PCR conditions, (3) production of tagged ITS amplicons, and (4) preparation of the multiplexed amplicon pool to be sequenced. The presented procedure is independent of the specific high-throughput sequencing technology used, which makes it highly versatile.

Do ectomycorrhizal exploration types reflect mycelial foraging strategies?

Ectomycorrhizal exploration types are commonly assumed to denote spatial foraging patterns and resource-related niches of extraradical mycelia. However, empirical evidence of the consistency of foraging strategies within exploration types is lacking. Here, we analysed ectomycorrhizal foraging patterns by incubating root-excluding ingrowth mesh bags filled with six different substrates in mature Picea abies forests. High-throughput sequencing was used to characterise ectomycorrhizal fungal communities in the mesh bags and on adjacent fine roots after one growing season. Contrary to expectations, many ectomycorrhizal genera of exploration types that are thought to produce little extraradical mycelium colonised ingrowth bags extensively, whereas genera commonly associated with ample mycelial production occurred sparsely in ingrowth bags relative to their abundance on roots. Previous assumptions about soil foraging patterns of exploration types do not seem to hold. Instead, we propose that variation in the proliferation of extraradical mycelium is related to intergeneric differences in mycelial longevity and the mobility of targeted resources.

Ericoid shrubs shape fungal communities and suppress organic matter decomposition in boreal forests.

Mycorrhizal fungi associated with boreal trees and ericaceous shrubs are central actors in organic matter (OM) accumulation through their belowground carbon allocation, their potential capacity to mine organic matter for nitrogen (N) and their ability to suppress saprotrophs. Yet, interactions between co-occurring ectomycorrhizal fungi (EMF), ericoid mycorrhizal fungi (ERI), and saprotrophs are poorly understood. We used a long-term (19 yr) plant functional group manipulation experiment with removals of tree roots, ericaceous shrubs and mosses and analysed the responses of different fungal guilds (assessed by metabarcoding) and their interactions in relation to OM quality (assessed by mid-infrared spectroscopy and nuclear magnetic resonance) and decomposition (litter mesh-bags) across a 5000-yr post-fire boreal forest chronosequence. We found that the removal of ericaceous shrubs and associated ERI changed the composition of EMF communities, with larger effects occurring at earlier stages of the chronosequence. Removal of shrubs was associated with enhanced N availability, litter decomposition and enrichment of the recalcitrant OM fraction. We conclude that increasing abundance of slow-growing ericaceous shrubs and the associated fungi contributes to increasing nutrient limitation, impaired decomposition and progressive OM accumulation in boreal forests, particularly towards later successional stages. These results are indicative of the contrasting roles of EMF and ERI in regulating belowground OM storage.

Transcriptomic markers of fungal growth, respiration and carbon-use efficiency.

Fungal metabolic carbon acquisition and its subsequent partitioning between biomass production and respiration, i.e. the carbon-use efficiency (CUE), are central parameters in biogeochemical modeling. However, current available techniques for estimating these parameters are all associated with practical and theoretical shortcomings, making assessments unreliable. Gene expression analyses hold the prospect of phenotype prediction by indirect means, providing new opportunities to obtain information about metabolic priorities. We cultured four different fungal isolates (Chalara longipes, Laccaria bicolor, Serpula lacrymans and Trichoderma harzianum) in liquid media with contrasting nitrogen availability and measured growth rates and respiration to calculate CUE. By relating gene expression markers to measured carbon fluxes, we identified genes coding for 1,3-ß-glucan synthase and 2-oxoglutarate dehydrogenase as suitable markers for growth and respiration, respectively, capturing both intraspecific variation as well as within-strain variation dependent on growth medium. A transcript index based on these markers correlated significantly with differences in CUE between the fungal isolates. Our study paves the way for the use of these markers to assess differences in growth, respiration and CUE in natural fungal communities, using metatranscriptomic or the RT-qPCR approach.

Climate warming dominates over plant genotype in shaping the seasonal trajectory of foliar fungal communities on oak.

Leaves interact with a wealth of microorganisms. Among these, fungi are highly diverse and are known to contribute to plant health, leaf senescence and early decomposition. However, patterns and drivers of the seasonal dynamics of foliar fungal communities are poorly understood. We used a multifactorial experiment to investigate the influence of warming and tree genotype on the foliar fungal community on the pedunculate oak Quercus robur across one growing season. Fungal species richness increased, evenness tended to decrease, and community composition strongly shifted during the growing season. Yeasts increased in relative abundance as the season progressed, while putative fungal pathogens decreased. Warming decreased species richness, reduced evenness and changed community composition, especially at the end of the growing season. Warming also negatively affected putative fungal pathogens. We only detected a minor imprint of tree genotype and warming × genotype interactions on species richness and community composition. Overall, our findings demonstrate that warming plays a larger role than plant genotype in shaping the seasonal dynamics of the foliar fungal community on oak. These warming-induced shifts in the foliar fungal community may have a pronounced impact on plant health, plant-fungal interactions and ecosystem functions.

A group of ectomycorrhizal fungi restricts organic matter accumulation in boreal forest.

Boreal forest soils are important global carbon sinks, with significant storage in the organic topsoil. Decomposition of these stocks requires oxidative enzymes, uniquely produced by fungi. Across Swedish boreal forests, we found that local carbon storage in the organic topsoil was 33% lower in the presence of a group of closely related species of ectomycorrhizal fungi - Cortinarius acutus s.l.. This observation challenges the prevailing view that ectomycorrhizal fungi generally act to increase carbon storage in soils but supports the idea that certain ectomycorrhizal fungi can complement free-living decomposers, maintaining organic matter turnover, nutrient cycling and tree productivity under nutrient-poor conditions. The indication that a narrow group of fungi may exert a major influence on carbon cycling questions the prevailing dogma of functional redundancy among microbial decomposers. Cortinarius acutus s.l. responds negatively to stand-replacing disturbance, and associated population declines are likely to increase soil carbon sequestration while impeding long-term nutrient cycling.

A tipping point in carbon storage when forest expands into tundra is related to mycorrhizal recycling of nitrogen.

Tundra ecosystems are global belowground sinks for atmospheric CO2 . Ongoing warming-induced encroachment by shrubs and trees risks turning this sink into a CO2 source, resulting in a positive feedback on climate warming. To advance mechanistic understanding of how shifts in mycorrhizal types affect long-term carbon (C) and nitrogen (N) stocks, we studied small-scale soil depth profiles of fungal communities and C-N dynamics across a subarctic-alpine forest-heath vegetation gradient. Belowground organic stocks decreased abruptly at the transition from heath to forest, linked to the presence of certain tree-associated ectomycorrhizal fungi that contribute to decomposition when mining N from organic matter. In contrast, ericoid mycorrhizal plants and fungi were associated with organic matter accumulation and slow decomposition. If climatic controls on arctic-alpine forest lines are relaxed, increased decomposition will likely outbalance increased plant productivity, decreasing the overall C sink capacity of displaced tundra.

The Taxon Hypothesis Paradigm-On the Unambiguous Detection and Communication of Taxa.

Here, we describe the taxon hypothesis (TH) paradigm, which covers the construction, identification, and communication of taxa as datasets. Defining taxa as datasets of individuals and their traits will make taxon identification and most importantly communication of taxa precise and reproducible. This will allow datasets with standardized and atomized traits to be used digitally in identification pipelines and communicated through persistent identifiers. Such datasets are particularly useful in the context of formally undescribed or even physically undiscovered species if data such as sequences from samples of environmental DNA (eDNA) are available. Implementing the TH paradigm will to some extent remove the impediment to hastily discover and formally describe all extant species in that the TH paradigm allows discovery and communication of new species and other taxa also in the absence of formal descriptions. The TH datasets can be connected to a taxonomic backbone providing access to the vast information associated with the tree of life. In parallel to the description of the TH paradigm, we demonstrate how it is implemented in the UNITE digital taxon communication system. UNITE TH datasets include rich data on individuals and their rDNA ITS sequences. These datasets are equipped with digital object identifiers (DOI) that serve to fix their identity in our communication. All datasets are also connected to a GBIF taxonomic backbone. Researchers processing their eDNA samples using UNITE datasets will, thus, be able to publish their findings as taxon occurrences in the GBIF data portal. UNITE species hypothesis (species level THs) datasets are increasingly utilized in taxon identification pipelines and even formally undescribed species can be identified and communicated by using UNITE. The TH paradigm seeks to achieve unambiguous, unique, and traceable communication of taxa and their properties at any level of the tree of life. It offers a rapid way to discover and communicate undescribed species in identification pipelines and data portals before they are lost to the sixth mass extinction.

Higher host plant specialization of root-associated endophytes than mycorrhizal fungi along an arctic elevational gradient.

How community-level specialization differs among groups of organisms, and changes along environmental gradients, is fundamental to understanding the mechanisms influencing ecological communities. In this paper, we investigate the specialization of root-associated fungi for plant species, asking whether the level of specialization varies with elevation. For this, we applied DNA barcoding based on the ITS region to root samples of five plant species equivalently sampled along an elevational gradient at a high arctic site. To assess whether the level of specialization changed with elevation and whether the observed patterns varied between mycorrhizal and endophytic fungi, we applied a joint species distribution modeling approach. Our results show that host plant specialization is not environmentally constrained in arctic root-associated fungal communities, since there was no evidence for changing specialization with elevation, even if the composition of root-associated fungal communities changed substantially. However, the level of specialization for particular plant species differed among fungal groups, root-associated endophytic fungal communities being highly specialized on particular host species, and mycorrhizal fungi showing almost no signs of specialization. Our results suggest that plant identity affects associated mycorrhizal and endophytic fungi differently, highlighting the need of considering both endophytic and mycorrhizal fungi when studying specialization in root-associated fungal communities.

Accounting for environmental variation in co-occurrence modelling reveals the importance of positive interactions in root-associated fungal communities.

Understanding the role of interspecific interactions in shaping ecological communities is one of the central goals in community ecology. In fungal communities, measuring interspecific interactions directly is challenging because these communities are composed of large numbers of species, many of which are unculturable. An indirect way of assessing the role of interspecific interactions in determining community structure is to identify the species co-occurrences that are not constrained by environmental conditions. In this study, we investigated co-occurrences among root-associated fungi, asking whether fungi co-occur more or less strongly than expected based on the environmental conditions and the host plant species examined. We generated molecular data on root-associated fungi of five plant species evenly sampled along an elevational gradient at a high arctic site. We analysed the data using a joint species distribution modelling approach that allowed us to identify those co-occurrences that could be explained by the environmental conditions and the host plant species, as well as those co-occurrences that remained unexplained and thus more probably reflect interactive associations. Our results indicate that not only negative but also positive interactions play an important role in shaping microbial communities in arctic plant roots. In particular, we found that mycorrhizal fungi are especially prone to positively co-occur with other fungal species. Our results bring new understanding to the structure of arctic interaction networks by suggesting that interactions among root-associated fungi are predominantly positive.

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.

Rhizosphere allocation by canopy-forming species dominates soil CO2 efflux in a subarctic landscape.

In arctic ecosystems, climate change has increased plant productivity. As arctic carbon (C) stocks predominantly are located belowground, the effects of greater plant productivity on soil C storage will significantly determine the net sink/source potential of these ecosystems, but vegetation controls on soil CO2 efflux remain poorly resolved. In order to identify the role of canopy-forming species in belowground C dynamics, we conducted a girdling experiment with plots distributed across 1 km2 of treeline birch (Betula pubescens) forest and willow (Salix lapponum) patches in northern Sweden and quantified the contribution of canopy vegetation to soil CO2 fluxes and belowground productivity. Girdling birches reduced total soil CO2 efflux in the peak growing season by 53%, which is double the expected amount, given that trees contribute only half of the total leaf area in the forest. Root and mycorrhizal mycelial production also decreased substantially. At peak season, willow shrubs contributed 38% to soil CO2 efflux in their patches. Our findings indicate that C, recently fixed by trees and tall shrubs, makes a substantial contribution to soil respiration. It is critically important that these processes are taken into consideration in the context of a greening arctic because productivity and ecosystem C sequestration are not synonymous.

Fungal ecological strategies reflected in gene transcription - a case study of two litter decomposers.

Microbial communities interplay with their environment through their functional traits that can be a response or an effect on the environment. Here, we explore how a functional trait-the decomposition of organic matter, can be addressed based on genetic markers and how the expression of these markers reflect ecological strategies of two fungal litter decomposer Gymnopus androsaceus and Chalara longipes. We sequenced the genomes of these two fungi, as well as their transcriptomes at different steps of Pinus sylvestris needles decomposition in microcosms. Our results highlighted that if the gene content of the two species could indicate similar potential decomposition abilities, the expression levels of specific gene families belonging to the glycoside hydrolase category reflected contrasting ecological strategies. Actually, C. longipes, the weaker decomposer in this experiment, turned out to have a high content of genes involved in cell wall polysaccharides decomposition but low expression levels, reflecting a versatile ecology compare to the more competitive G. androsaceus with high expression levels of keystone functional genes. Thus, we established that sequential expression of genes coding for different components of the decomposer machinery indicated adaptation to chemical changes in the substrate as decomposition progressed.

Nitrogen dynamics of decomposing Scots pine needle litter depends on colonizing fungal species.

In boreal ecosystems plant production is often limited by low availability of nitrogen. Nitrogen retention in below-ground organic pools plays an important role in restricting recirculation to plants and thereby hampers forest production. Saprotrophic fungi are commonly assigned to different decomposer strategies, but how these relate to nitrogen cycling remains to be understood. Decomposition of Scots pine needle litter was studied in axenic microcosms with the ligninolytic litter decomposing basidiomycete Gymnopus androsaceus or the stress tolerant ascomycete Chalara longipes. Changes in chemical composition were followed by 13C CP/MAS NMR spectroscopy and nitrogen dynamics was assessed by the addition of a 15N tracer. Decomposition by C. longipes resulted in nitrogen retention in non-hydrolysable organic matter, enriched in aromatic and alkylic compounds, whereas the ligninolytic G. androsaceus was able to access this pool, counteracting nitrogen retention. Our observations suggest that differences in decomposing strategies between fungal species play an important role in regulating nitrogen retention and release during litter decomposition, implying that fungal community composition may impact nitrogen cycling at the ecosystem level.