Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Open Access

Peer-reviewed

Research Article

Prevailing Negative Soil Biota Effect and No Evidence for Local Adaptation in a Widespread Eurasian Grass

Abstract

Background

Soil biota effects are increasingly accepted as an important driver of the abundance and distribution of plants. While biogeographical studies on alien invasive plant species have indicated coevolution with soil biota in their native distribution range, it is unknown whether adaptation to soil biota varies among populations within the native distribution range. The question of local adaptation between plants and their soil biota has important implications for conservation of biodiversity and may justify the use of seed material from local provenances in restoration campaigns.

Methodology/Principal Findings

We studied soil biota effects in ten populations of the steppe grass Stipa capillata from two distinct regions, Europe and Asia. We tested for local adaptation at two different scales, both within (ca. 10–80 km) and between (ca. 3300 km) regions, using a reciprocal inoculation experiment in the greenhouse for nine months. Generally, negative soil biota effects were consistent. However, we did not find evidence for local adaptation: both within and between regions, growth of plants in their ‘home soil’ was not significantly larger relative to that in soil from other, more distant, populations.

Conclusions/Significance

Our study suggests that negative soil biota effects can prevail in different parts of a plant species' range. Absence of local adaptation points to the possibility of similar rhizosphere biota composition across populations and regions, sufficient gene flow to prevent coevolution, selection in favor of plasticity, or functional redundancy among different soil biota. From the point of view of plant - soil biota interactions, our findings indicate that the current practice of using seeds exclusively from local provenances in ecosystem restoration campaigns may not be justified.

Citation: Wagner V, Antunes PM, Ristow M, Lechner U, Hensen I (2011) Prevailing Negative Soil Biota Effect and No Evidence for Local Adaptation in a Widespread Eurasian Grass. PLoS ONE 6(3): e17580. https://doi.org/10.1371/journal.pone.0017580

Editor: Justin Wright, Duke University, United States of America

Received: October 19, 2010; Accepted: February 9, 2011; Published: March 29, 2011

Copyright: © 2011 Wagner et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Funding: This study was made possible through a grant from the German Academic Exchange Service, DAAD, to V.W. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

Introduction

Soil biota effects are increasingly recognized as an important factor driving the abundance and distribution of plants [1]–[3]. These reciprocal interactions can result in different functional outcomes; from neutral, mutualistic to parasitic [4], [5]. By influencing seedling germination, plant growth and fitness, soil biota effects are an important factor controlling community composition, species diversity, and successional dynamics [6]–[9].

There is accumulating evidence that the functional outcomes of plant and soil biotic interactions within a habitat derive from coevolved relationships. This is supported by studies on invasive species, in which plant performance was compared in soil from within (‘home soil’) or outside the distribution range (‘away soil’). Several species were shown to be suppressed by their native ‘home soil’ and to be released from its negative effects when grown in soil from the introduced range, e.g. Centaurea maculosa [10], Ammophila arenaria [11], Carpobrotus edulis and Carpobrotus×cf. acinaciformis [12], and Tragopogon dubius [13]. Exotic plants can also be limited in their naturalization if suitable symbionts are not available in the introduced range, as documented for Cytisus scoparius [14], and Pinaceae species [15].

While these studies indicate that plants have coevolved with their native soil biota, it is less clear to what extent local adaptation varies among populations within the native range of a plant species [4], [16]–[19]. Evidence for local adaptation to soil biota was found in the case of Trifolium repens and Rhizobium bacteria [20], Amphicarpaea bracteata and Bradyrhizobium bacteria [21], Allium vineale and Uniola paniculata to arbuscular mycorrhizal fungi [22], [23]. Contrasting results have been reported for Acacia species and Rhizobium bacteria [24], Pinus and Rhizopogon fungi [25], Ammophila arenaria and nematodes [26]. Most of these studies were however limited because they did not consider the hypothesis that adaptation might be manifested at different spatial scales [27]–[29]. Generally, these studies were confined to specific organism groups whose effects may be overridden by more important soil organism groups under field conditions [13]. Furthermore, they tested for direct interactions between pairs of species, leaving out putative additive or net effects in interactions with soil biota.

The question of local adaptation between plants and their soil biota has far-reaching implications for studies on plant-soil biota interactions and land management. In restoration campaigns, for example, much effort is spent in the common practice of using seed material from local provenance [30]. If local adaptation between plants and soil biota exists, this would be one more justification in favor of this effort [31].

In this study, we tested whether plants are locally adapted to their soil biota by carrying out a reciprocal transplant experiment in the greenhouse; a recommended approach to examine local adaptation [32]. We used plant and soil material from populations located in two different regions of a species' range, and addressed our question at two spatial levels, within (ca. 10–80 km) and between (ca. 3300 km) two regions. We chose the steppe grass Stipa capillata L. as a model species because S. capillata is an important element of Eurasian dry grasslands [33] and grasses have been shown to be sensitive to soil biota [2]. Furthermore, this species has a wide distribution range in Eurasia, which allows testing for local adaptation at a large scale [34]. We used seed and soil material from populations in the centre of its distribution range, in Asia, and from the range periphery, in Europe. Peripheral populations were spatially and genetically more isolated than those in the central range (V. Wagner, W. Durka, I. Hensen, unpublished data), increasing the possibility for local coevolution with soil biota. Furthermore, as the two regions were 3300 km apart and differed climatically [34] we hypothesized local adaptation to be stronger between than within regions.

