Mycorrhizal Association Better Predicts Tree Effects on Soil Than Leaf Habit

Predicting the ecological consequences of shifting forest composition requires an understanding of how tree species alter the soils beneath them. However, given the huge number of tree species on earth, it is critical to assess the extent to which easily observable traits facilitate the prediction of soil biogeochemistry. Both aboveground and belowground traits are important drivers of forest function. The objective of this study was to assess the relative importance of leaf habit (evergreen or deciduous) and root mycorrhizal association (arbuscular [AM] or ectomycorrhizal ECM]) on biogeochemistry. The relative importance of these two traits for carbon (C) and nitrogen (N) dynamics have proven difficult to disentangle as most deciduous tree species associate with AM fungi and most evergreen tree species associate with ECM fungi. Using planted, single-species forestry plots at The Morton Arboretum (DuPage County, IL), we found that mycorrhizal association better predicted tree effects on soil than leaf habit. While both leaf habit and mycorrhizal association drove variation in soil C:N ratios, mycorrhizal association alone predicted differences in soil pH, available N pools, extracellular enzyme activities, and C and N cycling rates. Surprisingly, ammonium concentrations were higher and net nitrification rates were faster in ECM plots than in AM plots. Greater N-degrading enzyme activities in ECM plots and intrinsically high soil pH across the garden likely drove more rapid N cycling in ECM plots. Overall, this study supports the inclusion of mycorrhizal association in terrestrial biosphere models, but suggests that the effects of mycorrhizal association on N dynamics may be site-dependent.


