Can polyploidy confer invasive plants with a wider climatic tolerance? A test using Solidago canadensis

Abstract Polyploidy can cause variation in plant functional traits and thereby generate individuals that can adapt to fluctuating environments and exploit new environments. However, few empirical studies have tested for an association between ploidy level and climatic tolerance of invasive cytotypes relative to conspecific native‐range cytotypes. Here, we used an invasive plant Solidago canadensis to test whether invasive populations had a higher proportion of polyploids, greater height and stem‐base diameter, and occupied a wider range of climatic conditions than conspecific native‐range populations. We also tested whether the invasive populations had overcome genetic founder effects. We sampled a total of 80 populations in parts of the invaded range in China and native range in North America for in situ measurements of plant height and stem‐base diameter in the field and for population genetic and cytotype analyses. To examine climatic correlates, we augmented our field‐sampled data with occurrence records obtained from Global Biodiversity Information Facility. All, except one, of the populations that we sampled in China occurred in a humid subtropical climate. In contrast, the North American populations occurred in humid continental, humid subtropical, and semi‐arid climatic zones. All populations of S. canadensis in China were purely hexaploid, while the North American populations were diploid, tetraploid, and hexaploid. The invasive hexaploids were significantly taller and had a larger stem‐base diameter than native hexaploids. Native hexaploids were significantly taller and had larger stem‐base diameter than native diploids. Climatic correlate assessment found that invasive and native populations occupied different climatic envelopes, with invasive populations occurring in warmer and less seasonal climates than native populations. However, there was no significant correlation between ploidy level and climatic envelope of S. canadensis. Molecular phylogeography data suggest reduced genetic founder effects in the invaded range. Overall, these results suggest that polyploidy does not influence S. canadensis climatic tolerance.


| INTRODUC TI ON
Invasions by exotic plant species pose substantial threats to native biodiversity and alter ecosystem processes across a broad range of environments (McGeoch et al., 2010;Pyšek et al., 2012). To improve effectiveness of invasive plant management, it is necessary to understand the mechanisms that underpin invasive spread of exotic plants (Dawson, Burslem, & Hulme, 2009). Species with high variation in functional traits (morphological, physiological, and phenological) can adapt to different stress levels and increase the uptake of limiting resources across a broad range of climatic conditions (Nicotra et al., 2010). Thus, plants that possess higher mean values of, and variation in, functional traits in the invaded ranges than in their native ranges may have an enhanced capacity to colonize a wide range of climatic conditions Richardson & Pyšek, 2006). Therefore, an understanding of the factors that determine variation in functional trait values in invasive plants is key to predicting invasion risks (Leffler, James, Monaco, & Sheley, 2014;te Beest, Esler, & Richardson, 2015).
It also remains unresolved whether a pronounced shift in cytotype distribution that has been reported in several invasive plant species is due to stochastic genetic founder events or adaptive evolutionary processes. Population genetic theory predicts that invasive species are likely to experience genetic founder effects, which may reduce their ability to adapt to novel environmental conditions (Dlugosch & Parker, 2008). Founder effects can cause a reduction in genetic diversity because the introduced populations are often a subset of the original, larger populations (Dlugosch & Parker, 2008).
Nevertheless, studies have found evidence for local adaptation in several invasive species (Oduor, Leimu, & Kleunen, 2016), which implies that the species have overcome genetic founder effects.
Invasive species may overcome (or at least minimize) founder effects through multiple introduction events (i.e., introductions of diverse genetic lineages from different source populations in the native range) or a single introduction event from a native-range source population having high genetic variation (Genton, Shykoff, & Giraud, 2005;Le Roux et al., 2011;Meyerson & Cronin, 2013;. However, few studies have used cytobiogeographical and molecular phylogeographical approaches to determine whether invasive cytotypes represent a full complement of genetic diversity that occurs in the native range. Here, we used the invasive plant S. canadensis to test a prediction that polyploidy can enable invasive plant species to colonize a wider range of climatic conditions. Specifically, we combined a cytobiogeographical approach, estimates of functional traits and climatic correlates, and molecular phylogeography methods to test whether invasive populations of S. canadensis: 1) had a higher proportion of polyploids, 2) were larger, 3) had a wider climatic tolerance, and 4) had overcome founder effects, relative to native populations of S. canadensis.

