Are Invasive Genotypes Superior ? An Experimental Approach Using Native and Invasive Genotypes of the Invasive Grass Phalaris Arundinacea

The admixture and recombination of individuals from the native range into a new range may lead to the production of invasive genotypes that have higher fitness and wider climatic tolerances than the native genotypes. In this paper, we compare the survival and growth of native EU and invasive NA genotypes when planted back into the native EU range near where the EU genotypes were collected. We test this hypothesis using the invasive wetland grass Phalaris arundinacea. If invasive genotypes have evolved to have higher survival and growth, then they should outperform the native EU genotypes under field conditions that are better suited to the EU genotypes. Individual plants of the wetland grass, Phalaris arundinacea collected from native Europe (Czech Republic (CZ) and France (FR)) and North America (Vermont (VT) and North Carolina (NC)) were planted into common gardens in Trebon, Czech Republic (49.0042 ̊N, 14.7721 ̊E) and Moussac, France (43.9808 ̊N, 4.2241 ̊E). Invasive genotypes from North Carolina (NC) survived as well or better than native genotypes in both the Trebon and Moussac garden. Additionally, invasive NC genotypes suffered higher herbivore damage than native genotypes but their growth and survival were not significantly different than genotypes from the other regions. A companion field experiment that simulated biomass removal through grazing indicated that invasive NC genotypes recovered faster following grazing than genotypes from other regions. Our results suggest that not all invasive genotypes are superior and regional differences in aggressiveness between invasive genotypes are as great as differences between individuHow to cite this paper: Molofsky, J., Collins, A.R., Imbert, E., Bitinas, T. and Lavergne, S. (2017) Are Invasive Genotypes Superior? An Experimental Approach Using Native and Invasive Genotypes of the Invasive Grass Phalaris arundinacea. Open Journal of Ecology, 7, 125-139. https://doi.org/10.4236/oje.2017.72010 Received: December 1, 2016 Accepted: February 6, 2017 Published: February 9, 2017 Copyright © 2017 by authors and Scientific Research Publishing Inc. This work is licensed under the Creative Commons Attribution International License (CC BY 4.0). http://creativecommons.org/licenses/by/4.0/


