Soil Nutrient Content Influences the Abundance of Soil Microbes but Not Plant Biomass at the Small-Scale

Small-scale heterogeneity of abiotic and biotic factors is expected to play a crucial role in species coexistence. It is known that plants are able to concentrate their root biomass into areas with high nutrient content and also acquire nutrients via symbiotic microorganisms such as arbuscular mycorrhizal (AM) fungi. At the same time, little is known about the small-scale distribution of soil nutrients, microbes and plant biomass occurring in the same area. We examined small-scale temporal and spatial variation as well as covariation of soil nutrients, microbial biomass (using soil fatty acid biomarker content) and above- and belowground biomass of herbaceous plants in a natural herb-rich boreonemoral spruce forest. The abundance of AM fungi and bacteria decreased during the plant growing season while soil nutrient content rather increased. The abundance of all microbes studied also varied in space and was affected by soil nutrient content. In particular, the abundance of AM fungi was negatively related to soil phosphorus and positively influenced by soil nitrogen content. Neither shoot nor root biomass of herbaceous plants showed any significant relationship with variation in soil nutrient content or the abundance of soil microbes. Our study suggests that plants can compensate for low soil phosphorus concentration via interactions with soil microbes, most probably due to a more efficient symbiosis with AM fungi. This compensation results in relatively constant plant biomass despite variation in soil phosphorous content and in the abundance of AM fungi. Hence, it is crucial to consider both soil nutrient content and the abundance of soil microbes when exploring the mechanisms driving vegetation patterns.


Introduction
Resource heterogeneity in space and time is a ubiquitous phenomenon in natural ecosystems and plays a fundamental role in determining plant community dynamics [1,2].Small-scale resource heterogeneity can be especially important in determining plant productivity since plants can concentrate their root biomass in areas of high nutrient content [3][4][5].However, plant growth can also be strongly influenced by soil microorganisms [6].Symbiotic microorganisms, such as mycorrhizal fungi, contribute to plant nutrient uptake, thereby affecting plant growth and biomass allocation [7].Similarly, bacteria play an important role in increasing nutrient availability for plants [6].Despite increasing recognition of the roles played by soil nutrients and especially soil microbes in determining vegetation patterns, information about how soil nutrients, microbes and plant biomass covary in natural ecosystems remains limited.
The spatial distribution of soil microbes is heterogeneous at various scales [8,9] and can be closely related to soil fertility [7,[10][11][12].Less is known about how temporal variation in soil nutrient content [1] influences soil microbial communities, though the abundance of fungi and bacteria has been shown to display significant temporal variation, with maximum abundance in spring and autumn [13].Indeed, the abundance of soil microbes is expected to be high in spring and autumn when nutrient availability is high and decline as the plant growth season proceeds and nutrient availability subsequently decreases [14].However more evidence is needed to better understand these dynamics and the mechanisms driving them.
Studies have found significant small-scale patchiness in the abundance in AM fungi, ectomycorrhizal fungi, other fungi and bacteria [15][16][17] but knowledge about the simultaneous dynamics of different microbial groups remains limited.Both positive and negative associations in the abundance of co-occuring microbes have been reported: in saprotrophic fungi and bacteria [18] and AM fungi and bacteria [19,20].Positive interactions may be the result of one group providing substrate for the other [19], while negative interactions may indicate competition for nutrients [20].A recent study demonstrated that high N:P ratio in soil can provide an advantage, in the form of plant assimilated carbon, to AM fungi, in competition with non-mycorrhizal fungi that also inhabit plant roots [21].
Mycorrhizal fungi contribute to plant nutrient acquisition from the soil [7] and their abundance has been found to decrease with increasing nutrient availability [10,11,22].Hence, mycorrhizal fungi and other soil microbes are expected to play an important role in determining plant productivity in ecosystems with low soil nutrients [6].However, it is unknown whether the effect of soil microbes on plant productivity remains dependent on soil nutrient content at small scales.A recent experimental study did not find any change in the effect of AM fungi on plant productivity in conditions where local soil fertility varied [23], but that study did not consider the effect of microbes besides AM fungi.
We simultaneously measured the temporal and spatial variation of soil nutrient content, soil microbial biomass and plant biomass in a natural herb-rich boreonemoral spruce forest.We focused on a scale that corresponds to the approximate size of adult plant individuals in the understory (i.e.small scale, 15 cm 2 ).The study site is characterized by heterogeneous light and soil conditions, patchy understory vegetation [24], and remarkably taxon rich AM fungal communities in the roots of herbaceous species [25,26].We ask: i) Does soil nutrient content and microbial abundance vary during the first half of the growing season when plant growth is rapid?
ii) Are there negative associations in the small-scale spatial distribution of soil microbial groups?
iii) Is the spatial distribution of soil microbes and plant biomass related to soil nutrient content?

