Positive citation bias and overinterpreted results lead to misinformation on common mycorrhizal networks in forests.

A common mycorrhizal network (CMN) is formed when mycorrhizal fungal hyphae connect the roots of multiple plants of the same or different species belowground. Recently, CMNs have captured the interest of broad audiences, especially with respect to forest function and management. We are concerned, however, that recent claims in the popular media about CMNs in forests are disconnected from evidence, and that bias towards citing positive effects of CMNs has developed in the scientific literature. We first evaluated the evidence supporting three common claims. The claims that CMNs are widespread in forests and that resources are transferred through CMNs to increase seedling performance are insufficiently supported because results from field studies vary too widely, have alternative explanations or are too limited to support generalizations. The claim that mature trees preferentially send resources and defence signals to offspring through CMNs has no peer-reviewed, published evidence. We next examined how the results from CMN research are cited and found that unsupported claims have doubled in the past 25 years; a bias towards citing positive effects may obscure our understanding of the structure and function of CMNs in forests. We conclude that knowledge on CMNs is presently too sparse and unsettled to inform forest management.

Seasonal and annual dynamics of western Canadian boreal forest plant communities: A legacy data set spanning four decades.

As boreal forests rapidly warm due to anthropogenic climate change, long-term baseline community data are needed to effectively characterize the corresponding ecological changes that are occurring in these forests. The combined seasonal dynamics (SEADYN) and annual dynamics (ANNDYN) data set, which documents the vegetative changes in boreal forests during the snow-free period, is one such source of baseline community data. These data were collected by George H. La Roi and colleagues in Alberta, Canada from 1980 to 2015 within permanent sampling plots established in the Hondo-Slave Lake area (eight stands; 1980-2015) in central Alberta and the Athabasca Oil Sands (AOS) region (17 stands; 1981-1984) near Fort McMurray in northeastern Alberta. Various data were collected, with temporal and spatial coverage differing by data set. These data sets include, but are not limited to, cover of each identified vascular plant and bryoid (moss, liverwort, and lichen) species; forest mensuration; forest litter production; and soil temperature and moisture. Notably, permanent sampling plots were set up as a grid, which will facilitate analyses of spatial relations. These data can be used to analyze long-term changes in seasonal dynamics and succession within boreal forest communities and serve as a baseline for comparison with future forest conditions in unmanaged, managed, and reclaimed forests. Data are released under a CC-BY license; please cite this data paper when using the data for analyses.

Choices on sampling, sequencing, and analyzing DNA influence the estimation of community composition of plant fungal symbionts.

Plant root symbionts, namely mycorrhizal fungi, can be characterized using a variety of methods, but most of these rely on DNA. While Sanger sequencing still fulfills particular research objectives, next-generation sequencing currently dominates the field, thus understanding how the two methods differ is important for identifying both opportunities and limitations to characterizing fungal communities. In addition to testing sequencing methods, we also examined how roots and soils may yield different fungal communities and how disturbance may affect those differences. We sequenced DNA from ectomycorrhizal fungi colonizing roots of Pinus banksiana and found that operational taxonomic unit richness was higher, and compositional variance lower, for Illumina MiSeq-sequenced communities compared to Sanger-sequenced communities. We also found that fungal communities associated with roots were distinct in composition compared to those associated with soils and, moreover, that soil-associated fungi were more clustered in composition than those of roots. Finally, we found community dissimilarity between roots and soils was insensitive to disturbance; however, rarefying read counts had a sizeable influence on trends in fungal richness. Although interest in mycorrhizal communities is typically focused on the abiotic and biotic filters sorting fungal species, our study shows that the choice of methods to sample, sequence, and analyze DNA can also influence the estimation of community composition.

Long-term nitrogen addition does not sustain host tree stem radial growth but doubles the abundance of high-biomass ectomycorrhizal fungi.