Materials and Methods

Study Species

Stipa capillata L. (Poaceae) is a perennial tussock grass that grows in dry grasslands on nutrient poor, sandy to loamy soil. Its native range covers large areas of Eurasia and, in the core, in the steppes of Russia and Kazakhstan (Figure 1a), it is one of the most common plant species [33]. However, at its north-western range periphery in Central Europe (Figure 1a) S. capillata is rare, isolated and red-listed [34]. The species is recognized by its inflorescence, bearing single stalked florets with a long and naked awn. Flowers are known to be facultative cleistogamous [35]. Caryopses ripen and are dispersed in animal fur or by wind in late August to October. Mycorrhizal colonization in roots of S. capillata was reported from Russia by Mukhin and Betekhtina [36] and from China by Shi et al. [37].

thumbnail
Download:
Figure 1. Map of the study localities.

A) Overview of the study regions within the native distribution range of Stipa capillata (grey color). Detailed map of the study regions in (B) Europe and (C) Asia. Population localities are marked as black stars.

https://doi.org/10.1371/journal.pone.0017580.g001

Study Sites

We collected seeds and soil from five populations within each of the two regions (i.e., Europe and Asia) that were approximately 3300 km apart. European populations, in the lower Oder basin, in Brandenburg, Germany (Figure 1B) were 9.7–82 km apart from each other (mean: 46.6 km). Plants were located in grassland fragments on slopes that were separated by fields, settlements and forests. Climate in this region is mild, with annual precipitation of 540 mm (interpolated climatic data based on [38]), cool summers (18°C mean temperature) and warm winters (−1°C). Soils in European populations consisted of nutrient poor sandy brown earths on glacial and periglacial sand deposits (see Table S1 for soil chemical properties). Asian populations, in the Naurzum district, Kostanayskaya Oblast, Kazakhstan (Figure 1C), were 9.6–69.6 km apart from each other (mean: 39.1 km) and covered flat plains with few settlements and forests in between. Climate in this region is dry with 260 mm annual precipitation, warm during the summer (22°C) and cold in the winter (−13°C). Soils were comprised of nutrient poor sandy chestnut soils on sand and loam sediments (see Table S1 for soil chemical properties).

Collection of Soil and Seeds

Soil and seeds were collected in September 2008 within 10×10 m plots; each considered as representing a population. For each population, we collected ca. 100 g soil with a metal corer (diameter: 2 cm, length: 20 cm) from underneath six randomly chosen tussocks, respectively, thus obtaining a total of approximately 600 g of soil for each population. The soil corer was thoroughly washed with water before collecting soil from another population. We collected and pooled seeds from all individuals within each plot. Soil was air-dried and transported overland in closed plastic bags to the University of Halle-Wittenberg, Germany. The soil was sieved (2 mm mesh), roots were cut into 1 cm long pieces and this substrate was refrigerated at 5°C before further processing. Seeds were also stored at 5°C.

General Experimental Setup

For our greenhouse experiment, we used a sterilized (autoclaved at 121°C three times for 1 h) 1∶4 sand∶loam mixture as background soil. By using a sterile background soil, we minimized the presence of contaminant soil biota in our pots and were able to provide equal soil abiotic growing conditions for our study plants. Seeds were surface sterilized (2 min in 50% sodium hypochlorite bleach, 2 min in 70% alcohol) to eliminate adherent soil biota. We pre-germinated seeds on heat sterilized (36 h, 200°C) sand and planted individual seedlings into 450 ml cone-tainer pots (Stuewe and Sons Inc., Corvallis, OR, USA). Similar cone-tainer pots were used in subsequent experiments. To avoid cross-contamination over the course of our experiment, we sterilized all equipment by flame or by soaking in bleach for one hour.

The experiment was conducted at the greenhouse facilities of the University of Halle-Wittenberg, from March to October 2009. Before conducting Experiments 1–3, we trained soil to build up local soil biota in the sterile background soil. By choosing this approach we avoided maternal effects of the soil biota and equalized abiotic differences among soil samples [39]. For this training stage, we planted ten seedlings from each of the ten populations, with one seedling per pot containing a 9∶1 mixture of sterilized background and field-collected soil from their own population. We watered plants once or twice a week individually using equal amounts of deionized water, being cautious to avoid leaching. After 13 weeks, we removed the above ground biomass and mixed the trained soil including chopped roots into one bulk sample for each respective population. Trained soil was later used as inoculum (15% per pot) aimed at testing general soil biota effects and local adaptation in plants. Using whole soil as inoculum is a common tool in soil biota studies (see [10], [13], [40] for similar approaches) given the technical challenge of teasing apart the effects of a large variety of different organisms in the soil [3], [41]. Moreover, this strategy allows studying overall soil biota effects on plants similar to what plants experience in the field [2], [13]. Using a low inoculum dosage excluded the confounding effects of abiotic soil properties, a tactic similarly used by Bever [40].

Experiment 1: General Soil Biota Effect