INTRODUCTION
Anthropogenic disturbances-including climate change, air pollution, invasive forest pests, and harvesting regimes-are altering the distribution of tree species and the composition of forest communities, with consequences for ecosystem functioning (Crowley et al., 2016;Crowley and Lovett, 2017). For instance, increases in temperature and precipitation have shifted the northern hardwood-boreal ecotone upward in elevation (Beckage et al., 2008), and climate-induced shifts in temperate forest community composition drive changes in soil biogeochemistry (Rollinson et al., 2012;McDaniel et al., 2013McDaniel et al., , 2014. As such, it is critical to characterize the effects of trees on soils in order to better predict the consequences of changing forest composition. Tree species differ in their influences on soil carbon (C) and nutrient cycling, in part, due to differences in the quantity and chemistry of their organic matter inputs (e.g., leaf and root litters and root exudates). However, given that there are over 60,000 known tree species (Beech et al., 2017), a species-level approach is impractical for both data collection and modeling applications. Instead, plant functional groups, which reflect suites of traits correlated through ecoevolutionary trade-offs (Reich, 2014), offer a practical solution for predicting consequences of tree distribution and forest composition shifts.
Leaf habit (evergreen or deciduous) and mycorrhizal association (arbuscular mycorrhizal [AM] or ectomycorrhizal [ECM]) are two traits commonly used to bin tree species into functional groups. Indeed, both leaf habit and mycorrhizal association have proven useful for predicting range shifts in response to climate change; deciduous forests are predicted to expand in response to climate change while the trailing range edge of AM trees constricts more rapidly that of ECM trees in response to climate change (Cramer et al., 2001;Lankau et al., 2015). Both classification schemes also distinguish tree nutritional strategies, with consequences for soil biogeochemistry. For instance, compared to deciduous trees, evergreen trees exhibit greater nitrogen (N) use efficiency and lower litter decomposition rates owing to greater litter recalcitrance (e.g., higher litter lignin:N ratios) (Aerts, 1995;Scott and Binkley, 1997;Silver and Miya, 2001;Reich and Oleksyn, 2004). Subsequently, soils under evergreen trees commonly have lower pH, wider C:N, and slower N cycling rates than those under deciduous trees (Binkley, 1995;Vesterdal et al., 2008;Cools et al., 2014;Augusto et al., 2015;Mueller et al., 2016). Similarly, due to a combination of litter chemistry and mycorrhizal nutrient acquisition strategies, stands dominated by AM-associated trees (AM stands) and stands dominated by ECM-associated trees (ECM stands) also have distinct soil pH, C:N, and N "economies" (Read, 1991;Cornelissen et al., 2001;Phillips et al., 2013;Averill et al., 2014;Lin et al., 2017;Zhu et al., 2018;Keller and Phillips, 2019). Specifically, relative to those under ECM-associated trees, soils under AM-associated trees tend to have higher pH, narrower C:N, greater inorganic N availability, and faster N mineralization and nitrification rates (Phillips et al., 2013;Averill et al., 2014;Lin et al., 2017;Zhu et al., 2018).
While the effects of both leaf habit and mycorrhizal association on soil processes have been considered separately, evaluating their relative importance has proven challenging due, in part, to correlations among traits in commonly studied clades. For instance, in a recent analysis of fine root decomposition rates of AM and ECM trees, 60% of broadleaf ECM species were from the order Fagales and 90% of the conifers were ECM . Few studies account for this phylogenetic autocorrelation (but see Averill et al., 2019). Similarly, most previous studies on leaf habit effects have represented evergreen trees with Pineaceae (Augusto et al., 2015). Pinaceae are also associated with ECM fungi, while most angiosperms associate with AM fungi (Averill et al., 2019), making it difficult to identify whether leaf habit or mycorrhizal association drives differences between these groups. Disentangling leaf habit and mycorrhizal association effects is even more challenging in soil biogeochemistry studies given that dominant ECM species in the wild are primarily from two clades-Fagales and Pinaceae-and AM deciduous trees are often represented by Acer Keller and Phillips, 2019). Considering these biases, extra care must be taken when interpreting the effects of leaf habit and mycorrhizal association on plant traits and ecosystem functions (summarized in Tedersoo and Bahram, 2019). Capturing broad phylogenetic breadth is paramount for robustly evaluating the impacts of mycorrhizal association and leaf habit on soil biogeochemistry.
Additionally, teasing apart tree effects on soils vs. soils as a driver of tree community composition has remained a challenge. Hans Jenny, in his seminal work, referred to this conundrum as "a real bugberrer" (Jenny, 1941). For example, sites with less soil N have high concentrations of evergreen trees, implying that deciduous trees out-compete evergreen trees in sites with high soil N (Binkley and Giardina, 1998;Lusk et al., 2003). Similarily, sites with less soil N tend to contain larger concentrations of trees associating with ECM fungi than with AM fungi, suggesting that AM-associated trees are also superior competitors in Nrich sites Jo et al., 2019). In contrast, in planted common garden experiments, tree effects on soil are isolated from tree colonization preferences. Common garden studies occasionally reveal unexpected patterns (Binkley, 1995;Mueller et al., 2016). In fact, in a meta-analysis of leaf habit effects on soil N mineralization, Mueller et al. (2016) found higher rates of N mineralization under deciduous trees compared to evergreen trees, but only when naturally occurring stands were included; when their analysis was limited to common garden studies, they found no differences in N mineralization rates between deciduous and evergreen stands. Thus, an open question remains-to what extent do trees with different leaf habits and/or mycorrhizal associations create sites with distinct soil properties? In this study, we used nearly 100-year-old tree plantations spanning a broad phylogenetic breadth to evaluate the relative importance of mycorrhizal association and leaf habit on soil C and N dynamics.