| Study species
Solidago canadensis L. (Asteraceae) is a rhizomatous clonal perennial forb native to North America and has invaded China and other regions (Europe, Australia, New Zealand, and Japan) (Abhilasha, Quintana, Vivanco, & Joshi, 2008;Jin, Yuan, Gao, Oduor, & Li, 2020;Xu & Qiang, 2011). Solidago canadensis reproduces both clonally and sexually (Li et al., 2017). The plant is self-incompatible and produces small wind-dispersed seeds for long-distance dispersal (Li et al., 2017). However, the plant can spread clonally at a local scale (Li et al., 2017). In China, S. canadensis was first introduced to Zhejiang, Jiangsu, and Shanghai provinces in the early 1930s as an ornamental garden species (Xu et al., 2012;Zhao et al., 2015). It escaped from cultivation and experienced explosive spread across eastern China after 1980s, and is now distributed widely across eastern China where it abounds in disturbed habitats, including agricultural lands, roadsides, railways, and city suburbs (Dong, Yu, & He, 2015;Zhao et al., 2015). The species has diploid, tetraploid, and hexaploid cytotypes in the native range, but only hexaploids have been found in China (Li et al., 2017). Previous studies (Xu et al., 2012;Zhao et al., 2015) suggest that the Chinese populations were founded through multiple introduction events.

| Sampling and functional trait measurement
To identify occurrence locations of S. canadensis in China and North America, we conducted field surveys in summer and autumn for the period 2013-2016. We sampled a total of 2,051 S. canadensis individuals in 19 native and 61 invasive populations. The number of individuals sampled per population ranged from five to 33. The native range in North America was represented by 577 individuals, while the invaded range in China was represented by 1,474 individuals (Appendix 1). Any two sample populations were at least 5 km apart.
In each population, we sampled at random plants that occurred within a radius of 10-20 m to minimize a potential bias due to environmental heterogeneity within a population (e.g., soil moisture and nutrients). To avoid sampling plants from the same genet within the same radius, we sampled plants that were at least two meters apart from each other. For each S. canadensis individual, we took in situ measurements of stem-base diameter and plant height (in cm). These functional traits are generally linked to ecological strategy axes of plants. As a plant developmental stage (i.e., phenophase) can influence the values of these traits, we took measurements only from mature plants that had set seeds and were therefore presumably not undergoing further development in height and stem diameter. We omitted leaf measurements because an individual S. canadensis plant can have hundreds of leaves (especially in the non-native range), and hence, finding the largest leaf for measurement can be extremely difficult. In the native range, the ecological niches of S. canadensis and those of its congeners S. gigantea and S. altissima overlap and the three species can coexist (Benelli et al., 2019;Weber, 1998).
However, S. canadensis prefers loose and drier soils than the congeners, and hence, it occurs mostly near urban areas, roadsides, and railways, while the congeners occur mainly on riverside and embankments (Benelli et al., 2019). Therefore, within a locality, we sampled only the preferred habitats of S. canadensis to avoid accidental sampling of the congeners. Moreover, we conducted a test of phylogenetic relatedness among the populations that we sampled in the field to determine whether they all belonged to S. canadensis (see the section on population genetics and phylogenetic analyses below). All the sampled populations in China occurred in a humid subtropical climate except one that occurred in a subtropical highland climate.

| Determination of ploidy level
We took fresh leaf samples from the same S. canadensis individuals used for functional trait measurements described above. The leaf materials were held in self-sealing plastic bags containing silica and transported to the laboratory for ploidy level determination. We applied a modified flow cytometry method by Suda and Trávníček (2006)