Global change has altered nitrogen availability in boreal forest soils. As ectomycorrhizal fungi play critical ecological functions, shifts in their abundance and community composition must be considered in the response of forests to changes in nitrogen availability. Furthermore, ectomycorrhizas are symbiotic, so the response of ectomycorrhizal fungi to nitrogen cannot be understood in isolation of their plant partners. Most previous studies, however, neglect to measure the response of host trees to nitrogen addition simultaneously with that of fungal communities. In addition to being one-sided, most of these studies have also been conducted in coniferous forests. Deciduous and "dual-mycorrhizal" tree species, namely those that form ecto- and arbuscular mycorrhizas, have received little attention despite being widespread in the boreal forest. We applied nitrogen (30 kg ha-1  year-1 ) for 13 years to stands dominated by aspen (Populus tremuloides Michx.) and hypothesized that tree stem radial growth would increase, ectomycorrhizal fungal biomass would decrease, ectomycorrhizal fungal community composition would shift, and the abundance of arbuscular mycorrhizal (AM) fungi would increase. Nitrogen addition initially increased stem radial growth of aspen, but it was not sustained at the time we characterized their mycorrhizas. After 13 years, the abundance of fungi possessing extramatrical hyphae, or "high-biomass" ectomycorrhizas, doubled. No changes occurred in ectomycorrhizal and AM fungal community composition, or in ecto- and AM abundance measured as root colonization. This dual-mycorrhizal tree species did not shift away from ectomycorrhizal fungal dominance with long-term nitrogen input. The unexpected increase in high-biomass ectomycorrhizal fungi with nitrogen addition may be due to increased carbon allocation to their fungal partners by growth-limited trees. Given the focus on conifers in past studies, reconciling results of plant-mycorrhizal fungal relationships in stands of deciduous trees may demand a broader view on the impacts of nitrogen addition on the structure and function of boreal forests.

Assessing the dual-mycorrhizal status of a widespread tree species as a model for studies on stand biogeochemistry.

Viewing plant species by their mycorrhizal type has explained a range of ecosystem processes. However, mycorrhizal type is confounded with plant phylogeny and the environments in which mycorrhizal partners occur. To circumvent these confounding effects, "dual-mycorrhizal" plant species may be potential models for testing the influence of mycorrhizal type on stand biogeochemistry. To assess their use as models, duality in mycorrhizas within a single host species must be confirmed and factors underlying their variation understood. We surveyed roots, soils, and leaves of mature aspen (Populus tremuloides) across 27 stands in western Canada spanning two biomes: boreal forest and parklands. Aspen roots were mostly ectomycorrhizal with sporadic and rare occurrences of arbuscular mycorrhizas. We further tested whether a climate moisture index predicted abundance of ectomycorrhizal roots (number of ectomycorrhizal root tips m-1 root length) surveyed at two depths (0-20 cm and 20-40 cm) and found that ectomycorrhizal root abundance in subsoils (20-40 cm) was positively related to the index. We subsequently examined the relationships between ectomycorrhizal root abundance, leaf traits, and slow and fast pools of soil organic carbon and nitrogen. The ratio of leaf lignin:N, but not its components, increased along with ectomycorrhizal root abundance in subsoils. Soil carbon and nitrogen pools were independent of ectomycorrhizal root abundance. Our results suggest that (1) categorizing aspen as dual-mycorrhizal may overstate the functional importance of arbuscular mycorrhizas in this species and life stage, (2) water availability influences ectomycorrhizal root abundance, and (3) ectomycorrhizal root abundance coincides with leaf quality.

Changes in soil fungal community composition depend on functional group and forest disturbance type.

Disturbances have altered community dynamics in boreal forests with unknown consequences for belowground ecological processes. Soil fungi are particularly sensitive to such disturbances; however, the individual response of fungal guilds to different disturbance types is poorly understood. Here, we profiled soil fungal communities in lodgepole pine forests following a bark beetle outbreak, wildfire, clear-cut logging, and salvage-logging. Using Illumina MiSeq to sequence ITS1 and SSU rDNA, we characterized communities of ectomycorrhizal, arbuscular mycorrhizal, saprotrophic, and pathogenic fungi in sites representing each disturbance type paired with intact forests. We also quantified soil fungal biomass by measuring ergosterol. Abiotic disturbances changed the community composition of ectomycorrhizal fungi and shifted the dominance from ectomycorrhizal to saprotrophic fungi compared to intact forests. The disruption of the soil organic layer with disturbances correlated with the decline of ectomycorrhizal and the increase of arbuscular mycorrhizal fungi. Wildfire changed the community composition of pathogenic fungi but did not affect their proportion and diversity. Fungal biomass declined with disturbances that disrupted the forest floor. Our results suggest that the disruption of the forest floor with disturbances, and the changes in C and nutrient dynamics it may promote, structure the fungal community with implications for fungal biomass-C.