The first treatment aimed to test the general effect of soil biota on plant growth in all populations. We grew plants in sterilized background soil either with soil inoculum from their own population (‘home soil’) or without soil inoculum (‘control soil’). A total of 200 pots were used: 2 treatments×5 populations×2 regions×10 replicates. All pots were distributed randomly on the greenhouse bench using 28 supports that were rotated every third week. To avoid contamination, pots were arranged in a chess-like manner within each support and watered individually to avoid splash transfer. Initially, two seedlings were planted. After three weeks, we removed the smaller seedling and started measuring the maximum height of the remaining plant at three week intervals to estimate growth rate. After 19 weeks, we harvested the shoot and root biomass. Aboveground material was dried in an oven at 80°C for 48 h. Roots were rinsed, dried in silica gel and stored in plastic bags with silica gel. For each plant, we measured dry weight (shoot and root). Furthermore, we quantified mycorrhizal colonization, including arbuscules, in roots of five plants per population that were grown in their ‘home soil’ or in ‘control soil’ according to the method by Trouvelot et al. [42].

Experiment 2: Local Adaptation Within Regions

The purpose of the second treatment was to test for local adaptation within the two regions, respectively. We compared fully crossed plant - soil inoculum combinations from different source populations within the regions by growing plants in increasingly distant ‘away soils’. This approach allowed testing for local adaptation in a clinal manner by evaluating plant performance as a function of the continuous distance between plant and soil inoculum origin [25]. We used plants grown in ‘home soil’ from Experiment 1 for the distance of 0 km and grew an additional number of 400 plants: 4 treatments ( = 4 additional distances)×5 populations×2 regions×10 replicates. Planting of seedlings, handling of the experiment and harvest were as in Experiment 1.

Experiment 3: Local Adaptation Between Regions

In a third treatment, we tested adaptation between Europe and Asia, by comparing growth of plants in ‘home soil’ and ‘away soil’ from the other region. In the latter case, we did not fully cross all plant and soil populations but for each plant population chose randomly (without replacement) soil from one population from the other region as soil inoculum. We used the following pairs: EU1/KZ3, EU2/KZ6, EU4/KZ7, EU5/KZ2, and EU6/KZ1. We used plants grown in ‘home soil’ from Experiment 1 and compared them to 100 additional ‘away soil’ pots: 5 populations×2 regions×10 replicates. Planting of seedlings, handling of the experiment and harvest were as in Experiment 1.

Statistical Analysis

Soil biota effects can be analyzed in two ways, using either a net soil biota effect index or original dry weight values. Because both approaches have advantages and disadvantages [39], we used both methods in our analysis.

First, we performed an analysis with a net soil biota effect index similar to the index of relative interaction intensity in Armas et al. [43] as (XNS−XS)/(XNS+XS), with XNS being the dry weight of a plant in non-sterile soil and XS the dry weight of a plant in sterile soil. As we did not use paired sterile and non-sterile replicates, we employed a bootstrap procedure to estimate the net soil biota effect [44]. For that purpose, we calculated mean XNS and mean XS based on 10 bootstrapped samples, respectively, in the boot package [45] in R [46]. These values were incorporated to calculate a net soil biota effect. We repeated this calculation 999 times. We used a t-test to compare the general effect of soil biota between regions (Experiment 1). Local adaptation within regions was tested with linear mixed-effects models using the nlme package in R [47], with mean net soil biota effect as a response variable, distance and region as fixed predictor variables (Experiment 2). As plant and soil populations were crossed in a fully reciprocal manner we were able to analyze local adaptation in plants growing in soil from increasingly distant populations, and in plants of different origins growing in a soil population [48]. Thus we incorporated plant or soil population as random factors, respectively, and tested whether their interaction with distance improved the model fit [49]. To analyze local adaptation between regions (Experiment 3), we used a linear generalized least squares model, with mean net soil biota effect as a response variable, region and soil treatment as predictor variables. Using linear mixed-effects models, we tested whether plant population and the interaction between plant population and soil treatment as random effects improved the model fit.

Secondly, we analyzed original dry weight per plant (sum of shoot and root weight) with linear mixed-effects models. Throughout all analyses, we used soil treatment, region and supports (block) as fixed effects. If the variable block and its interactions were not found to significantly improve the model they were removed from the final model. We used soil treatment as a categorical factor when analyzing the general effects of soil (‘home soil’ vs. ‘control soil’, Experiment 1) and local adaptation between regions (‘home soil’ vs. ‘away soil’ from the other region, Experiment 3). In the case of local adaptation within regions (Experiment 2), plant and soil populations were fully reciprocal among populations, so that distance between plant and soil populations (‘away soil’) was used as a continuous predictor and local adaptation was tested from two perspectives (see above). Model selection followed the methodology proposed by Zuur et al. [49]. Frequency of mycorrhizal fungi in the root system was calculated in the program MycoCalc according to the method proposed by Trouvelot et al. [42].

We used repeated measures analysis to analyze the influence of soil treatment and region on plant height during the first nine weeks of the experiment. We chose this time frame because the linear plant height – time relationship became asymptotic after this date across all treatments. We used generalized least squares and linear mixed-effects models in the package nlme, as implemented in the program R, to construct our models. In the case of local adaptation within regions (Experiment 2), we treated distance between soil and plant origin again as a continuous variable and analyzed data from two perspectives (see above). To account for possible temporal pseudoreplication and variance heterogeneity through time, we tested models with different correlation and variance-covariance structures by inspecting the AIC [50]. The autoregressive heterogeneous variance-covariance structure had the best model fit accounting for temporal correlation among individual measurements and different variances at each time point. The covariable block was not found to significantly improve the model and was omitted from the analysis.