| Statistical analyses
To test whether S. canadensis plants from the invaded range had significantly greater height and stem-base diameter than S. canadensis plants from the native range, we fitted linear mixed-effect models.
In the models, a trait value was treated as a dependent variable, while S. canadensis range (invaded vs. native) was treated as a fixedeffect independent variable. Population identity of S. canadensis was treated as a random-effect factor and nested within range. To test whether polyploidy in general conferred higher mean trait values, we also fitted linear mixed-effect models to test whether the functional trait values differed between hexaploid and diploid individuals in the native range. Because only hexaploids were found in the invaded range, this latter test was not possible for invasive plants.
In the test, ploidy level (polyploid vs. diploid) was treated as a fixedeffect independent variable, while S. canadensis population identity was treated as a random-effect independent variable. We did additional analyses to test for the potential confounding effect of climate on trait expression by comparing traits of invasive plants with those of native plants from the same climatic zone only. As nearly all invasive plants were sampled in a humid subtropical climate, we used a subset of plants from this climatic zone only. We did separate comparisons for groups of plants that occurred in populations with mixed-ploidy and for groups of plants that occurred in populations with one ploidy type. Moreover, we compared traits of diploid plants with those of hexaploid plants that occurred in mixed-ploidy populations only to exclude the potential confounding effect of pure and mixed-ploidy populations occupying habitats with contrasting ecological conditions. In all the analyses, we fitted the models using maximum likelihood with the lme function in the nlme package (Pinheiro, Bates, DebRoy, & Sarkar, 2007) in r v 3.5.2 (R Core Team, 2018). We used likelihood ratio tests to assess the significance of each factor by comparing a model without the factor with a full model.

| Assessing climatic tolerance
We described the climatic space occupied by S. canadensis with 19 variables with a 30 arc-sec resolution (~1 km) derived from the WorldClim database (Hijmans, Cameron, Parra, Jones, & Jarvis, 2005). After standardization of the climatic variables (mean equal to zero and standard deviation equal to 1), we reduced the dimensionality of the climatic hyperspace to two axes (PC1 and PC2) based on a PCA (Broennimann et al., 2012). In addition to our 80 field-sampled records of S. canadensis in China and North America, we collected information on occurrences of S. canadensis provided by the Global Biodiversity Information Facility (GBIF) database (http://www.gbif.org/) in order to assess more comprehensively the climatic tolerance of the species. We retained only the GBIF occurrence records in the native (United States, Canada, and Mexico) and invaded (Australia, New Zealand, Japan, and China) ranges where the species is most widespread. Occurrences flagged as invalid (e.g., in the oceans) or with doubtful coordinates were removed as well as those separated by less than 30 arc-seconds to avoid pseudoreplication. We then extracted values of PC1 and PC2 for each of the 2,214 retained occurrence locations (i.e., our 80 field-sampled records and 2,134 GBIF records, which included 1,624 records in the native range and 510 records in the invaded range).
To test whether climatic space occupied by S. canadensis differed significantly between the invaded and native ranges, we calculated the kernel-smoothed density of occurrence of the two ranges based on the density of the different combinations of climatic conditions available to each range (Broennimann et al., 2012). The overlap between the climatic envelope of the two ranges was estimated with the Schoener's D metric, which ranges from 0 (no envelope overlap) to 1 (identical envelopes) (Schoener, 1968). We then tested for envelope equivalency and the less restrictive hypothesis of envelope similarity (comparing the observed values of envelope overlap to the 95th percentile density of the simulated values) (Warren, Glor, & Turelli, 2008). To achieve this, the D value between the two ranges was compared to a null distribution of D values computed between simulated envelopes built through randomization procedures (10,000 randomly selected "background" points generated 100 times). This was computed with the r package "ecospat" (Di Cola et al., 2017).
As a D value does not allow a determination of whether an absence of envelope overlap results from the lack of overlap for one or two PCA axes, we used the function "niceOverPlot" implemented in R for easier interpretation. To test whether S. canadensis inhabited a wider range of climatic conditions in its invaded range than in the native range, we performed an F-test of equality of variances (homoscedasticity) along the two PCA axes. We also used the same method to compare the climatic zones occupied by diploid and hexaploid populations of S. canadensis in its native range. Finally, we compared D values and envelope equivalencies and similarities between the following groups of S. canadensis populations: 1) North American versus Chinese populations that we had sampled in the field and with the same ploidy level, 2) North American versus Chinese populations that we had sampled in the field and without controlling for variation in ploidy level, 3) a global data set (i.e., the 80 populations that we had sampled in the field and GBIF records) of North America and Mexico grouped together versus China, and 4) a global data set of North America and Mexico grouped together versus all occurrence records in Australia, New Zealand, Japan, and China. Comparisons 1 and 2 tested the potential effect of variation in ploidy, while comparisons 3 and 4 tested whether our sampling was representative of the invaded and native ranges. For comparisons 3 and 4, it was not possible to control for variation in ploidy level because information on ploidy level was missing in the GBIF database.