Introduction
An introduced species may behave invasively in the new range because it already possesses traits that confer invasiveness [1] [2], or alternatively, it may evolve invasiveness in situ in the new range [3] [4] [5] [6]. Moreover, the introduction history of a species can affect the likelihood that it will become invasive [7] [8]. While many plant species suffer a genetic bottleneck when introduced in the new range [9], multiple introductions of a plant species may make it more likely that a species becomes invasive in part because multiple introductions can inflate genetic diversity [8] [10] [11]. In addition, multiple introductions may result in the admixture of genomes that have never come into contact with each other creating novel, invasive genotypes that may express different traits and enhanced fitness [11] [12] [13] [14].
A successful introduction can also depend upon the relationship between the introduced individuals and the new environment. Introduced individuals may have different climatic tolerances than their native counterparts [15] and/or wider climatic tolerances [16] [17] and this may contribute to their success. Testing whether invasive genotypes have a wider climatic tolerance than native genotypes requires planting clones of both known native and invasive genotypes in their home climate and in a different climate to test whether invasive genotypes have greater ability to survive and grow under a new climatic condition.
The increased growth of invasive genotypes found under controlled conditions may not be observed under field conditions [18]. Disentangling whether the reason is due to the traits of the introduced individuals or environmental factors in the field is difficult because traits favored in the new range may incur a cost in the native range. For example, reduced herbivore loads in the new range may have selected for genotypes that have reallocated resources from defense to growth resulting in invasive genotypes having faster growth than native ones [2] [19] [20]. But, if these genotypes are transported back to the native region in which herbivores are more abundant, they may experience higher levels of herbivore damage than their native counterparts and this can negate any increases in growth. Yet, if invasive genotypes are superior, then invasive genotypes may still experience greater growth even if they are preferentially preyed upon. Documenting higher growth following herbivore damage is difficult in a common garden because of variability in herbivore prevalence and local conditions. To document such a pattern requires simulating herbivore damage experimentally [21] [22].
In this paper, we examine the performance of invasive genotypes when compared to that of native genotypes in the native's own range in the invasive wetland plant, Phalaris arundinacea. The invasive grass Phalaris arundinacea is a good model system to address the issues of the emergence of novel and superior invasive genotypes [23]. Invasive genotypes of Phalaris arundinacea have been shown to be the product of multiple introductions and subsequent admixture [8] and in common greenhouse conditions, invasive genotypes were shown to have faster growth rate, to be taller, have more tillers and greater final biomass [8].
Additionally, invasive genotypes were also shown to have a smaller genome size [24] and had no consistent differences in genetic architecture [25]. However, in order to test whether these invasive genotypes are superior requires that native and invasive genotypes be compared in the native range to determine if the invasive genotypes created in the new range are superior (in survival, growth and/or reproduction) to native genotypes even under conditions in which native genotypes have evolved [26] [27] [28]. In this way, we can determine if the superior performance of invasive genotypes results in overall greater performance under conditions in which the native genotypes should be favored. In addition, by planting genotypes in a field common garden, we can assess in situ herbivore damage to determine if invasive genotypes suffer more herbivore damage than native ones [29] [30]. However, greater herbivore damage may not necessarily result in reduced growth and/or biomass as greater herbivore damage may be compensated for by a faster growth rate [21]. We test this idea in an experimental common garden in which the same genotypes are subject to biweekly biomass removal to simulate grazing by large herbivores. Finally, by planting genotypes collected from northern and southern populations into both northern and southern gardens, we can determine if differences in performance are due to genotypes having adapted to a similar climate rather than native and invasive differences.

Methods
Phalaris arundinacea is a C-3 perennial grass, native to wetlands and wet meadow habitats in Europe and Asia [23]. Considered highly invasive in the midwest and eastern United States, it was originally introduced repeatedly as a wet forage grass and has also been included as seed in conservation mixtures used for soil stabilization [23]. In Vermont, Phalaris has been used as a forage crop in wet pastures for over 50 years but in North Carolina it has been primarily introduced in conservation mix for ditch stabilization (Molofsky personal communication).  (Table 1). Northern and southern populations were chosen to sample different climatic regions of Phalaris arundinacea's range [8]. Thus, our sampling protocol allowed for a comparison of how genotypes from the native (CZ, FR) and invasive (VT, NC) range respond when planted outside their home climate (i.e. NC, FR in Trebon, Czech Republic and VT, CZ in Moussac, France).
Genotypes were identified through allozyme analysis [8] and a subset were selected for experimentation. We selected 18 invasive genotypes (9 VT and 9 NC) and 18 native genotypes (12 CZ and 6 FR). Fewer French genotypes were chosen because one French population contained only hexaploids, while all other populations contained tetraploids; we thus, eliminated hexaploid individuals from our study. Chosen genotypes were transplanted into pots in the University of Vermont greenhouse, where they were maintained, and then sequentially propagated prior to experimentation to remove any maternal environmental effects. Selected tillers of the chosen genotypes were placed into greenhouse flats, placed on their side and allowed to produce replicate tillers. In this way, we created identical copies of all genotypes used in the experiments.    Table 2) but the differences were not related to the native/invasive dichotomy; rather the invasive region had both the highest (NC) and the lowest (VT) survivorship (see Results section).
In each garden, we transplanted nine replicates of each of 36 genotypes (i.e. a total of 324 plants per field site). The nine replicate 5 m × 5 m blocks were specifically chosen to encompass the natural variation present in our field sites. Individual rhizomes were planted into the existing vegetation with minimal disturbance to the surrounding soil and vegetation following the methods outlined We conducted an experiment to determine genotypic and regional differences