Site Description
This study was conducted in the forestry plots, a planted experiment similar to a common garden, at The Morton Arboretum in Lisle, Illinois (41.81 • N, 88.05 • W). The region has a continental climate with temperatures as low as −6 • C in January and up to 22 • C in July, with 800-1,000 mm mean annual precipitation. Soils found here are deep and moderatelyto poorly-drained Alfisols formed from a thin layer of loess (0.3-1 m) underlain by glacial till and Mollisols formed from alluvium. The major soil series are Ashkum, Beecher, and Ozaukee silt loams and Sawmill silty clay loam (Soil Survey Staff and Natural Resources Conservation Service United States Department of Agriculture, 2019). Prior to European settlement, The Morton Arboretum consisted mainly of burr and white oak savannah in addition to prairie and oak-dominated deciduous forest (Bowles et al., 1994). Following settlement, "high grading" or selective cutting of large diameter trees occurred, followed by some clear cutting in 1917 (Wilhelm, 1991).
Beginning in 1922, monoculture plots ranging from half an acre to four acres (0.2-1.6 hectares) in size were established with the goal of demonstrating the practical value of reforestation in the Midwestern U.S. By extension, the plots were also established to test and study "all the timber trees of the world which might come under consideration for reforestation purposes in this part of the country" (Morton Arboretum Staff, 1929). Currently, 38 of the initial plots remain intact and viable for examining tree effects on soil properties ( Table 1). Plots are distributed throughout the 1,700 acre (688 hectare) property of The Morton Arboretum ( Figure 1A). Plots with only 3-4 stems remaining were deemed too small and excluded from our study; plots were also excluded if the understories consisted of frequently mowed turf grass or mulch. These remaining plots include 28 species and range in size from 265 m 2 to 3.4 acres. Living trees in the plots were planted between 1922 and 1948 according to accession records (http://bol.mortonarb.org); library records and published bulletins indicate that most plots were planted in the 1920's. While five plots contain no accessioned trees and two contain accessioned trees with FIGURE 1 | Locations of forestry plots at The Morton Arboretum. Trees varying in mycorrhizal association (AM or ECM) and leaf habit (deciduous or evergreen) (A) were planted as monoculture plots on Alfisols (Ozaukee and Beecher soil series) and Mollisols (Sawmill or Ashkum soil series) (B). In (A), pink, AM plots; blue, ECM plots; light hues, deciduous plots; dark hues, evergreen plots. In (B), dark orange, Ozaukee series (Alfisol); light orange, Beecher series (Alfisol); dark green, Sawmill series (Mollisol); light green, Ashkum series (Mollisol); blue, water; gray, other soil series.
unknown planting dates, three of these plots were planted in the 1920's according to published bulletins, and one additional plot is clearly marked on a 1942 Evergreen trail map (Watts, 1942), indicating it was planted in the 1930's, at the latest. Only four plots have no accession or other published planting records. As such, it is likely that all plots were at least 70 years old at the time of sampling with most between 85 and 95 years old. One exception is a Robinia pseudoacacia plot that has no accessioned trees or planting records. Tree cores indicate the oldest trees in this plot were established in 1978 (Christine Rollinson, personal communication).
Tree species included in our study span broad phylogenetic and geographical-origin space with all combinations of mycorrhizal associations and leaf habits represented ( Table 1). As such, the species present in our study deviate from those commonly examined in other studies. There are 11 AM deciduous, 6 AM evergreen, 7 ECM deciduous, and 14 ECM evergreen plots. Mycorrhizal association assignments were made following Tedersoo and Brundrett (2017), Brundrett and Tedersoo (2019), Akhmetzhanova et al. (2012), and Wang and Qiu (2006). Species-level mycorrhizal associations were not available for Taxus cuspidata or Cladrastis kentukea, so mycorrhizal assignments were based on genus (Taxus cuspidata) and family (Cladrastis kentukea) information. All evergreen tree species are gymnosperms, and most deciduous tree species are angiosperms. Two Larix species are deciduous gymnosperms.
All plots were planted as monocultures of individual tree species except for two Picea plots that contain a mixture of Picea abies, Picea glauca, and Picea pungins. Though soil series and orders vary among the plots ( Figure 1B; Table 1), there are no systematic differences in soil properties underlying trees with varying leaf habits or mycorrhizal associations (X 2 (df = 1, N = 38) ≥ 0.173; P ≥ 0.356).

Soil Sampling and Processing
To evaluate the effects of tree traits on soil biogeochemistry, we collected soil samples four times during the 2018 growing season: in early May, mid-June, mid-August, and late October. Thirtytwo of the 38 plots were sampled in May, 37 were sampled in June, and all 38 were sampled in August and October. After removing the litter layer, we collected four cores from the top 15 cm of the mineral soil from each plot with a 5 cm diameter stainless steel soil corer. Organic horizons in these soils are thin or absent due to relatively high soil pH and nutrient status. Coring sites were purposefully located in the center of a triangle consisting of three dominant stems of the species of interest. The four coring sites were distributed evenly throughout a plot, and we took care to avoid edge effects; cores were collected at least 5 m from the edge of a given plot with edge borders defined by the outermost trees in a plot.
Upon returning to the laboratory, we separated cores into the upper 0-5 cm and lower 5-15 cm layers of soil. We aggregated layers for each plot and sieved them through a 2 mm sieve to homogenize the soil and remove rocks and roots. We stored sieved soils at 4 • C overnight before extracting C and N pools, conducting C mineralization assays, and initiating N cycling incubations. We conducted enzyme assays and chloroform fumigations within 12 months of collection. We stored soils for enzyme analyses at −80 • C prior to analyses; soils for microbial biomass fumigations at −20 • C prior to extraction; and all solutions at −20 • C prior to analyses.

Sample Analyses
We measured the water and organic matter content, as well as pH, of our samples using standard laboratory techniques. We dried subsamples of soil (5 g) at 105 • C for 24 h to determine gravimetric soil moisture. We subsequently measured organic matter (OM) content by ashing dry soils in a muffle furnace at 450 • C for 16 h. We measured soil pH in 0.01 M CaCl 2 solutions (5 g soil in 40 mL of solution) using a benchtop electrode pH meter (Orion 5 Star, Thermo Scientific, Beverly, Massachusetts).
To assess tree effects on soil C and N pools, we quantified total C and N, extractable inorganic N, extractable organic C and total N, and microbial biomass C and N concentrations. To measure total C and N, we dried soil samples at 55 • C and ground them to fine powders. We used the dry combustion method to measure total C and N concentrations (Vario El III, Elementar, Lengenselbold, Germany). To evaluate inorganic N pools, we extracted inorganic N (ammonium and nitrate) from 4 g soil samples with 2M KCl. We quantified ammonium concentrations using the salicylate-nitroprusside method (Sims et al., 1995) and measured absorbance at 660 nm on a microtiter plate reader (Synergy HTX, Biotek, Winooski, VT). We quantified nitrate concentrations using the VCl 3 /Griess method (Hood-Nowotny et al., 2010) and measured absorbance at 540 nm on a microtiter plate reader. Total inorganic N is the sum of ammonium and nitrate. We extracted organic C and total N from 10 g soil samples with 0.5M K 2 SO 4 and quantified total organic C and total N concentrations in extracts with high-temperature oxidation and chemiluminescence detection (TOC-L and TNM-L TOC-TN analyzer, Shimadzu, Kyoto, Japan). We calculated extractable organic N as the difference between total extractable N and total inorganic N. We determined soil microbial biomass C and N concentrations by quantifying changes in 0.5M K 2 SO 4extractable pools of C and N after 4 days of chloroform fumigation (Vance et al., 1987). We adjusted our microbial biomass C values to reflect an extractability of 45% (Beck et al., 1997) and our microbial biomass N values to reflect an extractability of 54% (Brookes et al., 1985).
To determine C mineralization rates, we incubated 5 g soil samples for 3 h at 23 • C and quantified C mineralization rates as the change in CO 2 concentrations over the course of the incubation (LI-6200, Li-Cor Incorporated, Lincoln, Nebraska). We determined net N mineralization rates by quantifying the change in 2M KCl-extractable pools of ammonium and nitrate in 4 g subsamples after an aerobic 14 day laboratory incubation at 23 • C. We measured net nitrification rates by quantifying the change in nitrate over the same time period. Ammonium and nitrate concentrations were quantified as above.

Statistical Analyses
We used mixed linear models to evaluate the effects of mycorrhizal association, leaf habit, and their interaction on soil biogeochemical properties. Our mixed linear models included mycorrhizal association (AM or ECM), leaf habit (evergreen or deciduous), and their interaction as fixed effects, plot ID as a random effect, and the above variables of interest as dependent variables to determine the impact of tree traits on soil properties. We ran separate sets of models for the two soil layers (0-5 cm and 5-15 cm). In addition to examining the above variables scaled by soil dry weights, we also evaluated the effects of mycorrhizal association, leaf habit, and their interaction on total soil C:N, microbial biomass C:N and the proportion of extractable N in inorganic forms. We also assessed tree trait effects on several enzyme ratios (BG:NAG, BG:AP, and NAG:AP). Enzyme ratios integrate relative microbial C and nutrient demand and availability, often varying with mycorrhizal association (Midgley and Phillips, 2016;Cheeke et al., 2017). We conducted all statistical analyses using the lmerTest package for mixed linear models in R (Kuznetsova et al., 2017;R Core Team, 2019). Results were considered significant if they had P-values < 0.05.

RESULTS
Across all analyses, there were no statistically significant interactions between mycorrhizal association and leaf habit (P ≥ 0.075). In addition, mycorrhizal association effects on soil properties in the 5-15 cm soil layer were generally consistent with those in the 0-5 cm soil layer or dampened (Tables 2, 3). Hence, we largely highlight the results from the main effects of mycorrhizal association and leaf habit on soil biogeochemistry in the 0-5 cm soil layer. Both mycorrhizal association and leaf habit drove variation in basic soil properties. While soil organic matter and water content did not vary between AM and ECM or deciduous and evergreen plots (Figure 2A), both mycorrhizal association and leaf habit drove differences in soil C:N-soil C:N was wider in ECM plots than in AM plots and wider in evergreen plots than in deciduous plots ( Figure 2C). While soil C content tended to be higher in evergreen plots than in deciduous plots (P = 0.074), neither %N nor %C clearly drove C:N patterns. However, mycorrhizal association alone predicted differences in soil pH-pH was higher in AM plots than in ECM plots ( Figure 2B). These pH patterns also manifested in the 5-15 cm soil layer. Differences in extractable N pools were largely limited to ammonium concentrations in AM and ECM plots. Ammonium concentrations were higher in ECM plots than in AM plots, though there were no significant differences in nitrate concentrations between AM and ECM plots ( Figure 3A). On average, nitrate concentrations were nearly two times greater than ammonium concentrations; as such, inorganic N concentrations were not significantly different but tended to be higher in ECM plots than in AM plots (P = 0.072). Organic N and total extractable N concentrations also exhibited little difference between AM and ECM plots. Total extractable N concentrations tended to be higher in ECM plots than in AM plots (P = 0.067), and both organic and extractable N concentrations tended to be higher in evergreen plots than in deciduous plots (TN: P = 0.056; DON: P = 0.050). Overall, the proportion of inorganic N in the extractable N pool did not vary with plot mycorrhizal association or leaf habit; inorganic N was 46 ± 2% of the total extractable N pool across all plots (Figure 3B).
Similarly, neither mycorrhizal association nor leaf habit drove variation in extractable organic C, microbial biomass C and microbial biomass N pools, or in microbial biomass C:N ratios, though microbial biomass C:N tended to be wider in ECM plots than in AM plots (P = 0.068), consistent with total soil C:N. However, in the 5-15 cm soil layer, microbial biomass C and N were both higher in AM plots than in ECM plots; we found no effects of mycorrhizal association on microbial biomass C:N in the 5-15 cm soil layer.
Extracellular enzyme activities and their ratios exclusively varied with plot mycorrhizal association. While BG and oxidative enzyme activities did not vary between AM and ECM plots (Figure 4A), both NAG and AP activities were greater in ECM plots than in AM plots (Figures 4B,C). As such, BG:NAG and BG:AP ratios were wider in AM plots than in ECM plots. However, we found no differences in NAG:AP ratios between AM and ECM plots. In the 5-15 cm soil layer, BG activities were higher in AM plots than in ECM plots, which drove wider BG:AP ratios in AM plots relative to ECM plots, consistent with BG:AP ratios in the 0-5 cm soil layer. Both C and N flux rates varied with plot mycorrhizal association, but not leaf habit (Figure 5). While N mineralization rates only tended to be faster in ECM plots relative to AM plots (P = 0.063), nitrification rates were faster in ECM plots than in AM plots. On average, nitrification rates were also faster than N mineralization rates, indicating that both initially available ammonium and N mineralized during the incubation were largely nitrified ( Figure 5A). Carbon mineralization rates were also faster in ECM plots than in AM plots ( Figure 5B).

DISCUSSION
Overall, mycorrhizal association was a better predictor of tree trait-driven differences in soil C and N dynamics than leaf habit. Specifically, soils in ECM plots had larger ammonium pools, greater N-degrading enzyme activities, and faster C and N flux rates than soils from AM plots. Additionally, soils in AM plots had higher pH than those in ECM plots. In contrast, leaf habit alone did not predict any variable. Soil C:N was wider in evergreen plots than in deciduous plots, but soil C:N was also wider in ECM plots compared to AM plots. Global syntheses have found that mycorrhizal association is a more important trait for predicting tree tissue chemistry and decomposition rate than leaf habit Averill et al., 2019;Keller and Phillips, 2019;See et al., 2019). Here, we show that mycorrhizal association is also more important than leaf habit for predicting soil biogeochemistry. However, our study also complements a growing body of literature demonstrating that mycorrhizal effects on soil processes are not uniform across sites (e.g., Craig et al., 2019). Given the recent push to incorporate mycorrhizal associations into terrestrial biosphere and global climate models Sulman et al., 2017Sulman et al., , 2019, it is imperative for underlying assumptions to be critically considered when inferring model implications.