P values of random effects were calculated by a likelihood ratio test and by comparing the model with and without the random effect [49].

Results

General Soil Biota Effect (Experiment 1)

Net soil biota effect was negative in the majority of populations in Europe and Asia (Fig. 2). In each region, plants growing in ‘control soil’ produced significantly more dry weight than plants growing in ‘home soil’ (Fig. S1). Individuals from Europe produced generally more dry weight than those from Asia, but not for the ‘control soil’ (soil treatment * region, Fig. S1). Plants in ‘control soil’ also grew faster over the course of the first nine weeks compared to plants in ‘home soil’ (Time×soil treatment, Table S3). Roots of plants inoculated with ‘home soil’ had a significantly larger frequency of mycorrhizal colonization than those grown in ‘control soil’ (mean = 27.6% vs. 1.3%; Kruskal-Wallis Χ = 52.3, degrees of freedom (d.f.) = 1, P<0.001).

thumbnail
Download:
Figure 2. General soil biota effect in populations from Europe and Asia (Experiment 1).

Bars represent estimated mean net soil biota effect, calculated as (XNS−XS)/(XNS+XS), with XNS being the dry weight of a plant in home soil and XS the dry weight of a plant in sterile control soil. As XNS and XS were not paired during the experiment, we used mean XS and XNS values based on 10 bootstrapped samples when calculating the net soil biota effect and repeated this procedure 999 times to estimate the mean net soil biota effect. Stars indicate whether in 95% (***) or 90% (**) of bootstrap runs effect values were different from zero. Mean net soil biota effect was not significantly different between Europe and Asia (t-test, t = 0.74, d.f. = 8, P = 0.480).

https://doi.org/10.1371/journal.pone.0017580.g002

Local Adaptation Within Regions (Experiment 2)

The net soil biota effect on plants did not change significantly with geographic distance, in both regions (Fig. 3, Table 1A). When comparing how plants from different populations performed in a given soil population, there was similarly no effect of geographic distance, region and their interaction (Table 1B). Dry weight (Fig. S2, Table S2) and growth rate (time×soil treatment, Table S4) did similarly not change significantly in response to geographic distance. Individuals from Europe produced generally more dry weight than those from Asia but this effect did not change with geographic distance (geographic distance×region, Table S2).

thumbnail
Download:
Figure 3. Estimated mean net soil biota effect as a function of geographic distance and region (Experiment 2).

Points represent estimated mean net soil biota effect for Europe (solid) and Asia (open). Lines (solid for Europe and dashed for Asia) show the fitted linear mixed-effects model as predicted by the fixed effects geographic distance, region and their interaction for plants growing in soil from increasingly distant populations (random effects not shown). Geographic distance, region and their interaction were not found to significantly influence mean net soil biota effect (see Table 1A for details of the statistical test).

https://doi.org/10.1371/journal.pone.0017580.g003

thumbnail
Download:
Table 1. Linear mixed-effects models to test for local adaptation within regions (Experiment 2).

https://doi.org/10.1371/journal.pone.0017580.t001

Local Adaptation Between Regions (Experiment 3)

There was no significant difference in estimated mean net soil biota effect when plants from Europe and Asia were grown in their ‘home soil’ or in ‘away soil’ from the other region (Fig. 4). Similarly, we did not find a significant difference in dry weight (Figure S3) and growth rate (Table S5) between plants grown in ‘home soil’ or ‘away soil’ from the other region.

thumbnail
Download:
Figure 4. Mean net soil biota effect of plant populations in Europe and Asia grown in ‘home soil’ and ‘away soil’ from the other region (Experiment 3).

Region, soil treatment and their interaction were not found to significantly influence mean net soil biota effect, as shown by a linear model: FRegion = 0.01, d.f. = 1,16, P = 0.912, FSoil treatment = 0.11, d.f. = 1,16, P = 0.749, FRegion×soil treatment = 0.53, d.f. = 1,16, P = 0.475. Including plant population or the interaction of plant population with soil treatment as random effects in a linear mixed-effects model did not significantly improve the fit of the model.

https://doi.org/10.1371/journal.pone.0017580.g004

Discussion

General Soil Biota Effect

Our study shows that Stipa capillata is controlled by negative soil biota effects: mean net soil biota effect was predominantly negative, dry weight and growth rate decreased significantly when plants were grown in their ‘home soil’ as compared to the sterile ‘control soil’. Thus, our study is consistent with a growing body of literature that reported negative soil biota effects on plants, especially grasses [2], [5], [8], [51]. Furthermore, our study is one of the first to show that negative soil biota effects can prevail in different parts of a species' range. A recent study by van Grunsven et al. [13] came to the same conclusion and showed that Tragopogon pratensis experienced negative soil biota effects throughout its range in Europe.

It is unlikely that our observed differences were caused by a higher proportion of background soil in control pots and thus by a higher nutrient content that arises from the process of sterilization. We minimized this putative confounding factor by keeping inoculum dosage low and pre-training our soil inoculum [40]. Thus, nutrients should have been similar across treatments. Furthermore, since inoculum dosage was low, it is unlikely that allelochemicals played a role in diminishing plant performance.