The effects of ectomycorrhizal fungal networks on seedling establishment are contingent on species and severity of overstorey mortality.

For tree seedlings in boreal forests, ectomycorrhizal (EM) fungal networks may promote, while root competition may impede establishment. Thus, disruption to EM fungal networks may decrease seedling establishment owing to the loss of positive interactions among neighbors. Widespread tree mortality can disrupt EM networks, but it is not clear whether seedling establishment will be limited by the loss of positive interactions or increased by the loss of negative interactions with surrounding roots. Depending upon the relative influence of these mechanisms, widespread tree mortality may have complicated consequences on seedling establishment, and in turn, the composition of future forests. To discern between these possible outcomes and the drivers of seedling establishment, we determined the relative importance of EM fungal networks, root presence, and the bulk soil on the establishment of lodgepole pine and white spruce seedlings along a gradient of beetle-induced tree mortality. We manipulated seedling contact with EM fungal networks and roots through the use of mesh-fabric cylinders installed in soils of lodgepole pine forests experiencing a range of overstorey tree mortality caused by mountain pine beetle. Lodgepole pine seedling survival was higher with access to EM fungal networks in undisturbed pine forests in comparison with that in beetle-killed stands. That is, overstorey tree mortality shifted fungal networks from being a benefit to a cost on seedling survival. In contrast, overstorey tree mortality did not change the relative strength of EM fungal networks, root presence and the bulk soil on survival and biomass of white spruce seedlings. Furthermore, the relative influence of EM fungal networks, root presence, and bulk soils on foliar N and P concentrations was highly contingent on seedling species and overstorey tree mortality. Our results highlight that following large-scale insect outbreak, soil-mediated processes can enable differential population growth of two common conifer species, which may result in species replacement in the future.

Density-dependent processes fluctuate over 50 years in an ecotone forest.

Spatial patterns can inform us of forest recruitment, mortality, and tree interactions through time and disturbance. Specifically, successional trajectories of self-thinning and heterospecific negative density dependence can be interpreted from the spatial arrangement of forest stems. We conducted a 50-year spatial analysis of a forest undergoing succession at the ecotone of the southwestern Canadian boreal forest. The forest progressed from early to late sere and experienced repeated severe droughts, forest tent caterpillar outbreaks (Malacosoma disstria), as well as the outbreak of bark beetles. Cumulatively, the forest lost 70% of stems due to natural succession and a combination of disturbance events. Here, we describe spatial patterns displaying signals of successional self-thinning, responses to disturbance, and changes in patterns of density dependence across 50 years. Forest succession and disturbance events resulted in fluctuating patterns of density-dependent mortality and recruitment that persisted into late seral stages. The combined effects of conspecific and heterospecific density-dependent effects on mortality and recruitment resulted in near-spatial equilibrium over the study period. However, the strength and direction of these demographic and spatial processes varied in response with time and disturbance severity. The outbreak of forest tent caterpillar, pronounced drought, and bark beetles combined to reduce stand aggregation and promote a spatial equilibrium. Density-dependent processes of competition and facilitation changed in strength and direction with succession of the plot and in combination with disturbance. Together these results reinforce the importance of successional stage and disturbance to spatial patterns.

Expanding and testing fluorescent amplified fragment length polymorphisms for identifying roots of boreal forest plant species.

PREMISE OF THE STUDY: Identifying roots to species is challenging, but is a common problem in ecology. Fluorescent amplified fragment length polymorphisms (FAFLPs) can distinguish species within a mixed sample, are high throughput, and are inexpensive. To broaden the use of this tool across ecosystems, unique size profiles must be established for species, and its limits identified. METHODS: Fragments of three noncoding cpDNA regions were used to create size profiles for 193 species common to the western Canadian boreal forest. We compared detection success among congeners using FAFLPs and Sanger sequencing of the trnL intron. We also simulated and experimentally created communities to test the influence of species richness, cpDNA regions used, and extraction/amplification biases on detection success. RESULTS: Of the 193 species, 54% had unique size profiles. This value decreased when fewer cpDNA regions were used. In simulated communities, ambiguous species identifications were positively related to the species richness of the community. In mock communities, some species evaded detection owing to poor extraction or amplification. Sequencing did not increase detection success compared to FAFLPs for a subset of 24 species across nine genera. DISCUSSION: We recommend FAFLPs are best suited to confirm rather than discover species occurring belowground.