Soil C and N Dynamics Vary With Mycorrhizal Association
For one soil property-soil C:N-AM deciduous trees and ECM evergreen trees represented end members along a spectrum. The widest soil C:N was found under ECM evergreen trees while the narrowest soil C:N was found under AM deciduous trees. This finding is consistent with a myriad of other studies finding wider C:N in ECM soils than in AM soils (i.e., Averill et al., 2014;Lin et al., 2017;Zhu et al., 2018) and in evergreen soils than in deciduous soils (Vesterdal et al., 2008;Cools et al., 2014;Augusto et al., 2015). Two mechanisms are commonly attributed to C:N differences among stands: the "Gadgil" effect and differences in plant tissue C:N Gadgil, 1971, 1975;Averill et al., 2014;Cools et al., 2014;Averill and Hawkes, 2016). However, the Gadgil effect-a build-up of soil C in ECM soils due to competition between ECM and saprotrophic fungi for resources-is minimized in sites rich in inorganic N (Fernandez et al., 2019;Maaroufi et al., 2019;Smith and Wan, 2019), such as ours. Furthermore, we found no evidence of increased C concentrations in ECM plots relative to AM plots. As such, Gadgil effects are not likely driving differences in C:N ratios in our study. Instead, differences in soil C:N likely reflect litter C:N. Deciduous leaf litter tends to have a narrower C:N ratio than evergreen leaf litter (Cornwell et al., 2008;Augusto et al., 2015) and ECM leaf litter tends to have a narrower C:N ratio than AM leaf litter (Keller and Phillips, 2019). Thus, our soil C:N ratios are likely correlated with the C:N of the leaf litter than falls on it (Cools et al., 2014). Carbon dynamics exclusively varied with mycorrhizal association and largely aligned with patterns found at other sites. ECM soils from the 0-5 cm layer exhibited faster C mineralization rates than AM soils, and AM soils from the 5-15 cm layer had greater BG activities than ECM soils. In contrast, mycorrhizal association did not drive variation in oxidative enzyme activities. Similarly, across a series of sites in southern Indiana, soils from ECM stands had faster C mineralization rates and lower BG activities than AM stands, and there were no consistent differences in oxidative enzyme activities (Brzostek et al., 2015;Midgley and Phillips, 2016;Cheeke et al., 2017). While other recent studies have detected faster heterotrophic respiration and soil respiration rates in AM soils or in AM-dominated stands compared to ECM soils and stands (Taylor et al., 2016;Wang and Wang, 2018;Lang et al., 2019), these discrepancies are likely due to differences in methodologies and calculations. For instance, while Taylor et al. (2016) found that heterotrophic respiration (per unit C) was faster from AM soils than from ECM soils, AM soils also had lower C concentrations. Additionally, in situ soil respiration rate measurements encompass both root and heterotrophic respiration; faster respiration rates from AM stands were largely driven by higher soil temperatures in a recent study at the Hubbard Brook Experimental Forest (Lang et al., 2019), while in our lab incubations, temperature was held constant. Overall, our C dynamic patterns are not unusual, though drivers of these patterns and their relationships with soil respiration remain to be explored.
Based on previous studies, we expected soil NAG activities to be higher, and soil inorganic N pools and fluxes to be lower, in ECM plots relative to AM plots-in other words, ECM plots would have an "organic N economy" (Phillips et al., 2013). While soil NAG activities were higher in ECM plots than in AM plots, inorganic N pools (ammonium, in particular) and nitrification rates were actually larger and faster in ECM plots than in AM plots. In contrast, N mineralization and nitrification rates are commonly slower in ECM stands than in AM stands (Phillips and Fahey, 2006;Phillips et al., 2013;Midgley and Phillips, 2016;Lin et al., 2017). Faster inorganic N cycling rates and availability in ECM plots are especially surprising given the wider C:N and lower pH soils found in ECM plots compared to AM plots; N mineralization and nitrification rates are often negatively correlated with C:N and positively and positively correlated with pH (Scott and Binkley, 1997;Finzi et al., 1998a,b;Gundersen et al., 1998;Lovett et al., 2004;Hobbie et al., 2006). Perhaps most surprisingly, we found no evidence of differences in the N "economies" of AM and ECM plots; the relative abundance of soil inorganic N was similar between AM and ECM plots.

