Predicting Plant-Soil Feedback in the Field: Meta-Analysis Reveals That Competition and Environmental Stress Differentially Influence PSF

Past research on plant-soil feedbacks (PSF), largely undertaken in highly controlled greenhouse conditions, has established that plant species differentially alter abiotic and biotic soil conditions that in turn affect growth of other conspecific and heterospecific individuals in that soil. Yet, whether feedbacks under controlled greenhouse conditions reflect feedbacks in natural environments where plants are exposed to a range of abiotic and biotic pressures is still unresolved. To address how environmental context affects PSF, we conducted a meta-analysis of previously published studies that examined plant growth responses to multiple forms of competition, stress, and disturbance across various PSF methodology. We asked the following questions: (1) Can competition, stress, and disturbance alter the direction and/or strength of PSF? (2) Do particular types of competition, stress, or disturbance affect the direction and/or strength of PSF more than others? and (3) Do methods of conducting PSF research (i.e., greenhouse vs. field experiments and whether the source of soil inoculum conditioning is from the field vs. greenhouse) affect plant growth responses to PSF or competition, stress, and disturbance, or their interactions? We discovered four patterns that may be predictive of what future PSF studies conducted under more realistic conditions might reveal. First, relatively little is known about how PSF responds to environmental stress and disturbance compared to plant-plant competition. Second, specific types of competition enhanced negative effects of soil microbes on plant growth, and specific environmental stressors enhanced positive effects of soil microbes on plant growth. Third, whether PSF experiments are conducted in the field or greenhouse can change plant growth responses. And, fourth, how the soil conditioning phase is conducted can change plant growth responses. With more detail than previously shown, these results confirm that environmental context writ large can change plant growth responses in PSF experiments. These data should aid theory and predictions for conservation and restoration applications by showing the relative importance of competition, stress, and disturbance in PSF studies over time. Lastly, these data demonstrate how variation in experimental methods can alter interpretation and conclusions of PSF studies.

Past research on plant-soil feedbacks (PSF), largely undertaken in highly controlled greenhouse conditions, has established that plant species differentially alter abiotic and biotic soil conditions that in turn affect growth of other conspecific and heterospecific individuals in that soil. Yet, whether feedbacks under controlled greenhouse conditions reflect feedbacks in natural environments where plants are exposed to a range of abiotic and biotic pressures is still unresolved. To address how environmental context affects PSF, we conducted a meta-analysis of previously published studies that examined plant growth responses to multiple forms of competition, stress, and disturbance across various PSF methodology. We asked the following questions: (1) Can competition, stress, and disturbance alter the direction and/or strength of PSF? (2) Do particular types of competition, stress, or disturbance affect the direction and/or strength of PSF more than others? and (3) Do methods of conducting PSF research (i.e., greenhouse vs. field experiments and whether the source of soil inoculum conditioning is from the field vs. greenhouse) affect plant growth responses to PSF or competition, stress, and disturbance, or their interactions? We discovered four patterns that may be predictive of what future PSF studies conducted under more realistic conditions might reveal. First, relatively little is known about how PSF responds to environmental stress and disturbance compared to plant-plant competition. Second, specific types of competition enhanced negative effects of soil microbes on plant growth, and specific environmental stressors enhanced positive effects of soil microbes on plant growth. Third, whether PSF experiments are conducted in the field or greenhouse can change plant growth responses. And, fourth, how the soil conditioning phase is conducted can change plant growth responses. With more detail than previously shown, these results confirm that environmental context writ large can change plant growth responses in PSF

