Introduction {#s1}
============

Many exotic organisms become much more abundant and have greater impact on other species in their non-native than in their native ranges [@pone.0082281-Callaway1], [@pone.0082281-Vil1]. Many of these invasive species have been shown to evolve different trait expression in their non-native ranges [@pone.0082281-Ellstrand1]--[@pone.0082281-Ridenour1]. These changes can be due to adaptation, genetic drift, hybridization and/or founder effects [@pone.0082281-Bossdorf1] and have the potential to contribute substantially to invasion. However, whether these changes cause invasive success is speculative. An opportunity to explore the causal link between morphological changes and invasion is to compare shifts expressed by exotic invasive species to those of exotic congeners that naturalize in their new habitats without becoming unusually abundant or having strong impacts [@pone.0082281-Williamson1], [@pone.0082281-Hierro1].

Comparing the traits of exotic species that differ in their ability to invade may help to understand the mechanisms that promote invasion [@pone.0082281-Nijs1]--[@pone.0082281-Milbau1]; but there have been fewer studies focusing on differences among exotic invasive and exotic non-invasive species [@pone.0082281-Thbaud1]--[@pone.0082281-Pyeck1]. Thus the combination of both approaches: 1) the study of exotic invasive and exotic non-invasive species, and 2) studying them in both their native and non-native ranges, has a great deal of potential to shed light on traits that might be important for invasive success and, perhaps more importantly, which adaptive trait shifts between native and non-native ranges may contribute the most to an exotic species evolving in a way to become more invasive [@pone.0082281-VanKleunen1].

A substantial body of literature shows a strong and general tendency for plants from populations in their non-native ranges to increase in size, germination rate, and reproductive output when compared with their native ranges [@pone.0082281-Maron1]--[@pone.0082281-Bossdorf1], [@pone.0082281-Blossey1]--[@pone.0082281-Jakobs1]. In turn, a less common response is the loss of herbivore defensive capacity by plants in non-native ranges [@pone.0082281-Daehler1]--[@pone.0082281-Maron2], but when this occurs it is interpreted in the context of tradeoffs (the hypothesis of evolution of increased competitive ability, EICA) [@pone.0082281-Blossey1]. These evolutionary tradeoffs provide a major hypothesis for how exotic species might transform into invaders; however, to our knowledge there have been only a few works comparing inter-regional trait shifts among invasive and non-invasive congeners [@pone.0082281-Graebner1], [@pone.0082281-Montesinos1].

We compared trait shifts for a set of three very closely related species: *Centaurea solstitialis*, *C. calcitrapa*, and *C. sulphurea*. All three species have highly overlapping ranges both within their native range of Spain and within their non-native range of California. To explore trait changes for these three congeners, we grew plants in common greenhouse conditions and asked the following questions: 1) is the invasive *C. solstitialis* inherently larger, more fecund, or better defended than the two non-invasive congeners, 2) have important trait shifts between the native and non-native ranges occurred only for the invasive *C. solstitialis*?