Drivers of Mycorrhizal Effects on N Dynamics
What underlies these surprising N dynamics? While ECM fungi and their hosts may have co-evolved in low N environments (Arnoldi et al., 2019;Lu and Hedin, 2019), the characteristics that enable them to persist in low N environments (i.e., N acquisition abilities and N use efficiency) may enhance N availability in high N environments, as we found in this study. In fact, in a uniform N environment, ECM trees had faster root exudation rates than AM trees (Liese et al., 2018), indicating that faster exudation rates from ECM trees are not simply a phenotypic response to low N environments. ECM trees and fungi enhance N-rich compound degradation via root exudate priming and Ndegrading enzyme production (Read and Perez-Moreno, 2003;Talbot et al., 2008;Yin et al., 2014;Lin et al., 2018), and ECM trees take up relatively more organic N than and have higher N use efficiency than AM trees (Liese et al., 2018;Zhang et al., 2018). Thus, despite generally high soil N availability, ECM fungi and trees still produce high concentrations of Ndegrading enzymes in our ECM plots, and the ECM trees may continue to be efficient in their N use. These mechanisms account for the build-up of inorganic N, ammonium in particular, in ECM plots.
Faster nitrification rates in ECM plots relative to AM plots are the likely product of higher ammonium concentrations in ECM plots and the generally high soil pH at our site. As is commonly found across other sites, AM plots had higher soil pH than ECM plots (Phillips et al., 2013;Taylor et al., 2016;Chen et al., 2018), likely owing to the production of organic acids by roots and litter-degrading microbes. Nitrification is often suppressed in lower pH sites (Ste-Marie and Paré, 1999), and this may contribute to slower ECM nitrification rates found in other studies comparing AM and ECM stands (Phillips et al., 2013;Midgley and Phillips, 2016). However, the soil pH at our site was generally high−5.6 in ECM sites and 6.3 in AM sites, in the 0-5 cm soil layer. As such, despite the acidifying effects of ECM trees and fungi, soil pH did not mediate nitrification rates. Given the abundant ammonium and high pH under ECM trees at our site, subsequently higher nitrification rates in ECM soils are to be expected.
Stepping back, these seemingly surprising N dynamics are not particularly unusual. For instance, a study from Whitehall Forest in Georgia found no differences in inorganic N pools under AM and ECM trees (Taylor et al., 2016). Their pH ranges and averages were similar to ours (ECM = 5.6; AM = 6.1), and their nitrate pools were an order of magnitude larger than ammonium pools. Similarly, under naturally-occurring trees in northeast China, soils under ECM trees tended to have faster N mineralization rates and higher NAG enzyme activities relative to those under AM trees (Chen et al., 2018). Again, ECM trees had soils with an average pH of ∼5.5 compared to AM trees that had an average soil pH ∼6. In contrast, examples of AM and ECM stands with distinct N "economies" are largely from the same series of forests in southern Indiana, which have inherently low pH (AM = 5.2 and ECM = 4.3 at Moores Creek; [Phillips et al., 2013]). In a cross-site study, nitrate decreased as ECM dominance increased in low pH sites, but this relationship did not persist at higher pH sites (Craig et al., 2019). Similarly, in a study from Denmark that evaluated the effects of several tree species on soil biogeochemistry, there were no detectable effects of tree species on soil N dynamics at a higher pH site with an agricultural land-use history while at two previously forested sites with lower pH, there were clear "spruce effects" (Ribbons et al., 2018). In short, tree species effects, and AM and ECM tree effects, in particular, are modulated by site conditions, and inherent soil pH may be a particularly important factor mediating inorganic N pools and fluxes.

