Introduction

Phosphorus (P) availability can be a fundamental constraint to primary productivity in many ecosystems, especially older soils where total P typically decreases and the ratio of organic P (Po) to inorganic P (Pi) increases [1]. P-depleted soils can also be found with podzolization under wetter climates, where precipitation encourages the solubilization of primary mineral P, accelerates conversion of Pi to Po by microorganisms, and increases the formation of secondary mineral P, thereby reducing the availability of orthophosphate for uptake by roots and microbes [2]. P deficiencies are also starting to occur in European and North American temperate forests affected by N pollution and climate change [3, 4].

Members of the Pinaceae are obligately symbiotic with ectomycorrhizal (EM) fungi [5]. In these trees, virtually all of the fine absorptive roots form ectomycorrhizas; hence, the EM roots can be considered the nutrient-absorbing organs of these plants and the ability of the fungi to absorb nutrients will influence plant nutrient status [6,7,8]. EM fungi use several strategies to enhance mobilization of both inorganic and organic P [9,10,11]. Of particular importance for P mobilization is extracellular hydrolytic enzymes (exoenzymes) that adhere to the surfaces of mantle and extramatrical (i.e., external to the mantle) hyphae of ectomycorrhizas [8, 12,13,14]. Two key types of exoenzymes involved in the direct hydrolysis of P from phosphoesters of soil organic matter are acid phosphomonoesterases (APMs) and phosphodiesterases [15, 16]. While phosphodiesterase activities of ectomycorrhizas have been little studied, increases in phosphomonoesterase activity of litter colonized by EM fungi have been associated with loss of P from the litter and accumulation of P in EM seedlings [17, 18]. Other fungal exoenzymes, such as cellobiohydrolase, β-glucosidase (GU), β-glucuronidase, xylosidase (X), and laccase (LAC), play an indirect role in P acquisition by degrading cell walls in plant and fungal necromass, thereby providing access to the nutrients within [8, 19]. EM fungal species differ considerably in their expression of these enzymatic P-acquisition traits [20,21,22]; consequently, the species and diversity of EM fungi colonizing a plant influences its P content (e.g., [6, 23]).

EM fungal communities are highly diverse, with community composition filtered by deterministic processes (host and soil environment), along with some stochastic (dispersal) influences [24,25,26]. Evidence is accumulating that compatibility between soil chemistry and EM fungal nutrient acquisition traits, including extracellular enzyme activities and ion flux rates, is part of the environmental filter for EM fungal community assembly. This has been observed in forests with different host species [27], and over successional [28, 29] and nitrogen (N) gradients [30]. Such congruence between functional traits and local edaphic conditions can arise through several mechanisms: species selection, adaptation of fungal populations to site conditions over time, or phenotypic plasticity (i.e., acclimation; [27]). Site-specific EM fungal communities possessing traits that increase access to growth-limiting nutrients are expected to contribute to the overall fitness of the plant partner and, hence, the symbiosis [31]. Indeed, the ability of seedling genotypes to survive on a specific soil is mediated, at least in part, by being able to associate with the EM fungi adapted to those soils [32]. If tree genotypes are not able to associate with the locally adapted EM fungi, their growth may be impaired relative to local tree genotypes [33]. As more forest regions become P-limited, the availability of EM fungal communities adapted to accessing specific form of soil P will be important for survival and/or migration of trees; however, little is known about the natural variation in P-acquisition strategies exhibited by EM fungal communities as sources of P change across landscapes.

Although we have focussed on response to soil P levels alone, a potentially more comprehensive evaluation of EM fungal response to nutrient deficiencies would be through resource stoichiometry, particularly the simultaneous constraints of N and P (expressed as an N:P ratio) on metabolism and growth [34]. The balance of N and P limitations is particularly relevant to the principle of optimal allocation, which suggests that mycorrhizal fungi will preferentially allocate energy toward acquiring the resources that are most limited in supply [31, 35, 36]. This principle has been widely established experimentally for decomposer biota such as bacteria and saprotrophic fungi [37]. For example, microbes will maximize N-acquiring enzymes under low N:P resources such as leaf litter with limited N mineralization rates, but then shift to P-acquiring enzymes under high N:P litter where N supply is less constraining [38]. The sensitivity of exoenzyme production to both N and P supply, rather than P alone, could therefore broaden our perspectives on stand nutrition and the rationale behind nutrient acquisition strategies of mycorrhizal fungi.