Our approach of using field collected whole soil as inoculum enabled testing for soil biota effects on plants as similar to field conditions as possible. It was not our goal to pinpoint the exact components of the soil biota responsible for the observed negative effect. However, our results call for future, more mechanistic, investigations to identify rhizosphere pathogens of grasses and their specific roles. Generally, many different groups of soil organisms can have a negative effect on plants belowground, including bacteria, protozoa, and nematodes [41]. Nevertheless, attempts to narrow down general soil organism groups responsible for negative effects have implicated fungi [5] and oomycetes in particular for being important [52].

Local Adaptation Within and Between Regions

We found no evidence for local adaptation of S. capillata to soil biota. Plants grown in their ‘home soil’ and in ‘away soil’ did not experience a significantly different net soil biota effect index, dry weight, or growth rate. Contrary to our hypothesis, no local adaptation was detected at a larger scale, when plants were reciprocally transplanted between the two regions. This result is striking given that the study regions were 3000 km apart and differ drastically in macroclimatic properties that were likely to have favored different soil biota.

The lack of evidence for local adaptation in our experiment may be explained by several non-mutually exclusive causes. First, it is possible that soil biota composition was similar within and between the regions and, as a consequence, selection strength for locally adapted plants weak [53], [54]. Second, even if soil biota was differentiated across different localities [55], [56], such differences did not necessarily have to result in significant plant growth effects. Soil biota is compromised of a high diversity of species and strains [57], many of which might be functionally redundant [58]. Furthermore, past temporal variation in soil biota composition or even moderate migration rates could have enabled plants to respond to a variety of different soil organisms [32], [59]. Alternatively, it has been suggested that in species rich communities, coevolutionary responses are not governed by single pairwise interactions as much as by the diffuse interactions of multiple species [60], [61]. Thus, putative effects by some soil organisms on plant performance could have been overridden by other species [62]. From the plant perspective, gene flow among populations could have counteracted selection to local soil biota [32], [53], [54]. Although peripheral populations in Europe were genetically more differentiated than those in Asia (V. Wagner, W. Durka, I. Hensen, unpublished data) it is possible that even diminished gene flow was strong enough to counteract selection in this region. In addition, if genotype×genotype interactions are constrained by the environment (genotype×genotype×environment, [27], [32], [41], [63]), then experimental greenhouse conditions might have masked signs of local adaptation in our study. Therefore, future studies on local adaptation to soil biota should explore the role of environmental factors.

Plants can be locally adapted to a variety of environmental factors, including soil organisms [32]. Local adaptation is a significant consideration for land managers, for example when plant material is to be introduced to a target site. However, little has been done to test whether plant material needs to be derived from the same locality to perform best with soil biota at a restoration site. Our study provides a first indication that plants may not consistently be adapted to their local soil biota, so that this factor could potentially be precluded in management decisions. However, further studies are needed with other equally ecologically relevant plant species that take into consideration plant-plant competitive interactions [64], [65] as well as multiple abiotic factors.

Conclusions

We found persistent negative soil biota effects in a widespread Eurasian grassland species, in two climatically contrasting parts of its native distribution range. Contrary to our hypothesis, we did not find evidence for local adaptation of Stipa capillata to soil biota, neither at the scale of 10–80 km, nor at the scale of 3300 km. Very little is known about coevolution and the roles of pathogen rhizosphere communities on their host plants and aboveground plant communities. Moreover, there is a need for further research that explores the role of environmental factors in soil biota effects, especially in a community context. Nevertheless, based on our data, when specifically considering possible effects of soil biota, managers should take into account that local adaptation to soil biota may not necessarily be present in plants.

Supporting Information

Figure S1.

Dry weight of plants from populations in A) Europe and B) Asia grown in ‘home soil’ and in sterile ‘control soil’ (Experiment 1). A linear mixed-effects model showed that soil treatment and region had a significant effect on dry weight: FSoil treatment = 25.62, degrees of freedom (d.f.) = 1,187, P<0.001, FRegion = 5.63, d.f. = 1,8, P = 0.045, FSoil treatment×region = 0.10, d.f. = 1,187, P = 0.748. Plant population was used as a random effect (P<0.001).

https://doi.org/10.1371/journal.pone.0017580.s001

(TIF)

Figure S2.

Dry weight of plants as a function of geographic distance and region (Experiment 2). Points represent dry weight of plants for Europe and Asia, respectively. Lines show the fitted linear mixed-effects model as predicted by the fixed effects geographic distance, region and their interaction for plants growing in soil from increasingly distant populations. Plant population was used as a random effect (lines not shown). Distance, region and their interaction were not found to significantly affect dry weight per plant (see Table S2 for details of the statistical test).

https://doi.org/10.1371/journal.pone.0017580.s002

(TIF)

Figure S3.

Dry weight of plants from populations in A) Europe and B) Asia grown in ‘home soil’ and ‘away soil’ from the other region (Experiment 3). A linear mixed-effects model showed that soil treatment, region and their interaction did not significantly influence dry weight production in plants: FSoil treatment = 0.13, d.f. = 1,184, P = 0.718, FRegion = 4.8, d.f. = 1,8, P = 0.059, FSoil treatment×region = 0.73, d.f. = 1,184, P = 0.393, FBlock = 4.0, d.f. = 1,184, P = 0.048. Plant population×soil treatment was used as a random effect (P<0.001).

