Induced plasticity alters responses to conspecific interactions in seedlings of a perennial grass

Plants can interact with different individuals in their lifetime which may lead to plastic response that affect performance. If conspecific interactions are altered through previous plastic responses that could affect stabilizing niche mechanisms, in which conspecifics compete more intensely to promote diversity and coexistence. Here, I show interactions between Pascopyrum smithii and conspecifics resulted in largely canalized traits, whereas P. smithii with an invasive grass, Bromus tectorum resulted in plastic responses for root mass (p = 0.02), shoot mass (p < 0.0001), root mass fraction (p = 0.003) and plant height (p < 0.0001). A subset of individuals transplanted from these two interaction treatments which were moved with new, same aged conspecifics showed that previous interactions led to differing trait relationships: increases in the number of leaves for the interspecific-induced plants were related to increases in non-focal leaf production, whereas increases in the number of leaves for the intraspecific-induced plants were related to decreases in the non-focal plants (R2 = 0.52, p = 0.006). These results suggest that previous intraspecific interactions intensify conspecific competition and stabilize subsequent interactions with conspecifics by imposing greater competition, and that invasive-interspecific interactions can weaken stabilizing niche mechanisms, thus negatively influencing species coexistence.

Plant plastic responses include morphological, physiological, and behavioral changes that can influence resource capture, fitness, and survival [1][2][3][4] . Plants respond plastically to both biotic and abiotic factors, for example, water stress induces greater root length which coincides with greater water uptake 5 and maize competing with Sorghum reduced growth rate and root production in maize 6 . The identity of the interacting neighbor matters; evidence suggests that species in their invasive range induce plastic responses in co-occurring native species, and their presence results in transgenerational and adaptive signals in the native plant [7][8][9][10] . Additionally, Bromus tectorum, is highly competitive with root system traits primed for rapid resource capture 11,12 which influenced greater allocation to root mass and fine roots by the native Poa secunda following competition 13 , demonstrating the impact invasive species interactions have on plastic responses of native species.
In addition to induced plastic responses, a dearth of evidence shows that biotic-induced plastic responses influence subsequent interactions [14][15][16] . One example showed that leaf damage that simulated herbivory in clipped sagebrush plants induced production of volatile compounds in nearby tomato plants which reduced herbivory for the tomato plants 17 . Plants, particularly perennial plants, can interact with different plant species in their lifetime due to immigration, exposure to disturbance, or changes in native and invasive species cover, that can influence subsequent interactions. Additionally, native species out planted in restoration settings may experience interactions in nursery settings that could affect plant performance when out planted. Alexander et al. 18 transplanted plant communities to a warmer climate and showed that this resulted in intensified competition with new competitors. While this provides support that plant-plant interactions can change when meeting new neighbors, most plant studies that evaluate the impact of biotic-induced plasticity involve impacts from organisms of another trophic level (e.g., insects 17 ). This reveals a knowledge gap in our understanding of plant-plant induced changes to subsequent plant interactions.
Importantly, there are key relationships that can have ramifications for plant communities if altered by induced plasticity. Specifically, the strength and direction of conspecific interactions are important to species coexistence for stabilizing niche mechanisms 19,20 , in which population self-limiting mechanisms important to coexistence and diversity are supported by negative conspecific interactions. Therefore, changes to conspecific interactions have the potential to affect stabilizing mechanisms and species coexistence. Furthermore, fitness inequalities [19][20][21] are also important to species coexistence and they are trait differences that lead to inferior and  22 showed via simulations under competition-colonization trade-off-based species coexistence mechanism, that plasticity broadened the number of species, environments, and competitive contexts that resulted in coexistence. Plasticity also helped to maintain competitive ability across different water levels and strengthened stabilizing differences of some traits 23 . Generally, plasticity can contribute to complementary resource capture and trait differentiation in communities that may contribute to niche differentiation 24,25 , and thus plasticity is crucial to consider to advance our understanding of plasticity effects on species coexistence and community dynamics. This can be achieved by testing the impacts of previous interactions on subsequent interactions via studies that directly manipulate plasticity and quantify the effects on interactions 26 .
To address these knowledge gaps on plant-plant induced changes and whether this can affect conspecific interactions, I empirically determine whether intraspecific or interspecific interactions induce plastic responses, and whether these induced responses impact subsequent intraspecific interactions-interactions important to species coexistence. I performed transplant interaction experiments in a greenhouse and used above-and belowground traits and morphology to quantify plasticity and interaction outcomes. I used a grass native to the Western United States, Pascopyrum smithii, and the co-occurring invasive grass, Bromus tectorum, which has been shown to induce plastic responses in native grasses 13 . Then I assessed whether plasticity was induced and whether it led to differing subsequent intraspecific interactions that underlie stabilizing niche mechanisms. I hypothesized that (1) intraspecific and interspecific interactions will induce plastic responses in P. smithii, (2) intraspecific-induced and interspecific-induced focal P. smithii plants will differ in traits and morphology with new conspecifics, and (3) intraspecific-induced and interspecific-induced P. smithii will differently affect new conspecific neighbors.