Statistical Analyses
We performed separate analyses on each of our gardens to compare the performance of individuals within the garden but not between gardens. To assess survival differences among individuals, we used an ordinal logistic regression [32].
The effects in the model included block, range, region nested within range, pop-ulation nested within region and range, and genotype nested within population, region and range.
In addition to survival, we assessed overall plant performance for three plant traits (tiller number (log transformed), leaf number and stem height (log-transformed). We analyzed leaf number and stem height separately because they were weakly correlated (less than 0.49) and each trait may represent a different growth strategy. In Phalaris arundinacea, flowering only occurs once a plant reaches a threshold height (~30 cm Collins pers. observation). Thus, height represents an important feature of the plants that can influence competitive ability (for light) and reproductive ability.
We performed a mixed ANOVA with range and region within range as the main effects and block, population nested within range and genotypes nested within population region, range as random factors.
We also recorded the presence (1) or absence of herbivory on each plant in the Moussac garden. Data on the presence/absence of herbivory was analyzed using a Chi-Square test.
In the simulated grazing experiment, we analyzed how each plant recovered following biomass removal. Total accumulated removed biomass was analyzed using a general linear model with range and region within range as main effects.
The biomass data was square root transformed to achieve normality.

Results
The  Figure   1); however, despite the overall low survival, we still found significant differences in regional survival ( population and genotype were all highly significant (Table 2). However, average survivorship did differ between the northern and southern genotypes and cut across invasive and native and regional differences (Figure 1). The CZ and VT genotypes had lower survivorship (58% and 55% respectively) than the FR and NC genotypes (78% and 80% respectively).
In both the Trebon and Moussac gardens, tiller number was low (approximately 3 and 2, in Trebon and Moussac, respectively). The small number of tillers per plant made it unlikely that our main effects of range or region within range were significant and in fact, neither of our main effects were significant in either garden (Table 3(a) and Table 3(b)). However, in Trebon, the random  For the trait of plant height, in the Trebon garden, differences were not explained by differences in either range or region within range (Table 3( Our experimental grazing experiment provides evidence for how plants recover following major herbivore damage such as grazing. We predicted that invasive genotypes would recover faster than native genotypes and, in fact, we did find that invasive genotypes regrew more biomass following biomass removal than native genotypes (SS 1,102 = 96.34, F = 4.78, p = 0.0312). Although regional differences were not significant in the full model, multiple comparison tests indicate that NC genotype produced the most cumulative total biomass and FR genotypes the least (Figure 3).

Discussion
Introduced individuals can become invasive in a new range through a process of  selection for greater growth and higher fitness [1]. This selection process can occur through genotypic sorting post introduction or evolution in situ may select for traits that can enhance performance, promoting invasiveness in the new range [5]. The introduction of new genetic variants into a new range can also allow genomes that have never been in contact with each other to hybridize, potentially resulting in new combinations with different traits that have higher fitness than the original native counterparts [11] [14]. Introduced species that have been introduced multiple times for economic or conservation uses such as Phalaris arundinacea, may have already been selected prior to introduction for traits that confer advantages such as higher growth rate, greater biomass production, and greater ability to withstand herbivore damage [8] [23].
To understand whether introduced individuals have acquired traits that make them more invasive, one needs to compare native and invasive individuals when grown together under environmental conditions for which the natives has been selected but the introduced individuals have not. For the invasive and native genotypes studied here, we find only minimal differences between native and invasive genotypes; the only significant native/invasive difference was leaf number in one garden. Thus, for the individuals in this study, invasive individuals were not superior to their native counterparts. In fact, population differences may predominate over any native/invasive difference. In a common garden study on invasive Solidago canadensis in Europe, van Kluenen and Schmidt 2003 [33] found no consistent native/invasive differences among the 9 invasive European and 10 native United States populations when planted into a common garden in the invasive range. In our study, genotypes from NC outperformed other genotypes but these differences may not necessarily be related to their invasive status. Invasive VT genotypes enjoyed no such advantages, even when grown in similar climatic conditions such as Trebon, Czech Republic.
Invasive phenotypes may only be expressed in the correct environmental context [34]. For the same genotypes used in this study, Molofsky and Collins 2015 [18] found that invasive genotypes only outperformed native genotypes under ideal growing conditions but not under more stressful growing conditions. In this study, differences in survival among genotypes occurred but were not linked to invasion status. Rather in the northern garden in Trebon, the invasive genotypes had both the highest (NC genotypes) and the lowest (VT) survivorship. In the southern garden in Moussac, the southern genotypes (NC and FR) had higher survival than the northern genotypes (VT and CZ) suggesting that Phalaris genotypes had adapted to higher temperatures. Heat tolerance in this cool season grass is likely to be a stronger selective force than cold tolerance since plants can become dormant during the winter months.
Several plant traits such as faster growth rate, taller and greater leaf production are all assumed to be higher in invasive individuals [35]. In the common gardens here, no such trend was found. Plant height was primarily environmentally determined in Trebon (44% of the variation) and only correlated with plant identity in Moussac (22% of the variation). Tiller number, a measure of vegetative spread in a clonal plant was also not correlated with native/invasive dichotomy and was primarily environmentally determined. For leaf number, we did find a native/invasive dichotomy, but the differences were only present in one garden and too slight to be of biological significance.
Invasive populations have been postulated to have greater growth at the cost of reduction in defense [19] [36] [37]. If invasive individuals have reallocated resources away from defense and into growth, then we would expect that inva-sive individuals would be preferentially attacked in the native range. In common garden studies of Silene latifolia in the native European range, Wolfe et al. 2004 [37] showed that although invasive individuals had greater herbivore and fungal damage than natives in the native European garden, the invasive populations still outperformed native populations. Similarly, in a study on the invasive liana kudzu (Peuraria montana var. lobata) where invasive and native populations were planted into the native range, invasive individuals suffered higher herbivore damage than the native individuals, but produced greater biomass [30].
Here we find a similar result; although NC genotypes are the most heavily damaged, they survive as well as the native locally adapted FR genotypes. A companion experiment documented the response of native and invasive genotypes to experimental grazing and found that NC genotypes produced the greatest biomass and FR genotypes produced the least. Brodersen et al. 2008 [38] showed that FR genotypes have higher photosynthetic rates but slower growth rates in a common greenhouse study suggesting that increased carbon gained through photosynthesis is being allocated to defense.
Introduction history can have a large influence on the subsequent success of an introduced species [11]. Phalaris has been used as a forage crop in wet pastures for over 50 years in Vermont but in North Carolina it has been primarily introduced in conservation mix for ditch stabilization (Molofsky pers com). We therefore expected that the genotypes used for feed to have greater growth rates and greater ability to recover following grazing. Surprisingly, we find that NC genotypes have higher growth rates following grazing than genotypes from VT or the native range. However, the superior performance of the NC genotypes may be somewhat arbitrary and not related to their invasive status.
Our results have implications for the continued spread of invasive plant species. Firstly, if invasive plant species have wider climatic tolerances then native species than invasive species may spread at the expense of native ones. In our study, only the NC genotypes had greater climatic tolerances as evidence by their relatively high survival in both gardens as compared to the VT genotypes. Secondly, although biological control can be an important component to controlling invasive plant species, our results show that some invasive genotypes may be able to compensate with higher growth rates in the face of even extreme grazing pressure and may be able to maintain their fitness in the field even when preferentially preyed upon. Thirdly, an important conclusion of this research is the wide range of responses from presumably invasive genotypes when planted into field settings. Thus, studies that look for traits associated with invasiveness may find them in "ideal" settings but these same traits may not be expressed under more stressful field conditions. Finally, our study emphasizes that not all introduced populations are equally invasive and that it is important to examine multiple populations of introduced genotypes before drawing conclusions about whether invasive individuals outperform native ones [39]. Therefore, any conclusions about invasive behavior must be experimentally verified and tested in a range of habitats and for multiple populations.