- 1School of Life and Environmental Sciences, Sydney Institute of Agriculture, The University of Sydney, Camden, NSW, Australia
- 2School of Agriculture, Food and Wine, The University of Adelaide, Adelaide, SA, Australia
Plants spend a high proportion of their photosynthetically fixed carbon (C) belowground to support mycorrhizal associations in return for nutrients, but this C expenditure may decrease with increased soil nutrient availability. In this study, we assessed how the effects of nitrogen (N) fertiliser on specific root respiration (SRR) varied among mycorrhizal type (Myco type). We conducted a multi-level meta-analysis across 1,600 observations from 32 publications. SRR increased in ectomycorrhizal (ECM) plants with more than 100 kg N ha−1 applied, did not change in arbuscular mycorrhizal (AM) and non-mycorrhizal (NM) plants, but increased in plants with a dual mycorrhizal association in response to N fertilisation. Our results suggest that high N availability (>100 kg N ha−1) could disadvantage the growth of ECM plants because of increased C costs associated with maintaining higher root N concentrations, while the insensitivity in SRR by AM plants to N fertilisation may be because AM fungi are more important for phosphorus (P) uptake.
Introduction
Plants allocate up to 50% of photosynthetically fixed carbon (C) to root biomass and rhizodeposition (Pausch and Kuzyakov, 2018). In addition, plants spend a significant portion of fixed C to support mycorrhizal fungi in exchange for soil-derived nutrients (Smith et al., 2009); most of this C expenditure is associated with root and hyphal respiration (Hughes et al., 2008). Up to 86% of all flowering plants can form symbiotic associations with mycorrhizal fungi (Brundrett, 2009; van der Heijden et al., 2015), where associations with arbuscular mycorrhizal (AM) and ectomycorrhizal (ECM) fungi are the most widespread type of mycorrhiza (Brundrett, 2009). The C cost associated with forming and maintaining different types of mycorrhizas remains highly uncertain, with estimates as low as 4% and as high as 25% of the net primary production of a plant allocated to the fungi (Hobbie, 2006; Johnson and Gehring, 2007; Stuart and Plett, 2020). Plants associated with AM and ECM fungi differ in their nutrient economies and related biogeochemical transformations of C and nutrients (Phillips et al., 2013). Thus, it can, therefore, be expected that variation in nutrient availability will also affect root respiratory C costs of AM and ECM plants, with consequences for their productivity and abundance in terrestrial ecosystems.
In general, the mycorrhizal contribution to plant nitrogen (N) acquisition is higher for ECM than for AM fungi, while AM fungi tend to have a more important role in phosphorus (P) uptake than ECM fungi (van der Heijden et al., 2015). Whereas, ECM plants also tend to dominate in organic soils with relatively high C:N ratios (Averill et al., 2014; Jo et al., 2019; Soudzilovskaia et al., 2019), AM plants tend to dominate in mineral soils with lower C:N ratios and higher net N mineralisation rates (Johnson, 2010; Lin et al., 2017; Jo et al., 2019). ECM plants have an advantage in organic soils, because ECM fungi are able to decompose soil organic matter and acquire organic N via the production of specialised enzymes (Hobbie et al., 2013). However, because the production of these C enzymes is expensive (Soudzilovskaia et al., 2015), ECM symbiosis of C may also be costly to the plant under these conditions. Unlike ECM, AM fungi have no specialised enzymes for organic N acquisition (Terrer et al., 2016; Jansa et al., 2019) and are not able to decompose soil organic matter (SOM) directly, but they can accelerate the decomposition indirectly by translocating plant-derived C to soil microorganisms to stimulate their activity (Hodge and Fitter, 2010; Taylor et al., 2016).
Nitrogen availability can affect the plant-fungi relationship by changing the balance between costs and benefits of forming mycorrhizal associations (Phillips et al., 2013). Reliance of plants on mycorrhizal symbiosis reduces when N availability increases and is more easily accessible to plants; therefore, plants may reduce their C investment to mycorrhizas with increased N fertilisation (Högberg et al., 2010; Thirkell et al., 2019). N fertilisation further can induce P limitation, which could shift C expenditure from acquiring N to acquiring P (Johnson, 2010; Treseder et al., 2018). The effects of N availability on mycorrhizas have mostly been assessed through changes in colonisation percentage per length of root (for AM) or root tips (for ECM), showing an overall decrease of 15% with N fertilisation across different studies, in a meta-analysis (Treseder, 2004), but there is large variation among studies. Indeed, AM colonisation could increase with N fertilisation because AM fungi also have a substantial N demand, thereby competing with plants for N (Puschel et al., 2016; Wang et al., 2018), or because of a shift towards increased plant P demand (Corkidi et al., 2002).
Given that N fertilisation affects mycorrhizal colonisation, this should result in changes in C costs for plants to maintain mycorrhizal associations. Indeed, specific root respiration (SRR) rate, which is defined as the root respiration rate per unit of root biomass (Pregitzer et al., 2007), tends to be greater for roots in association with mycorrhizal fungi, because mycorrhizal respiration also contributes to SRR (Hughes et al., 2008). SRR can be directly measured on excised roots (Makita et al., 2012; Wang and Liu, 2014), which in most cases severs external hyphae, and, therefore, most likely only includes the contribution of intraradical hyphae (for AM) and other mycorrhizal structures in or on the surface of roots (for AM and ECM). When root respiration is indirectly measured as the difference in soil respiration rates between soils with and without live roots present (Ding et al., 2010; Tu et al., 2013), it not only includes respiration of external hyphae but also respiration from recent rhizodeposits.
In theory, increased N availability with N fertilisation reduces plant reliance on mycorrhizal symbiosis, which may result in a decline of mycorrhizal colonisation and, thence, a lower SRR rate. However, evidence for this is sparse; most N fertilisation studies have concentrated on examining changes in root respiration and have not accounted for concurrent changes in root biomass. For instance, root respiration in an ECM boreal forest increased at a low level of N fertilisation (20 kg N ha−1) but decreased at a high level of N (100 kg N ha−1) (Hasselquist et al., 2012), while chronic N fertilisation had no effect on root respiration in northern hardwood forests dominated by AM trees (Burton et al., 2004). It is difficult to pinpoint the cause of the variability in responses in root respiration, which could be due to the level of N fertiliser applied and due to the variable N fertilisation effects on root biomass (so that fertilisation effects on root respiration may not be the same as SRR), on mycorrhizal colonisation rates of ECM and AM, and/or on root N concentrations. Indeed, increased root N concentrations are linked to higher SRR rates (Reich et al., 2008), but increased N in roots can be a result of either direct root uptake of N or mycorrhizal support.
While mycorrhizas represent a significant sink for plant photoassimilates, and the amount of C invested is likely to differ with Myco type and N supply, the impacts of Myco type and N fertilisation on SRR are less clear. Therefore, our meta-analysis included an investigation of the effects of N fertilisation on SRR to assess N fertiliser effects on C loss per unit of root biomass for plants associated with different types of mycorrhizas. The effects of N fertilisation on soil C fluxes and their components (heterotrophic and autotrophic respiration) have been reviewed (Zhou et al., 2014), but to our knowledge, this is the first study synthesising N fertilisation effects on SRR by plants varying in mycorrhizal symbiosis type through a meta-analysis. We focused on SRR (i.e., respiration per unit of root biomass) rather than total root respiration as a measure of plant C costs because it accounts for potential confounding effects of N fertilisation on root biomass. Environmental conditions such as soil moisture (Soil moist) and temperature can also affect SRR mediated by N fertilisation (Zhang et al., 2014). We, therefore, used multi-model inference to identify the most important predictors on SRR rates in response to N fertilisation; predictors include Soil moist and soil temperature (Soil temp), mycorrhizal symbiosis types [AM, ECM, AM+ECM, and non-mycorrhizal (NM)], root respiration measurement method, and N fertiliser amount. To examine potential confounding effects of plant growth forms from that of mycorrhizal association, we compared N fertilisation effects on the variation of SRR between woody and herbaceous plants. We further assessed the effect of N fertilisation on root N concentration for the same studies that also had information on SRR. We hypothesised that N fertilisation would reduce SRR in ECM plants because of reduced reliance on the symbiotic relationship with ECM fungi but that N fertilisation would not affect SRR in plants associated with AM because these plants still rely on the symbiotic relationship with AM fungi for P uptake. Furthermore, given the lack of symbiosis in NM plants, we hypothesise that N fertilisation would have no effect on SRR in NM plants.
Methods
Literature Survey and Inclusion Criteria
We collected data related to SRR from published literature studies by searching in the Scopus and ScienceDirect databases in June 2021. We used the advanced search engine and searched for research articles published using the following keywords: (1) “specific root respiration,” (2) “root respiration” AND “root biomass,” (3) “autotrophic respiration” AND “root biomass” OR “fertilisation,” (4) “rhizosphere respiration” AND “root biomass” OR “fertilisation,” and (5) “below-ground respiration” AND “root biomass” OR “fertilisation.” After excluding conference study papers, we found a total of 1,051 study papers. We applied the following criteria to be included in our meta-analysis: studies should have an N fertilisation treatment (without or in combination with other treatments) and should either report SRR directly or report root biomass and root respiration (or autotrophic respiration) allowing us to calculate SRR. We considered three different methods for calculating SRR: (1) based on the measurements of root respiration on excised roots (“root excavation”), (2) based on the difference in respiration between trenched and non-trenched plots (“trenching”), and (3) based on the difference in respiration between planted and unplanted or root-free/bare soil treatments (“unplanted”). For the trenching and unplanted methods, root biomass measurements in non-trenched and planted treatments also needed to be reported to calculate SRR. We excluded publications that were conducted in (1) hydroponic systems, (2) pots filled with sand only, (3) experiments that used oxygen (O2) consumption to estimate respiration, and (4) C isotope pulse labelling experiments that were unable to quantify root respiration. Multiple respiration measurements with time were included as separate observations. We found a total of 32 study papers that fit these criteria (Supplementary Table 1) and retrieved a total of 1,600 observations from these study papers (1,559 observations from field experiments and 41 observations from pot experiments, Supplementary Information).
Data Collection
We used Plot Digitizer software Ver. 2.6.8 to extract the data from figures. We also extracted Soil moist, temperature, and N concentrations in roots when reported. Soil moist units were reported differently in different publications including volumetric, gravimetric, and water-filled pore space (WFPS). We converted all Soil moist data to gravimetric moisture using the bulk density when reported (including five study papers), or otherwise, we used a bulk density of 1.2 g cm−3 (another five study papers) (Bache et al., 2008). For each plant species, we collected information on mycorrhizal association type, including AM association, ECM association, AM+ECM (dual association), and non-symbiotic (NM) association based on plant species known to form an association with different types of mycorrhizas (Brundrett, 2009; Cosme et al., 2018; Teste et al., 2020). We considered AM+ECM association for those observations from mixed forests containing species associated with both AM and ECM and where respiration measurements were obtained by the trenching method. We further recorded whether plants were herbaceous or woody. Fertiliser application amount, experimental type (field vs. pot experiment), and method used for respiration measurement were also recorded. The fertiliser application rate was grouped into low ( ≤ 100 kg N ha−1) and high fertilisation rates (> 100 kg N ha−1).
Calculation and Statistical Analysis
Data were back-transformed when they were reported in a natural log. For those study papers that reported SE, we calculated the SD (σ) by multiplying the SE with the root square of the sample size (N). In cases where root respiration was calculated based on the difference between the mean total (non-trenched, planted) and heterotrophic (trenched, unplanted) respiration, the SD for root respiration was calculated using the pooled Equation (1) (Kaltenbach, 2012):
where σtr and n are the SD and sample size for total respiration, and σhr and m are the SD and sample size for heterotrophic respiration, respectively.
The Taylor expansion Equation (2) was used to calculate the SD for SRR when it was derived from the mean of root respiration divided by the mean of root biomass (Kendall, 1987):
where σrr/rb is the SD of SRR, μrr and σrr are the mean and SD for root respiration, and μrb and σrb are the mean and SD for root biomass, respectively.
Finally, for those means with no SD reported, the SD based on the coefficient of variance (CV) averaged across all other observations where SDs were reported was calculated.
The effect size of N fertilisation (LnRR) on SRR and root N concentration and their variance of the effect size Var(LnRR) (Borenstein et al., 2009; Lajeunesse, 2011) were calculated using Equations (3) and (4), respectively:
where μt and μc are the mean of SRR or root N concentration under fertilisation and control treatments, and σ and N are the SDs and sample size of μt and μc, respectively. An inverse variance method including two variance components (see below) was used to calculate the weight of individual studies (Veroniki et al., 2016). The restricted maximum likelihood (REML) method was used to estimate the between-study variance τ2 (residual heterogeneity), and a Q-profile method was applied to calculate the 95% CIs (Veroniki et al., 2016).
Mixed effect models were used for meta-estimates and meta-regressions to examine the association of categorical and continuous moderators (or predictor as explanatory variables) with the magnitude of the effect size (Terrer et al., 2019; Rubio-Aparicio et al., 2020). We examined how effect sizes differed among categorical factors of Myco type (AM, ECM, AM+ECM, and NM, both at fertilisation levels ≤ 100 and > 100 kg N ha−1), experimental type (field and pot experiments), method (root excavation, trenching, and unplanted), and plant growth form (herbaceous and woody), and continuous factors of Soil moist, Soil temp, and N fertilisation level. Effect sizes and their corresponding variances were calculated in R using the function “escalc” from the package metafor (Viechtbauer, 2010). Because the corresponding effect size estimates for studies with multiple observations (multiple measurements over time) are likely to be correlated (Gleser and Olkin, 2009; Lajeunesse, 2011; Nakagawa and Santos, 2012), we generated a variance-covariance (VCV) matrix and conducted a multi-level mixed effect model that accounts for the non-independency among effect sizes (Noble et al., 2017; Midolo et al., 2019). We used the function “rma.mv” from the package metafor where observation ID was nested within the study as a random effect (second variance σ2 component), and a VCV matrix was used for weighing the effect sizes (Noble et al., 2017; Midolo et al., 2019; Terrer et al., 2021).
A multi-model inference was further conducted by modelling all possible predictor (moderator) combinations to examine which predictor provided the best fit with N fertilisation effect sizes on SRR, based on the sum of Akaike weights (AICc) and p > 0.05 (Burnham, 2002; Gurka, 2006). We further conducted a meta-regression relating LnRR of SRR to N fertilisation with root N concentration in fertilised treatments.
All analyses were done in R, version 4.0.0 (R Development Core Team, 2020). We used the package metafor to conduct the multi-level meta-analysis and meta-regressions (Viechtbauer, 2010) and glmulti for model selection and multi-model inference (Calcagno and de Mazancourt, 2010).
Results
Across 1,600 observations from 32 publications, SRR was not affected by N fertilisation in field or pot experiments (Table 1). Because the majority of observations came from field experiments (1,559 observations), the remainder of our analysis is focused on these studies. Given the large variation in effect sizes among observations, we conducted a multi-model inference to find out the most important moderators on effect sizes of SRR. We included Myco type, Soil moist, Soil temp, N fertiliser level ( ≤ 100 kg N ha−1 and > 100 kg N ha−1), and measurement method. Results of our model selection showed that Myco type and N fertiliser level were the two most important moderators explaining variation in N fertilisation effect sizes, with both moderators showing a sum of AICc larger than the cut-off value 0.8 (Figure 1).
Table 1. Effect sizes (LnRR) of specific root respiration (SRR) and root N concentration to N fertilisation for all observations and separated by field and pot experiments.
Figure 1. Sum of Akaike weights (AICc) from the multi-model inference for predictors including mycorrhizal type (Myco type), N fertilisation levels (N level), soil temperature (Soil temp), gravimetric soil moisture (Soil moist), Specific root respiration measurement method on N fertilisation effect sizes of specific root respiration. Cut-off value (blue line) is set at 0.8 differentiating between important and less important predictors. N.B. only observations from field experiments are included.
We found no significant variation for SRR to N fertilisation among different respiration measurement methods (Table 2) and not between woody and herbaceous plants (Table 3).
Table 2. Effect sizes (LnRR) of SRR to N fertilisation for different respiration measurement methods.
Among the Myco types, SRR increased in plants with dual association (AM+ECM) but was not affected in plants associated with AM only, ECM only, and NM plants only (NM), although NM plants showed a higher variation in the effect size compared to the other Myco types (Figure 2).
Figure 2. Effect sizes (LnRR) of specific root respiration (SRR) to N fertilisation in plants with no mycorrhizal associations (NM), with ectomycorrhizal associations (ECM) only, arbuscular mycorrhizal associations (AM) only, or AM+ECM associations. Numbers in brackets represent the number of observations, and error bars represent the 95% CIs. **indicates significance at p < 0.01. Only observations from field experiments are included.
We then tested whether the effect size of SRR in plants associated with different types of mycorrhizas varied between low ( ≤ 100 kg N ha−1) and high N (> 100 kg N ha−1) fertilisation levels (Figure 3). SRR significantly increased in plants with AM+ECM under low N fertilisation (p = 0.01, Figure 3A). Under high N fertilisation, SRR increased in ECM plants (Figure 3B), while it showed a non-significant increase in plants associated with AM+ECM (but note the low number of observations for AM+ECM association with high variation in effect sizes in response to high N fertilisation). SRRs in AM and NM plants were not affected by either low or high N fertilisation (also note the relative low number of observations and high variation in effect sizes in response to both low and high N fertilisation for NM plants, Figure 3).
Figure 3. Effect sizes (LnRR) of SRR in plants with no mycorrhizal associations (NM), with ECM associations only, AM associations only, or AM+ECM associations, for (A) low N fertilisation levels (N ≤ 100 kg ha−1) and (B) high N fertilisation levels (N > 100 kg ha−1). Numbers in brackets represent the number of observations and error bars represent the 95% CIs. ** and * indicate significance at p < 0.01 and p < 0.05, respectively. Only observations from field experiments are included.
We examined whether differences among Myco type in SRR effect sizes to N fertilisation could be explained by their differences in root N concentration. As expected, root N concentration increased with N fertilisation across all Myco types in field experiments (p = 0.03) but not in pot experiments (Table 1). Plants associated with ECM showed the largest increase in root N concentration with N fertilisation and with smaller increases in root N concentration of AM and NM plants (Figure 4A), while there were no observations for AM+ECM plants. Due to the low number of observations, we were unable to test whether there were differences in the effect sizes of root N concentration between low and high N fertilisation levels. A meta-regression relating SRR affected by N fertilisation to root N concentration in fertilised treatments revealed a positive relationship explaining 44% of the variation (p < 0.001, r2 = 0.44, Figure 4B).
Figure 4. Effect sizes (LnRR) of root N concentration to N fertilisation (A) in association with different mycorrhizal symbiosis and (B) meta-regression relating LnRR of SRR to N fertilisation with root N concentration in fertilised treatments (NM, non-mycorrhizal—red symbols; ECM, ectomycorrhizal—blue symbols; and AM, arbuscular mycorrhizal—grey symbols). ***, **, and *indicate significance at p < 0.0001, p < 0.01, and p < 0.05, respectively. Shade band around the regression fit represents 95% CIs. Only observations from field experiments are included.
Discussion
The mycorrhizal type was the most important predictor of N fertilisation effect sizes on SRR, with no changes for AM, ECM, and NM plants and no increases for AM+ECM plants. N fertilisation level was the second most important moderator of SRR effect sizes. In particular, when both factors were examined together, SRR increased in ECM plants at high N fertilisation rates (>100 kg N ha−1) and for AM+ECM plants at low N fertilisation rates ( ≤ 100 kg N ha−1). While N fertilisation effects on SRR can be due to several factors, such as fertiliser-induced changes in root architecture and morphology (López-Bucio et al., 2003), root N concentrations (Reich et al., 2008), and mycorrhizal colonisation (Nilsson and Wallander, 2003; Treseder, 2004; Emmanuel et al., 2012), in this study, we focus on explaining these effects from a C cost and N benefit perspective of mycorrhizal symbiosis (Phillips et al., 2013).
Measurement method did not play an important role on SRR (Table 2), despite the fact that respiration measured by “trenching” and “unplanted” methods would also include external mycelial respiration (and possibly decomposition of rhizodeposition), while the “root excavation” method would only include internal mycelial respiration in SRR. This suggests that external mycelial respiration and decomposition of rhizodeposition were little affected by N fertilisation (Nilsson and Wallander, 2003; Hobbie, 2008; Wang et al., 2021). Further, we found no variation in SRR to N fertilisation between woody and herbaceous plants (Table 3), indicating that our results were not confounded by plant growth form. Interestingly, Soil temp and moist did not contribute significantly to explaining the variance in SRR effect sizes. Although Soil temp and moisture have been shown as effective factors controlling SRR (Atkin et al., 2000; Zhou et al., 2016), our results highlight that N fertilisation plays a more important role than Soil temp and moist on root respiration (Chen et al., 2019).
From the C cost and N benefit perspective, increased N availability (i.e., with increased N fertilisation) can decrease plant C cost as a result of direct N uptake so that mycorrhizal symbiosis may become too C expensive. A decrease in C cost may particularly occur in ECM plants. For instance, the excretion of extracellular enzymes by ECM fungi to decompose organic matter and release N for plant use (Pritsch and Garbaye, 2011; Nicolás et al., 2019) is C expensive, and it may be more C cost-effective for those plants to take up N directly from the soil when the availability of N increases. Furthermore, the extraradical mycelial networks that ECM fungi build to forage for nutrients in nutrient-poor soils (Chen et al., 2019) require a substantial amount of plant C (Soudzilovskaia et al., 2015) and may not be a C-efficient way for plants to acquire nutrients. Consequently, a reduction in SRR, due to reduced formation of mycorrhizal associations with increased N availability (Nilsson and Wallander, 2003; Treseder, 2004; Emmanuel et al., 2012), is expected as a plant C cost management strategy for ECM plants. To our surprise, we did not observe this but, instead, observed an increase in SRR in ECM plants at high levels of N fertilisation.
There are several possibilities that could explain our results. For instance, high rates of N fertilisation may not have decreased mycorrhizal colonisation and hyphal growth of ECM fungi and associated C cost. Although this is not supported by some studies in ECM trees where N fertilisation caused a decrease in mycelial production (Nilsson and Wallander, 2003; Ekblad et al., 2016), in a meta-analysis Treseder (2004) observed no significant effect of N fertilisation on ECM colonisation. It is also possible that N fertilisation may have changed the mycorrhizal communities (Avis et al., 2003; Parrent et al., 2006) towards species that require more C from their plant host, thereby increasing SRR. It is further possible that N fertilisation increased SRR in ECM plants because of an associated increase in root N concentration. Our results showed that when plants respired more C per unit of root biomass, they also took up more N with increased N fertilisation. The positive relationship between root N concentration and SRR suggests that higher N uptake results in greater C cost, possibly as a result of increased metabolic activities related to nutrient uptake, assimilation, and transportation (Reich et al., 2008; Burton et al., 2012). Therefore, any reductions in C costs due to lower mycorrhizal colonisation may have been outweighed by the increased C costs associated with the higher root N concentrations through increased direct uptake, which resulted in increased SRR at high N fertilisation rates (>100 kg N ha−1) (see conceptual diagram in Figure 5).
Figure 5. Conceptual model illustrating the influence of soil N availability on SRR and root N concentration for plants in association with (A) AM and (B) ECM fungi. With increased soil N availability, both ECM and AM plants spend less C to support mycorrhiza resulting in reduced SRR associated with mycorrhizal symbiosis (a dark blue component of SRR) but respire more C due to an increase in root N concentration (a light blue component of SRR). However, the increase in plants C costs associated with higher root N concentration outweighs the reduction in C to support ECM fungi, while it counterbalances the reduction in C to support AM fungi. Furthermore, in AM plants, C costs for mycorrhizal symbiosis may remain relatively high at high soil N availability, because AM symbiosis is more important for P uptake. We refrained from drawing models for NM and AM+ECM plants due to limited data.
Unlike ECM, our results showed that N fertilisation did not affect SRR in AM plants, under low or high N (Figure 3). One possible explanation is that AM colonisation is less affected by N availability and is more important for P supply so that the fungal contribution to SRR may not be affected by N availability (Hughes et al., 2008). However, AM colonisation on average decreased with N fertilisation by 15% in a meta-analysis (Treseder, 2004), although there was large variation among studies, possibly due to large variation in available N in control treatments. An alternative explanation may be that a reduction in mycorrhizal colonisation with N fertilisation may have been counterbalanced by an increase in root N concentration on SRR (Figure 5). Therefore, in the case of AM plants, reduced C costs via reduced AM colonisation may have been offset by the increased C costs associated with greater plant N uptake (Reich et al., 2008), resulting in no net change in SRR in response to N fertilisation.
It might be expected that SRR declines with N fertilisation in plants with the dual mycorrhizal association (AM+ECM) because supporting both types of mycorrhizal association could be more C expensive for plants than supporting one type. However, our results showed that SRR increased with N fertilisation (significant at N ≤ 100 kg ha−1) in plants with the dual association (Figure 3A). Considering the different capabilities of AM and ECM fungi to acquire nutrients from organic and inorganic sources, they may play a complementary role in N and P uptake (Teste et al., 2020). Increased inorganic N availability with N fertilisation could shift the dominance of mycorrhizal association from ECM acquiring organic forms of N to AM fungi acquiring inorganic N (Phillips et al., 2013), without causing a reduction in mycorrhizal colonisation, but where fertilisation-induced enhancement in root N concentrations would increase respiratory C cost. N fertilisation can also induce a shift from N to P limitation in plants, which could increase plant C cost, given that soil P availability is usually much lower and less mobile compared to N.
Specific root respiration in NM plants was not affected by N fertilisation. If the effects of N fertilisation on SRR are controlled by mycorrhizal associations, then this would be expected. Conversely, if N fertilisation also increased root N concentration, then this would increase SRR. However, the large variation in the effect size of SRR in NM plants among studies precludes making any generalisations. NM plants can be categorised into early colonisers that have no morphological root adaptations and thrive in disturbed or relatively fertile soils, and plants with strong morphological root adaptations, such as cluster or dauciform roots that help plants with P acquisition in severely P-impoverished soils (Lambers and Teste, 2013). These two groups grow in soils that contrast in nutrient availability, and this could explain the large variation in N responses for NM plants, although we did not have enough observations to examine this for low and high N fertilisation levels.
We found that SRR increased in ECM plants supplied with more than 100 kg N ha−1 but not in AM and NM plants, while root N concentrations also increased the most in ECM plants in response to N fertilisation. These results provide insights into why ECM plants tend to dominate in organic or high C:N ratio soils and AM plants in mineral or low C:N ratio soils (Taylor et al., 2016; Jo et al., 2019; Soudzilovskaia et al., 2019). N acquisition through ECM association seems less C expensive for plants when soil N availability is low. This frequently occurs in high C:N ratio soils where most of the soil N is in organic form. Under these conditions, ECM plants may have an advantage over AM plants in acquiring N from organic sources via excretion of extracellular enzymes (Pritsch and Garbaye, 2011; Sulman et al., 2017). With increased inorganic N availability caused by N fertilisation, root N concentrations increase, and at some point this may become too C expensive and no longer economical for plants to support ECM associations. In contrast, supporting AM associations when soil N availability is high may still be beneficial to plants. Although C costs could increase because of an increase in root N concentration, other nutrients like P may become more limiting to plant growth and AM symbiosis may be advantageous (Jo et al., 2019), given that AM fungi play an important role in P acquisition (van der Heijden et al., 2015) and are believed to be less C costly for plants than the ECM symbiosis (Viertelhauzen, 2013). SRR may, therefore, be a key variable in explaining responses of plant growth and their dependence on mycorrhizal associations with N fertilisation and increased atmospheric N deposition. Our meta-analysis further highlights the need for future research studies addressing mechanisms (e.g., biochemical and metabolic pathways) causing variation in SRR and plant N uptake associated with different Myco types in response to N fertilisation.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.
Author Contributions
FD conceived the research idea and participated in the evaluation of the results. BB collected the data, conducted the meta-analysis, performed the model-selection and multi-model inference, and drafted the manuscript. All authors contributed to editing the manuscript.
Funding
This research was supported with a Postgraduate Research Scholarship funded by the Australian Research Council Discovery (Grant No. DP190102262) and the Thomas Lawrance Pawlett Postgraduate Research Scholarship from the Sydney Institute of Agriculture.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2021.711720/full#supplementary-material
References
Atkin, O. K., Edwards, E. J., and Loveys, B. R. (2000). Response of root respiration to changes in temperature and its relevance to global warming. New Phytol. 147, 141–154. doi: 10.1046/j.1469-8137.2000.00683.x
Averill, C., Turner, B. L., and Finzi, A. C. (2014). Mycorrhiza-mediated competition between plants and decomposers drives soil carbon storage. Nature 505, 543–545. doi: 10.1038/nature12901
Avis, P. G., McLaughlin, D. J., Dentinger, B. C., and Reich, P. B. (2003). Long-term increase in nitrogen supply alters above- and below-ground ectomycorrhizal communities and increases the dominance of Russula spp. in a temperate oak savanna. New Phytol. 160, 239–253. doi: 10.1046/j.1469-8137.2003.00865.x
Bache, B. W., Chesworth, W., Chesworth, W., Gessa, C., and Lewis, D. T. (2008). “Bulk density,” in Encyclopedia of Soil Science, ed W. Chesworth (Dordrecht: Springer Netherlands). doi: 10.1007/978-1-4020-3995-9_80
Borenstein, M., Hedges, L. V., Higgins, J. P. T., and Rothstein, H. R. (2009). Introduction to meta-analysis. John Wiley & Sons, Ltd doi: 10.1002/9780470743386
Brundrett, M. C. (2009). Mycorrhizal associations and other means of nutrition of vascular plants: understanding the global diversity of host plants by resolving conflicting information and developing reliable means of diagnosis. Plant Soil. 320, 37–77. doi: 10.1007/s11104-008-9877-9
Burnham, K. P. (2002). Model Selection and Multimodel Inference: A Practical Information-Theoretic Approach. 2nd edn. New York, NY: Springer Netherlands.
Burton, A. J., Jarvey, J. C., Jarvi, M. P., Zak, D. R., and Pregitzer, K. S. (2012). Chronic N deposition alters root respiration-tissue N relationship in northern hardwood forests. Glob. Chang. Biol. 18, 258–266. doi: 10.1111/j.1365-2486.2011.02527.x
Burton, A. J., Pregitzer, K. S., Crawford, J. N., Zogg, G. P., and Zak, D. R. (2004). Simulated chronic NO3– deposition reduces soil respiration in northern hardwood forests. Glob. Chang. Biol. 10, 1080–1091. doi: 10.1111/j.1365-2486.2004.00737.x
Calcagno, V., and de Mazancourt, C. (2010). glmulti: An R package for easy automated model selection with (generalized) linear models. J. Stat. Softw. 1:12. doi: 10.18637/jss.v034.i12
Chen, F., Yan, G., Xing, Y., Zhang, J., Wang, Q., Wang, H., et al. (2019). Effects of N addition and precipitation reduction on soil respiration and its components in a temperate forest. Agric For Meteorol. 271, 336–345. doi: 10.1016/j.agrformet.2019.03.021
Corkidi, L., Rowland, D. L., Johnson, N. C., and Allen, E. B. (2002). Nitrogen fertilization alters the functioning of arbuscular mycorrhizas at two semiarid grasslands. Plant Soil. 240, 299–310. doi: 10.1023/A:1015792204633
Cosme, M., Fernández, I., Van der Heijden, M. G. A., and Pieterse, C. M. J. (2018). Non-mycorrhizal plants: The exceptions that prove the rule. Trends Plant Sci. 23, 577–587. doi: 10.1016/j.tplants.2018.04.004
Ding, W., Yu, H., Cai, Z., Han, F., and Xu, Z. (2010). Responses of soil respiration to N fertilization in a loamy soil under maize cultivation. Geoderma. 155, 381–389. doi: 10.1016/j.geoderma.2009.12.023
Ekblad, A., Mikusinska, A., Ågren, G. I., Menichetti, L., Wallander, H., Vilgalys, R., et al. (2016). Production and turnover of ectomycorrhizal extramatrical mycelial biomass and necromass under elevated CO2 and nitrogen fertilization. New Phytol. 211, 874–885. doi: 10.1111/nph.13961
Emmanuel, B., Fagbola, O., and Osonubi, O. (2012). Influence of fertiliser application on the occurrence and colonisation of arbuscular mycorrhizal fungi (AMF) under maize/Centrosema and sole maize systems. Soil Res. 50, 76–81. doi: 10.1071/SR11254
Gleser, L. J., and Olkin, I. (2009). “Stochastically dependent effect sizes,” in The Handbook of Research Synthesis and Meta-Analysis. 2nd ed, eds H. Cooper, L. V. Hedges, J. C. Valentine (New York, NY: Russell Sage Foundation).
Gurka, M. J. (2006). Selecting the best linear mixed model under REML. Am. Stat. 60, 19–26. doi: 10.1198/000313006X90396
Hasselquist, N. J., Metcalfe, D. B., and Högberg, P. (2012). Contrasting effects of low and high nitrogen additions on soil CO2 flux components and ectomycorrhizal fungal sporocarp production in a boreal forest. Glob. Chang. Biol. 18, 3596–3605. doi: 10.1111/gcb.12001
Hobbie, E. A. (2006). Carbon allocation to ectomycorrhizal fungi correlates with belowground allocation in culture studies. Ecology. 87, 563–569. doi: 10.1890/05-0755
Hobbie, E. A., Ouimette, A. P., Schuur, E. A. G., Kierstead, D., Trappe, J. M., Bendiksen, K., et al. (2013). Radiocarbon evidence for the mining of organic nitrogen from soil by mycorrhizal fungi. Biogeochemistry 114, 381–389. doi: 10.1007/s10533-012-9779-z
Hobbie, S. E. (2008). Nitrogen effects on decomposition: a five-year experiment in eight temperate sites. Ecology 89, 2633–2644. doi: 10.1890/07-1119.1
Hodge, A., and Fitter, A. H. (2010). Substantial nitrogen acquisition by arbuscular mycorrhizal fungi from organic material has implications for N cycling. PNAS 107, 13754–13759. doi: 10.1073/pnas.1005874107
Högberg, M. N., Briones, M. J., Keel, S. G., Metcalfe, D. B., Campbell, C., Midwood, A. J., et al. (2010). Quantification of effects of season and nitrogen supply on tree below-ground carbon transfer to ectomycorrhizal fungi and other soil organisms in a boreal pine forest. New Phytol. 187, 485–493. doi: 10.1111/j.1469-8137.2010.03274.x
Hughes, J. K., Hodge, A., Fitter, A. H., and Atkin, O. K. (2008). Mycorrhizal respiration: implications for global scaling relationships. Trends Plant Sci. 13, 583–588. doi: 10.1016/j.tplants.2008.08.010
Jansa, J., Forczek, S. T., Rozmoš, M., Püschel, D., Bukovská, P., and Hršelová, H. (2019). Arbuscular mycorrhiza and soil organic nitrogen: network of players and interactions. Chem. Biol. Technol. Agric. 6:10. doi: 10.1186/s40538-019-0147-2
Jo, I., Fei, S., Oswalt, C., Domke, G., and Phillips, R. (2019). Shifts in dominant tree mycorrhizal associations in response to anthropogenic impacts. Sci. Adv. 5, 1–8. doi: 10.1126/sciadv.aav6358
Johnson, N. C. (2010). Resource stoichiometry elucidates the structure and function of arbuscular mycorrhizas across scales. New Phytol. 185, 631–647. doi: 10.1111/j.1469-8137.2009.03110.x
Johnson, N. C., and Gehring, C. A. (2007). “Mycorrhizas: symbiotic mediators of rhizosphere and ecosystem processes,” in The Rhizosphere, eds Z. G. Cardon and J. L. Whitbeck (Burlington: Academic Press). doi: 10.1016/B978-012088775-0/50006-9
Kaltenbach, H.-M. (2012). A Concise Guide to Statistics. SpringerBriefs in Statistics. Heidelberg: Springer. doi: 10.1007/978-3-642-23502-3
Lajeunesse, M. J. (2011). On the meta-analysis of response ratios for studies with correlated and multi-group designs. Ecology 92, 2049–2055. doi: 10.1890/11-0423.1
Lambers, H., and Teste, F. P. (2013). Interactions between arbuscular mycorrhizal and non-mycorrhizal plants: do non-mycorrhizal species at both extremes of nutrient availability play the same game? Plant Cell Environ. 36, 1911–1915. doi: 10.1111/pce.12117
Lin, G., McCormack, M. L., Ma, C., and Guo, D. (2017). Similar below-ground carbon cycling dynamics but contrasting modes of nitrogen cycling between arbuscular mycorrhizal and ectomycorrhizal forests. New Phytol. 213, 1440–1451. doi: 10.1111/nph.14206
López-Bucio J. Cruz-Ramirez A. and, Herrera-Estrella, L. (2003). The role of nutrient availability in regulating root architecture. Curr. Opin. Plant Biol. 6, 280–287. doi: 10.1016/S1369-5266(03)00035-9
Makita, N., Kosugi, Y., Dannoura, M., Takanashi, S., Niiyama, K., Kassim, A. R., et al. (2012). Patterns of root respiration rates and morphological traits in 13 tree species in a tropical forest. Tree Physiol. 32, 303–312. doi: 10.1093/treephys/tps008
Midolo, G., De Frenne, P., Hölzel, N., and Wellstein, C. (2019). Global patterns of intraspecific leaf trait responses to elevation. Glob. Chang. Biol. 25, 2485–2498. doi: 10.1111/gcb.14646
Nakagawa, S., and Santos, E. S. A. (2012). Methodological issues and advances in biological meta-analysis. Evol. Ecol. 26, 1253–1274. doi: 10.1007/s10682-012-9555-5
Nicolás, C., Martin-Bertelsen, T., Floudas, D., Bentzer, J., Smits, M., Johansson, T., et al. (2019). The soil organic matter decomposition mechanisms in ectomycorrhizal fungi are tuned for liberating soil organic nitrogen. ISME J. 13, 977–988. doi: 10.1038/s41396-018-0331-6
Nilsson, L. O., and Wallander, H. (2003). Production of external mycelium by ectomycorrhizal fungi in a norway spruce forest was reduced in response to nitrogen fertilization. New Phytol. 158, 409–416. doi: 10.1046/j.1469-8137.2003.00728.x
Noble, D. W. A., Lagisz, M., O'Dea, R. E., and Nakagawa, S. (2017). Nonindependence and sensitivity analyses in ecological and evolutionary meta-analyses. Mol. Ecol. 26, 2410–2425. doi: 10.1111/mec.14031
Parrent, J. L., Morris, W. F., and Vilgalys, R. (2006). CO2-enrichment and nutrient availability alter ectomycorrhizal fungal communities. Ecology 87, 2278–2287. doi: 10.1890/0012-9658(2006)87[2278:CANAAE]2.0.CO;2
Pausch, J., and Kuzyakov, Y. (2018). Carbon input by roots into the soil: Quantification of rhizodeposition from root to ecosystem scale. Glob. Chang. Biol. 24, 1–12. doi: 10.1111/gcb.13850
Phillips, R. P., Brzostek, E., and Midgley, M. G. (2013). The mycorrhizal-associated nutrient economy: a new framework for predicting carbon-nutrient couplings in temperate forests. New Phytol. 199, 41–51. doi: 10.1111/nph.12221
Pregitzer, K. S., Zak, D. R., Loya, W. M., Karberg, N. J., King, J. S., and Burton, A. J. (2007). “The contribution of root - rhizosphere interactions to biogeochemical cycles in a changing world,” in The Rhizosphere, eds Z. G. Cardon and J. L. Whitbeck (Burlington: Academic Press). doi: 10.1016/B978-012088775-0/50009-4
Pritsch, K., and Garbaye, J. (2011). Enzyme secretion by ECM fungi and exploitation of mineral nutrients from soil organic matter. Ann. For. Sci. 68, 25–32. doi: 10.1007/s13595-010-0004-8
Puschel, D., Janouskova, M., Hujslova, M., Slavikova, R., Gryndlerova, H., and Jansa, J. (2016). Plant-fungus competition for nitrogen erases mycorrhizal growth benefits of Andropogon gerardii under limited nitrogen supply. Ecol. Evol. 6, 4332–4346. doi: 10.1002/ece3.2207
R Development Core Team (2020). R: A Language and Environment for Statistical Computing. Vienna: The R Foundation for Statistical Computing.
Reich, P. B., Tjoelker, M. G., Pregitzer, K. S., Wright, I. J., Oleksyn, J., and Machado, J.-L. (2008). Scaling of respiration to nitrogen in leaves, stems and roots of higher land plants. Ecol. Lett. 11, 793–801. doi: 10.1111/j.1461-0248.2008.01185.x
Rubio-Aparicio, M., López-López, J. A., Viechtbauer, W., Marín-Martínez, F., Botella, J., and Sánchez-Meca, J. (2020). Testing categorical moderators in mixed-effects meta-analysis in the presence of heteroscedasticity. J. Exp. Educ. 88, 288–310. doi: 10.1080/00220973.2018.1561404
Smith, F. A., Grace, E. J., and Smith, S. E. (2009). More than a carbon economy: nutrient trade and ecological sustainability in facultative arbuscular mycorrhizal symbioses. New Phytol. 182, 347–358. doi: 10.1111/j.1469-8137.2008.02753.x
Soudzilovskaia, N. A., van Bodegom, P. M., Terrer, C., Zelfde, M., and McCallum, I. (2019). Global mycorrhizal plant distribution linked to terrestrial carbon stocks. Nat. Commun. 10:5077. doi: 10.1038/s41467-019-13019-2
Soudzilovskaia, N. A., van der Heijden, M. G. A., Cornelissen, J. H. C., Makarov, M. I., Onipchenko, V. G., Maslov, M. N., et al. (2015). Quantitative assessment of the differential impacts of arbuscular and ectomycorrhiza on soil carbon cycling. New Phytol. 208, 280–293. doi: 10.1111/nph.13447
Stuart, E. K., and Plett, K. L. (2020). Digging deeper: In search of the mechanisms of carbon and nitrogen exchange in ectomycorrhizal symbioses. Front. Plant Sci. 10:1658. doi: 10.3389/fpls.2019.01658
Sulman, B. N., Brzostek, E. R., Medici, C., Shevliakova, E., Menge, D. N. L., and Phillips, R. P. (2017). Feedbacks between plant N demand and rhizosphere priming depend on type of mycorrhizal association. Ecol. Lett. 20, 1043–1053. doi: 10.1111/ele.12802
Taylor, M. K., Lankau, R. A., and Wurzburger, N. (2016). Mycorrhizal associations of trees have different indirect effects on organic matter decomposition. J. Ecol. 104, 1576–1584. doi: 10.1111/1365-2745.12629
Terrer, C., Jackson, R. B., Prentice, I. C., Keenan, T. F., Kaiser, C., Vicca, S., et al. (2019). Nitrogen and phosphorus constrain the CO2 fertilization of global plant biomass. Nat. Clim. Change. 9, 684–689. doi: 10.1038/s41558-019-0545-2
Terrer, C., Phillips, R. P., Hungate, B. A., Rosende, J., Pett-Ridge, J., Craig, M. E., et al. (2021). A trade-off between plant and soil carbon storage under elevated CO2. Nature 591, 599–603. doi: 10.1038/s41586-021-03306-8
Terrer, C., Vicca, S., Hungate, B. A., Phillips, R. P., and Prentice, I. C. (2016). Mycorrhizal association as a primary control of the CO2 fertilization effect. Science 353, 72–74. doi: 10.1126/science.aaf4610
Teste, F. P., Jones, M. D., and Dickie, I. A. (2020). Dual-mycorrhizal plants: their ecology and relevance. New Phytol. 225, 1835–1851. doi: 10.1111/nph.16190
Thirkell, T., Cameron, D., and Hodge, A. (2019). Contrasting nitrogen fertilisation rates alter mycorrhizal contribution to barley nutrition in a field trial. Front. Plant Sci. 10:1312. doi: 10.3389/fpls.2019.01312
Treseder, K. K. (2004). A meta-analysis of mycorrhizal responses to nitrogen, phosphorus, and atmospheric CO2 in field studies. New Phytol. 164, 347–355. doi: 10.1111/j.1469-8137.2004.01159.x
Treseder, K. K., Allen, E. B., Egerton-Warburton, L. M., Hart, M. M., Klironomos, J. N., Maherali, H., et al. (2018). Arbuscular mycorrhizal fungi as mediators of ecosystem responses to nitrogen deposition: A trait-based predictive framework. J. Ecol. 106, 480–489. doi: 10.1111/1365-2745.12919
Tu, L., Hu, T., Zhang, J., and Li, X. (2013). Nitrogen addition stimulates different components of soil respiration in a subtropical bamboo ecosystem. Soil Biol. Biochem. 58, 255–264. doi: 10.1016/j.soilbio.2012.12.005
van der Heijden, M. G. A., Martin, F. M., Selosse, M.-A., and Sanders, I. R. (2015). Mycorrhizal ecology and evolution: the past, the present, and the future. New Phytol. 205, 1406–1423. doi: 10.1111/nph.13288
Veroniki, A. A., Jackson, D., Viechtbauer, W., Bender, R., Bowden, J., Knapp, G., et al. (2016). Methods to estimate the between-study variance and its uncertainty in meta-analysis. Res. Synth. Methods 7, 55–79. doi: 10.1002/jrsm.1164
Viechtbauer, W. (2010). Conducting meta-analyses in R with the metafor package. J. Stat. Softw. 36, 1–48. doi: 10.18637/jss.v036.i03
Viertelhauzen, L. (2013). Arbuscular vs Ectomycorrhizal fungi: Limiting or Stimulating Tree Growth With N Deposition? MSc thesis, Utrecht University, Utrecht.
Wang, G., and Liu, F. (2014). Carbon allocation of Chinese pine seedlings along a nitrogen addition gradient. For Ecol Manag. 334, 114–121. doi: 10.1016/j.foreco.2014.09.004
Wang, R., Bicharanloo, B., Shirvan, M. B., Cavagnaro, T. R., Jiang, Y., Keitel, C., et al. (2021). A novel 13C pulse-labelling method to quantify the contribution of rhizodeposits to soil respiration in a grassland exposed to drought and nitrogen addition. New Phytol. 230, 857–866. doi: 10.1111/nph.17118
Wang, X.-X., Wang, X., Sun, Y., Cheng, Y., Liu, S., Chen, X., et al. (2018). Arbuscular mycorrhizal fungi negatively affect nitrogen acquisition and grain yield of maize in a N deficient soil. Front. Microbiol. 9:418. doi: 10.3389/fmicb.2018.00418
Zhang, C., Niu, D., Hall, S. J., Wen, H., Li, X., Fu, H., et al. (2014). Effects of simulated nitrogen deposition on soil respiration components and their temperature sensitivities in a semiarid grassland. Soil Biol. Biochem. 75, 113–123. doi: 10.1016/j.soilbio.2014.04.013
Zhou, L., Zhou, X., Shao, J., Nie, Y., He, Y., Jiang, L., et al. (2016). Interactive effects of global change factors on soil respiration and its components: a meta-analysis. Glob. Chang. Biol. 22, 3157–3169. doi: 10.1111/gcb.13253
Keywords: association, carbon cost, expenditure, meta-regression, multi-level, multi-model inference, symbiosis, uptake
Citation: Bicharanloo B, Cavagnaro TR, Keitel C and Dijkstra FA (2021) Nitrogen Fertilisation Increases Specific Root Respiration in Ectomycorrhizal but Not in Arbuscular Mycorrhizal Plants: A Meta-Analysis. Front. Plant Sci. 12:711720. doi: 10.3389/fpls.2021.711720
Received: 19 May 2021; Accepted: 12 July 2021;
Published: 06 August 2021.
Edited by:
Sergio Saia, University of Pisa, ItalyReviewed by:
Néstor Fernández Del-Saz, University of Concepcion, ChileQiang-Sheng Wu, Yangtze University, China
Copyright © 2021 Bicharanloo, Cavagnaro, Keitel and Dijkstra. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Bahareh Bicharanloo, YmFoYXJlaC5iaWNoYXJhbmxvbyYjeDAwMDQwO3N5ZG5leS5lZHUuYXU=