Plant-Soil Feedbacks of Plantago lanceolata in the Field Depend on Plant Origin and Herbivory

Soil biota involved in plant-soil feedbacks (PSFs) have an impact on the growth of plant individuals. So far, studies investigating the role of soil-biota mediated PSFs in plant performance were mostly performed in greenhouses and focused predominantly on species differences, whereas the contribution of PSFs to plant performance under field conditions and intraspecific variation in PSFs among plant populations remain poorly investigated. Here, we performed a PSF pot experiment under field conditions to study intraspecific variation in plant responses to soil biota. We used seeds from multiple seed families of Plantago lanceolata L. together with Plantago-conditioned soils from contrasting habitats (three non-fertilized pastures vs. three fertilized mown pastures) to test whether plants show a positive or negative response to their parental soil biota. We furthermore tested whether these PSFs depend on microclimate and insect herbivory. To this end, we reciprocally transplanted plants and their soil biota between the two habitat types and excluded aboveground herbivores from half of the plants, respectively. When grown without herbivores, plants from both habitat types showed similar and neutral PSFs independently of the transplant site. In contrast, in the presence of herbivores, PSFs for plants from non-fertilized pastures were negative in both habitats (i.e. plants performed better when they grew with foreign soil biota), whereas PSFs for plants from fertilized mown-pastures remained neutral. Our results suggest that soil biota alone might only play a minor role for performance of P. lanceolata and that the outcome of soil-biota mediated PSFs is modulated by effects of herbivores in different habitats.