Materials and methods
Study materials and process. I conducted two experiments in a greenhouse (  www.nature.com/scientificreports/ The focal plant P. smithii is a perennial rhizomatous bunchgrass native to the western United States and the non-focal plant is B. tectorum, a non-native invasive annual grass species. Seeds of P. smithii were purchased in July 2014 (Central Milling Wheatland, UT, USA), and B. tectorum seeds were wild collected from Grand County, UT (38.998°N, − 110.177°W) in June 2013. Seed collection was permitted by the director of the Bonderman Field Station at Rio Mesa and research performed conformed to the Chicago Botanic Garden Ethics Code.
I germinated 500 seeds (germination determined via radicle emergence). Germination timing varied, but I planted all interacting germinants with similar sized emerged radicles to each pot simultaneously to minimize size differences.
Experiments one and two ran from 6 Jun. 2017 to 18 Sept. 2017 with an average of 21 °C/ 16 °C day/night with lighting supplementing ambient lighting from outdoors. I used 20-7.6 cm × 7.6 cm × 15.2 cm rectangular pots (Stuewe & Sons, OR, USA) per treatment for each experiment with fine commercial sand. I watered plants three times weekly, rotated trays weekly, and added Murashige-Skoog (Sigma Aldrich, MO, USA) half strength nutrients (4 g/L) at week two of both experiments. Experiment one: interactions to induce plasticity. The treatments applied were intraspecific interactions with three interacting P. smithii seedlings in each pot (n = 15); interspecific interactions with one P. smithii seedling surrounded by two B. tectorum seedlings in each pot (n = 18); and one P. smithii growing alone (n = 14) serving as the control group (Fig. 1a). Plants were evaluated at the seedling stage as seedlings are more sensitive to competition 27,28 and can exhibit plastic responses in above-and belowground traits 29 , thus plastic responses can be induced at an early life stage. After 4 weeks of growth, plants were randomly designated to either be harvested for data collection or transplanted to new pots with two P. smithii neighbors for experiment two.
At harvest, I gently separated the plants at the root systems under water. Next, I measured the plant height then scanned each focal plant using an Epson expression 10000XL scanner (Epson, CA, USA) and counted the number of leaves and root tips from the scanned images. I recorded root and shoot mass after drying in an oven at 95 °C for 5 days and calculated the functional trait root mass fraction (RMF: root mass/ plant mass) (data; 30 ).

Experiment two: induced plasticity impacts on intraspecific interactions. Focal plants were
transplanted using spatulas from the pots used in the first study after watering the pots and minimizing focal plant root disturbance (Fig. 1b, Supplementary Figure S1). I transplanted focal seedlings to randomized pots with two P. smithii plants of the same age (4 weeks old). Pascopyrum smithii plants in the recipient pots were propagated and grown using the same procedure as experiment one. This resulted in two treatments: (1) intraspecific-induced plants transplanted with two same-age intraspecific plants (n = 14), and (2) interspecificinduced plants transplanted with two same-age intraspecific plants (n = 14) (Fig. 1b), and an undisturbed control treatment in which one P. smithii grew alone for 8 weeks (n = 12; not used for analytical comparisons as the plants were not transplanted or disturbed as in the treatments). I washed plants of sand after harvesting and I counted the cumulative number of rhizomes of all three plants in each pot because the root systems were highly tangled together and with mesh used to keep the sand in place. I scanned the plants then counted the number of leaves for the focal and non-focal plants and dried plants in a drying oven at 95 °C for 5 days prior to weighing focal plant aboveground biomass.
Statistical analyses. All analyses were performed in R version 4.0.2 as well as figures 31 . I applied a square root transformation to shoot mass and RMF data to meet assumptions of normality for experiment one. In experiment one, to test whether plant morphology differs between the control group and both intraspecific and interspecific interaction treatments, I contrasted treatment means of the focal plant with the control group to assess if plasticity was induced. I used linear models with treatment as a categorical predictor variable, and shoot mass, plant height, root mass, and RMF as response variables in separate models. I used a generalized linear model with a Poisson error distribution as the count data distribution 32

Results
Plant morphology differed between the control group and treatment groups for some measures when intraspecific neighbors (P. smithii conspecifics) and interspecific neighbors (P. smithii and B. tectorum) interacted. Root mass, shoot mass, and plant height were greatest in the interspecific treatment and similar between the control and intraspecific treatment, however RMF was lowest in the interspecific treatment and similar between the control and intraspecific treatment, and root tip count and number of leaves did not differ by treatment. Root mass contrast marginally differed by treatment compared to the control (R 2 = 0.15, F 2,36 = 3.15, p = 0.055, Fig. 2a) and between the control and interspecific interaction (p = 0.02), but not for the intraspecific interaction (p = 0.4). Shoot mass contrast differed by treatment compared to the control (R 2 = 0.4, F 2,36 = 12.14, p < 0.0001, Fig. 2b) and between the control group for intraspecific interaction (p = 0.04) and the interspecific interaction (p < 0.0001). Root mass fraction differed by treatment compared to the control (R 2 = 0.22, F 2,36 = 5.1, p = 0.01, Fig. 1c) and between the control and the interspecific (p = 0.003), but not for the intraspecific interaction (p = 0.2). Plant height differed by treatment compared to the control (R 2 = 0.35, F 2,36 = 9.3, p = 0.0006, Fig. 1d) and between the control and the interspecific (p < 0.0001) and the intraspecific interaction (p = 0.05). Root tip count and number of leaves did not differ by treatment compared to the control (p = 0.9 and p = 0.8, respectively 30 ). When intraspecific-induced and interspecific-induced P. smithii interacted with new conspecific neighbors, shoot mass and cumulative rhizome count did not differ between the interspecific-induced and intraspecificinduced focal plants and pots, but number of leaves differed between treatments whereby, interspecific-induced plants produced more leaves. Shoot mass (R 2 = 0.03, F 1,16 = 0.5, p = 0.5, Fig. 3a), and cumulative rhizome count did not differ between treatments (R 2 = 0.03, F 1,16 = 0.56, p = 0.5, Fig. 3b), however the number of leaves differed between the two treatments (R 2 = 0.62, F 1,16 = 0.5, p = 0.0001, Fig. 3c). Leaf trait hierarchy was higher in interspecific-induced (x̅ = 0.05) than intraspecific-induced (x̅ = − 0.57) treatments but did not statistically differ between treatments (R 2 = 0.07, F 2,14 = 1.1, p = 0.3, data not shown 30 ). However, the induced focal plants impact on non-focal plants differed (R 2 = 0.52, F 2,14 = 7.5, p = 0.006, Fig. 4) in which the relationship between induced focal plant leaf production negatively affected average number of leaves produced by the non-focal plants in the intraspecific-induced treatments (y = − 0.5x + 5.9), whereas there was a positive relationship in the interspecific induced treatment (y = 0.3x + 2.8).  www.nature.com/scientificreports/ Interspecific interactions with Bromus tectorum induced plastic responses in Pascopyrum smithii plants, as evidenced by increases in traits and biomass responses for root mass, shoot mass, and plant height, but not for root mass fraction (RMF). Positive interactions with B. tectorum is surprising given the suppressive impacts of B. tectorum on other native grass competitors 13,34 , but similar outcomes have been seen between P. smithii and B. tectorum 35 . Another interpretation could be of strong competition between these plants in which P. smithii increased allocation as a competitive response to reduce potential suppression and performance losses under the presence of B. tectorum. Knowing the response of B. tectorum would provide more insight into this relationship, but measures on B. tectorum were not recorded. Plastic response in plant height to competition may have in-part been due to facilitation or due to the need to improve access to light for photosynthesis, coupled with  www.nature.com/scientificreports/ shade avoidance responses to neighboring B. tectorum or acted as a cue indicating the presence of a neighbor 36 . Plant height is an important trait that influences light capture for photosynthesis and plays an important role in ecological strategy 37 . Shade avoidance responses can include increases in plant height that aid in access to light and can in-turn impose competition on neighbors 38 , but the responses of B. tectorum and non-focal P. smithii are unknown in this study. Traits and morphology following intraspecific interactions, on the other hand, were largely canalized in which root mass and shoot mass did not differ from the control group. RMF and plant height differed from the control indicating plastic responses in these traits. Plastic response in the functional trait, RMF, has ecological significance: RMF can increase strongly to stress 39 , such as water and nutrient stress (a putative stressor in competition) and to competition 6,13,40 . Changes in RMF and plant height are adaptive response as they improve persistence and establishment under the new environment 41 . This may provide benefits in the long-term for interacting P. smithii conspecifics under intense competition. Furthermore, for species coexistence, Bennett et al. 42 showed that stronger intraspecific competition resulted when root allocation was greater and more intense intraspecific competition should directly influence and strengthen stabilizing niche mechanisms 26 .

This study addressed whether interactions induced plastic responses and impact subsequent interactions and shows that interspecific interactions with
Root tip count and number of leaves were canalized for both intraspecific and interspecific interactions. This is surprising for root tips because lateral root tip proliferation responds highly plastically to competition and resources stress 43 , for instance, Lolium perenne experienced an eightfold increase in lateral roots (a major constituent of root tip count) under drought stress 43 .
In experiment two, leaf count was suppressed in intraspecific-induced relative to the interspecific-induced treatment indicating that differences in performance between the focal plants persisted in interactions with new conspecific neighbors. The previous interaction environment also produced a different relationship between the focal plant and non-focal neighbor's leaf production indicating that interspecific-induced and intraspecificinduced plants differed in responses to trait hierarchy with the new neighbors. Trait hierarchies are often explored between species 33,42 , however trait differences exist below the species level 34,44,45 and influence coexistence 46 . For example, Abbott and Stachowicz 46 showed that increasing root mass distance in pairwise competing genotypes of eelgrass (Zostera marina) decreased the probability of genotype coexistence. Generally, the results of the present study suggest that previous intraspecific interactions act to stabilize subsequent interactions. This is achieved by leading to more intense conspecific competition that may result in self-limiting population dynamics that act under a strong source of competition from conspecifics 19,20 . This is corroborated by results from a meta-analysis of phenomenological studies of intraspecific and interspecific interactions 47 . This study presents novel evidence that intraspecific interactions can be weakened by previous interactions with interspecific invasive plants and may play a role in weakening stabilizing niche mechanisms. This invasive plant-mediated change to conspecific interactions may provide some insight into invasive species-mediated population decline if the increasing facilitation results in positive population growth that is not stable under low densities. However, longer term data would be required to address changes to population dynamics. This study originally aimed to also test the impacts of intraspecific-induced and interspecific-induced plants on both new conspecific and new heterospecific neighbors. However, high mortality of the heterospecific plants precluded this assessment. To further understand both intraspecific and interspecific components of stabilizing niche mechanisms and plasticity, future work should also assess effects onto new conspecifics and non-invasive heterospecific plants.
Future work that aims to elucidate the influence of plant-plant induced plastic responses on subsequent interactions should explore study implications in a less controlled environment to assess the influence of environmental heterogeneity, though trait assessments in controlled settings have relevant relationships to performance in the field 48 . This study provides further support to previous research showing that the identity of associates growing with plants prior to out-planting can impact that plant's performance and interactions when transplanted to a new environment 49 . Plants grown as plugs for commercial purposes are usually grown with many intraspecific neighbors 50 including maternal plants grown for seed harvesting. Additionally, plastic responses to intraspecific neighbors compared to interspecific neighbors can also be apparent over multiple generations (e.g., transgenerational plasticity) 51 , so these factors should be considered when growing plants for restoration practices to ensure expected plant performance and meeting restoration goals. Lastly, plant life stage can influence the outcomes of species interactions 52 . For instance, the intensity of conspecific competition for the cycad Dioon sonorense decreased with age, whereas heterospecific competition increased with age 52 . Seedlings, like the ones used here, are thought to be more sensitive to competition 27,28 , suggesting that the current study on seedlings is informative to plant interaction knowledge base. Transplant studies have been instrumental in elucidating the effects of novel competitors and environments 18 , but excavation can impact plant growth. However, because all treatments were excavated in the same manner, this allows conclusions to be drawn based on the interaction treatments imposed. Future studies should explicitly incorporate plant interactions at different life stages in transplant studies and incorporate multiple interactions to determine consequences of plasticity in subsequent interactions.

Conclusions
This study provides evidence of how plasticity induced by interactions may affect stabilizing mechanisms by altering intraspecific interactions and has implications for its impact on species coexistence. I show that previous intraspecific interactions may function to further stabilize niche differences, whereas previous interspecific interactions with a non-native invasive species can weaken them, which has implications for better understanding community dynamics of plants that encounter new neighbors in their lifetime. Future research should assess induced plasticity in different abiotic and biotic contexts to provide more empirical evidence for the ecological consequences of plasticity.