Implications and Future Directions
Future studies should identify the penultimate drivers of our observed differences in soil C and N pools and fluxes by examining soil property relationships with continuous tree traits, such as leaf and root litter and root exudate quantity and chemistry (C:N, lignin:N, calcium content, etc.). Stand-level properties including tree productivity and standing biomass, light incidence to the forest floor, and understory community composition may also drive variable C and N dynamics. For instance, if trees or understory plants in our AM stands are more productive than those in our ECM stands, this could account for the lower concentrations of available ammonium in AM plots (Binkley and Giardina, 1998). Studies that have considered both tree species and functional groups (e.g., mycorrhizal association, leaf habit) have often found that variation within these groups is greater than variation between groups (i.e., Taylor et al., 2016;Chen et al., 2018;Sun et al., 2018). As such, identifying the continuous traits that best predict differences among species within these groups will be key to incorporating variation into terrestrial biosphere models. Root properties may be particularly important for soil biogeochemistry. However, while some studies have found differences in exudation rates and root decomposition between AM and ECM trees (Yin et al., 2014;Liese et al., 2018;See et al., 2019), others have not (Brzostek and Finzi, 2011;Sun et al., 2018). Furthermore, the effects of mycorrhizal fungi on C and N dynamics vary markedly across fungal lineages (Lindahl and Tunlid, 2015;Pellitier and Zak, 2018). General understanding of variation in root and fungal characteristics and their relationships with soil C and N dynamics remains limited; this topic clearly merits further investigation. Our study complements a growing body of research suggesting that mycorrhizal association is a more important trait than leaf habit; as such, mycorrhizal association is an important trait to incorporate into terrestrial biosphere and global climate models Sulman et al., 2017Sulman et al., , 2019. However, intrinsic soil properties may mediate the impacts of trees on soils (Callesen et al., 2013;Ribbons et al., 2018;Craig et al., 2019). While percent clay is often considered (e.g., Averill et al., 2014), we have argued here that soil pH may be important for predicting mycorrhizal effects on N cycling. Soil pH can be predicted by soil and climate (Slessarev et al., 2016), and, thus, could be incorporated into global models. This could assist in identifying where tree-driven changes in pH will or will not lead to differences in soil N dynamics. While mycorrhizal association is increasingly incorporated into global models Sulman et al., 2019), these are parameterized assuming ECM trees exist in low N environments. Here, we show that this may not be universally true, and we urge future modeling efforts to account for this context-dependency.
Broadly speaking, some mycorrhizal-associated plant traits and ecosystem properties seem to be very stable across sites and studies (e.g., leaf and litter N; soil C:N and pH), others may be context-dependent (leaf litter decomposition; soil inorganic N; rhizosphere effects), and others still may be driven by factors unrelated to tree mycorrhizal associations (root decomposition and other root traits). Our common garden-style design, which captures a broad phylogenetic range of tree species, may better reflect tree effects on soil properties than studies performed in natural settings. New common garden studies, designed to evaluate mycorrhizal effects across broad phylogenetic space and across sites with varying inherent soil properties, are needed to obtain a better picture of AM vs. ECM effects on soils, their interactions with other traits, the mechanisms driving these effects, and their context-dependency.

DATA AVAILABILITY STATEMENT
The datasets generated for this study are available on request to the corresponding author.

AUTHOR CONTRIBUTIONS
MM formulated the research questions and designed the research. MM and RS collected samples, analyzed samples, and conducted statistical analyses. MM wrote the first draft of the manuscript with significant intellectual and editing input from RS. All authors contributed to manuscript revision, and read and approved the final submission.