| DNA extraction, PCR, and sequencing
To estimate gene flow in S. canadensis from North America to China and phylogenetic relationships among the 80 invasive and native populations that we sampled in the field, we sequenced two chloro- All the 536 sequences were aligned using MEGA v 10.0.5 software (Kumar, Stecher, Li, Knyaz, & Tamura, 2018) and then trimmed to 782 bp, which included insertions and deletions.

| Population genetic diversity and structure analyses
To test whether patterns of genetic diversity differed between S. canadensis populations in the invaded (China) and native (North America) ranges, we used sequences of the chloroplast spacer regions (psbA-trnH and trnL-F) to estimate the number of polymorphic sites, the total number of haplotypes, haplotypic diversity (H d ), and nucleotide diversity (n) using ARLEQUIN v 3.5.2.2 (Excoffier & Lischer, 2010). These diversity indices were computed for each population separately. We then fitted linear mixed-effect models (in r v 3.5.2) to test whether S. canadensis populations from the invaded range differed significantly in mean values of the diversity indices from S. canadensis populations in the native range. In the models, the diversity indices were treated as dependent variables, while S. canadensis range (invaded vs. native) was treated as a fixed-effect independent variable. Population identity of S. canadensis was treated as a random-effect independent variable and nested within range.
To test whether S. canadensis populations in the invaded range had a different genetic structure than S. canadensis populations in the native range, we used the same sequences above to perform hierarchical analysis of molecular variance (AMOVA) (Excoffier, Smouse, & Quattro, 1992). The AMOVA was performed using the global data set to compare genetic structuring between invaded and native ranges, and hierarchically for the invaded and native ranges separately using individual populations. The hierarchical AMOVA divides the total genetic variance into components due to interindividual differences within a population and interpopulation differences within a range. Genetic differentiation between ranges and among populations within ranges was compared by the fixation index (F ST ) and tested by AMOVA in ARLEQUIN v 3.5.2.2 (Excoffier & Lischer, 2010). Significance of genetic differentiation was tested by 1,000 random permutations. The minority cytotype exclusion hypothesis predicts that there could be a restricted gene flow between polyploid and progenitor diploid populations (Levin, 1975), which could influence genetic structure of populations that contain different cytotypes. Therefore, we also performed AMOVA based only on purely hexaploid populations that occurred in the invaded range and mixed-ploidy populations in the native range that contained at least one hexaploid cytotype.