https://doi.org/10.1371/journal.pone.0017580.s003

(TIF)

Table S1.

Chemical properties of field collected soil and sterilized background soil used in the experiment.

https://doi.org/10.1371/journal.pone.0017580.s004

(DOC)

Table S2.

Linear mixed-effects models to test for local adaptation within regions (Experiment 2). Dry weight was used as a function of the fixed effects geographic distance, region (and block). Geographic distance denotes the distance between the population of the plant and the soil inoculum. Models were constructed in two ways: A) from the perspective of plant populations growing in increasingly distant soil populations, by including plant population as a random effect, and B) from the perspective of soil populations in which plants from increasingly distant populations were grown, by including soil population as a random effect. P values≤0.01 are in bold.

https://doi.org/10.1371/journal.pone.0017580.s005

(DOC)

Table S3.

Repeated measures analysis of plant growth during Experiment 1. Plant height was measured on individuals every third week within a nine week period and used as a response variable of the fixed effects time, soil treatment (‘home soil’, ‘control soil’), region and their interaction. Including plant population or the interaction between plant population and time as random effects did not significantly improve the fit of the model. P values≤0.01 are in bold.

https://doi.org/10.1371/journal.pone.0017580.s006

(DOC)

Table S4.

Repeated measures analysis of plant growth during Experiment 2. Plant height was measured on individuals at the beginning and every third week of a nine week period and used as a response variable of the fixed effects time, geographic distance, and their interaction. Geographic distance denotes the distance between the population of the plant and the soil inoculum. We used population and individual nested within population as random factors. Models were constructed in two ways: A) from the perspective of plant populations growing in increasingly distant soil populations, by including plant population as a random effect, and B) from the perspective of soil populations in which plants from increasingly distant populations were grown, by including soil population as a random effect. P values≤0.01 are in bold.

https://doi.org/10.1371/journal.pone.0017580.s007

(DOC)

Table S5.

Repeated measures analysis of plant growth during Experiment 3. Plant height was measured on individuals at the beginning and every third week of a nine week period and used as a response variable of the fixed effects time, soil treatment (‘home soil’, ‘away soil’ from the other region), region and their interaction. Population and individual nested within population were used as a random intercept. P values≤0.01 are in bold.

https://doi.org/10.1371/journal.pone.0017580.s008

(DOC)

Acknowledgments

We would like to thank the NGO Aq Tyrna and Hans Pfestorf for support with field work in Kazakhstan and Europe, respectively, Christine Voigt for assistance during soil analysis and the experiment, Eva Bremer, Heidi Hirsch, Sandra Nähring, Marco Seifert, Sabrina Träger and Heike Zimmermann for help with greenhouse work. We are also thankful to anonymous reviewers who helped to improve the manuscript.

Author Contributions

Conceived and designed the experiments: VW PA MR IH UL. Performed the experiments: VW PA MR UL. Analyzed the data: VW. Contributed reagents/materials/analysis tools: VW IH UL. Wrote the paper: VW PA IH MR UL.