Here we report on P-acquisition exoenzyme activities associated with ectomycorrhizas along a podzolization gradient of coastal forests across Vancouver Island, where there has been a considerable depletion of orthophosphate and primary mineral P, alongside an accumulation of Po [7]. Our goal was to test whether EM fungal P-acquisition strategies changed along this soil P gradient, while holding host tree species and tree age constant. We measured the activity of a series of exoenzymes found on the surfaces of EM root tips of coastal Douglas-fir (Pseudotsuga menzesii var. menzesii) stands and relate the exoenzyme concentrations to soil P, including inorganic P as well as organic fractions revealed by 31P NMR. We hypothesized that activities of all exoenzymes would be higher in strongly podzolized soils, where Po contributes the majority of total P. In addition, we tested whether exoenzyme activities were predominantly related to soil P availability or whether a combined index with N (as soil or host N:P) would better reflect allocation strategies by EM fungi. Last, we explored whether the highest levels of phosphatase activities would be associated with EM fungal species restricted in their distribution to the podzolized, low-P soils. Results of our study will provide further insights into the manner and extent EM fungal communities of a single tree host have responded to declining soil P availability during soil development.

Methods

Site descriptions

We selected ten locations across an ~100-km extent of southern Vancouver Island (British Columbia, Canada) that encompassed a podzolization gradient driven by an orographic rainshadow, as described in ref. [7]. The sites were selected to capture a gradient in soil N and P availability, ranging from upland mesotrophic stands in both dry (referred to as “Upland Brunisol,” N = 4) and wet (“Upland Podzol,” N = 4) coastal forests, along with moist, nutrient-rich soils in wet coastal forests (“Lowland Podzol,” N = 4). We established one plot (25 × 25 m in size) at each site, with the exception of Branch 167 and 247, where the complex topography enabled us to sample both an Upland and Lowland Podzol plot (~400 m apart at both sites). All locations had second-growth P. menzesii stands that ranged in age from 40 to 60 years.

Soil and foliar attributes

Soil and foliar data were first reported in ref. [7], with the additions here of Bray P and 31P NMR spectroscopy conducted later on the same soil samples. Total soil P (Pt) was the sum of inorganic (Pi) and organic phosphorus (Po) concentrations as determined by the ignition method using a sulfuric acid (0.5 M) digest [39]. P in soil solution, or weakly adsorbed to soil colloids, was determined by a Bray P extraction [39]. One representative forest floor sample (the partly and well-decomposed horizons of the surface organic layer [40]) from two plots of each soil type (Upland Brunisol, Upland Podzol, Lowland Podzol) was selected for 31P NMR analysis following [41]. Only forest floor samples were analyzed because mineral soil samples did not have high enough P concentrations to detect peaks. Briefly, samples were shaken in 0.5-M NaOH and 0.1-M disodium EDTA solution for 4 h, centrifuged, filtered to 6 µm, freeze-dried, and re-dissolved in minimum amount of deuterated solvent. Dimethyl methyl phosphonic acid was used as an internal standard. Proton-decoupled spectra were obtained on a Varian 400-MHz spectrometer and processed using MestreNova v12.0 software (Supplementary Method 1).

Sampling of EM root tips

Sampling of EM roots took place over 3 weeks, from late May to mid-June 2018. At the start of each week, we sampled four plots that included one or two replicates of each soil type. Seven soil cores, ~15-cm wide and 15-cm deep, were collected randomly from each plot using a long knife and pruning clippers. The soil cores included both forest floor and mineral soil substrates, and many EM roots were found near the interface of these two substrates. Soil cores were transported on ice to the lab for overnight storage at 4 °C. Exoenzyme assays took place over the subsequent 3 days of each sample week.

Microplate exoenzyme activity assays