| Testing for population demographic expansion
Changes in demographic history can influence the frequency of alleles, the distribution of mutations, and the coalescent times of gene copies (Zhang, Edwards, Kang, & Fuller, 2012). Hence, we inferred the effects of past demographic expansion on the current genetic variation in invaded and native ranges using Tajima's D and Fu's F S neutrality tests (Fu, 1997;Tajima, 1989) and pairwise mismatch distributions between haplotypes (Rogers & Harpending, 1992). The neutrality tests were performed with DnaSP v. 6 (Rozas et al., 2017), while the mismatch distribution analyses, which test patterns of nucleotide variation against a null model expected under a sudden population expansion, were implemented in ARLEQUIN v 3.5.2.2 (Excoffier & Lischer, 2010) using 10,000 bootstrap replicates. In the neutrality tests, examination of deviation from neutrality was based on 1,000 coalescent simulations. Nonstatistical difference from zero rejects the null hypothesis of neutral evolution. Significant negative values of Tajima's D and Fu's F S indicate an excess of young or rare alleles in the genealogy, which suggest recent population expansion or purifying selection (Fu, 1997;Tajima, 1989). On the other hand, significant positive values indicate processes such as recent population bottlenecks or balancing selection (Fu, 1997;Tajima, 1989). A mismatch distribution analysis calculates various population parameters, such as the raggedness index (r), with significantly ragged populations having stable demographic histories and non-significantly ragged populations having sudden population growth (Rogers & Harpending, 1992). For data distributed according to a sudden population expansion model (i.e., non-significant r indices), the mismatch distribution analysis computes parameters of demographic expansion such as moment estimators of time to the expansion (Tau = 2Tu), effective population size before expansion (Theta0, θ0), and effective population size after expansion (Theta1, θ1) between the observed and expected mismatches (Rogers & Harpending, 1992). Departure from a model of sudden expansion was tested for each population by summing the squared differences (SSD) between observed and estimated mismatch distribution. To further infer the demographic history of S. canadensis, we determined genealogical relationships among haplotypes using the statistical parsimony algorithm implemented in ARLEQUIN v 3.5.2.2 (Excoffier & Lischer, 2010). A starshaped genealogy indicates a lack of (or limited) geographic structure and suggests a signature of rapid population expansion (Slatkin & Hudson, 1991). Genealogical relationships among the haplotypes were inferred using a combined data set of invasive and native haplotypes and separately for the two sets of haplotypes.

| Analysis of phylogenetic relatedness
To test whether all the 80 populations of S. canadensis that we sampled in the invaded and native ranges belonged to the same species (i.e., S. canadensis), we performed analyses of phylogenetic relatedness using sequences of 34 haplotypes that were derived from 536 individual sequences of the two chloroplast spacer regions (psbA-trnH and trnL-F). We applied three methods, namely the Maximum Parsimony, Maximum Likelihood, and Neighbor-Joining.

| Functional traits
Invasive hexaploids were significantly taller and had a greater stembase diameter than native hexaploids (Figure 2a & b). Similar results were found when invasive hexaploids were compared with native hexaploids from the same climatic zone (i.e., humid subtropical) (Appendix 3a & b). Native hexaploids had greater mean values of the two traits than native diploids (Figure 2c & d). Similar results were found when native hexaploids from mixed-ploidy populations were compared with native diploids from mixed-ploidy populations (Appendix 3c & d).