INTRODUCTION
For a better understanding of the interplay between ecosystem dynamics and evolution, we not only need to understand how species interact, but also how populations evolve in response to environmental drivers (Van Nuland et al., 2016;Ware et al., 2019a). The concept of ecoevolutionary feedbacks looks at the interaction and adaptation between plants and their biotic and abiotic environment (Bailey et al., 2009). So-called plant-soil feedbacks (PSFs; Bever et al., 1997) are an ideal system to investigate such eco-evolutionary dynamics, as soil microbiota have short life-cycles and can quickly adapt while plant fitness is very closely linked to its interaction with soil biota (TerHorst and Zee, 2016). PSFs influence plant biomass production and performance and thus competitiveness and coexistence of species in plant communities (Klironomos, 2002;Van der Putten et al., 2013;Heinze et al., 2015a). On top of that, other environmental players, such as herbivory or habitat, can alter PSF responses (Heinze and Joshi, 2018;Pfennigwerth et al., 2018;Ware et al., 2019b). Simultaneously, PSFs may depend on the plant genotype, and different plant populations may evolve diverging PSF responses (Bukowski and Petermann, 2014;Wagg et al., 2015;Luo et al., 2016;Allen et al., 2018;Hawkins and Crawford, 2018). Given the ecological impact of PSFs on plant communities and the substantial intraspecific variation in PSFs, there is a need for an integrated approach that considers evolutionary-and ecological-scale processes simultaneously in order to improve our understanding of PSFs.
To measure PSF effects, experiments typically compare plant growth on "home" soils (i.e., self-cultivated) to plant growth on "away" soils (i.e., non-self-cultivated; Bever et al., 1997;Brinkman et al., 2010;Van der Putten et al., 2013). PSFs are positive when plant growth is greater on home than on away soils and negative when plant growth is greater on away than on home soils (Bever, 1994). Soil biota play an important role in PSFs (Ehrenfeld et al., 2005;Brinkman et al., 2010). Negative PSFs mostly result from the harmful effects of soil pathogens and herbivores such as fungal pathogens, insect larvae or nematodes, whereas positive PSFs follow from symbioses with mycorrhizal fungi and other microorganisms such as decomposers involved in nutrient cycling (Van der Putten et al., 2016). Across the literature PSFs are predominantly negative (e.g.,  suggesting that pathogenic or parasitic soil biota accumulate over time in "home" soils. Subsequently, these negative PSFs would lead to a decrease in competitive ability (Kardol et al., 2007) and would cause local rarity (Klironomos, 2002) or even extinction in local communities when competing with species that gain positive PSFs from soil biota. Similar to species in a community, populations within a species can respond differently to their soil biota as a result of evolutionary processes (Bukowski and Petermann, 2014;Wagg et al., 2015).
Agriculturally used grasslands provide an ideal system to study intraspecific variation in PSFs. Agricultural land use affects environmental conditions belowground by influencing abiotic soil properties (Alt et al., 2011;Birkhofer et al., 2012), which in turn influences the composition of soil biota (Herold et al., 2014). Plant roots, representing a large part of a plant's biomass (Yang et al., 2010), are exposed to these soil biota and interact with them. Research on such interactions has shown that land-use mediated changes in soil biota affect plant performance and coexistence (Heinze et al., 2015a,b). Plants also alter abiotic soil properties and soil biota through litter production, nutrient-uptake and exudation processes and hence create soil-legacy effects (Ehrenfeld et al., 2005;Van der Putten et al., 2013). Such initial plant-induced changes of soil properties can then influence the establishment and growth of subsequently establishing plants (Bever, 1994;Bever et al., 1997;Kardol et al., 2007). Interestingly, long-term land-use management may elicit evolutionary responses in soil biota, plants and their interaction, which has been well investigated for plants (Turesson, 1922;Warwick and Briggs, 1979;Silvertown et al., 2006). Moreover, additional factors, such as herbivores and habitat, can be influenced by land use as well (Gardiner and Hassall, 2009;Gossner et al., 2014;Simons et al., 2014) and may have evolutionary implications for PSFs. Therefore, including multiple environmental variables may help to understand the eco-evolutionary dynamics of PSFs under realistic environmental conditions. The majority of PSF experiments have been performed in greenhouse conditions that fail to place PSFs in the context of environmental conditions that are likely to affect the role of plant growth in communities (Heinze et al., 2016). Relative to PSFs measured in greenhouse experiments, measurements of PSFs under field conditions have been found to vary mainly due to the diverse abiotic and biotic interactions that plants and soils receive under natural field conditions (Casper et al., 2008;Heinze et al., 2016). Hence, there is a need to perform PSF experiments under field conditions ( Van der Putten et al., 2016), where biotic interactions, such as herbivory, are present. Under such natural conditions recent research found, for example, that aboveground insect herbivory modulated the outcome of PSFs (Heinze and Joshi, 2018) and that the PSF response increased with the intensity of herbivory (Heinze et al., 2019).
We performed a reciprocal pot transplant experiment in the field and investigated intraspecific variation in biotic PSFs (i.e., using inoculated standardized soils) in the widespread plant species, Plantago lanceolata L. We included six P. lanceolata populations, three each from two contrasting habitats with contrasting land-use intensity (non-fertilized pastures vs. fertilized mown pastures). In this experiment we also tested the effect of abiotic habitat factors by reciprocally transplanting plants into the two habitat types, and the effect of herbivore presence versus absence, using a herbivory-exclosure treatment, on the outcome of intraspecific biotic PSFs.
Specifically, we asked: 1. Do P. lanceolata populations differ in biotic PSFs between the two contrasting plant origins? 2. Are biotic PSFs of P. lanceolata affected by abovegroundinsect herbivory and habitat?
3. Does herbivory damage to plants vary with respect to plant and soil origin? 4. What is the relative impact of soil biota (i.e., biotic PSFs) compared to aboveground-insect herbivory and habitat on the overall biomass variation?

Study Site and Species
To test whether populations of P. lanceolata from two habitat types of contrasting land use show PSFs to local soil biota and how PSFs are affected by herbivory and abiotic habitat factors, we performed an experiment under field conditions in the UNESCO Biosphere Reserve Schorfheide-Chorin embedded in a glacial landscape in the lowlands of north-eastern Germany ( Figure S1A). The experiment was conducted within the framework of the Biodiversity Exploratories, a large-scale and long-term project investigating the effects of land use on biodiversity in Germany (Fischer et al., 2010). The two contrasting habitat types-non-fertilized pastures and fertilized mown pastures-were chosen, because the factors fertilization and mowing were shown to affect the composition of soil biota (Herold et al., 2014) as well as aboveground phenotypic trait differentiation of plants (Völler et al., 2013) on the sites of the Biodiversity Exploratories. For each habitat type, three sites were chosen resulting in three pairs of non-fertilized pasture vs. fertilized mown pasture (mean geographic distance between pairs: 1.3 km, Figure S1B).
The perennial herbaceous plant species, Plantago lanceolata L., was chosen as a model organism because it is widespread in both habitat types (Joshi et al., 2001). P. lanceolata is an outbreeding (Ross, 1973) and mostly wind-pollinated (Clifford, 1962) grassland species. In previous studies P. lanceolata has been shown to exhibit strong trait differentiation and adaptation to specific land use practices (e.g. Wolff and Delden, 1987;Van Tienderen and van der Toorn, 1991;Joshi et al., 2001). For instance, in a reciprocal transplant experiment using three populations from sites with contrasting land use, fitness was always higher in the home site and fitness-related traits such as flowering time varied (Van Tienderen and van der Toorn, 1991).

Seed and Soil Collection and Preparation
Seeds of P. lanceolata were collected in summer 2017 in all six sites from five randomly chosen individual plants, which were growing at least 1 m apart from each other. Seeds collected from a single plant individual were considered a seed family. To avoid microbial contaminations, seeds were surface-sterilized using 7% sodium hypochlorite (NaClO) solution (Heinze et al., 2017). Afterwards seeds were washed with autoclaved water (20 min, 121 • C), germinated in sterilized (5-times in 24 h; 20 min, 121 • C) sand (grain size: 2 mm; Brun and Böhm, Potsdam, Germany) in petri dishes (diameter: 9 cm) and placed in sterile plastic chambers (32 cm × 50 cm × 14 cm; Meyer; Germany) in a greenhouse at the University of Potsdam.
To investigate the effect of local soil biota on plants from different origins we used species-specific field-conditioned soils of P. lanceolata in accordance with the "natural experiment" approach (Kulmatiski and Kardol, 2008). Previous results from the Biodiversity Exploratories have shown that the land-use treatments affected abiotic soil properties (Alt et al., 2011;Birkhofer et al., 2012) and the composition of soil biota (Herold et al., 2014) on the chosen sites. In spring 2018 on the respective sites, similar to Heinze and Joshi (2018), soil material was collected from below P. lanceolata individuals located in the center of larger patches (diameter >20 cm) of P. lanceolata to make sure that the soil was conditioned by our target species and over longer time periods in both habitats. Following Brandt et al. (2014) we collected rhizosphere soil and soil directly adjacent to the rhizosphere from 20 individuals per site. After sampling, the soil was stored for 3 weeks at 4 • C in the dark until use in the experiment.

PSF Experiment
Seeds and soils of plants from six study sites were used in an experiment to investigate intraspecific variation in PSFs and how effects of land use modulate PSF responses. Within the three habitat pairs (i.e., a non-fertilized pasture and a fertilized mown pasture), rhizosphere soils from the same habitat served as "home" soil, whereas soils from the contrasting habitat served as "away" soil. As we were interested in site-specific PSFs and not within-site variation in PSFs we mixed the 20 individual soil samples to one bulk sample per site. Although this mixing procedure has been criticized for its potential to increase the likelihood of falsely detecting PSFs by decreasing variance in plant responses among individual soil samples (Reinhart and Rinella, 2016) this procedure was appropriate for our specific research question. Furthermore, several studies reported that soil mixtures produce similar PSFs compared to independent soil samples and suggest that soil handling methods should be dependent on specific research questions and feasibility (e.g., Cahill et al., 2016;Kulmatiski, 2016;Gundale et al., 2019;Teste et al., 2019). To reduce potential differences in soil nutrient availability among soils and to focus on effects of soil biota we used the collected rhizosphere soils as inoculum (10%) into an autoclaved soil:sand mixture (Brinkman et al., 2010). The soil:sand mixture consisted of a 1:1 mixture of sieved (mesh size: 7 mm) field soil collected from a meadow at the field site of the University of Potsdam (N52 • 24 ′ 29.76 ′′ , E13 • 1 ′ 13.74 ′′ , Brandenburg, Germany) and purchased sand (grain size: 2 mm; Brun and Böhm; Potsdam, Germany).
Pots (Deepots D25: volume 0.41 L; height 25 cm; diameter 5 cm; Stuewe and Sons; USA) were prepared with an autoclaved fleece strip (6 cm × 25 cm) covering 10 cm of the pots' inside and extending 15 cm below the pot to enable continuous watering from below. The pots were subsequently filled with the inoculated soils. To limit cross-contamination of soil biota between the pots, each pot was placed in a separate plastic cup (volume 0.3 L; height 15.2 cm; diameter 5.9 cm) and received an additional layer (1 cm) of sterilized sand on top.
In May 2018, 1-week old seedlings were transplanted into the prepared pots. The planting scheme followed a reciprocal transplant design: per habitat pair, seedlings of both habitat types were planted in their home and away soils. Each plant origin × soil origin × herbivory treatment (see below and Figure 1) combination was replicated 10 times in each site, by including two offspring each from the five randomly chosen seed families per plant origin.
After planting, seedlings were transported to a protected outdoor location on the field site of the University of Potsdam and were allowed to acclimatize for 2 weeks.

Herbivory Treatment
To investigate the impact of aboveground insect herbivory on PSFs of the different plant origins we performed a herbivoryexclusion treatment in accordance with Heinze and Joshi (2018) and Heinze et al. (2019). This herbivore-exclusion treatment was applied on all six experimental sites. In each site we established two plots (120 cm × 160 cm) that were spaced 80 cm apart. Each plot was equipped with a cage (length 160 cm × width 120 cm × height 100 cm) that was either completely covered with fly mesh (mesh size: 1.3 mm; Meyer; Germany) or only shaded by fly mesh on top. The fully covered cages excluded herbivorous insects (MacDonald and Kotanen, 2010), whereas the shaded cages allowed aboveground herbivorous insects to reach the experimental plants while providing the same levels of shade and precipitation as the cage treatment (Heinze and Joshi, 2018). In the fully covered plots, we removed the aboveground vegetation with the top 5 cm of the soil to exclude non-developed aboveground-herbivorous insects whose eggs might be attached to plants or buried in the soil. Furthermore, in the fully covered plots the fly screen was buried into the soil and one side was prepared as a door to water the plants. The fully covered plots that excluded aboveground herbivorous insects (>1.3 mm) are referred to as "without herbivory, " whereas the shaded plots are referred to as "with herbivory" treatment throughout the manuscript.

Experimental Set-Up and Measurements
After 2 weeks of acclimatization in early-June, the prepared pots were brought to the six prepared sites. At each site, plants from both habitat types were grown in their home and away soil in the presence vs. absence of aboveground herbivory (Figure 1) resulting in a total of 480 experimental plants (3 habitat pairs × 2 sites × 2 plant origins × 2 soil origins × 2 herbivory treatments × 5 seed families × 2 replicates).
As we wanted to exclude direct competition between experimental plants and neighboring plants and as we were interested in the effects of aboveground invertebrates (excluding slugs) the experimental plants were placed in boxes (78 cm × 50 cm × 30 cm). The two replicates per treatment combination were divided between the two boxes at each plot and placed at random positions within those boxes. In that way each box contained one replicate for every treatment combination, either in the with herbivory treatment or the without herbivory treatment. At all sites the experimental set-up was fenced (3 m × 3 m) to prevent damages by cattle.
During the experiment, plants were watered every third week. After 12 weeks in the field, damage by aboveground chewing insects without further discrimination of feeding guilds was visually estimated for whole plants. Afterwards shoot biomass was harvested and roots were washed. Root and shoot biomass were dried to constant weight (minimum 48 h at 80 • C) and kept in the drying oven until it was weighed.

Statistical Analysis
To compare plant performance on their home vs. away soil, PSFs were calculated for each seed family individually. PSF values were calculated using total plant biomass (shoot + root). To obtain a quantitative measure of PSFs, where positive values show positive PSF and negative values negative PSF, respectively, the log ratio of home vs. away biomass was calculated (Brinkman et al., 2010). Home biomass was the biomass of a certain seed family on its own soil and away biomass was the biomass of the same seed family on the soil from the contrasting site.
Due to the random distribution of the pots in the boxes a block-wise analysis was not appropriate. Therefore, in accordance with Heinze et al. (2016) we calculated PSFs using all possible combinations of seed family-specific biomass production of one seed family in home soil in comparison to biomass production of this seed family in away soils (Equation 1).
For each replicate of the six plant origins, the home/away ratio in multiple comparisons with all 10 replicates in the same away group was calculated.
The PSF values were used in a linear mixed-effects model, testing for effects of the experimental treatments. Since the factor soil origin was already incorporated in our calculation of the PSFs, the factors plant origin (non-fertilized pasture vs. fertilized mown pasture), habitat (idem) and herbivory (with vs. without) and all possible interactions were included as fixed factors whereas habitat pair was included as a random factor.
To test for the effect of the herbivory treatment, a linear mixed-effects model was performed using the estimated percentage of feeding damage on aboveground biomass as response variable against all the factors included in our design, i.e., plant origin, habitat, soil origin and herbivory, and all possible two-way interactions. We included habitat pair and seed family as random factors in this model.
To determine the amount of variance in total plant biomass explained by the experimental treatments, a variance component analysis was performed using a linear mixed-effects model in which the same explanatory factors were included as in the model for feeding damage, but now all factors were treated as random factors.
The residuals of all models were checked for normality of distribution and homogeneity of variance. Data analysis was performed using R version 3.4.2 (R Core Team, 2017) and the lme4 package (Bates et al., 2014).

Biotic PSFs Between Contrasting Plant Origins in the Absence or Presence of Aboveground Insect Herbivory
In our PSF model, plant origin as well as its interaction with the herbivory treatment had a significant influence on PSFs, while other factors remained non-significant (Table 1).
To explore the environmental interactions in more detail, we performed two separate models, one using data for plants excluded from herbivory and the other on data of plants exposed to natural herbivory.
When grown without herbivores, biotic PSFs of plants from both habitat types were neutral, independently of transplant site (Table 2; Figure 2A). In the presence of aboveground insect herbivores, however, the plants from non-fertilized pastures showed negative biotic PSFs in both habitat types, whereas the The model tested for all the factors included in our study design, except for soil origin which is already incorporated in the calculation for the PSF quantification. Those are plant origin, habitat (non-fertilized pastures vs. fertilized mown pastures) and herbivory (with vs. without). PSF was calculated using log(biomass on home soil / biomass on away soil). Significant P-values (P < 0.05) are indicated in bold typeface.  Figure 2B).
Although some feeding damage was observed in the withoutherbivory plots, the percentage of aboveground feeding damage was significantly higher in the with-herbivory plots [F (340, 7) = 13.79; P = 0.0002; Table S1; Figure 3], but only in the fertilized mown pastures (Figure 3). Damage in the with-and without-herbivory plots was similar in the non-fertilized pastures (Figure 3). These different herbivory effects had, however, no significant effect on plant biomass in this experiment (see below and Table 3).