Exoenzyme activity assays (Supplementary Method 2) were modified from protocols of refs. [20, 42, 43] and the Kennedy lab website (https://cbs.umn.edu/kennedy-lab/protocols). Starting the day after field sampling, the soil cores of each plot were examined and fresh, vigorous clusters of EM root tips of the same morphology were selected and excised to proceed with exoenzyme activity assays. Distinct EM morphotypes were taken from each soil sample to ensure a wide representation of EM fungal species per plot, although some common morphotypes (e.g., Cenococcum) were selected from multiple soil samples within a plot. We retrieved between 10 and 16 morphotype clusters per plot, for a total of 157 clusters across all 12 plots.

Surface debris and soil particles were carefully removed from the EM root tips under a dissecting scope. Emanating hyphae were kept intact as much as possible during cleaning for those mantles with short-, medium- and long-distance exploration types [44]. Five individual EM root tips from each cluster were removed for the exoenzyme assay. The exoenzyme activities were assayed in series on each root tip (rinsing the roots between assays) in the following order: xylosidase (X), N-acetyl glucosaminidase (NAG), β-GU, APM, acid phosphodiesterase (APD), and LAC. The exoenzymes can be categorized as nutrient acquisition enzymes (NAG, APM, and APD), which hydrolytically cleave the target nutrient from the organic substrate, and lignocellulosic enzymes (X, GU, and LAC), which act to break down plant and fungal cell walls, in turn releasing the cellular contents for further action of nutrient acquisition enzymes. All the exoenzymes assayed are known to be cell-wall bound [45], and thus, their activities are not dependent on a continued supply of fixed carbon from the plant. The activities recorded in the lab should be considered as an index of activities present at the time of sample collection, including those of microbial communities living on the surface of the mycorrhizal tips. The APD assay was modified from ref. [43] for use on root tips rather than soil, with volumes for all reagents adjusted for 96-well microplates. Following the fluorescence or absorbance assays, the root tips were then transferred into a clear 96-well microplate and scanned (Epson V800) for an estimate of projected area using WinRhizo software (Regent Instruments Inc., Québec, Canada) (Supplementary Method 2).

Molecular analysis

For each cluster that was assayed for exoenzymatic activity, 5–10 fresh EM tips were retained at −80 °C for amplification of fungal DNA following a CTAB buffer extraction protocol (Supplementary Method 3). Fungal DNA was amplified using the universal ITS5/ITS4 or the basidiomycete-specific ITS1F/LR21 primer pairs. Forward and reverse sequences were aligned using Sequencher 4.7 (Gene Codes Corporation, Ann Arbor, MI USA). Sequences were BLAST-searched through the UNITE sequence database to assign Species Hypothesis based upon a >98.5% similarity criteria [46].

Data analysis

Prior to analysis, the exoenzymatic activities for all five root tips per cluster were averaged. Soil nutrient content (0–20-cm mineral soil plus forest floor) was determined using bulk density, coarse fragment content, and forest floor depths as outlined in ref. [7]. Relationships between exoenzymatic activities and soil P characteristics (Pi, Po, and Bray P as kg ha−1), soil C:N (molar), soil N:Pt (molar), foliar N (%), foliar P (%), and foliar N:P ratio (molar) were evaluated by linear regression using lm() function in R 3.4.1 [47]. We used all the exoenzymatic data in the analysis, regardless of whether the fungal partner was identified molecularly, as we considered morphological classification sufficient to determine mycorrhization. Before regressions were conducted, exoenzymatic activities were averaged by plot and then log-transformed in order to meet the normality criteria for linear regressions. Regression results for Pi and Bray P were virtually identical so for brevity only Pi outputs are presented here. A nonmetric multidimensional (NMDS) ordination of community composition between the three soil types, based on species incidence (Bray–Curtis dissimilarity), was undertaken with metaMDS() function and function ggplot2() [47].

Results

Soil P status

Upland Brunisols were well-endowed with readily available P, as indicated by Bray P concentrations ~25 times greater (average 83 [SE 14] mg kg−1 in mineral soil) than Podzols (average 3.2 [SE 0.6] mg kg−1; Table 1). A similar contrast was found with inorganic P (Pi) concentrations, which averaged 760 [SE 101] mg kg−1 in mineral soil for Brunisols and only 98 [SE 23] mg kg−1 across Podzol plots. Organic P concentrations (Po), by contrast, were generally highest across lowland Podzols, especially in mineral soils and forest floors (Table 1). Data on the P content (kg ha−1) of the upper profiles of these soils were reported previously, along with foliar N% and P% [7].

Table 1 Average soil phosphorus concentrations of forest floor and mineral soil for each plot, including Bray P, inorganic P (Pi) and organic P (Po) (SE in brackets).
Full size table

In addition to differences in total Po, NMR spectra revealed considerable differences in the chemical make-up of the organic P among soil types (Supplementary Fig. 1). Podzols had approximately twofold higher phosphodiester to orthophosphate ratios when compared to Brunisols, while phosphomonoester ratios were lowest in Upland Podzols and highest in Lowland Podzols (Table 2). Phosphonates were reasonably abundant in Lowland Podzols but were not detected in Brunisols (Table 2).

Table 2 Integration values for all NMR peaks when orthophosphate was set to 1.00.
Full size table

Exoenzyme activity in relation to soil nutrient status

Contrary to our hypothesis, none of the six exoenzyme activities measured on EM root-tip surfaces was related to soil Po content (kg ha−1 of the upper profile; Table 3). However, APD activity had a highly significant (p = 0.002, adj. R2 = 0.59), negative association with Pi content (Fig. 1a), along with a similar but more marginal decline in APM activity with increasing Pi (Table 3). The relationship of APD activity to soil Pi was logarithmic, with a marked increase at Pi contents <150 kg ha−1 (equivalent to <5 kg ha−1 Bray P), for a fourfold increase (0.03–0.12 μmol mm−2 min−1) in activity along the gradient from Brunisols to Podzols (Fig. 1a). The corresponding increase in APM activity was ~2.5-fold (0.15–0.35 μmol mm−2 min−1), although the concentrations of APM were 3–4× greater than APD overall. There were no significant changes in exoenzyme activity with soil N availability (as soil C:N ratio) over the podzolization gradient, but APD and β-GU showed a positive trend with soil N:Pt (Table 3 and Fig. 1b).

Fig. 1: Linear regressions between average acid phosphodiesterase (APD) activity (log-transformed) and soil properties.
figure 1

APD activity in relation to a soil inorganic P content, b soil N:Pt molar ratio (N = 12), and c ratio of diester to orthophosphate concentration (N = 6). (SE of plots as error bars, SE of regressions as shaded area).

Full size image
Table 3 Statistical significance of exoenzyme activities (μmol mm−2 min−1) in linear regressions with soil inorganic P (Pi) and organic P (Po) content (kg ha−1), along with soil C:N and N:Pt ratios (enzymes values log-transformed, p values < 0.05 in bold).
Full size table

While exoenzyme activity displayed no significant relationships with total soil Po content, specific components of the organic P identified by NMR appeared to be important. Specifically, we noted a strong, positive relationship (p = 0.002, adj. R2 = 0.90) between APD activity and the ratio of phosphodiesters to orthophosphate across the subset of sites analyzed (Fig. 1c). The relationship between APM activity and the relative concentration of phosphodiesters was also significant but weaker (data not shown; p = 0.048; adj. R2 = 0.58).

Exoenzyme activity in relation to foliar nutrient concentrations and stoichiometry

Foliar attributes, particularly P% and N:P ratio, had stronger relationships than the soil variables with exoenzyme activities (Table 4). We found clear declines in four of the six exoenzyme activities with increasing foliar P, but only one marginal gain (LAC) in relation to foliar N (Table 4 and Fig. 2). The extent of P deficiencies in relation to N, as measured by foliar N:P, was a more precise predictor than P% of exoenzyme activities, with significant positive correlations for all six exoenzymes. The relationships for foliar N:P were particularly robust (adj. R2 > 0.5) for APD, APM, and LAC activity (Table 4 and Fig. 2b, f). When aligned with increasing foliar N:P, the net gains in LAC, β-GU, and X exoenzyme activity ranged from two- to fivefold (Fig. 2d, f).

Table 4 Statistical significance of exoenzyme activities (μmol mm−2 min−1) in linear regressions with foliar nitrogen, phosphorus, and N:P ratio (enzyme values log-transformed, p values < 0.05 in bold).
Full size table
Fig. 2: Linear regressions of exoenzyme activity (log-transformed) with foliar attributes.
figure 2

First column depicts foliar P% and second column foliar N:P (molar ratio) for a, b acid phosphomonoesterase (APM) and phosphodiesterase (APD); c, d xylosidase (X) and b-glucosidase (GU); and e, f laccace (LAC) and N-acetylglucosaminidase (NAG). (N = 12, SE of plots as error bars, SE of regressions as shaded area).

Full size image

Exoenzyme activity in relation to EM fungal species

Root tips examined from the 12 plots were virtually 100% EM. Of 157 EM fungal clusters sampled, 130 were successfully sequenced (83%). Overall, we identified 31 EM fungal species across the Brunisol plots, 36 species in Upland Podzols, and 35 species on Lowland Podzols (Supplementary Table 1). The EM fungal communities of the Podzol and Brunisol soil types clustered separately in the NMDS ordination, with some overlap between Upland and Lowland Podzol plots (Fig. 3).

Fig. 3: Nonmetric multidimensional scaling ordination of EM communities (species incidence) between the three soil types.
figure 3

Ordination of EM communities from 12 plots, individual EM species represented by smaller points colored according to predominant soil type.

Full size image

A small subset of species, or multiple species from diverse genera (e.g., Cortinarius), was sampled frequently enough (>3) to illustrate some contrasts in exoenzyme activity among ectomycorrhizae formed by different fungal taxa (Fig. 4). The abundant ectomycorrhizae found only on Podzol sites, such as those formed by Lactarius substriatus, L. subviscidus, and Russula brevipes var. acrior had notably high expression of all three key exoenzymes (APM, APD, LAC). By contrast, Russula xerampelina, Cenococcum geophilum, and Piloderma olivaceum mycorrhizae, which were found almost exclusively on Brunisols, had much lower APD and LAC activity, while some genera, such as Cortinarius, were found across much of the edaphic gradient but featured little capacity for P-acquiring exoenzymes (Fig. 4).

Fig. 4: Average exoenzyme activity of acid phosphomonoesterase, acid phosphodiesterase, and laccase for select fungal taxa.
figure 4

Number of observations per fungal taxa in brackets, SE as bars.

Full size image

Discussion

The diversity of nutrient acquisition traits among EM fungal species has been posited to explain niche breadths of EM plants [48]. One challenging niche for the growth and establishment of trees is P-deficient soil, which can occur as a result of natural weathering as well as in landscapes where P imbalances have been induced by N deposition or climate change [6, 31]. In a naturally occurring gradient of P availability across southern Vancouver Island, we found evidence that EM fungal communities use different P-acquisition strategies depending on soil type, and that those strategies are tightly coupled with host foliar nutrition and stoichiometry. Specifically, EM fungal communities on Douglas-fir roots in P-deficient podzolized soils had much higher (two- to fivefold) activities of direct (APM, APD) and indirect (LAC, X, β-GU) P-acquiring exoenzymes than EM fungal communities in P-sufficient Brunisols. The large increases in exoenzyme activity by site-specific EM fungal communities and their associated microbiome [49] emphasizes the strong plant-soil feedbacks operating in these coastal ecosystems.

Even though phosphatase activities were highest in Podzols, with their high Po:Pi ratios, activities were more closely related to measures of inorganic than organic P in soil. This is consistent with other studies on ectomycorrhizas, where a negative relationship between inorganic P and EM mantle surface APM has been found ([50, 51], but see ref. [52]), indicating the likelihood of product (i.e., orthophosphate) inhibition [51]. Furthermore, in studies where Po has been applied to pure fungal cultures, soils, or EM plants [53,54,55,56], only one study found any relationship between Po and APM activity, in that case a negative one [57]. The lack of any relationship between phosphatases and total soil Po in our study and others may be because some forms of the organic P do not act as a substrate for these exoenzymes. For example, phosphonates, which were relatively abundant in the lowland Podzols, are generally considered very resistant to microbial degradation [58]. In addition, a large proportion of the Po was likely phytate [59], which, although it would have been grouped with monoesters in our NMR analysis, may not be used effectively by EM fungi [60]. When we quantified the relative amounts of diesters using 31P NMR, we were able to show a strong positive correlation with both APD and APM, suggesting that these exoenzymes were being produced in response to the increase of their substrates within the soil organic matter at Podzol sites. Phosphomonoesters are one of the products of phosphodiesterase activity.

Of the two types of phosphatase activity measured, APD appeared to be especially important in the response of ectomycorrhizas to low Pi. While the production of APD has been described as an important P-acquisition strategy of mycorrhizal fungi [10, 61], and other soil microbes [62, 63], it is often not included in EM enzymatic studies (but see refs. [64, 65]). Like APM, soil APD is integral to P cycling [15]. Phosphodiesters are considered to be important for the replenishment of Pi in soils because their phosphate groups are shielded from ionic interactions, making them less likely to be adsorbed to soil particles and more accessible to enzymatic cleavage [66, 67]. The activity of APD can be the rate-limiting step in Po turnover of some soils [68] because, unlike the P in monoesters, the release of P from soil phosphodiesters appears to be limited by enzyme availability [16]. In our study, the net increase in APD activity was twice that of APM across the edaphic gradient, emphasizing the focused allocation of resources to this exoenzyme.

The strong relationship we uncovered between EM exoenzyme activity and foliar N:P ratio was a key finding of this study and illustrates succinctly how the mycorrhizas likely participated in an optimal allocation strategy [31]. According to this model, plants and mycorrhizal fungi allocate biomass and energy to strategies that will acquire the more limiting nutrient(s), because organisms function optimally when stoichiometry is maintained. Consequently, at N:P foliar ratios below a threshold (generally 14–16) [69], resources would be expected to be focused on N acquisition, but above this threshold the limitations from P would equal and then surpass that of N. In this way the energy devoted by the EM fungi to exoenzyme allocation was more precisely governed by P deficiency in relation to N constraints, rather than P availability alone. The N:P ratios of EM fruiting bodies from these plots were nearly identical to P. menzesii foliar N:P [7], so the exoenzyme allocation strategy documented here reflects the stoichiometry of both the fungal symbiont and the tree, as would be expected if the host/symbiont relationship is viewed as a “holobiont” [70]. Foliar N:P, it should be noted, is a more expedient measure of resource availability than soil pools, especially Po, but can also be a somewhat dynamic property that is affected by, for example, stand development stage [71] or interannual climatic variation [72].

The prominent EM fungal species producing the highest quantities of P-acquiring exoenzymes (e.g., Lactarius substriatus, Russula brevipes var. acrior, Lactarius subviscidus) were found only on the strongly podzolized soils of the west coast of Vancouver Island, both in the sampling of this study and in previous, extensive EM fungal surveys of eastern Vancouver Island [73]. The wholescale species turnover within EM fungal communities likely selected for species that are adapted to exploiting Po, as indicated by enzymatic traits, and thereby able to outcompete other EM and saprotrophic fungi [29, 74]. Because of the dependence of P. menzesii on ectomycorrhizas for acquiring soil nutrients, these traits would also maximize fitness of the host in response to P scarcity, with better adapted EM fungi possibly favored by the host tree through the targeted provision of carbon [75], or by downregulation of plant defenses during root colonization [76]. Through competition with other fungi and/or support by the tree, assembly of the Podzol EM fungal communities could have involved selection of either fungal species [77] or individual fungi within the population with appropriate P-acquisition traits, or acclimation of existing fungi to the low Pi and high Po through upregulation of phosphatase genes. Given the high β diversity across this landscape, it appears less likely that capacity for physiological plasticity by fungal individuals or intraspecific variations were drivers of EM fungal community assembly here [27]. Some controlled, manipulative studies (e.g., reciprocal seedling transplants across soil types; [78]) could be used to test these capacities more thoroughly.

EM fungi are known to structure bacterial communities on their mantle surfaces and in soil [79]. The contribution of mantle-associated bacteria and archaea to exoenzyme activities measured in our study is not known, but incorporating antibiotics into assays have indicated a minor effect [80]. In addition, some of the variation in exoenzyme activities could reflect in part the morphology of the mantle, particularly for medium- and long-distance exploration types if extramatrical hyphae were lost during root preparation [44]. However, this potential bias should not have affected the average EM community response in exoenzyme activity since we found a mix of contact, short- and medium-distance mantles across all study plots. Furthermore, the considerable differences in exoenzyme activity among taxa within an exploration type, such as the contact mantle of Russula xerampelina and Lactarius substriatus, are consistent with a dominant effect of fungal species on these key traits [20,21,22, 41, 81].

Conclusions

A podzolization gradient across a short transect of temperate coastal forests provided an ideal field test for differences in P-acquisition strategies among EM fungal communities of P. menzesii. Our results establish a significant role for APD in the acclimation of EM fungi to low Pi soils. The close link in exoenzyme activity with host stoichiometry highlighted the importance of assessing P deficiency in relation to N constraints, rather than P alone, as would be consistent with an optimal allocation strategy. The EM fungal species with the highest expression of key exoenzymes were found only on Podzols and underscore how resource constraints, such as sharp declines in Pi with soil podzolization, likely drive species turnover and the assembly of distinct, well-adapted EM fungal communities. Large increases in many key exoenzymes by site-specific EM fungi illustrate how a diverse symbiotic community greatly enhances the ability of a single tree host to accommodate the severe edaphic constraints inherent to these wet coastal landscapes.