Role of urban ectomycorrhizal fungi in improving the tolerance of lodgepole pine (Pinus contorta) seedlings to salt stress.

With large forested urban areas, the city of Edmonton, Alberta, Canada, faces high annual costs of replacing trees injured by deicing salts that are commonly used for winter road maintenance. Ectomycorrhizal fungi form symbiotic associations with tree roots that allow trees to tolerate the detrimental effects of polluted soils. Here, we examined mycorrhizal colonization of Pinus contorta by germinating seeds in soils collected from different locations: (1) two urban areas within the city of Edmonton, and (2) an intact pine forest just outside Edmonton. We then tested the responses of seedlings to 0-, 60-, and 90-mM NaCl. Our results showed lower abundance and diversity of ectomycorrhizal fungi in seedlings colonized with the urban soils compared to those from the pine forest soil. However, when subsequently exposed to NaCl treatments, only seedlings inoculated with one of the urban soils containing fungi from the genera Tuber, Suillus, and Wilcoxina, showed reduced shoot Na accumulation and higher growth rates. Our results indicate that local ectomycorrhizal fungi that are adapted to challenging urban sites may offer a potential suitable source for inoculum for conifer trees designated for plating in polluted urban environments.

Additive or synergistic? Early ectomycorrhizal fungal community response to mixed tree plantings in boreal forest reclamation.

Ectomycorrhizal fungi are an important component to ecosystem function in the boreal forest. Underlying factors influencing fungal community composition and richness, such as host identity and soil type have been studied, but interactions between these factors have been less explored. Furthermore, mixed-species stands may have additive or synergistic effects on ectomycorrhizal fungi species richness, but this effect is challenging to test on natural sites due to difficulty in finding monospecific and mixed-species stands with similar site conditions and history. Forest reclamation areas can provide an opportunity to explore some of these fundamental questions, as site conditions and history are often known and managed, with the added benefit that knowledge emerging from these studies can be used to evaluate the recovery of degraded forest landscapes. Here, we compared the richness and composition of ectomycorrhizal fungi in young single- and mixed-species stands established on a reclamation area designed to inform strategies to restore upland boreal forests disturbed by oil sands mining. Seedlings of three host tree species (Populus tremuloides, Pinus banksiana, Picea glauca) were planted in single- and mixed-species stands on three different salvaged soils (forest floor material, peat, subsoil). After four growing seasons, there was no difference in total ectomycorrhizal fungi species richness and composition in mixed- versus combined single-species stands indicating that an additive effect of host tree species prevailed early in development. However, there were compositional shifts in fungal communities across both the host tree species and the salvaged soil type, with soil type being the strongest driver.

Accounting for local adaptation in ectomycorrhizas: a call to track geographical origin of plants, fungi, and soils in experiments.

Local adaptation, the differential success of genotypes in their native versus foreign environments, can influence ecological and evolutionary processes, yet its importance is difficult to estimate because it has not been widely studied, particularly in the context of interspecific interactions. Interactions between ectomycorrhizal (EM) fungi and their host plants could serve as model system for investigations of local adaptation because they are widespread and affect plant responses to both biotic and abiotic selection pressures. Furthermore, because EM fungi cycle nutrients and mediate energy flow into food webs, their local adaptation may be critical in sustaining ecological function. Despite their ecological importance and an extensive literature on their relationships with plants, the vast majority of experiments on EM symbioses fail to report critical information needed to assess local adaptation: the geographic origin of the plant, fungal inocula, and soil substrate used in the experiment. These omissions limit the utility of such studies and restrict our understanding of EM ecology and evolution. Here, we illustrate the potential importance of local adaptation in EM relationships and call for consistent reporting of the geographic origin of plant, soil, and fungi as an important step towards a better understanding of the ecology and evolution of EM symbioses.

Change in soil fungal community structure driven by a decline in ectomycorrhizal fungi following a mountain pine beetle (Dendroctonus ponderosae) outbreak.