Relative Impact of Treatments on Plant Performance
Overall, the factors -and their interactions-included in the study design together explained less than 10% of the variation in biomass of the experimental plants ( Table 3). Out of these factors the origin of soil biota accounted for most of the variance in biomass production (4.9%). The interaction between plant origin and habitat accounted for 2.0% of the variance whereas all other factors and their interactions explained less than 1% (Table 3).

PSFs Without Aboveground Herbivores
We found a main effect of plant origin on PSF as well as an interaction between plant origin and herbivory. When focusing on the results when natural aboveground insect herbivory was excluded, i.e., the sole effect of soil biota, no differences in intraspecific PSFs for populations of P. lanceolata were observed. Although only scarcely treated in the literature, there are examples of studies that found intraspecific variation in PSFs (Wagg et al., 2015;Luo et al., 2016;Eck et al., 2019). For example, Smith et al. (2012) found seedling survival of Populus angustifolia to be 2.5 times higher in soil biota from the rhizosphere of their mother tree, but neutral PSFs as well have been found (Bukowski and Petermann, 2014;Allen et al., 2018). However, our study did not test for individual effects of maternal soil, but for landuse effects (i.e., fertilization etc.) on soil biota and how different intensities of land use might influence PSFs in a population. Fertilization is known to influence the composition of soil biota (Herold et al., 2014) and thus the relative abundance of fungi and bacteria in soils (van der Heijden et al., 2008). Although we used soils from non-fertilized vs. fertilized sites in our experiment, we did not observe differences in PSFs when aboveground insect herbivores were excluded. Either the contrasting land-use regimes did not create as much differentiation in local soil biota as expected or there were indeed differences in local soil biota, but the plants respond neutrally to them. We do not know the effects of soil biota per se, since a soil inoculum was used without an unconditioned or sterilized control soil to test against the effect of soil biota.
Although the inoculation technique is optimal to investigate effects of soil biota (e.g., biotic PSFs; Brinkman et al., 2010) there are also several disadvantages. For example, inoculation into sterilized background soil might dilute the soil biota to an extent where effects, actually present in nature, are not detected anymore (Brinkman et al., 2010). Furthermore, since populations of bacteria grow much faster than those of fungi, the latter would decrease in their ratio when adding the inoculum to standardized soil. This effect could further minimize the effect of fungi in our experimental set-up. However, in contrast to other studies working on intraspecific PSFs using experimental conditioning, we used field-conditioned soil according to the "natural experiment" approach (Kulmatiski and Kardol, 2008) that is suggested to better reflect the soil communities in the field (Brinkman et al., 2010). Furthermore, with regard to our questions, the approach appeared favorable, as land use was shown to create substantial and long-lasting changes to the microbial soil community (Jangid et al., 2011;Dupouey et al., 2013) and thus soil legacies that cannot be changed via short-term conditioning with the targeted plants under artificial experimental conditions (Kulmatiski and Beard, 2011).
In general, experimental growth conditions -controlled conditions in greenhouses or climate chambers vs. field conditions-were shown to affect plant physiology (Ferreira et al., 1995) and growth (Mishra et al., 2012) as well as the development and composition of soil biota (Ge et al., 2011). In the field, harsh abiotic conditions (e.g., wind; Whitehead and Luti, 1962;Nicoll and Ray, 1996), fluctuating light regimes (Hunt and Burnett, 1973) as well as temperature (thermoperiodicity; Went, 1974) may increase biomass allocation to roots which are suggested to modulate interactions between plants and soil biota (Bardgett et al., 1998). Hence, our field experiment testing for intraspecific PSF effects fulfilled the need to consider PSFs under more natural conditions (Van der Putten et al., 2016).
The absence of a response to local soil biota under natural conditions -but without aboveground insect herbivory-does not imply that other environmental factors do not affect the performance of plants (see below). Further studies on PSFs might incorporate other important factors affecting plant growth and soil biota such as abiotic soil conditions that are influenced by land use (Alt et al., 2011;Birkhofer et al., 2012). However, we included one important environmental factor-aboveground insect herbivory-in our experiment that was shown to affect the outcome of PSFs under field conditions (Heinze and Joshi, 2018), depending on the intensity of herbivory (Heinze et al., 2019).