References

  1. 1. Reynolds HL, Packer A, Bever JD, Clay K (2003) Grassroots ecology: Plant-microbe-soil interactions as drivers of plant community structure and dynamics. Ecology 84: 2281–2291.
  2. 2. Kulmatiski A, Beard KH, Stevens JR, Cobbold SM (2008) Plant-soil feedbacks: a meta-analytical review. Ecology Letters 11: 980–992.
  3. 3. Bever JD, Dickie IA, Facelli E, Facelli JM, Klironomos J, et al. (2010) Rooting theories of plant community ecology in microbial interactions. Trends in Ecology and Evolution 25: 468–478.
  4. 4. Klironomos JN (2003) Variation in plant response to native and exotic arbuscular mycorrhizal fungi. Ecology 84: 2292–2301.
  5. 5. Petermann JS, Fergus AJF, Turnbull LA, Schmid B (2008) Janzen-Connell effects are widespread and strong enough to maintain diversity in grasslands. Ecology 89: 2399–2406.
  6. 6. van der Putten WH, Van Dijk C, Peters BAM (1993) Plant-specific soil-borne diseases contribute to succession in foredune vegetation. Nature 362: 53–56.
  7. 7. van der Heijden MGA, Boller T, Wiemken A, Sanders IR (1998) Different arbuscular mycorrhizal fungal species are potential determinants of plant community structure. Ecology 79: 2082–2091.
  8. 8. Packer A, Clay K (2000) Soil pathogens and spatial patterns of seedling mortality in a temperate tree. Nature 404: 278–281.
  9. 9. van der Heijden MGA (2004) Arbuscular mycorrhizal fungi as support systems for seedling establishment in grassland. Ecology Letters 7: 293–303.
  10. 10. Callaway RM, Thelen GC, Rodriguez A, Holben WE (2004) Soil biota and exotic plant invasion. Nature 427: 731–733.
  11. 11. Knevel IC, Lans T, Menting FBJ, Hertling UM, van der Putten WH (2004) Release from native root herbivores and biotic resistance by soil pathogens in a new habitat both affect the alien Ammophila arenaria in South Africa. Oecologia 141: 502–510.
  12. 12. van Grunsven RHA, Bos F, Ripley BS, Suehs CM, Veenendaal EM (2009) Release from soil pathogens plays an important role in the success of invasive Carpobrotus in the Mediterranean. South African Journal of Botany 75: 172–175.
  13. 13. van Grunsven RHA, van der Putten WH, Bezemer TM, Berendse F, Veenendaal EM (2010) Plant-soil interactions in the expansion and native range of a poleward shifting plant species. Global Change Biology 16: 380–385.
  14. 14. Parker MA, Malek W, Parker IM (2006) Growth of an invasive legume is symbiont limited in newly occupied habitats. Diversity and Distributions 12: 563–571.
  15. 15. Nuñez MA, Horton TR, Simberloff D (2009) Lack of belowground mutualisms hinders Pinaceae invasions. Ecology 90: 2352–2359.
  16. 16. Linhart YB (1991) Disease, parasitism and herbivory: Multidimensional challenges in plant evolution. Trends in Ecology and Evolution 6: 392–396.
  17. 17. Parker MA (1999) Mutualism in metapopulations of legumes and rhizobia. American Naturalist 153: S48–S60.
  18. 18. Ehrenfeld JG, Ravit B, Elgersma K (2005) Feedback in the plant-soil system. Annual Review of Environment and Resources 30: 75–115.
  19. 19. Macel M, Lawson CS, Mortimer SR, Šmilauerova M, Bischoff A, et al. (2007) Climate vs. soil factors in local adaptation of two common plant species. Ecology 88: 424–433.
  20. 20. Svenning MM, Junttila O, Solheim B (1991) Symbiotic growth of indigenous white clover (Trifolium repens) with local Rhizobium leguminosarum biovar trifolii. Physiologia Plantarum 83: 381–389.
  21. 21. Parker MA (1995) Plant fitness variation caused by different mutualist genotypes. Ecology 76: 1525–1535.
  22. 22. Ronsheim ML, Anderson SE (2001) Population-level specificity in the plant-mycorrhizae association alters intraspecific interactions among neighboring plants. Oecologia 128: 77–84.
  23. 23. Sylvia DM, Alagely AK, Kane ME, Philman NL (2003) Compatible host/mycorrhizal fungus combinations for micropropagated sea oats. I. Field sampling and greenhouse evaluations. Mycorrhiza 13: 177–183.
  24. 24. Burdon JJ, Gibson AH, Searle SD, Woods MJ, Brockwell J (1999) Variation in the effectiveness of symbiotic associations between native rhizobia and temperate Australian Acacia: within-species interactions. Journal of Applied Ecology 36: 398–408.
  25. 25. Hoeksema JD, Thompson JN (2007) Geographic structure in a widespread plant-mycorrhizal interaction: pines and false truffles. Journal of Evolutionary Biology 20: 1148–1163.
  26. 26. de la Peña E, Bonte D, Moens M (2009) Evidence of population differentiation in the dune grass Ammophila arenaria and its associated root-feeding nematodes. Plant and Soil 324: 307–316.
  27. 27. Thompson JN (2005) The geographic mosaic of coevolution. Chicago: University of Chicago Press. 443 p.
    • 28. Becker U, Colling G, Dostal P, Jakobsson A, Matthies D (2006) Local adaptation in the monocarpic perennial Carlina vulgaris at different spatial scales across Europe. Oecologia 150: 506–518.
    • 29. Bischoff A, Crémieux L, Smilauerova M, Lawson CS, Mortimer SR, et al. (2006) Detecting local adaptation in widespread grassland species - the importance of scale and local plant community. Journal of Ecology 94: 1130–1142.
    • 30. Gray AJ (2002) The evolutionary context: a species perspective. In: Perrow MR, Davy AJ, editors. Handbook of ecological restoration. Volume 1. Principles of restoration. Cambridge, UK: Cambridge University Press. pp. 66–82.
      • 31. Thrall PH, Slattery JF, Broadhurst LM, Bickford S (2007) Geographic patterns of symbiont abundance and adaptation in native Australian Acacia-rhizobia interactions. Journal of Ecology 95: 1110–1122.
      • 32. Kawecki TJ, Ebert D (2004) Conceptual issues in local adaptation. Ecology Letters 7: 1225–1241.
      • 33. Lavrenko EM (1970) Provinzial'noe razdelenie Prichernomorsko-Kazakhstanskoy podoblasti stepnoy oblasti Evrazii. [The division of the Black Sea-Kazakhstan subregion of the Eurasian steppe region]. Botanicheskiy Zhurnal 55: 609–625.
      • 34. Wagner V, von Wehrden H, Wesche K, Fedulin A, Sidorova T, et al. Similar performance in central and range-edge populations of a Eurasian steppe grass under different climate and soil pH regimes. Ecography. In press.
      • 35. Ponomarev AN (1961) Kleystogamiya u kovyley. [Cleistogamy in feather grasses]. Botanicheskiy Zhurnal 46: 1229–1236.
      • 36. Mukhin VA, Betekhtina AA (2006) Adaptive significance of endomycorrhizas for herbaceous plants. Russian Journal of Ecology 37: 1–6.
      • 37. Shi ZY, Chen ZC, Zhang LY, Feng G, Christie P, et al. (2007) Diversity and zonal distribution of arbuscular mycorrhizal fungi on the northern slopes of the Tianshan Mountains. Science in China Series D: Earth Sciences 50: Supp. I135–141.
      • 38. Hijmans RJ, Cameron SE, Parra JL, Jones PG, Jarvis A (2005) Very high resolution interpolated climate surfaces for global land areas. International Journal of Climatology 25: 1965–1978.
      • 39. Brinkman EP, van der Putten WH, Bakker E-J, Verhoeven KJF (2010) Plant-soil feedback: experimental approaches, statistical analyses and ecological interpretations. Journal of Ecology 98: 1063–1073.
      • 40. Bever JD (1994) Feedback between plants and their soil communities in an old field community. Ecology 75: 1965–1977.
      • 41. Katan J (2002) Interactions of roots with soil-borne pathogens. In: Waisel Y, Eshel A, Kafkafi U, editors. Plant roots: the hidden half. New York: Marcel Dekker. pp. 949–959.
        • 42. Trouvelot A, Kough JL, Gianinazzi-Pearson V (1986) Mesure du taux de mycorhization VA d'un système radiculaire. Recherche de méthodes d'estimation ayant une signification fonctionnelle. In: Gianinazzi-Pearson V, Gianinazzi S, editors. Physiological and genetical aspects of mycorrhizae: Proceedings of the 1st European Symposium on Mycorrhizae. Versailles: Service des Publications INRA. pp. 217–221.
          • 43. Armas C, Ordiales R, Pugnaire FI (2004) Measuring plant interactions: A new comparative index. Ecology 85: 2682–2686.
          • 44. Carvalho LM, Antunes PM, Martins-Loução MA, Klironomos JN (2010) Disturbance influences the outcome of plant-soil biota interactions in the invasive Acacia longifolia and in native species. Oikos 119: 1172–1180.
          • 45. Canty A, Ripley B (2010) boot: Bootstrap R (S-Plus) Functions. R package version 1.2-43.
            • 46. R Development Core Team (2010) R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing. Version. 2.12.0.
              • 47. Pinheiro J, Bates D, DebRoy S, Sarkar D, R Development Core Team (2010) nlme: Linear and nonlinear mixed effects models. R package version 3.1-97.
                • 48. Leimu R, Fischer M (2008) A meta-analysis of local adaptation in plants. PLoS ONE 3: e4010.
                • 49. Zuur AF, Ieno EN, Walker NJ, Saveliev AA, Smith GM (2009) Mixed effects models and extensions in ecology with R. New York: Springer. 574 p.
                  • 50. Pinheiro JC, Bates DM (2000) Mixed-effects models in S and S-Plus. New York: Springer. 528 p.
                    • 51. Bell T, Freckleton RP, Lewis OT (2006) Plant pathogens drive density-dependent seedling mortality in a tropical tree. Ecology Letters 9: 569–574.
                    • 52. Reinhart KO, Royo AA, van der Putten WH, Clay K (2005) Soil feedback and pathogen activity in Prunus serotina throughout its native range. Journal of Ecology 93: 890–898.
                    • 53. Linhart YB, Grant MC (1996) Evolutionary significance of local genetic differentiation in plants. Annual Review of Ecology and Systematics 27: 237–277.
                    • 54. García-Ramos G, Kirkpatrick M (1997) Genetic models of adaptation and gene flow in peripheral populations. Evolution 51: 21–28.
                    • 55. Fierer N, Jackson RB (2006) The diversity and biogeography of soil bacterial communities. Proceedings of the National Academy of Sciences, USA 103: 626–631.
                    • 56. Reinhart KO, Clay K (2009) Spatial variation in soil-borne disease dynamics of a temperate tree, Prunus serotina. Ecology 90: 2984–2993.
                    • 57. Roesch LFW, Fulthorpe RR, Riva A, Casella G, Hadwin AKM, et al. (2007) Pyrosequencing enumerates and contrasts soil microbial diversity. The ISME Journal 1: 283–290.
                    • 58. Allison SD, Martiny JBH (2008) Resistance, resilience, and redundancy in microbial communities. Proceedings of the National Academy of Sciences, USA 105: 11512–11519.
                    • 59. Sultan SE, Spencer HG (2002) Metapopulation structure favors plasticity over local adaptation. The American Naturalist 160: 271–283.
                    • 60. Fox LR (1981) Defense and dynamics in plant-herbivore systems. American Zoologist 21: 853–864.
                    • 61. Dyer AR, Rice KJ (1997) Intraspecific and diffuse competition: The response of Nassella pulchra in a California grassland. Ecological Applications 7: 484–492.
                    • 62. Stinchcombe JR, Rausher MD (2001) Diffuse selection on resistance to deer herbivory in the ivyleaf morning glory, Ipomoea hederacea. American Naturalist 158: 376–388.
                    • 63. Antunes PM, Koch AM, Morton JB, Rillig MC, Klironomos JN (2011) Evidence for functional divergence in arbuscular mycorrhizal fungi from contrasting climatic origins. New Phytologist 189: 507–514.
                    • 64. McCarthy-Neumann S, Kobe RK (2010) Conspecific and heterospecific plant-soil feedbacks influence survivorship and growth of temperate tree seedlings. Journal of Ecology 98: 408–418.
                    • 65. McCarthy-Neumann S, Kobe RK (2010) Conspecific plant-soil feedbacks reduce survivorship and growth of tropical tree seedlings. Journal of Ecology 98: 396–407.