INTRODUCTION
The past 25 years of research on plant-soil feedback (PSF), largely under highly controlled conditions, has established that plant species differentially alter abiotic and biotic soil conditions that in turn affect growth of other conspecific and heterospecific individuals in that soil (Bever, 1994;Bever et al., 1997;Wardle et al., 2004;Ehrenfeld et al., 2005;Kulmatiski and Kardol, 2008;Mangan et al., 2010;Putten et al., 2013). The broad conclusion from this research is that positive and negative PSF can shape community composition and ecosystem functioning by driving patterns of plant diversity, succession, and invasion (Bever et al., 1997;Mills and Bever, 1998;Klironomos, 2002;Reynolds et al., 2003;Mangan et al., 2010;Johnson et al., 2012;Van Der Heijden et al., 2018). Yet most PSF experiments are rarely conducted under the environmental conditions where large-scale ecological and evolutionary mechanisms occur, that is, in the field. Various findings indicate differences in PSF between greenhouse and field experiments (Putten et al., 2016;Schittko et al., 2016;Florianová and Münzbergová, 2018;Heinze and Joshi, 2018;Kivlin et al., 2018;Forero et al., 2019) which suggest that inferences from greenhouse studies may not accurately represent how PSF functions on the landscape (Kulmatiski and Kardol, 2008;Putten et al., 2013Putten et al., , 2016Smith-Ramesh and Reynolds, 2017;Crawford et al., 2019;De Long et al., 2019). While greenhouse experiments have the advantage of maintaining control over non-focal variables and are crucial for confirming and falsifying mechanisms, recent field based PSF research has demonstrated that thoughtful field study designs can generate tractable and interesting results . By not accounting for the surrounding environmental matrix in which plant-soil interactions occur, it is difficult to accurately predict community composition and productivity based on the plant species that are either inhibited by negative PSF or those that persist by positive PSF.
Competition, stress, and disturbance are common environmental pressures that occur across the landscape. To cope with these pressures, plants can differentially allocate resources to growth, reproduction, or maintenance. Selection for each of these strategies depends on tradeoffs for traits that allow tolerance to competition, stress, or disturbance events. This phenomenon, known as the competition-stress-disturbance (C-S-D) hypothesis (Grime, 1977), predicts that competitive plants thrive in ecosystems with high competition because of evolved traits such as large size and extensive fine roots. Likewise, traits such as high stem and root storage of resources allow stresstolerant plants to thrive under highly stressful conditions, while disturbance-resilient plants thrive in frequently or intensely disturbed ecosystems due to evolved traits such as fast growth rates and fast reproduction. Changes in plant traits in response to these environmental pressures may subsequently impact plant-soil interactions by altering soil biota and chemistry. Although the C-S-D hypothesis is well-studied (Herms and Mattson, 1992;Reich, 2014;Rosado and de Mattos, 2017), it does not directly describe the role of plant-soil interactions in plant response to the biotic and abiotic environment (which C-S-D represents) that may enhance or impede a plant's ability to cope with competition, stress, or disturbance. For example, when accounting for soil microbial dynamics under well-watered vs. drought-stress conditions, a greenhouse study (Lau and Lennon, 2012) found that Brassica rapa fitness increased under drought when focal plants were grown with "dry-adapted" soil microorganisms. As such, examining PSF across environmental gradients of competition, stress, or disturbance will provide insight into plant-soil dynamics that influence plant species' persistence or decline under more realistic field conditions.
Competition is a prominent driver of plant fitness, community composition, and coexistence (Tilman, 1982;Callaway and Walker, 1997;Chesson, 2000;Aschehoug et al., 2016) and may influence the strength and direction of PSF in several ways. While competition reduces resource availability, a plant may benefit from its "home" microbiome and thus exhibit positive PSF if soil mutualists increase the availability of limiting resources (Bessler et al., 2012). Conversely, the phenomenon of reduced survival and density of conspecific seedlings near mature or parent conspecifics (Burkey, 1994;Bell et al., 2006;Mangan et al., 2010;Rolhauser et al., 2011;Reinhart et al., 2012;Comita et al., 2014) suggests that intraspecific competition may promote negative PSF due to accumulation of host-specific soil pathogens. The relative importance of competition and PSF on plant performance was examined in a recent meta-analysis (Lekberg et al., 2018). Distinguishing between both inter-vs. intraspecific and low vs. high-density competition and across multiple PSF treatments Lekberg et al. (2018) found antagonistic interactive effects of competition and PSF. In other words, the combined effect of competition and PSF was greater than the individual effects of PSF or competition alone, often leading to highly reduced plant growth and demonstrating how biotic interactions can alter PSF outcomes. Plant growth responses to PSF are likely also modified by plant responses to changes in the abiotic environment, although this was not tested by Lekberg et al. (2018).
As human-induced global change accelerates, stress (defined here as prolonged or continuous environmental pressures) and disturbance events defined here as sudden, temporally constrained changes in the environment) (Hillebrand and Kunze, 2020) are predicted to rapidly increase (Allen et al., 2010;Dai, 2013;Barbero et al., 2015). These types of environmental pressures will likely have additionally profound effects on plant distributions and plant-soil interactions. Increasing evidence suggests that plant response to environmental stressors including drought (Lau and Lennon, 2012;Vílchez et al., 2016;Kannenberg and Phillips, 2017), herbivory (Badri et al., 2013), and salt tolerance (Qin et al., 2016) can be mediated in part by interactions with soil biota.
As with competition, environmental stress, or disturbance can alter resource availability (e.g., soil nutrients, light, water). Fertilization is often considered a benefit to plant growth, however there is evidence that it can cause an imbalance in PSF. A surplus of nutrients from nutrient deposition could promote positive PSF if plant defense against pathogens is enhanced with high resource availability (Smith-Ramesh and Reynolds, 2017). Positive PSF may also occur under low-resource stress if plants rely on host-specific soil microbes for limiting nutrients (Reynolds et al., 2003;Revillini et al., 2016). Alternatively, nutrient inputs in excess of plant growth demands can shift the balance of plant-microbe interactions and prompt dissociations between plants and their home soil biota (Wallenda and Kottke, 1998;Treseder and Allen, 2002;Revillini et al., 2016). Nutrient inputs from fertilizer applications can lead to limitation of other micronutrients (Whalen et al., 2018) and changes in soil chemistry (Erisman et al., 2013) that modify a plant's association with the soil microbiome. Negative PSF could also occur if soil pathogens thrive under resource-rich conditions (Hersh et al., 2012;Spear et al., 2015).
Other types of environmental stressors like aboveground herbivory could promote negative PSF if grazed plants are less able to defend from soil pathogens that are more abundant in conspecific-conditioned soil (Smith-Ramesh and Reynolds, 2017). Negative PSF under herbivory may also occur if grazed plants reduce carbon allocation to roots and thus reduce the ability to support host-specific microbial mutualists (Smith-Ramesh and Reynolds, 2017). Alternatively, positive PSF could occur if grazed plants allocate growth to roots and thus support host-specific microbial mutualists (Smith-Ramesh and Reynolds, 2017). Disturbance events such as fire could promote either negative or positive PSF depending on how soil microbes respond to fire disturbance. Fire has generally been found to decrease microbial biomass and diversity (Dooley and Treseder, 2012;Pressler et al., 2019;Whitman et al., 2019) due to heatsterilization of microbes and loss of soil carbon valuable for soil microbes (Dooley and Treseder, 2012). Negative PSF could occur if abundance of host-specific mutualists declines from fire, whereas positive PSF could occur if fire reduces soil pathogen abundance. However, it is possible that increases in ash deposition following fire could stimulate microbial growth by increasing availability of soil inorganic nitrogen and alleviating nutrient limitation (Rau et al., 2008;Schafer and Mack, 2010;Dooley and Treseder, 2012). In this scenario, negative or positive PSF could occur depending on whether host-specific mutualists or pathogens thrive from nutrient availability.
Predicting PSF accurately on the landscape also requires examining the methods under which plant growth response is measured and the methods that are used to condition the soil microbiome. PSFs have been found to vary between greenhouse and field experiments. An experiment by Schittko et al. (2016) for example, found that the majority of plant species in an experiment exhibited positive PSF under controlled greenhouse conditions, but did not exhibit any significant PSF under natural field conditions. A meta-analysis of the PSF literature seems to corroborate this observation-experiments conducted in the greenhouse often produce larger effect sizes in PSF than those conducted in the field Forero et al., 2019). This broad pattern likely suggests that the true importance of PSF on the landscape is different than what is measured in greenhouse studies. Importantly, experiment location informs interpretation of the specific mechanisms of PSF. Field experiments testing PSF allow for understanding how a natural environment influences soil effects on plants, whereas PSF experiments conducted in the greenhouse do not provide this level of understanding (Smith-Ramesh and Reynolds, 2017).
Similarly, how the soil microbiome is conditioned can greatly influence plant growth. An important distinction between greenhouse-conditioned soil inoculum and field-conditioned soil inoculum is that greenhouse-conditioned soil contains primarily only microbes associated with the focal conditioning plant(s), whereas field-conditioned soil contains microbes associated with site-specific edaphic characteristics in addition to microbes associated with the focal conditioning plant(s). Subsequent differences in microbial composition can alter plant growth. For example, the perennial forb Centaurea maculosa exhibited negative PSF with field-conditioned soil inoculum and positive PSF with greenhouse-conditioned soil inoculum (Callaway et al., 2004). In this instance, greenhouse-conditioned soil could result in the accumulation of microbial symbionts or the loss of pathogens compared to field-conditioned soil . Similar to the role of experiment location, soil conditioning source can also inform interpretation of the specific mechanisms of PSF. Using soil conditioned by focal plants in the greenhouse limits understanding how the environmental context of a natural field environment influences plant effects on soil (Smith-Ramesh and Reynolds, 2017).
Several recent papers have specifically called for empirical studies that measure PSF under manipulations of abiotic and biotic variation (Putten et al., 2013(Putten et al., , 2016Smith-Ramesh and Reynolds, 2017;Lekberg et al., 2018;Crawford et al., 2019). However, no concerted effort has been made to assess the current extent and relative importance of PSF across environmental gradients of competition, stress, or disturbance. We expanded upon a recent effort (Lekberg et al., 2018) that analyzed the interaction between PSF and plant competition. Lekberg et al. found that PSF varied when plants were exposed to competitors or not. We build on that study by examining PSF when plants are exposed to environmental stress or disturbance, to leverage our understanding of plant C-S-D strategies in making predictions of PSF. We also investigated how PSF responses to C-S-D may vary in greenhouse and field conditions to further develop PSF predictive frameworks. Using data from 300 independent manipulations from 76 publications, we examined plant growth responses to multiple forms of competition, stress, and disturbance across a range of PSF methods to identify trends in plant-soil interactions across biotic and abiotic environments and experimental conditions. To understand the relative importance of environmental context as well as common methods for PSF, we addressed the following questions to better understand and predict what PSF outcomes might be expected under realistic field conditions: (1) Can competition, stress, and disturbance alter the direction and/or strength of PSF? (2) Do particular types of competition, stress, and disturbance affect the direction and/or strength of PSF more than others? (3) Do methods of conducting PSF research (i.e., greenhouse vs. field experiments and whether the source of soil conditioning inoculum is from the field vs. greenhouse) affect plant growth responses to PSF or competition, stress, and disturbance, or the interaction of both?

Paper Selection
We collected previously published, peer-reviewed data of PSF in response to experimental manipulations of competition from studies used in a recent meta-analysis (Lekberg et al., 2018) supplemented with a search conducted in Web of Science for more recent competition studies as well as stress or disturbance. We searched through April 2019 for publications that crossed PSF experiments with one of the environmental manipulations of competition, stress, or disturbance using the following search terms "plant-soil feedback * " OR "PSF" AND "competi * "; "plantsoil feedback * " OR "PSF" AND "stress" OR "drought" OR "herbiv * " OR "precipitation" OR "temperature" OR "salinity" OR "light" OR "fertiliz * " OR "enrichment"; "plant-soil feedback * " OR "PSF" AND "disturbance" OR "mining" OR "mine tailings" OR "wind" OR "hurricane" OR "tornado" OR "fire" OR "grazing" OR "agriculture." We identified additional publications not found in our initial Web of Science search by searching for studies that had cited publications included in our initial data set. Several publications were included twice in our data set because they measured PSF under multiple manipulations of competition, stress, or disturbance (Larios and Suding, 2015;Shivega and Aldrich-Wolfe, 2017;Yu et al., 2017;Hawkins and Crawford, 2018;Zhao et al., 2018), or under more than one stress or disturbance level (Heinze et al., 2016;Valliere and Allen, 2016).
We screened publications for studies that included (1) soil treatment methods indicative of a manipulative PSF experimental design (as detailed below), (2) plant growth responses to soil treatments, specifically, either aboveground biomass or plant height, (3) factorial design of PSF treatments crossed with some manipulation of either competition, stress, or disturbance, where the experiments were undertaken in the field or in the greenhouse and (4) measures of mean, error, and sample size for plant growth in all treatments. We did not have any criteria for the length of the study. We excluded one publication (Brandt et al., 2015) and experiments from several publications (Coykendall and Houseman, 2014;Maron et al., 2016;Zuppinger-Dingley et al., 2016) from the Lekberg et al. (2018) data set because of the absence of a full factorial design. We found nine additional PSF x competition publications that were not included in the Lekberg et al. (2018) meta-analysis (de la Peña et al., 2010;Chen et al., 2012;Zhang et al., 2014;Chung and Rudgers, 2016;Bezemer et al., 2018;Hawkins and Crawford, 2018;Xue et al., 2018;Zhao et al., 2018;Lozano et al., 2019). Soil feedback manipulations for PSF experimental design were conducted in three ways in the included publications: (1) soil conditioned from the focal species ("home" soil) or from a heterospecific species ("away" soil) (hereafter, home-conditioned vs. away-conditioned), (2) live soil (majority of soil biota active; hereafter referred to as "active soil") or sterilized soil (majority of soil biota absent inactive; hereafter referred to as "inactive soil"), or (3) soil untreated or treated with fungicide (non-fungicide treated soil hereafter referred to as "diverse soil"). All studies that included treatments of untreated vs. fungicide soil or treatments of live vs. sterilized soil used soil that had been conditioned by plants. Publications that measured competition included one or two types of competition treatments, either the number of plants differed between treatments and the focal plant grew alone or with other plants (here referred to as "alone-together"). In other instances, competition was quantified as equal number of plants between treatments and the focal plant was exposed to either intra-or interspecific competition (here referred to as "interspecific-intraspecific"). The soil treatments were then added factorially, either directly as field-conditioned inocula (at various quantities and according to treatments) or as a second phase conditioned soil from a two-part design in which plants were initially grown in pots in a greenhouse conditioning phase using a field soil inoculum. Soil conditioned from the first phase was then used as the inoculum in the PSF phase of the experiment (see  for details of these designs).

Data Collection
Mean values and measures of plant growth were collected from text and tables in the main publication and/or Supplemental Information. We used GraphClick (Arizona Software) to extract mean and standard error values from figures when raw data was not provided. If not provided, standard deviations were back calculated from standard errors and sample sizes (SD = SE × √ n). In cases where data was not clearly available in the publication, we contacted the authors. We excluded two publications (Medina-Roldán et al., 2012;Kaisermann et al., 2017) for which we received no response from authors, and therefore could not include these studies. Some studies measured the performance of multiple focal species, and thus included multiple experiments. Some studies contained multiple trials within an experiment in which a focal species was examined under multiple treatments (i.e., multiple home/away soils, multiple competitors, or multiple levels of stress or disturbance). For each record in our dataset, we recorded the type of soil feedback manipulation (as described above) and environmental manipulation (competition manipulation described above). Stress manipulations consisted of drought, fertilization (as a representation of nutrient deposition), grazing/herbivory, shade (light availability), mining, and temperature. Disturbance manipulations consisted of fire and tornado. We differentiated stress manipulations as those that represented prolonged or continuous environmental pressures experienced by focal plant(s) and disturbance manipulations as those that represent sudden, temporally constrained changes in the environment (Hillebrand and Kunze, 2020).
We recorded duration of the feedback experiment and multiple descriptors of the focal plant. We also recorded the location of where the soil was collected, as well as the source of inoculum conditioning phase (i.e., collected directly from the field = conditioned in the field, or collected from a training phase in pots under controlled conditions = conditioned in the greenhouse). Field-conditioned soil represents microbiota associated with specific plant species in the field as well as microbiota associated with the edaphic conditions of that site (i.e., pH, nutrient levels, soil moisture), whereas greenhouseconditioned soil represents mostly microbiota associated with the plants that were used to condition the inoculum. For each record we also recorded experiment location. Experiment location was defined as greenhouse if the plant growth response to soil phase was conducted in the greenhouse. Location was defined as field if the response phase was conducted in a natural field environment. Extracted data of all publications included in the dataset is available in the (Table S1).

Effect Size Calculations
We conducted an interaction meta-analysis to assess effects of PSF, competition, stress, and disturbance across PSF methods on plant growth using the relative interaction intensity (RII) as the effect size metric (Armas et al., 2004). We preferred this metric over the log response ratio -a widely used metric in ecological meta-analyses-because unlike the log response ratio, RII is bounded between −1 and 1, and therefore symmetrical around zero, and it can be calculated if plant growth is zero in control groups. Following Armas et al. (2004), RII is calculated as: We therefore calculated the effect of PSF manipulation as: And the effect of competition, stress, or disturbance (i.e., C-S-D) manipulation as: To quantify the combined effect of PSF with C-S-D, we followed the calculation in Kivlin et al. (2013) modified from Armas et al. (2004) for a two-factor RII. The interactive effect of soil feedback and competition, stress, or disturbance can be described as the effect of PSF when competition, stress, or disturbance is present compared to the effect of PSF when competition, stress, or disturbance is absent. This was calculated as: In these equations, Y is the mean plant growth for t = treatment or c = control for the competition, or stress, or disturbance treatment, and f = soil feedback imposed (away-conditioned soil, live soil, or non-fungicide treated soil) or n = no soil feedback imposed (home-conditioned soil, sterilized soil, or fungicidetreated soil). To calculate the 95% confidence intervals around the means for each record, variance was weighted by the sample size (n) and calculated using the standard deviation (s). For each record, we followed the calculations used in Kivlin et al. (2013) to calculate variance. Variance for the main effect of PSF was calculated as: , and the variance of the main effect of competition, stress or disturbance as:

Data Analysis
We used the rma.mv function in the metafor package (Viechtbauer, 2010) in R 3.5.2 for all analyses. In all analyses, we separated the dataset by competition studies, stress studies, and disturbance studies and thus ran competition, stress, and disturbance models separately. Prior to analyses to assess our questions, we first identified the individual effects of PSF and C-S-D by building separate multivariate mixed effects models using RII PSF and RII CSD as response variables and PSF Vi and CSD Vi as the variance. These models tested for the main effects of PSF and C-S-D averaged across soil feedback manipulations and across competition, stress, disturbance manipulations. Thus, we tested for differences in plant growth in response to types of soil feedback and types of competition, stress, disturbance by including soil manipulation and C-S-D manipulation as moderators in separate models. A moderator in the rma.mv function is analogous to a fixed effect in an ANOVA model (Viechtbauer, 2010) and allows the model to calculate the effect size of specific levels of a factor.
To address the first question of whether competition, stress, and/or disturbance alters the direction and/or strength of PSF, we built multivariate mixed effects models using RII Int as the response variable and as Int Vi as the variance to identify the main interactive effect of PSF x C-S-D (averaged across soil feedback manipulations and across competition, stress, disturbance manipulations). To address the second question of whether particular types of competition, stress, and disturbance affect the direction and/or strength of PSF more than others, we built separate multivariate mixed effects RII Int models that included C-S-D manipulation as a moderator. To address the third question of whether methods of conducting PSF research alter the individual effects of PSF, or C-S-D, or the interactive effect of PSF x C-S-D, we built additional RII Int mixed effects models using experiment location (field vs. greenhouse) and soil inoculum conditioning source (field-conditioned vs. greenhouseconditioned) as moderators.
We included plant species as a random effect in all models because species are not independent and past evolutionary history may affect plant response regardless of treatment (Wooliver et al., 2017). Including plant species as a random effect also allows us to make comparisons across studies. To further account for phylogenetic variation of plant growth responses, we replicated the analysis using plant family as a random effect. Results did not vary between the species and family analyses. While it is likely that the strength of PSF varies with plant ontogeny (Kardol et al., 2013), the included studies lacked sufficient replication in experiment duration to use duration as another random effect. For all models we performed post-hoc tests using the linearHypothesis function in the car package (Fox et al., 2013) to test whether effect sizes differed from one another. We report QM as the test statistic for moderator coefficients of the rma.mv models. We considered results to be significantly different from zero if α ≤0.05.

Distribution of Studies
Using the selection criteria, we identified 300 studies from 76 publications that measured plant growth in response to different soil PSF methods in a specific environmental context (competition, stress, or disturbance). Of these, 199 studies (43 publications) measured PSF with competition, 95 studies measured PSF with stress (34 publications), and 5 studies measured PSF with disturbance (2 publications). The majority of studies were conducted in the greenhouse and 86% of all studies used grasses or forbs as the focal plant species. Field experiments comprised only 9% of all studies (29 studies from 7 publications); 26 field studies focused on competition (6 publications), and 3 field studies focused on stress (1 publication).

Main Effects of PSF, Competition, Stress, and Disturbance on Plant Growth
In general, plant growth responses to the main effects of PSF, stress, and disturbance were neutral and the response to competition was negative ( Figure 1A). PSF had a generally neutral effect on plant growth ( Figure 1A; p = 0.28; Table S2, row 3). However, the direction of PSF varied by type of soil feedback manipulation (Figure 2A; QM = 28.43, p < 0.0001; Table S2, rows 5-7). Active (i.e., live) soil reduced plant growth compared to sterile soil (p < 0.0001; Table S2, row 6), and the effect of PSF was neutral in studies that tested home vs. away-conditioned soil (p = 0.28; Table S2, row 7) or in those that tested diverse vs. nondiverse soil (i.e., untreated vs. fungicide treated soil) (p = 0.97; Table S2, row 5). Post-hoc analysis revealed that plant growth was reduced more in live vs. sterile soil manipulations than in home vs. away soil manipulations (χ 2 = 26.84, p < 0.0001; Table S2, row 92) or untreated vs. fungicide manipulations (χ 2 = 3.97, p = 0.05; Table S2, row 90). PSF was similarly neutral in studies that tested home vs. away soil manipulations and untreated vs. fungicide manipulations (χ 2 = 0.20, p = 0.66; Table S2, row 91).
Plant growth responses to competition, stress, and disturbance were variable and depended on the type of competition, stress, or disturbance imposed. Competition in general significantly reduced plant growth ( Figure 1A; p < 0.0001; Table S2, row 32), and the effect of plants grown together with a competitor was nearly 150% greater than the effect of inter-vs. intraspecific competition on plant growth (Figure 2A; χ 2 = 51.78, p < 0.0001; Table S2, row 113). Generally, stress had marginally negative effects on plant growth (Figure 1A; p = 0.22; Table S2, row 56), and plant growth varied among different types of stress manipulations (QM = 94.94, p < 0.0001; Table S2, rows 58-63). As expected, drought (Figure 2A; p = 0.0001; Table S2 row 58) and shade (Figure 2A; p < 0.0001; Table S2, row 61) greatly reduced plant growth, and plant growth was reduced similarly under drought and shade (χ 2 = 2.66, p = 0.20; Table S2, row 132). Fertilization, on the other hand, increased plant growth (Figure 2A; p = 0.0003; Table S2, row 59). Increases in plant growth in response to fertilization were significantly different than reductions in plant growth in response to drought (χ 2 = 29.65, p < 0.0001;

Do the Effects of PSF and C-S-D Vary Between Methods Used to Conduct PSF Research?
Whether the experiment was conducted in the field or greenhouse altered the effects of competition and stress on plant growth, but not the effect of PSF (Figure 3A). There was no difference in the effects of away-conditioned vs. homeconditioned soils (χ 2 = 0.24, p = 0.62;  phase occurred in the field or greenhouse. Studies that used live and sterile soils as the soil feedback manipulation were only conducted in the greenhouse. Conversely, competition reduced plant growth when plants were grown in the field (Figure 3A; p < 0.001; Table S2, row 42) and the greenhouse (Figure 3A; p < 0.001; Table S2, row 43), but the effect was 138% stronger in field experiments than greenhouse experiments (χ 2 = 48.21, p < 0.001; Table S2, row 122). This trend was driven by alone/together competition studies because inter/intra-specific competition studies were only conducted in the greenhouse. Only grazing/herbivory stress was conducted in both greenhouse and field settings, and thus it was the only stressor we could analyze between experiment location. Grazing/herbivory stress enhanced plant growth in the greenhouse (Figure 3A; p < 0.001; Table S2, row 71), whereas its effects on plant growth were neutral in field experiments ( Figure 3A; p = 0.48; Table S2, row 70). Grazing/herbivory stress enhanced plant growth over 300% more in greenhouse experiments than in field experiments ( Figure 3A; χ 2 = 11.31, p < 0.001; Table S2, row 157). Disturbance studies were only tested in greenhouse conditions. The interactions between PSF and competition were neutral in both field and greenhouse experiments (Figure 3B; QM = 2.52, p = 0.28; Table S2, rows 181-182) and post-hoc analysis showed that these effects were similar between field and greenhouse experiments (χ 2 = 0.95, p = 0.33; Table S2, row 245). Contrary to the finding that the main effect of grazing/herbivory was only positive in greenhouse experiments, there was a synergistic effect of PSF and grazing/herbivory only in field experiments ( Figure 3B; p = 0.04; Table S2, row 202). In other words, away-conditioned, active or diverse soil enhanced plant growth under grazing/herbivory in field experiments. However, post-hoc analysis showed that the interaction of PSF and grazing/herbivory in field experiments was not significantly different than the neutral effect of PSF and grazing/herbivory in greenhouse experiments (χ 2 = 2.34, p = 0.13; Table S2, row 268). Disturbance interactions with PSFs were only tested in greenhouse conditions.
In general, soil inoculum conditioning source modified the effect of PSF on plant growth (Figure 4A; QM = 34.93, p < 0.0001; Table S2, rows 23-24). Away-conditioned, active, or diverse soil reduced plant growth when soil inoculum was conditioned in the field (Figure 4A; p < 0.0001; Table S2, row 23) and enhanced plant growth when soil inoculum was conditioned in the greenhouse (Figure 4A; p = 0.007; Table S2, row 24). The various soil feedback manipulations also responded differently depending on the source of the conditioning phase. Active soil reduced plant growth when conditioned in the field (p = 0.03; Table S2, row 26) but had no effect on plant growth when conditioned in the greenhouse (p = 0.40; Table S2, row 27). The effects of away-conditioned soil on plant growth did not vary between field-conditioned (p = 0.20; Table S2, row 29) or greenhouse-conditioned soil (p = 0.50; Table S2, row 30).
The interaction of PSF and competition was similarly neutral between experiments that used field-conditioned and greenhouse-conditioned inoculum ( Figure 4B; χ 2 = 0.04, p = 0.84; Table S2, row 247). In general, the interaction of PSF and environmental stress was slightly synergistic when soil inoculum was conditioned in the field (Figure 4B; QM = 4.92, p = 0.08; Table S2, row 205). This trend was likely driven by drought, fertilization, and grazing/herbivory studies in which the average effects of PSF and each of these stressors was positive, though not statistically significant. However, post-hoc analysis showed that the interaction of PSF and stress with field-conditioned soil inoculum was not significantly different than the neutral effect of PSF and stress with greenhouse-conditioned soil inoculum (χ 2 = 0.52, p = 0.47; Table S2, row 270). Conversely, the interaction of PSF and mining stress was slightly synergistic when soil was conditioned in the greenhouse (QM = 4.73, p = 0.09; Table S2, row 215). Away-conditioned, active, or diverse soil enhanced plant growth 300% more under mining stress when soil inoculum was conditioned in the greenhouse rather than conditioned in the field (χ 2 = 3.81, p = 0.05; Table S2, row 276). Studies that manipulated temperature were not included in this analysis as these studies only included greenhouse-conditioned soil.

DISCUSSION
In general, our meta-analysis revealed important patterns that may aid predictions of PSF under natural or field conditions where plant-plant competition and environmental stressors are common. The results of this meta-analysis suggest that soil microbes may generally reduce plant growth where plant-plant competition for resources occurs and may enhance plant growth under drought stress conditions. These results also highlight an important research gap in examining PSF under environmental disturbance. Additionally, this meta-analysis suggests that where PSF experiments are conducted (field or greenhouse) affects the outcome of the study-especially when plants are competing. The nature of the soil conditioning manipulation (live vs. sterile, home vs. away, or fungicide application) can change plant growth responses, and again this is most evident when plants are competing for limiting resources. Our meta-analysis corroborates and expands upon Lekberg et al. (2018) to confirm that environmental context writ large can alter plant growth responses in PSF experiments, at least in grass and forb species with which the majority of PSF studies have been conducted.

Plant Competition, and Specific Environmental Stressors and Alter Direction and Magnitude of PSF
Competition shapes resource availability (Tilman, 1982(Tilman, , 1990 which can subsequently affect PSF, and the data to date show that PSF is often more examined in response to plant-plant competition than to environmental stress and disturbance. While growing next to a plant compared to growing alone more greatly reduced plant growth than when plants were growing with conspecifics or heterospecifics, types of competition reduced plant growth differently depending on the type of soil feedback imposed. When focal plants were grown in awayconditioned soil, interspecific competition reduced plant growth more than intraspecific competition. This likely indicates that interspecific competition either enhanced negative effects of nonhost associated soil biota on plant growth or reduced beneficial effects of non-host associated soil biota on plant growth. When focal plants were grown in active or diverse soil, interspecific competition reduced plant growth to a similar degree to that of intraspecific competition. This could suggest that taxonomic composition of the soil microbiome is more important than microbial diversity for reducing plant growth under interspecific competition. Conversely, in alone vs. together competition studies, the effect of away-conditioned, active, and diverse soil was similar. This may suggest that just the presence of a soil microbiome is influential for inhibiting plant growth when a plant is in close proximity to a competitor. While the frequency of environmental stress and disturbance are increasingly occurring in terrestrial landscapes, these data from 76 publications suggest that stress and disturbance can have both negative and positive effects on plant growth individually (e.g., drought and shade reduce growth but fertilization increases growth) but when combined with PSF, lead to few synergistic or antagonistic effects. While the combined effect of PSF and environmental stress was generally synergistic, the trend was driven by drought stress studies. Away-conditioned, active, or diverse soil enhanced plant growth under drought conditions. Though the exact mechanisms to explain this trend remain unclear, soil microbes may facilitate plant growth under lowresource stress if plants rely on microbes for acquiring limiting nutrients (Reynolds et al., 2003;Revillini et al., 2016). Although this meta-analysis found that the combined effect of PSF and grazing/herbivory was on average synergistic but not statistically significant, grazing/herbivory could potentially induce benefits of the soil microbiome to plant growth if plants allocated resources to roots and subsequently supported microbial mutualists (Smith-Ramesh and Reynolds, 2017).
This meta-analysis found that the combined effects of PSF and environmental disturbance were neutral for both fire and tornado disturbance treatments. It is unclear, however, whether this is a true pattern across fire and tornado disturbance or across many types of environmental disturbance because this dataset only includes five studies from two publications (Nagendra and Peterson, 2016;Senior et al., 2018). Soil microbiomes and their interactions with plants are likely to be affected by drastic landscape changes brought about by environmental stress and disturbance, yet this meta-analysis shows that there is not yet enough data to predict how PSF responds to environmentally disruptive events. This gap of knowledge highlights the necessity of PSF experiments to incorporate manipulations that are representative of environmental variation under global climate change (Putten et al., 2016;Long et al., 2019).

Experiment Location and Soil Conditioning Source Affect the Outcomes of PSF Research
Comparisons of field vs. greenhouse studies showed that PSF did not vary by experiment location. Competition had a more negative effect when experiments were conducted in the field rather than in the greenhouse. However, the combined effects of PSF and competition were similarly neutral between field and greenhouse experiments. Grazing/herbivory only increased plant growth in greenhouse conditions. The combined effect of PSF and grazing/herbivory was synergistic in field conditions and neutral in greenhouse conditions, though further analyses indicated that the interactive effects of PSF and grazing/herbivory were not significantly different between field and greenhouse experiments.
It is important to note that only 9% of these studies were conducted in the field, suggesting that we have little inference for the strength or direction that environmental variation will demonstrate when interacting with PSF to affect plant growth. The few interactive effects shown here both support and contradict studies showing that PSF can differ between greenhouse and field experiments in different environments (Putten et al., 2016;Schittko et al., 2016;Florianová and Münzbergová, 2018;Heinze and Joshi, 2018;Kivlin et al., 2018) demonstrating how little is understood about the function of PSF in nature (Forero et al., 2019).
PSF varied by soil inoculum conditioning source. Awayconditioned, active, or diverse soil reduced plant growth if soil was conditioned in the field and enhanced plant growth if soil was conditioned in the greenhouse. The combined effect of PSF and competition was similarly neutral between studies that used field-conditioned and greenhouse-conditioned soil inoculum. The combined effect of PSF and stress was slightly synergistic when soil inoculum was conditioned in the field, though this was not significantly different than the neutral effect of PSF and stress with greenhouse-conditioned soil inoculum. The combined effect of PSF and mining stress was slightly synergistic when soil was conditioned in the greenhouse. Awayconditioned, active, or diverse soil enhanced plant growth 300% more under mining stress when soil inoculum was conditioned in the greenhouse rather than in the field. While it is difficult to identify detailed mechanisms, the trend from this metaanalysis of increased plant growth in greenhouse-conditioned soil relative to field-conditioned soil may be due to lower microbial diversity in greenhouse-conditioned soil. Specifically, greenhouse-conditioned soil may contain lower abundance of pathogens than field-conditioned soil (Callaway et al., 2004;. These results demonstrate that there are multiple methodological approaches that both allow us to infer how field studies may respond and to show how variation in methods can change interpretation of results. The finding that effects of PSF can differ between greenhouse and field experiments provides justification for pairing PSF experiments in the greenhouse with those conducted in the field, when possible (Smith-Ramesh and Reynolds, 2017). Additionally, that the source of inoculum for PSF experiments can also mediate interactive effects show that careful interpretation is required of studies that use each method (Smith-Ramesh and Reynolds, 2017). Studies that use field-conditioned inoculum for example, should be mindful to infer that focal plants are responding to microbes associated with the focal conditioning plant(s) in addition to microbes associated with site-specific soil characteristics. Focal plants in studies that use greenhouse-conditioned inoculum, on the other hand, are responding primarily to microbes associated with the focal conditioning plant(s). Moreover, it should be noted that 86% of the focal plants examined in this meta-analysis are grasses and forbs. While, overall, we did not see mean differences in effect sizes among plant functional groups (data not shown-see results in Table S2), this result is related to by the reduced power to detect an effect due to low sample sizes of trees and shrubs that have been studied to date.
An important distinction between greenhouse-conditioned soil inoculum and field-conditioned soil inoculum is that greenhouse-conditioned soil contains primarily only microbes associated with the focal conditioning plant(s), whereas fieldconditioned soil contains microbes associated with site-specific edaphic characteristics in addition to microbes associated with the focal conditioning plant(s).

CONCLUSIONS
Environmental context writ large can change plant growth responses in PSF experiments, a conclusion supported by our analyses comparing the interactive effects of PSF with a range of competition, stress and disturbance types relative to their individual effects. Our results have direct application in ecological conservation and restoration because we show that positive PSF is synergistically enhanced by drought and stress. Additionally, our study suggests that field experiments may yield different responses than greenhouse experiments. Data from these 76 studies show the need for more research on PSF across environmental stresses and disturbances (i.e., two-thirds of these studies were conducted on competition) and the need for increased representation of a wide diversity of focal plant species, because the majority of PSF studies were conducted with grasses and forbs. The lack of studies investigating true gradients of stress and disturbance (i.e., multiple experimental levels, rather than presence-absence of different types of stress or disturbance) indicate that we know very little about how PSF effects will respond to stress and disturbance on the landscape. Our metaanalysis enables future research into plant community dynamics in a changing world. There is, therefore, much empirical work to look forward to.

DATA AVAILABILITY STATEMENT
All datasets generated for this study are included in the article/Supplementary Material.

AUTHOR'S NOTE
This manuscript has been authored in part by UT-Battelle, LLC, under contract DE-AC05-00OR22725 with the US Department of Energy (DOE). The US government and the publisher, by accepting the article for publication, acknowledges that the US government retains a non-exclusive, paid-up, irrevocable, worldwide license to publish or reproduce the published form of this manuscript, or allow others to do so, for US government purposes. DOE will provide public access to these results of federally sponsored research in accordance with the DOE Public Access Plan (http://energy.gov/downloads/doe-publicaccess-plan).

AUTHOR CONTRIBUTIONS
KB, JM, SB, CL, LM, MP, JLS, and IW screened all potential papers and coded data for papers that were included in the meta-analysis, curated by KB and JM. KB and JM ran all initial analyses and made figures. SK guided data analysis. KB ran all final analyses and made all final figures. KB, JM, JAS, and JB coordinated writing, to which all authors contributed. Ideas were conceived of by all authors.

ACKNOWLEDGMENTS
We thank the organizers of this special feature and the University of Tennessee. Special thanks to individual authors of published studies who generously contributed requested data. We also thank the two reviewers who provided very beneficial feedback on this manuscript.