| Climatic envelopes of S. canadensis
The PCA revealed two principal components that cumulatively explained 74% of variation in the climatic data ( nor similar (p > 0.14) ( Table 1). In all cases, climatic envelopes clearly differentiated along the temperature gradient (PC2), with invasive populations preferring warmer and less seasonal climates (Figure 3).
The invasive populations exhibited a contraction in climatic tolerance (i.e., the ratio of the climate inhabited by native populations to that by invasive populations was significantly greater than one) (Table 1). However, the invasive hexaploids and native hexaploids inhabited similar climates (p = 0.92 along PC1 and p = 0.46 along PC2) (Table 1). An analysis of a pooled data set of the 80 populations that we surveyed in China and North America did not find a significant difference in the climatic conditions among the populations (p = 0.74 along PC1 and p = 0.75 along PC2) (Table 1). In addition, diploid and hexaploid native populations occurred in similar climatic conditions (D = 0.78; equivalency test; p = 0.02; Appendix 4).

| Correlations between plant traits and bioclimatic variables
Stem-base diameter had significant positive correlations with five (annual mean temperature, mean diurnal range, precipitation of the F I G U R E 2 Mean (±1SE) height and stem-base diameter of Solidago canadensis plants. Shown are comparisons between invasive hexaploids and native hexaploids (a & b), and between native diploids and native hexaploids (c & d)

| Genetic diversity of S. canadensis in the invaded and native ranges
We found a total of 34 haplotypes from 536 sequences ( (Table 2). However, mean haplotypic and nucleotide diversities were similar between the invaded and native populations (Table 2).

| Genetic structure and phylogenetic relatedness
We detected significant genetic structures at all hierarchical levels (Table 3). There was a low (7.27%) but significant genetic TA B L E 1 Results of climatic envelope analysis for invasive and native populations of Solidago canadensis. The comparison "Invasive versus native hexaploids" compares Chinese hexaploid with North American hexaploid populations; "Among samples" compares variation among our field-sampled populations in China and North America pooled together (80 populations); "China versus Native range" compares occurrences in China with North America and Mexico (our field-sampled sites and GBIF records); "Entire invaded range versus Native range" compares occurrences in entire invaded range in China, Australia, New Zealand, and Japan with North America and Mexico (our field-sampled sites and GBIF records). The last three comparisons did not control for a variation in ploidy level. The column "D" refers to the value of the Schoener's D metric, "Equivalency" to the P value of the equivalency test, "Similarity" to the P value of the similarity test, and "Breadth ratio" to the ratio between the environmental tolerance of native populations to that of invasive populations along the axes PC1 and PC2 (when a significant difference in environmental tolerance was found). For the breadth ratios, n.s. means that the ratio is non-statistically different (p > 0.05) from one  (Table 3). Nonetheless, the invasive populations exhibited higher within-population genetic variation (67.10%) than the native populations (47.68%) ( Table 3).

| D ISCUSS I ON
The present study used S. canadensis to test whether invasive populations had a higher proportion of polyploids, were significantly larger, and occupied a wider range of climatic conditions than conspecific native-range populations. Because biogeographic differences between invasive and native populations in cytotype frequency and functional traits could be due to stochastic founder events, the study also tested whether the invasive populations had overcome (or minimized) founder effects. We found that all invasive S. canadensis   canadensis was introduced to China as an ornamental garden plant (Xu et al., 2012;Zhao et al., 2015), it is possible gardeners preferentially introduced hexaploids that had some unique desirable traits like larger sizes. Third, because invasive plant introductions are F I G U R E 4 Climatic envelopes of Solidago canadensis in its native (orange) and invaded (blue) ranges. The following comparisons were made: (a) invasive hexaploid versus native hexaploid populations; (b) all invasive versus native populations that we sampled in the field without controlling for variation in ploidy; (c) all occurrences in North America and Mexico (our field-sampled sites and GBIF records) versus China; (d) all occurrences in North America and Mexico versus entire invaded range in China, Australia, New Zealand, and Japan. Histograms on top and on the right side of each figure give the standardized density distribution along the first (PC1) and second (PC2) principal components, respectively often stochastic events, strong founder events may have caused a loss of cytotype diversity in the invaded range (Mráz et al., 2014).

| Cytobiogeography and functional trait differences
This may occur particularly when there is a spatial separation of the different cytotypes in the native range, which increases the chance that only a subset of the cytotypic diversity is introduced (Kubatova et al., 2008;Lafuma, Balkwill, Imbert, Verlaque, & Maurice, 2003). However, the second and third scenarios seem less plausible because our molecular phylogeography data (Tables 2-3 and Appendix 8) and previous studies (Xu et al., 2012;Zhao et al., 2015) suggest that the Chinese populations have minimized founder effects through multiple introductions.
The present finding that invasive S. canadensis hexaploids had greater height and stem-base diameter than native hexaploids ( Figure 2) and that the two traits were significantly correlated with various bioclimatic factors negatively and positively (Appendix 5) suggests that divergence in the trait means could be the outcome of differential natural selection imposed by the climatic factors. Evolution of increased sizes has been associated with invasiveness of several plant species (e.g., Oduor et al., 2013;Whitney & Gabler, 2008;Oduor, Stift, & van Kleunen, 2015a;Oduor, van Kleunen, & Stift, 2017 ). Other studies that compared traits of invasive polyploids and conspecific native-range polyploids in common-garden settings produced mixed results. For instance, the invasive Oxalis pes-caprae tetraploids produced more aboveground biomass and offspring bulbs than native O. pes-caprae tetraploids (Tavares et al., 2019). Within S. gigantea, invasive tetraploids produced more biomass than native-range hexaploids (Nagy et al., 2018). In contrast, invasive tetraploid C. stoebe had similar length and number of florets, length of capitula and longest branch, leaf weight and leaf area as conspecific native tetraploids (Mráz, 2011). In Vicia cracca, invasive tetraploids had greater seedling height and aboveground biomass than native tetraploids; however, the two groups of plants produced similar numbers of seeds and branches (Líblová, Eliášová, & Münzbergová, 2017).
Biogeographical comparisons of traits that use small and nonrepresentative subsets of populations may not give a complete picture of the ecological and evolutionary processes that cause trait divergence between invaded and native ranges (Rosche et al., 2019). Moreover, comparisons of plant traits between ranges through in situ field measurements can give different results from common-garden experiments (Rosche et al., 2019). Therefore, more comprehensive studies that span the entire invaded and native ranges of S. canadensis and combine manipulative common-garden experiments with cytobiogeography and molecular phylogeography are required to more rigorously test whether divergence in ploidy levels and functional traits between ranges could be due to natural selection.

| No correlation between ploidy level and climatic envelope of S. canadensis
The invasive S. canadensis populations occupied a novel climatic envelope (i.e., there was no significant envelope overlap with native S. canadensis) and exhibited a contraction in climatic tolerance (i.e., occurred in warmer and less seasonal climates) relative to native S. canadensis populations independent of variation in ploidy level (

| Population demographic expansions in both ranges
Results of Tajima

| Signatures of reduced genetic founder effects in the invaded range
The finding that invasive populations had significantly lower mean number of polymorphic sites and haplotypes than native populations (

| All sampled populations belonged to S. canadensis
Although the ecological envelope of S. canadensis can overlap with that of its congeners like S. gigantea and S. altissima in the native range, our analysis of phylogenetic relatedness among the study populations did not reveal any clear pattern of clade formation (i.e., all populations were monophyletic; Appendices 9-11). Thus, it is highly likely that the invasive and native populations that we sampled in the field all belonged to S. canadensis.

| CON CLUS ION
Our results suggest that polyploidy does not influence climatic tolerance of S. canadensis in its invaded and native ranges. While a caveat exists for our study because our sampling of S. canadensis cytobiogeography did not cover its entire invaded and native ranges, the results instead indicate that unknown genetic and environmental factors may influence S. canadensis climatic tolerance. Our results provide a useful first step to comprehensively test whether polyploidy influences climatic envelope of S. canadensis. For S. canadensis, effects of polyploidy at molecular and nuclear levels have not been investigated and would be a first step in determining how genetic changes scale up to physiological processes and ultimately biogeographic distribution patterns. Therefore, future mechanistic studies should couple cytobiogeographical and molecular phylogeography approaches with functional genomic approaches that document gene expression patterns under different ecological conditions. In addition, experimental studies such as multisite commongarden and reciprocal transplants in both invaded and native ranges will be required.

ACK N OWLED G M ENT
We thank three anonymous reviewers and the associate editor for constructive comments on previous versions of the manuscript. We

CO N FLI C T O F I NTE R E S T
The authors declare that there is no conflict of interest.