Study site and sampling design
The study was carried out in a boreonemoral forest near Koeru, Estonia (58u589N; 26u039E).The climate of the study site is transitional between maritime and continental.The mean annual precipitation in the area is 665 mm and the mean annual air temperature is 5.0 uC (25.4 uC in January and 16.6 uC in July on average; Estonian Meteorological and Hydrological Institute, 2013).The study site is a 130 ha patch of Hepatica type forest on a calcareous cambisol (cf.[24]).The soil of the area is relatively uniform in terms of pH, ranging from 5.2-5.7,but exhibits some variation in nutrient content [27,28].The tree layer is dominated by Picea abies (L.) Karst.with individuals of Fraxinus excelsior L. and Acer platanoides L. also present.The shrub layer is dominated by Corylus avellana L. Common species in the herb layer are Oxalis acetosella L., Hepatica nobilis Schreb., Paris quadrifolia L. and Viola mirabilis L. The forest can be classified as relatively undisturbed old-growth forest, but clear cutting of approximately 1-2 ha patches has taken place repeatedly in parts of the study site.Detailed descriptions of the area are provided by Moora et al. [24] and Zobel et al. [27].This study complies with the laws of Estonia in which no permits are required to carry out research on public land, unless specified otherwise in legislation.No regulation applies to the Koeru site and no protected species were sampled during the study.
Three 1056105 cm plots (A, B, C) were positioned at 5 m intervals in an area of homogenous forest in terms of vegetation cover and light conditions as estimated by eye.All plots were divided into 49 (15615 cm) quadrats.During preliminary measurements conducted in summer 2007, five plant individuals of the five most common understory species (Oxalis acetosella, Hepatica nobilis, Fragaria vesca L., Galeobdolon luteum Huds., Viola mirabilis) were excavated to a depth of 15 cm, and the root systems of those plants measured with a ruler to determine the size of an average adult plant root system.Soil cores studied during the preliminary study harboured most of the roots in the upper 10 cm of the soil and the average size of plant roots systems was 15615 cm (data not presented).The size of quadrats used in the main experiment therefore corresponds to the average size of the root system of an adult plant (the average of the species listed above), and is expected to harbour the majority of the biomass of herbaceous plant roots and microbes present in the rhizosphere.Plot A was used to address temporal variation, and plots A, B and C were used to address spatial variation in soil nutrient content and soil microbial activity at small scales.Plot A was sampled twice, at the beginning (26.05.2008) and in the middle (21.07.2008) of the growing season.The small spatial scale of the study meant that it was not possible to address temporal variation in the abundance of soil nutrients and microbes over a longer time period, as additional sampling would have necessitated considerable disturbance.Soil samples (50 g and 10 g) were collected from the soil surface (depth 10 cm) at both sampling times for soil chemical, and fatty acid analysis respectively.Soil samples collected at beginning of the plant growing season (in plot A only) were excavated with a soil corer (2 cm diameter and 10 cm deep) at the corner of each quadrat (n = 49) in order to minimize disturbance to the plots.The hole was filled with soil taken from nearby, to minimize sampling effects on soil conditions.Collected soil was mixed (to avoid any effect of vertical structure on measured parameters) and subjected to soil chemical and fatty acid analysis (see below).Soil samples collected in the middle of the growing season (July) were collected during destructive harvesting (see below).One core was taken from each quadrat; in plot A, the locations where the first samples were collected were avoided.The soil corer used during destructive harvesting had a diameter of 15 cm and a depth of 10 cm.All soil cores were transported to the laboratory where, after the removal of soil samples (50 and 10 grams as before) from all of the soil cores (n = 147), the shoots and roots of herbaceous plant species were separated carefully by hand from the samples from plots B and C (n = 98).Plant biomass was dried to a constant weight at 70 uC, and weighed.As the vegetation in plot A was disturbed by the sampling at the beginning of the growing season, biomass was not collected from that plot.