Western North American landscapes are rapidly being transformed by forest die-off caused by mountain pine beetle (Dendroctonus ponderosae), with implications for plant and soil communities. The mechanisms that drive changes in soil community structure, particularly for the highly prevalent ectomycorrhizal fungi in pine forests, are complex and intertwined. Critical to enhancing understanding will be disentangling the relative importance of host tree mortality from changes in soil chemistry following tree death. Here, we used a recent bark beetle outbreak in lodgepole pine (Pinus contorta) forests of western Canada to test whether the effects of tree mortality altered the richness and composition of belowground fungal communities, including ectomycorrhizal and saprotrophic fungi. We also determined the effects of environmental factors (i.e. soil nutrients, moisture, and phenolics) and geographical distance, both of which can influence the richness and composition of soil fungi. The richness of both groups of soil fungi declined and the overall composition was altered by beetle-induced tree mortality. Soil nutrients, soil phenolics and geographical distance influenced the community structure of soil fungi; however, the relative importance of these factors differed between ectomycorrhizal and saprotrophic fungi. The independent effects of tree mortality, soil phenolics and geographical distance influenced the community composition of ectomycorrhizal fungi, while the community composition of saprotrophic fungi was weakly but significantly correlated with the geographical distance of plots. Taken together, our results indicate that both deterministic and stochastic processes structure soil fungal communities following landscape-scale insect outbreaks and reflect the independent roles tree mortality, soil chemistry and geographical distance play in regulating the community composition of soil fungi.

Stress differentially causes roots of tree seedlings to exude carbon.

How carbon (C) flows through plants into soils is poorly understood. Carbon exuded comes from a pool of non-structural carbohydrates (NSC) in roots. Simple models of diffusion across concentration gradients indicate that the more C in roots, the more C should be exuded from roots. However, the mechanisms underlying the accumulation and loss of C from roots may differ depending on the stress experienced by plants. Thus, stress type may influence exudation independent of NSC. We tested this hypothesis by examining the relationship between NSC in fine roots and exudation of organic C in aspen (Populus tremuloides Michx.) seedlings after exposure to shade, cold soils and drought in a controlled environment. Fine root concentrations of NSC varied by treatment. Mass-specific C exudation increased with increasing fine root sugar concentration in all treatments, but stress type affected exudation independently of sugar concentration. Seedlings exposed to cold soils exuded the most C on a per mass basis. Through 13C labeling, we also found that stressed seedlings allocated relatively more new C to exudates than roots compared with unstressed seedlings. Stress affects exudation of C via mechanisms other than changes in root carbohydrate availability.

Home-field advantage? evidence of local adaptation among plants, soil, and arbuscular mycorrhizal fungi through meta-analysis.

BACKGROUND: Local adaptation, the differential success of genotypes in their native versus foreign environment, arises from various evolutionary processes, but the importance of concurrent abiotic and biotic factors as drivers of local adaptation has only recently been investigated. Local adaptation to biotic interactions may be particularly important for plants, as they associate with microbial symbionts that can significantly affect their fitness and may enable rapid evolution. The arbuscular mycorrhizal (AM) symbiosis is ideal for investigations of local adaptation because it is globally widespread among most plant taxa and can significantly affect plant growth and fitness. Using meta-analysis on 1170 studies (from 139 papers), we investigated the potential for local adaptation to shape plant growth responses to arbuscular mycorrhizal inoculation. RESULTS: The magnitude and direction for mean effect size of mycorrhizal inoculation on host biomass depended on the geographic origin of the soil and symbiotic partners. Sympatric combinations of plants, AM fungi, and soil yielded large increases in host biomass compared to when all three components were allopatric. The origin of either the fungi or the plant relative to the soil was important for explaining the effect of AM inoculation on plant biomass. If plant and soil were sympatric but allopatric to the fungus, the positive effect of AM inoculation was much greater than when all three components were allopatric, suggesting potential local adaptation of the plant to the soil; however, if fungus and soil were sympatric (but allopatric to the plant) the effect of AM inoculation was indistinct from that of any allopatric combinations, indicating maladaptation of the fungus to the soil. CONCLUSIONS: This study underscores the potential to detect local adaptation for mycorrhizal relationships across a broad swath of the literature. Geographic origin of plants relative to the origin of AM fungal communities and soil is important for describing the effect of mycorrhizal inoculation on plant biomass, suggesting that local adaptation represents a powerful factor for the establishment of novel combinations of fungi, plants, and soils. These results highlight the need for subsequent investigations of local adaptation in the mycorrhizal symbiosis and emphasize the importance of routinely considering the origin of plant, soil, and fungal components.