Materials and Methods {#s2}
=====================

*Centaurea solstitialis*, *C. calcitrapa*, and *C. sulphurea*, are closely related within the *Jacea* group of the *Centaurea* phylogeny [@pone.0082281-GarciaJacas1]. These three species have overlapping distributions in their native ranges in Spain and also in their non-native ranges in California. *Centaurea solstitialis* has been introduced into California since at least 1824 [@pone.0082281-Maddox1]. Although Southern Europe is generally accepted as *C. solstitialis* native region, Prodan \[as cited in 31\] discussed that East Mediterranean and Caucasus are in fact the original source, while Western Mediterranean areas (like Spain) could have been colonized later. *C. calcitrapa\'s* introduction in California is thought to be about 1896 [@pone.0082281-Robbins1], [@pone.0082281-Pitcairn1] and for *C. sulphurea* at least 1923 [@pone.0082281-Muth1], [@pone.0082281-Barbe1]. While *C. solstitialis* is likely to have experienced multiple introductions into California [@pone.0082281-DiTomaso1]--[@pone.0082281-Widmer1] and recent genetic analyses [@pone.0082281-Dlugosch1] show that, when compared with other regions in the native range (Hungary, Romania, Turkey and Georgia), Spanish populations have a genetic structure which is most similar to that of California and Argentina, thus confirming previous papers pointing to the importance of Spanish origin introductions on Californian plants [@pone.0082281-Gerlach1], [@pone.0082281-Hierro2]. All populations studied occurred in or nearby roadsides in highly disturbed areas subjected to Mediterranean climate. Two of these neo-allopatric species, *C. solstitialis* and *C. calcitrap*a, occur over broad native (Southern Europe) and nonnative ranges (Americas, Australia), whereas *C. sulphurea* has a highly restricted native range in Spain and Morocco, and occurs as only a few populations in California [@pone.0082281-GerlachJr1], [@pone.0082281-Montesinos1]. *Centaurea solstitialis* is much less common in Spain than *C. calcitrapa* [@pone.0082281-GarciaJacas1] but *C. solstitialis* has become an aggressive invader in California, while *C. calcitrapa* has not. Based on herbarium records and our observations, *Centaurea sulphurea* is not common in either range ([www.gbif.org](http://www.gbif.org)). All three species inhabit the same ruderal habitats, are winter annuals (although *C. calcitrapa* can occasionally be bi-annual), form basal rosettes, and develop single bolting flower stems from the rosette. All three species also form large spines on their capitula, providing a common trait for which to compare allocation to defense.

In the summer of 2009, we collected seeds from fifteen different individuals from each of 45 different populations across the distributional range of the three species in Spain and California, USA (hereafter, "regions"; see [Table S1](#pone.0082281.s001){ref-type="supplementary-material"}). We sampled eight *C. solstitialis* populations from Spain and 11 populations from California, 10 *C. calcitrapa* populations from Spain and nine from California, and four *C. sulphurea* populations from Spain and three from California. No permissions were required for seed collection. The study species are considered weeds both in the native and non-native ranges and seed sampling was done on roadsides. We confirm that the field studies did not involve endangered or protected species. In January 2010, three seeds, randomly selected from each individual mother plant regardless of their achene type (pappus/non-pappus), were sown in each 2.2 L pot in a 50∶50 mix of 20--30 grit sand and local soil from natural grasslands near Missoula, Montana (total N = 675 plants). After germination, and before seedlings could experience any competitive effect on each other, exceeding seedlings were manually removed so that only one plant remained in each pot. Germination rates were sufficient to reach the intended number of replicates. All plants survived until the end of the experiment. A small proportion (3%) of *C. calcitrapa* plants did not reach reproductive state (it is occasionally biennial); such plants were evenly distributed among experimental groups; data from these plants was not used for the analyses. Plants were randomly mixed in a common garden greenhouse experiment with a temperature range of 10--35°C, watered every 1--2 days, and fertilized biweekly with 100 mL of 1.16 g L^−1^ Scotts Miracle-Gro (15∶30∶15 + micronutrients). Plants were grown until they flowered in a pollinator-excluded greenhouse, and we measured several variables during the germination and growth period: germination rate, rosette diameter, number of capitula, and spine length. Rosette diameter was measured, in mm, 30 and 90 days after sowing in order to calculate relative growth rates. We also manually cross pollinated plants by rubbing receptive capitula from two different individuals from the same natural population with each other. Since *C. solstitialis* is self-incompatible [@pone.0082281-Montesinos1], crosses were randomly made between individuals within populations in order to obtain a seed-set from one flower per individual. Plants were harvested after capitula maturation (July, 2010), which was similar for all three study species. Harvested plants were dried for 48h at 70°C and weighed.

Data were analyzed with the statistical software R 2.15.2 [@pone.0082281-R1] by means of linear mixed-effects models after Laird and Ware [@pone.0082281-Laird1] but allowing for nested random effects of population. Species and regions (Spain or California) were fixed factors for each individual test, and the interaction between the two fixed factors was also studied. Tukey post-hoc tests were used by using R package "multcomp" and accounting with population as a nested random factor. Variables were transformed for normality when necessary.

Results {#s3}
=======

For all populations and ranges combined, individuals of the non-invasive *C. calcitrapa* produced greater total mass than either of the other two species, *C. sosltitialis* was intermediate in total mass, and *C. sulphurea* was the smallest ([Fig. 1A](#pone-0082281-g001){ref-type="fig"}; *F~species~* = 15.46; df = 2,41 *P*\<0.001; [Tables S2](#pone.0082281.s002){ref-type="supplementary-material"}, [S3](#pone.0082281.s003){ref-type="supplementary-material"}). The rosette relative growth rates (RGR) of *C. solstitialis* and *C. calcitrapa* were greater than those of *C. sulphurea* ([Fig. 1B](#pone-0082281-g001){ref-type="fig"}; *F~species~* = 43.32; df = 2,39 *P*\<0.001).

![Total plant biomass (g; mean ± SE) (A); relative growth rate of rosette diameter (mm mm^−1^ day^−1^; mean ± SE) (B) for each species from each region.\
Different letters indicate statistically significant differences (*P*≤0.05).](pone.0082281.g001){#pone-0082281-g001}

The non-invasive *C. calcitrapa* produced the highest number of capitula per plant, followed by *C. solstitialis* and then by *C. sulphurea* ([Fig. 2A](#pone-0082281-g002){ref-type="fig"}; *F~species~* = 43.84; df = 2,42; *P\<*0.001; Tukey\'s post hoc *P\<*0.001 for all species pairs). The number of seeds per capitulum was higher for *C. sulphurea* while *C.calcitrapa* and *C. solstitialis* did not differ ([Fig. 2B](#pone-0082281-g002){ref-type="fig"}; *F~species~* = 78.14; df = 2,3212; *P\<*0.001; [Table 1](#pone-0082281-t001){ref-type="table"}). There were no significant differences among species for seed germination rates ([Fig. 2C](#pone-0082281-g002){ref-type="fig"}; *F* ~species~ = 1.85; df = 2,45; *P* = 0.168). *Centaurea solstitialis* seed-set is lower than previously reported [@pone.0082281-Maddox2], which could be due to natural variability within the species, or to a lower pollination efficiency of our manual treatment than that of natural insect pollinators. In any case, and since our tests compare data obtained from identical treatments, our results are informative at the comparative level.

![Number of capitula per plant (mean ± SE) (A); Number of seeds per capitulum (mean ± SE) (B); Germination percentage (mean ± SE) (C); Spine length (mm; mean ± SE) (D) for each species and region.\
Different letters indicate statistically significant differences (*P*≤0.05).](pone.0082281.g002){#pone-0082281-g002}

10.1371/journal.pone.0082281.t001

###### Percent increase in the mean trait value for plants from the non-native range, compared to mean of the trait from the native range.

![](pone.0082281.t001){#pone-0082281-t001-1}

  Trait                 *C. sulphurea*   *C. calcitrapa*   *C. solstitialis*
  -------------------- ---------------- ----------------- -------------------
  Biomass                 **8%^a^**        **19%^b^**           **10%**
  Rosette RGR               26%^a^         **36%^bc^**          −6%^c^
  Capitula per plant        1%^a^             3%^b^             −15%^c^
  Seeds per capitula       −20%^a^          −16% ^b^            22%^b^
  Germination rate       **123%^a^**       **24%^a^**         **63%^a^**
  Spine length              −2%^a^            2%^b^           **21% ^c^**

*P*≤0.05) and different letters indicate significant differences among species (*P*≤0.05). Numbers in bold indicate significance levels for inter-regional differences (

Across all populations and ranges *C. calcitrapa* had longer spines than either of the two congeners, *C*. sosltitialis was intermediate, and *C. sulphurea* produced the smallest spines ([Fig. 2D](#pone-0082281-g002){ref-type="fig"}; *F* ~species~ = 278.07; df = 2,39; *P*\<0.001; Tukey\'s post hoc *P\<*0.001 for all species pairs).

For all three species, individual plants grown from seed collected in the non-native ranges produced more total biomass than plants grown from seed collected in the native region (*F~region~* = 4.07; df = 1,41; *P* = 0.051), but the largest shift in this trait was expressed by the non-invasive *C. calcitrapa* ([Fig. 1A](#pone-0082281-g001){ref-type="fig"}, [Table 1](#pone-0082281-t001){ref-type="table"}). *Centaurea calcitrapa* was the only species to show an inter-regional shift in rosette RGR, with plants from California growing faster than plants from Spain ([Fig. 1B](#pone-0082281-g001){ref-type="fig"}, [Table 1](#pone-0082281-t001){ref-type="table"}; *F~region~* = 3.91; df = 1,39; *P* = 0.055, *F~species\ x\ region~* = 3.79; df = 2,39; *P* = 0.031; Tukey\'s post hoc *P* = 0.007). We found no differences in the number of capitula per plant between native and non-native regions for any of the three species ([Fig. 2A](#pone-0082281-g002){ref-type="fig"}; *F~region~* = 0.49; df = 1,41; *P* = 0.487). *Centaurea sulphurea* was the only species that showed a shift in the number of seeds per capitulum between regions, with Californian plants producing fewer seeds per capitulum than Spanish plants ([Fig. 2B](#pone-0082281-g002){ref-type="fig"}, [Table 1](#pone-0082281-t001){ref-type="table"}; *F~region~* = 78.13; df = 2,3212; *P* = 0.005; *F~species\ x\ region~* = 0.057). Seeds produced by plants of all species from the non-native range had higher germination rates, with the average increase across the three species of 61% when compared to the native region ([Fig. 2C](#pone-0082281-g002){ref-type="fig"}, [Table 1](#pone-0082281-t001){ref-type="table"}; *F~region~* = 22.90; df = 1,45; *P*\<0.001). *Centaurea sulphurea* experienced the highest change in germination rate between ranges followed by *C. solstitialis* (*P* = 0.006; *P* = 0.002). *Centaurea solstitialis* plants from California were the only species that demonstrated an increase in spine length in the non-native range, with Californian plants producing spines that were 21% longer than their conspecifics from the native range ([Fig. 2D](#pone-0082281-g002){ref-type="fig"}; *F~region~* = 8.86; df = 1,39; *P = *0.005; *F~species\ x\ region~* = 4.39; df = 2,39; *P* = 0.019).

Discussion {#s4}
==========

In general, traits related to growth, fecundity, and defense that we measured were not inherently greater in value for the exotic invasive *C. solstitialis* than for non-invasive congeners. More importantly, the length of spines on the calyces was the only trait for which *C. solstitialis* demonstrated a greater increase in the non-native range than the other congeners. Thus the most important interpretation of our results is that, for these species and traits, we found no evidence that the invasive species exhibited any stronger shifts in key traits than very similar but non-invasive exotic congeners.

Our results are consistent with a large body of literature demonstrating substantial differences in morphology or size for plants from populations in the native and non-native ranges. For example, an increase in size for plants in the non-native ranges has been widely documented and overall, invasive species tend to have higher values also for leaf-area allocation, shoot allocation, growth rate, physiology, and fitness [@pone.0082281-Ridenour1], [@pone.0082281-Bossdorf1], [@pone.0082281-VanKleunen1], [@pone.0082281-Blossey1], [@pone.0082281-Siemann1], [@pone.0082281-Buckley1], [@pone.0082281-Maddox2]--[@pone.0082281-Elton1]. Seed and seedling size has been demonstrated to differ between native and non-native ranges of *C. solstitialis*, but not for the non-invasive *C. calcitrapa* or *C. sulphurea* [@pone.0082281-Graebner1]. However seedling size only differed when *C. solstitialis* was grown in competition with the European native grass *Bromus hordeaceus*, when Californian individuals grew more than their Spanish counterparts, but not when in competition with the American native *Poa secunda*, thus showing that individuals from the non-native range present increased competitive ability under certain circumstances [@pone.0082281-Graebner1]. Germination rates have been tested much less, but Hierro *et al.* [@pone.0082281-Hierro2] found higher germination rates for seeds produced on *C. solstitialis* plants from the non-native range of California than plants from the native range of Turkey in a common garden; Ridenour *et al.* [@pone.0082281-Ridenour1] found that the mean germination rate of *C. stoebe* from non-native North American populations was 81% higher than that of native European populations. Kudoh *et al.* [@pone.0082281-Kudoh1] observed patterns consistent with adaptation to fall germination for invasive strains of *Cardamine hirsuta*. These results indicate that exotic plants may experience strong selective pressure in their non-native ranges and respond to this pressure rapidly; a rapid accumulation of ecological adaptations which has been found to lead to incipient degrees of reproductive isolation between native and non-native ranges of introduced species [@pone.0082281-Montesinos1]. However, our results for these *Centaurea* congeners suggest caution in interpreting these apparent adaptive responses as *drivers* of invasion. Both exotic invasive and non-invasive congeners exhibited substantial trait shifts, potentially very important for local adaptation but, in view of our results, not necessarily key to the dramatic biogeographic shifts in abundance and impact manifest by *C. solstitialis* but not *C. calcitrapa* or *C. sulphurea*.

Invasive species often grow much faster than the native species they exclude [@pone.0082281-Graebner1], [@pone.0082281-Pyeck1], [@pone.0082281-Cronk1]--[@pone.0082281-Burns1]. Graebner *et al*. [@pone.0082281-Graebner1] found that, when in competition with native grass species, the relative growth rate (RGR) of *C. solstitialis* biomass, from seed to seedling, was greater than that of *C. calcitrapa* or *C. sulphurea*, but there was no difference in this trait between the native and non-native ranges for any of the three congeners. Here, we also found that the overall RGR of *C. solstitialis* rosettes was 11% greater than that of *C. calcitrapa* and 42% greater than that of *C. sulphurea*, but that only *C. calcitrapa* showed evidence for evolving higher RGR in the non-native range. Perhaps these inherently rapid growth rates were critical for the early success of *C. solstitialis* relative to its congeners.

*Centaurea solstitialis* was intermediate in spine length, but the only species to produce significantly larger spines in California than in Spain. Increased spine length could simply be correlated with an increase in total biomass, but spine length augmented 21% in non-native populations; whereas biomass increased 10%. Given that *C. solstitialis* in California often occurs with domestic livestock, perhaps this strong generalist herbivore pressure is a more important selective force in the non-native range of California. Indeed, cows, sheep and goats have been described as effective grazers of *C. solstitialis* and have been effectively used as a control agent in California [@pone.0082281-Thomsen1], [@pone.0082281-Ditomaso1].

Time since introduction is an important variable when studying biological invasions, and time lags of more than 50 years between introduction and invasion are common [@pone.0082281-Daehler2], [@pone.0082281-Hobbs1]. Several multispecies studies found a positive relationship between the capacity to spread of exotic plants and the time since their introduction in the non-native areas [@pone.0082281-Pyek1]--[@pone.0082281-Cadotte1]. The cause of the lags is partly inherent to the dynamics of population growth and range expansion [@pone.0082281-Crooks1], [@pone.0082281-Taylor1] but frequently it is also related to the need to accumulate sufficient genetic diversity via repeated introductions; or to develop key adaptations to the non-native range which are crucial for invasive success [@pone.0082281-VanKleunen1], [@pone.0082281-Daehler2], [@pone.0082281-Crooks1]. This has generated a phenomena called "invasion debt", a hypothesis suggesting that past human activities involving the introduction of exotic species will have a future impact on ecosystems after a time lag since date of introduction [@pone.0082281-Essl1]--[@pone.0082281-Seabloom1]. The three species considered in this study differ in time since introduction, from 189 years for the invasive *C. solstitialis*, to 117 years and 90 years for the non-invasive *C. calcitrapa* and *C. sulphurea*, respectively. Interestingly, our results show that these three closely related species show a similar amount of trait-shifts, suggesting that trait-shifts can happen very rapidly after introduction although, visibly, they do not necessarily lead to invasive success at the same speed.

A number of studies have attempted to elucidate the key traits that drive invasions [@pone.0082281-Pyeck1], [@pone.0082281-Baker1]--[@pone.0082281-Rejmnek1]. These efforts however, have led to the conclusion that the important traits may differ among species, and even among sites for the same invasive species. For instance, Hierro et al. [@pone.0082281-Hierro2] studied germination patterns for *C. sosltitialis*, and found differences not only between native and non-native ranges, but also between populations from the non-native ranges of Argentina and California. Previous studies [@pone.0082281-Graebner1] with the species system used in the present work found competitive advantages related to seed and seedling size of *C. solstitialis* under certain circumstances; however, both studies considered different variables under different conditions, and the absence of significance for some of the studied variables does not preclude the possibility of some key variable not being considered.

Multispecies comparisons have been fruitful, but tend to compare invasive species with natives within the invaded region [@pone.0082281-Pyeck1], [@pone.0082281-Peuelas1], [@pone.0082281-Knapp1]. The trait differences detected so far between native and non-native species could be due in part more to their different geographic origins, and the likelihood of being transport than to actually advantageous traits [@pone.0082281-Knapp1]. Traits have been compared among several naturalized invasive and non-invasive species. For instance Muth and Pigliucci [@pone.0082281-Muth1] compared *Crepis* and *Centaurea* in their non-native range, and found that invasiveness corresponded well with species-specific trait interactions and introduction histories. This approach does not explain if the differentiating traits are acquired after introduction or were present in the native range, thus studies involving individuals from both the native and the non-native ranges are necessary [@pone.0082281-VanKleunen1]. Our results do not preclude that a unique and synergistic combination of local adaptations drive invasion success. However, since we found that non-invasive species show some of the same trait shifts between native and non-native regions as an invasive species, our results suggest caution in assuming that shifts in traits thought to be important to invasion actually cause the invasive success of a particular species.

Supporting Information {#s5}
======================

###### 

**Location of** ***Centaurea*** **populations for each of the species from each of the studied regions.** Latitude and longitude coordinates are datum WGS84.

(DOCX)

###### 

Click here for additional data file.

###### 

**Statistical results for the Linear Mixed Models for each studied trait and factor.** Significant differences are in bold.

(DOCX)

###### 

Click here for additional data file.

###### 

**Tukey\'s post-hoc p-values following Linear Mixed Models for species to species comparisons for each trait.** Significant differences are in bold, n/a: not applicable.

(DOCX)

###### 

Click here for additional data file.

Thanks to Margarita Zorrilla and to everyone at Callaway Lab for their help with lab work.

[^1]: **Competing Interests:**The authors have declared that no competing interests exist.

[^2]: Conceived and designed the experiments: DM DA RMC. Performed the experiments: DM RMC. Analyzed the data: YG DM. Contributed reagents/materials/analysis tools: AD DM RMC. Wrote the paper: YG RMC AD DM.