Soil chemical and fatty acid analysis
Soil samples for chemical analyses were kept at 280uC prior to analysis, in order to prevent the loss of soil nutrients.Analysis of soil nutrient content followed International Organization for Standardization (ISO) protocols.Specifically, determination of soil ammonium (NH 4 ) and nitrate (NO 3 ) content (mg/kg) was based on ISO 14256-2: 2005.Organic carbon (C) and total nitrogen (N) content (%) were determined based on ISO 10694: 1995 and ISO 13878: 1998, respectively.As soil NH 4 , NO 3 and C content were highly correlated with total N content, only the latter was used in further analysis.The Mehlich III procedure was used to determine soil phosphorus (P) and potassium (K) content (mg/kg): P was extracted by reaction with acetic acid and fluoride compounds; exchangeable K was extracted by the action of ammonium nitrate and nitric acid [29].Chemical analyses were performed at the laboratory of the Agricultural Research Centre in Saku, Harjumaa, Estonia.
Ester-linked fatty acid (ELFA) content in soil samples was used as a proxy for the biomass of major microbial groups: AM fungi, other fungi and bacteria.Fatty acid profiling is widely used to measure microbial biomass and coarse microbial community composition [10,[30][31][32][33][34][35].The abundance of soil microbes can also be estimated by quantifying their DNA by qPCR [35] but the advantage of fatty acid profiling is that it enables the biomass of different microbial groups to be determined simultaneously.Three grams of fresh soil (randomly taken from collected 10 g) was lyophilized immediately after separation and stored at room temperature.Phospholipids were converted to ester-linked fatty acids and extracted from the soil following the protocol of Schutter & Dick [36], which methylates only ester-linked and not free fatty acids.Briefly, 15 ml of 0.2 M KOH in methanol was added to 2 g of dry soil, mixed and incubated with shaking at 200 rpm at 37uC for 1 h.The sample was neutralized with 3 ml of 1 M acetic acid.Fatty acid methyl esters were extracted by adding 10 ml of hexane, the sample was mixed by vortexing and left to settle into two layers.The upper hexane layer was transferred to a clean tube and evaporated under a vacuum in a Jouan RC1022 Vacuum Concentrator.Fatty acids were re-suspended in 100 ml hexane and transferred to a GC vial for analysis using gas chromatography.
The ester-linked fatty acid peak areas were normalized against the internal standard and expressed as mg per g of dry soil.Fatty acid nomenclature follows Ranneklev & Ba ˚a ˚th [37].Ester-linked fatty acid 16:1v5c was used as a biomarker for AM fungi.This marker is widely used as an indicator for AM fungi [30,32], though the different fatty acid fractions of cells may be of different indicative value for AM fungi.Namely, NLFA (neutral lipid fatty acids, storage lipids) 16:1v5c occurs only in AM fungi, while PLFA (phospholipid fatty acids, structural lipids in cell membranes) 16:1v5c can also occur in some Gram-negative bacteria [38,39].NLFA is therefore proposed as a more suitable indicator of AM fungi in samples with potentially high-bacterial content, such as soils [39,40].However, ELFA 16:1v5c has been shown to correlate very well with AM fungal inoculation in an experiment containing both AM fungi and microbes from the AM fungal inocula, where ELFA 16:1v5c was absent in microbial washtreated, but AM fungal-free soils [41].Therefore, ELFA measurements used in this paper represent whole-cell fatty acids, i.e. they contain both the PLFA and NLFA fraction.Fatty acid 18:2v6,9 was used as a biomarker for fungi other than AM fungi, as it is a dominant fatty acid in most fungi [31] but is only found at very low levels in AM fungi [30].We calculated the ratio between the levels of fatty acids 16:1v5c and 18:2v6,9 (hereafter AM fungi: other fungi ratio) to estimate the relative abundance of AM fungi in the soil fungal community.Fatty acids i15:0, a15:0, 15:0, i16:0, 16:1v9c, br17:0, 10Me16:0, i17:0, a17:0, 17:0, cy17:0, br18:0, 10Me17:0, 18:1v7c, 10Me18:0, cy19:0 v 9, cy19:0 v 7 were used as bacterial biomarkers [31,42].

Statistical analyses
Linear Mixed-Effect Models (LMM) [43] were used to address temporal differences in soil nutrient content (N, P, K) and soil microbial biomass (fatty acid biomarkers for AM fungi, other fungi, AM fungi: other fungi ratio, bacteria).The abundance of other fungi and bacteria as well as the AM fungi: other fungi ratio  was log transformed prior to analysis to meet the assumptions of parametric tests.For each variable, a model using time (sampling time: May or July) as a fixed factor and quadrat identity as a random factor were fitted.The homogeneity assumption for model residuals was evaluated graphically [44].When any sign of heterogeneity was detected in the residuals, it was tested if keeping the variance fixed between the two time periods resulted in homogeneity for the residuals (i.e.including a fixed variance structure in the model; [44]).Models were fitted using Restricted Maximum Likelihood (REML) in order to compare models with different variance structures [44].Moran's I was used to determine the presence of spatial autocorrelation in the distribution of soil microbes.Moran's I calculates the degree of correlation between observations as a function of the spatial distance separating them [45].A correlogram is a graphical representation of the Moran's I statistic against distance classes, and provides information about the spatial autocorrelation of a variable and the size of patches.We constructed a spatial correlogram for each group of soil microbes (AM fungi, other fungi and bacteria) in each plot (A, B and C).Moran's I was evaluated at different distance classes (0 to 50 cm) and its significance was assessed by performing 9999 unrestricted permutations.Since the quadrats were distributed regularly, all correlograms were constructed using distance classes with increments of 0.075 m up to a maximum distance of 0.75 m.The first distance class (0.000-0.075) was discarded since it did not contain any pairs of quadrats.
The Modified t-Test for correlation [46] was used to study pairwise relationships between the abundance of soil microbes: AM fungi, other fungi and bacteria in each plot (A, B and C).This method evaluates the correlation between two spatial patterns using the Pearson's product-moment correlation, but corrects the significance of the correlation for the presence of spatial autocorrelation.
Generalized least square models (GLS) [43] were used to determine which factors best explain variations in soil microbial abundance (in plots A, B, C) and both above-and belowground biomass of herbaceous plants (in plots B, C).The set of predictor variables were different depending on the response variable considered.Specifically, in the models of soil microbial abundance (AM fungi, other fungi, bacteria and AM fungi: other fungi ratio), the influence of plot identity (A, B and C) and soil nutrients (N, P, K) were considered.In the models of above-and belowground biomass of herbaceous plants, the influence of plot identity, soil nutrients and microbial abundance were considered.In order to select the most parsimonious model for each variable, the protocol proposed by Zuur et al. [44] was followed.First, a model which included the full set of covariates as potential predictors was fitted for each variable.After that, the homogeneity of the residuals of the model was evaluated graphically.In addition, due to the nested structure of the sampling design, within plot spatial autocorrelation of residuals was studied.When spatial autocorrelation was present, the inclusion of correlation structures of different form were evaluated [47].Finally, models were fitted with Maximum Likelihood (ML) and the best set of predictor variables was selected using Akaike's Information Criterion (AIC).If several models showed similar AIC values (difference lower than two), the most parsimonious one, i.e. the one with fewer fixed components, was chosen [44].The final model was then refitted using REML to obtain estimates of factor effects and their significance.
LMM, GLS and modified t-Tests and correlograms were performed in the R environment [48], using the packages nlme [49], SpatialPack [50] and ncf [51].Maps of the spatial patterns of AM fungal abundance and soil P content were created using ArcGIS 9.3 software (ESRI Redlands, NY, USA).For each variable, the ''geostatistical analyst'' extension of ArcGIS 9.3 was used to generate semivariograms of density and create kriging maps.

Temporal variation
After the removal of outliers, data from 46 quadrats from plot A was subjected to analysis of temporal variation in soil nutrient content and soil microbial abundance.The linear mixed-effect models (LMMs) revealed significant temporal variation in the abundance of arbuscular mycorrhizal (AM) fungi, bacteria and soil K content, but not in the other variables (Table 1, Figure 1).Means of the measured variables are given in Table S1.The abundance of AM fungal and bacterial fatty acid biomarkers in the soil was higher in May than in July.Soil K content was lower at the beginning of the growing season.The model fitted to explain temporal variation in the abundance of AM fungi included a fixed variance structure as it substantially improved the model fit.

Spatial variation
After the outliers were removed 47, 46 and 44 quadrats from plots A, B and C respectively were subjected to analysis to test the spatial variation of soil nutrient content and the abundance of soil microbes.Limited evidence was found for the presence of spatial autocorrelation in the abundance of soil microbes studied in each plot, since only the abundance of bacteria showed significant autocorrelation (Figure S1).The abundance of AM fungi and soil bacteria tended to show positive spatial autocorrelation up to a distance of 20 cm, but this trend was significant only in the case of bacteria for plots A and C (Figure S1).The Modified t-Test for correlation revealed significant correlations between the abundance of AM fungi and other fungi, AM fungi and bacteria and between bacteria and other fungi (Table 2).
Generalized least square models revealed that plot identity was an important predictor of the abundance of AM fungi (Table 3); the mean abundance was significantly lower in plots A and B compared to plot C (Table 4).In addition, the abundance of AM fungi was significantly positively related to the content of soil N and negatively affected by soil P content (Table 3, Figure 2).Plot identity also played a significant role in predicting the abundance of other fungi (Table 3); the mean abundance was considerably higher in plot C compared to plots A and B (Table 4).In addition, the abundance of other fungi was significantly positively influenced by soil K (Table 3).The AM fungi: other fungi ratio was also predicted by plot identity (Table 3); the mean AM fungi: other fungi ratio differed significantly between all of the plots, being lowest in plot A and highest in plot C (Table 4).In addition, the AM fungi: other fungi ratio was negatively associated with soil P content (Table 3).The abundance of bacteria was the only variable that was best explained by a model where a spatial correlation structure was included, indicating small-scale autocorrelation in bacterial biomass.The abundance of bacteria was explained by plot identity (Table 3), as the mean differed significantly between all plots, being lowest in plot A and highest in plot C (Table 4).In addition, the abundance of bacteria was significantly positively influenced by soil K content (Table 3).
Of all the variables considered, the AM fungi: other fungi ratio was the only one with a significant and positive effect on the shoot biomass of herbaceous plants (F = 5.29; P = 0.02).However, the mean shoot biomass of herbaceous plants did not differ between plots (Table 4).None of the studied variables explained the root biomass of herbaceous plants, which also did not vary significantly between plots (Table 4).

Discussion
The small-scale abundance of soil microbes varied through time in the studied boreonemoral herb-rich spruce forest.We recorded a lower abundance of soil microbes in the middle compared to the beginning of the growing season, whereas no decrease in soil nutrient content occurred.Contrary to our expectations, the smallscale spatial abundance of different soil microbes, including bacteria, arbuscular mycorrhizal (AM) fungi and other fungi, was positively related, suggesting little or no importance of resource competition among these groups in our study system.However, the abundance of soil microbes was significantly influenced by soil nutrient content.In particular, the abundance of AM fungi was strongly negatively related to soil P content.Plant biomass on the other hand showed little variation.Based on these results, we suggest that plants potentially compensate for low soil P content with alternative resource acquisition mechanisms, for example via symbiosis with AM fungi, and therefore achieve equal biomass production despite variations in soil nutrient content and the abundance of soil microbes.Soil microbes (especially plant Table 4.The mean, standard deviation (SD) and the range of the soil nutrient content, ester-linked fatty acid (ELFA) biomarkers of AM fungi, other fungi, bacteria and vegetation characteristics measured at 1 m 2 scale.symbionts such as AM fungi) thus play a fundamental role in the relationship between soil nutrients and plant productivity.The abundance of soil microbes (especially AM fungi and bacteria) was lower in the middle of the plant growing season, which matches the findings of previous studies that have addressed the temporal variation of soil microbes at larger and smaller spatial scales [13,52].Active plant growth during the first half of the growing season did not bring about a decrease in soil nutrient content.Thus, the importance of soil nutrients in determining temporal changes in the abundance of soil microbes seem to be minor [52].Some previous studies have reported a contrasting pattern of AM fungal abundance during the plant growing season: an initial increase with a peak in the summer, followed by a decrease in the autumn [53,54].A potential explanation for this inconsistent variation of AM fungal abundance in soil can be the particular weather conditions of the sampling year, as soil AM fungal hyphal density is considerably affected by temperature and precipitation [54,55].Thus, data from long-term measurements with multiple sampling times during the growing season is needed in order to draw solid conclusions about the dynamics of soil microbes.
We found a strong positive association between the abundance of bacteria, AM fungi and other fungi, suggesting a limited role of resource competition among these microbial groups in a natural ecosystem.Several previous experimental studies have reported negative associations between the abundance of ectomycorrhizal fungi, AM fungi and bacteria [20,56], as well as between the abundance of ectomycorrhizal and saprotrophic fungi [57], and AM fungi and other fungi [21,58].On the other hand, studies conducted in natural conditions have demonstrated similar biomass dynamics of fungi other than AM fungi and bacteria during both the growing season [13] and over decades [59].These contrasting results highlight the importance of studies addressing the pattern of soil microbes in natural environments.The strong positive relationship between the abundance of bacteria, AM fungi and other fungi found in this study suggests that the use of slightly different resources by soil microbes promotes their coexistence.Indeed, AM fungi, which are expected not to be highly saprotrophic themselves, have been shown to promote the activity of saprotrophic microbes [60].In addition, fungal hyphae and spores can function as substrates for bacteria [19,61], and carbon released from dead bacteria can promote the growth of saprotrophic fungi [32].Also, living roots excrete carbon compounds [62], which can be consumed by saprotrophic soil microorganisms [13].In addition, AM fungi as well as ectomycorrhizal fungi have been shown to increase their growth in the presence of root exudates [63,64].
The abundance of all groups of soil microorganisms tested depended on the sampling location (plot identity), which confirms the small-scale spatial heterogeneity of soil microbes previously demonstrated [12].In addition, soil nutrients played an important role in determining the abundance of soil microorganisms.The negative effect of soil P content on the abundance of AM fungi has also been reported in earlier field fertilization experiments in grasslands and forests [65,66].As AM fungi can improve the acquisition of soil nutrients, especially P, by plants [7,67], they can play an important role when soil nutrient content is low.When soil nutrient levels increase, the cost to a plant of supporting AM fungi may outweigh the benefits, and can result in a decrease of AM fungal colonization of plant roots [7].However, soil AM fungal fatty acid content has also been reported to be positively or insignificantly related to soil P content [68,69].We recorded a positive effect of soil N on the abundance of AM fungi, a pattern that has also been registered at a large scale [35] and can be explained by the high N demand of AM fungi [70].In addition, AM fungal communities have been shown to be highly influenced by soil pH [71].However, Dumbrell et al. [71] focused on sites with high variation in soil pH (ranging from 4.5-8), whereas our study area is characterised by fairly uniform pH values [27,28].Therefore, the role of soil nutrients in determining the abundance of AM fungi in current study site could be expected to be more pronounced, while the effect of small-scale variation in soil pH remains to be fully considered in future studies.
Saprotrophic and ectomycorrhizal fungi as well as bacteria were positively influenced by soil K content.A recent study demonstrated an increase in bacterial biomass, and growth of extramatrical mycelia of ectomycorrhizal fungi in response to K addition in peatland soils with severe K deficiency [72].Our results suggest that K also has an important role to play in determining the abundance of soil microbes in soils where K content is not extremely low.On the other hand, the observed positive relationship can also reflect increased weathering and K influx resulting from active growth of bacteria and saprotrophic fungi [73,74].
Plants have been shown to concentrate their biomass into patches with high nutrient content [3][4][5], but in the current study, conducted in an old and stable ecosystem, we found no direct effect of nutrients on plant biomass.One explanation for the absence of the expected pattern is that plants are continuously foraging, while the pattern of biomass measured in studies such as this one reflects dynamics that have operated over numerous years.Greenhouse experiments have shown that nutrient poor patches contain higher plant biomass when next to nutrient rich patches [4,75] but those experiments were short-term.In the longer term, there will be natural turnover of nutrient rich patches (both temporal and spatial), and active foraging can therefore result in relatively uniform plant biomass.Secondly, the natural variation in nutrient content could be too low to generate distinct patterns.Indeed, understory plant biomass in this study area has not shown changes in response to induced variation in soil nutrient content either [23].Thirdly, the effect of soil nutrients can be mitigated by the influence of soil microbes such as AM fungi.AM fungi seemed to be the only group of soil microbes to affect understory plant growth, as the relative increase of AM fungi compared to other fungi resulted in increased plant shoot biomass.As AM fungi can greatly improve plant P acquisition [7,67], the presence of AM fungi can be more important for plant growth when soil P content is low and less beneficial for plants when soil nutrient content is high.In accordance with this, a decrease of AM fungal hyphal density in soil and AM colonization in plant roots has been reported in response to increased P availability [22,76]; whereas the nutrient content of plant tissues and plant biomass did not exhibit any change [22,68].Our study, based on data from a natural environment, demonstrates stable plant biomass in conditions of low abundance of AM fungi and high P content and vice versa.We support the suggestion of Beauregard et al. [68] that plants can compensate for low available nutrient concentration with more efficient symbiosis with soil microbes.Hence, it is crucial to consider both the soil nutrient content and the abundance of soil microbes when exploring the mechanisms driving vegetation patterns.Table S1 Temporal variation in the abundance of soil nutrients and microbes.Summary statistics of the soil nutrient content and ester-linked fatty acid (ELFA) biomarkers of arbuscular mycorrhizal (AM) fungi, other fungi and bacteria measured from plot A at the beginning (May) and in the middle (July) of the growing season.The mean, standard deviation (SD) and the range are presented for each variable.(PDF)

Figure 1 .
Figure 1.The abundance of soil microbes (A) and nutrients (B) in May (open bars) and July (closed bars).The abundance of arbuscular mycorrhizal (AM) fungi, other fungi and bacteria was estimated as the content of respective ester-linked fatty acid (y-axis on graph A, see Methods for details).Soil N concentration was measured in % (left y-axis of graph B) whereas P and K content were measured as mg/kg (right y-axis of graph B).Significant differences (Linear Mixed-Effect Models; p,0.05) in the measured parameters over time are marked with *. doi:10.1371/journal.pone.0091998.g001

Figure 2 .
Figure 2. Spatial pattern of soil phosphorus content and abundance of arbuscular mycorrhizal fungi (ester-linked fatty acid 16:1v5c in soil) in plots A, B and C. Maps were created by kriging, using data from soil samples collected from 15615 cm quadrats in each 1056105 cm plot.doi:10.1371/journal.pone.0091998.g002

Figure S1
Figure S1 Spatial correlograms of microbes.Spatial correlograms of arbuscular mycorrhizal fungi (a-c), other fungi (d-f), and bacteria (g-i) in plot A (a, d, g), B (b, e, h) and C (c, f, i).Distances with significant spatial autocorrelation are marked with solid circles.(PDF)

Table 1 .
Summary statistics of the Linear Mixed-Effect Models analysis fitted to study the influence of sampling time (May or July) on soil nutrient content and ester-linked fatty acid (ELFA) biomarkers of arbuscular mycorrhizal (AM) fungi, other fungi and bacteria measured from plot A.

Table 2 .
Results of the modified t-Test for correlation performed to study the relationship between the abundance of arbuscular mycorrhizal (AM) fungi, other fungi and bacteria in the three plots (A, B, C).

Table 3 .
Results of the generalized least square models performed to study the influence of environmental conditions on soil microbes.
For each soil microbial group, the coefficient associated with the explanatory variable is presented.In addition, the AIC weight of the regression model (i.e. its importance as compared to other models containing a different subset of explanatory variables) and the R 2 are shown.* P,0.05; ** P,0.01.doi:10.1371/journal.pone.0091998.t003