MycoDB, a global database of plant response to mycorrhizal fungi.

Plants form belowground associations with mycorrhizal fungi in one of the most common symbioses on Earth. However, few large-scale generalizations exist for the structure and function of mycorrhizal symbioses, as the nature of this relationship varies from mutualistic to parasitic and is largely context-dependent. We announce the public release of MycoDB, a database of 4,010 studies (from 438 unique publications) to aid in multi-factor meta-analyses elucidating the ecological and evolutionary context in which mycorrhizal fungi alter plant productivity. Over 10 years with nearly 80 collaborators, we compiled data on the response of plant biomass to mycorrhizal fungal inoculation, including meta-analysis metrics and 24 additional explanatory variables that describe the biotic and abiotic context of each study. We also include phylogenetic trees for all plants and fungi in the database. To our knowledge, MycoDB is the largest ecological meta-analysis database. We aim to share these data to highlight significant gaps in mycorrhizal research and encourage synthesis to explore the ecological and evolutionary generalities that govern mycorrhizal functioning in ecosystems.

Biases underlying species detection using fluorescent amplified-fragment length polymorphisms yielded from roots.

BACKGROUND: Roots of different plant species are typically morphologically indistinguishable. Of the DNA-based techniques, fluorescent amplified-fragment length polymorphisms (FAFLPs) are considered reliable, high throughput, inexpensive methods to identify roots from mixed species samples. False-negatives, however, are not uncommon and their underlying causes are poorly understood. We investigated several sources of potential biases originating in DNA extraction and amplification. Specifically, we examined the effects of sample storage, tissue, and species on DNA yield and purity, and the effects of DNA concentration and fragment size on amplification of three non-coding chloroplast regions (trnT-trnL intergenic spacer, trnL intron, and trnL-trnF intergenic spacer). RESULTS: We found that sample condition, tissue and species all affected DNA yield. A single freeze-thaw reduces DNA yield, DNA yield is less for roots than shoots, and species vary in the amount of DNA yielded from extractions. The effects of template DNA concentration, species identity, and their interaction on amplicon yield differed across the three chloroplast regions tested. We found that the effect of species identity on amplicon production was generally more pronounced than that of DNA concentration. Though these factors influenced DNA yield, they likely do not have a pronounced effect on detection success of fragments and only underscore the restriction on the use of FAFLPs for measuring species presence rather than their abundance. However, for two of the regions tested-the trnT-trnL intergenic spacer and the trnL intron-size-based fragment competition occurred and the likelihood of detection was higher for smaller than larger fragments. This result reveals a methodological bias when using FAFLPs. CONCLUSIONS: To avoid potential bias with the use of FAFLPs, we recommend users check for the disproportionate absence of species detected belowground versus aboveground as a function of fragment size, and explore other regions, aside from the trnT-trnL intergenic spacer and trnL intron, for amplification.

Ectomycorrhizal fungi mediate indirect effects of a bark beetle outbreak on secondary chemistry and establishment of pine seedlings.

Dendroctonus ponderosae has killed millions of Pinus contorta in western North America with subsequent effects on stand conditions, including changes in light intensity, needle deposition, and the composition of fungal community mutualists, namely ectomycorrhizal fungi. It is unknown whether these changes in stand conditions will have cascading consequences for the next generation of pine seedlings. To test for transgenerational cascades on pine seedlings, we tested the effects of fungal inoculum origin (beetle-killed or undisturbed stands), light intensity and litter (origin and presence) on seedling secondary chemistry and growth in a glasshouse. We also tracked survival of seedlings over two growing seasons in the same stands from which fungi and litter were collected. Fungal communities differed by inoculum origin. Seedlings grown with fungi collected from beetle-killed stands had lower monoterpene concentrations and fewer monoterpene compounds present compared with seedlings grown with fungi collected from undisturbed stands. Litter affected neither monoterpenes nor seedling growth. Seedling survival in the field was lower in beetle-killed than in undisturbed stands. We demonstrate that stand mortality caused by prior beetle attacks of mature pines have cascading effects on seedling secondary chemistry, growth and survival, probably mediated through effects on below-ground mutualisms.