PSFs With Aboveground Insect Herbivory
When experimental plants were exposed to natural aboveground insect herbivory, PSFs became negative for plants from the nonfertilized pastures but remained neutral for plants from the fertilized mown pastures. Such neutral to negative PSFs were often found in grasslands (Ehrenfeld et al., 2005;Van der Putten et al., 2016). In our study on intraspecific variation in PSF, the performance of P. lanceolata from non-fertilized pastures increased when growing on away soils, but only when insect herbivores were present. This result indicates that P. lanceolata harbors intraspecific variation in PSF between populations with contrasting land use. This is a remarkable result given that the six study sites were located in The percentage shows the amount of total variance explained.
a geographically small area and that gene flow is thus likely to be substantial among the study populations of this wind-pollinated outcrossing species. We thus propose that natural selection for differential PSF effects in the respective sites is stronger than the counteracting effects of gene flow. Simultaneously, it may be that gene flow allows the plant populations to evolve resistance against rapidly evolving antagonistic soil organisms. Both processes suggest a strong ecological significance for the observed intraspecific differences in PSF and implies that eco-evolutionary feedbacks play a role in shaping these grassland ecosystems. The mechanism behind the negative PSF effects remains unclear. Since aboveground insect herbivory obviously has a direct effect on shoot biomass, it is important to assess whether our results could be driven purely by shoot biomass loss. When testing for an effect of soil origin under the presence of insect herbivory on shoot and root biomass separately, we found that root biomass was more strongly affected than shoot biomass ( Figure S2). Moreover, aboveground herbivory itself significantly affected root biomass but not shoot biomass (Table S2). Since root biomass can only be indirectly influenced by aboveground herbivory, the observed negative PSFs on total biomass under the influence of herbivory is thus-at least in part-indirectly caused by aboveground herbivory and seems to involve changes in biomass allocation by the plant as the integrative result of both soil biota and insect herbivory.
Similarly, Mursinoff and Tack (2017) found that the response of P. lanceolata to local vs. foreign soil differed under the presence of a specialist leaf pathogen. Also Zhu et al. (2018) found similar effects in these plant-soil-herbivore interactions in P. lanceolata. Interestingly, some studies revealed that herbivory can change plant-soil-fungi interactions (Bardgett and Wardle, 2003;Bennett et al., 2009;Kostenko et al., 2012), while other studies demonstrated effects of plant-soil fungi interactions on herbivory (Fontana et al., 2009). From our results we cannot infer the exact mechanisms acting here, i.e., how the three parties are influencing each other. Further experiments are needed to elucidate the causal directions between herbivory and soil biota through P. lanceolata in our system and to understand the evolutionary drivers shaping this intraspecific variation in PSF. We speculate that the negative relationship between insect abundance and land-use intensity (Simons et al., 2017) as well as the negative relationship between symbiotic soil fungi and land-use intensity (de Vries et al., 2006) may act as selection pressures causing the observed intraspecific variation in PSF.
An important result is that insect herbivory was needed to elicit differential responses to soil origin between plants from different habitats. Although strictly speaking PSF solely involves plants and their soil biota and is ideally measured under controlled experimental conditions such as a greenhouse, it is clear that such approaches take plants and their associated soil biota out of their natural ecological context. By conducting PSF experiments under more natural conditions, we may be better able to assess the significance of various ecological drivers and their interactions, although the complexity of such experiments increases as well.

Relative Effect of Treatments on Plant Performance
With 4.9%, soil differences explained a small, but nevertheless the biggest part of total variance in our variance component analysis. Interestingly, differences in habitat only accounted for 0.7% of explained intraspecific variation in P. lanceolata. Assuming that land use is the primary driver of the differences in soil biota, this would imply that results of reciprocal transplant experiments may be driven substantially by soil biota rather than abiotic habitat effects. This is interesting given that data on the importance of climate vs. soil in plant local adaptation is scarce (Macel et al., 2007). However, it should be noted that the six study sites were located in a geographically small area and variation in microclimate among the habitats is therefore presumably small and solely influenced by land-use effects on biotic factors such as vegetation height. The interaction between habitat and plant origin accounted for 2% of the explained variation. Although such an interaction indicates variable responses of plants from different origins when transplanted to different habitats and may suggest local (mal-)adaptation (Kawecki and Ebert, 2004), this factor was not significant in our PSF model ( Table 1) and does not reflect a consistent pattern among habitat pairs.

CONCLUSION
Our experiment investigated the effects of soil biota on intraspecific variation among plant populations under natural conditions. We did not consider plants and their soils in isolation, but included their interaction with other biotic and abiotic factors. Our study therefore contributes to the growing literature considering eco-evolutionary feedbacks to explain how complex interactions between multiple parties influence ecological and evolutionary dynamics (Van Nuland et al., 2016;Ware et al., 2019a).
The most interesting finding of our study is that intraspecific differences between populations of P. lanceolata from sites with contrasting land-use intensity only became apparent when plants were exposed to insect herbivory. Thus, our study offers support to the theory that interactions between plants and soilmicrobiota can be mediated through aboveground-herbivores and the responses they induce in plants. Additionally, our study is among the first to show such complex interactions under field conditions.

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

AUTHOR CONTRIBUTIONS
AK, JJ, JS, and JH designed the study. LK, AK, JS, and JH collected seeds and soils, performed the experiment, and measured the data. LK, JS, and JH analyzed the data and wrote the first draft of the manuscript. LK, AK, JJ, OB, JS, and JH contributed to later versions of the manuscript.

ACKNOWLEDGMENTS
We would like to thank Nils Landmann from the ABP Landwirtschaft Neugrimnitz and colleagues as well as Torsten Meene, Markus Rubenbauer and Uta Schuhmacher for help in the field and organization. We also thank the manager of the Exploratory Schorfheide-Chorin, Miriam Teuscher, and all former managers for their work in maintaining the plot and project infrastructure, Christiane Fischer for giving support through the central office, Andreas Ostrowski for managing the central data base, and Markus Fischer, Eduard Linsenmair, Dominik Hessenmöller, Daniel Prati, Ingo Schöning, François Buscot, Ernst-Detlef Schulze, Wolfgang W. Weisser and the late Elisabeth Kalko for their role in setting up the Biodiversity Exploratories project. Field work permits were issued by the responsible state environmental offices of Brandenburg.