Low genetic diversity contrasts with high phenotypic variability in heptaploid Spartina densiflora populations invading the Pacific coast of North America

4992 | Ecology and Evolution. 2018;8:4992–5007. www.ecolevol.org Received: 29 September 2017 | Revised: 6 March 2018 | Accepted: 15 March 2018 DOI: 10.1002/ece3.4063 ORIGINAL RESEARCH Low genetic diversity contrasts with high phenotypic variability in heptaploid Spartina densiflora populations invading the Pacific coast of North America Jesús M. Castillo1 | Blanca Gallego-Tévar1 | Enrique Figueroa1 | Brenda J. Grewell2 | Dominique Vallet3 | Hélène Rousseau3 | Jean Keller3 | Oscar Lima3 | Stéphane Dréano4 | Armel Salmon3 | Malika Aïnouche3 This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited. © 2018 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd. 1Departamento de Biología Vegetal y Ecología, Universidad de Sevilla, Sevilla, Spain 2Department of Plant Sciences MS-4, USDAARS Exotic & Invasive Weeds Research Unit, University of California, Davis, California 3UMR CNRS 6553 ECOBIO, Université Rennes 1, Rennes, France 4Faculté de Médecine, Institut de génétique et Développement de Rennes (IGDR), UMR6290, CNRS, Université de Rennes1, Rennes, France Correspondence Jesús M. Castillo, Departamento de Biología Vegetal y Ecología, Universidad de Sevilla, Sevilla, Spain. Email: [email protected] Abstract Species can respond to environmental pressures through genetic and epigenetic changes and through phenotypic plasticity, but few studies have evaluated the relationships between genetic differentiation and phenotypic plasticity of plant species along changing environmental conditions throughout wide latitudinal ranges. We studied inter- and intrapopulation genetic diversity (using simple sequence repeats and chloroplast DNA sequencing) and inter- and intrapopulation phenotypic variability of 33 plant traits (using field and commongarden measurements) for five populations of the invasive cordgrass Spartina densiflora Brongn. along the Pacific coast of North America from San Francisco Bay to Vancouver Island. Studied populations showed very low genetic diversity, high levels of phenotypic variability when growing in contrasted environments and high intrapopulation phenotypic variability for many plant traits. This intrapopulation phenotypic variability was especially high, irrespective of environmental conditions, for those traits showing also high phenotypic plasticity. Withinpopulation variation represented 84% of the total genetic variation coinciding with certain individual plants keeping consistent responses for three plant traits (chlorophyll b and carotenoid contents, and dead shoot biomass) in the field and in commongarden conditions. These populations have most likely undergone genetic bottleneck since their introduction from South America; multiple introductions are unknown but possible as the population from Vancouver Island was the most recent and one of the most genetically diverse. S. densiflora appears as a species that would not be very affected itself by climate change and sealevel rise as it can disperse, establish, and acclimate to contrasted environments along wide latitudinal ranges. KEYWORDS climate change, cordgrass, invasive species, microsatellites, phenotypic plasticity | 4993 CASTILLO eT AL. 1 | INTRODUCTION Invasive plant species can respond to environmental pressures in novel environments through genetic change at individual level (local adaptation), gene expression regulation and epigenetic processes (Richards, 2012), and phenotypic plasticity in response to different environmental conditions (Fuller et al., 2010; Miner, Sultan, Morgan, Padilla, & Relyea, 2005; Stomp et al., 2008). In addition, novel genotypes resulting from hybridization between native and invasive plant species may facilitate rapid evolution (Strong & Ayres, 2013). These different evolutionary processes are closely related as phenotypic plasticity may itself evolve through genetic changes (Crispo, 2008) and may contribute to the evolution of traits that support increased niche breadth of invasive species as plasticity facilitates persistence of individuals facing novel environmental pressures (Drenovsky et al., 2012). Phenotypic plasticity is adaptive, affecting the performance and reproductive success of individuals, being a trait itself, and, therefore, is subjected to evolution (Matesanz, Gianoli, & Valladares, 2010). Additionally, intrapopulation trait variation may modulate the effects of genetic diversity on populations (Hughes, 2014). Moreover, polyploidization, and especially hybridization, have been related to high levels of genetic diversity, epigenetic variation, and phenotypic plasticity (Braeutigam et al., 2013; Meyerson et al., 2016) and they are considered important in shaping the geographic range of plant species (Weiss- Schneeweiss, Emadzade, Jang, & Schneeweiss, 2013). In this context, few studies have evaluated the relationships between genetic differentiation and phenotypic plasticity of invasive plant species to novel, changing environments along wide latitudinal ranges. These studies have shown that species with low genetic variation can colonize wide latitudinal ranges due to high plasticity levels, sometimes related to polyploidy and hybridization, and their spread into dissimilar environments can lead to genetic differentiation (Counts, 1993; De Kort et al., 2016; Meyerson et al., 2016; Quilot- Turion et al., 2013; Zhao, Yang, Xi, Gao, & Sun, 2012). For example, the high level of epigenetic changes recorded in invasive Spartina anglica (Parisod et al., 2009; Salmon, Ainouche, & Wendel, 2005) has been related to its high levels of morphological plasticity and its large ecological amplitude (Thompson, McNeilly, & Gay, 1991), contrasting with low interindividual genetic diversity (Ainouche, Baumel, & Salmon, 2004; Ayres & Strong, 2001; Baumel, Ainouche, & Levasseur, 2001). In contrast, the cosmopolitan glycophyte Phragmites australis (Cav.) Trin. ex Steud. shows high genetic diversity, with some lineages behaving as invasive out of their native range, at the same time that their high phenotypic plasticity improves the responses to global change factors (Eller et al., 2017). Recently, Liu, Strong, Pennings, and Zhang (2017) have showed that both phenotypic plasticity and rapid local adaptation can contribute to the spread of an invasive species, being both processes interrelated and depending on environmental conditions changing with latitude. The knowledge of these dialectic relationships between genetic and epigenetic changes and trait expression is especially important for making decisions about conservation of endangered species and for management of invasive species in a scenario of globally changing environment that is modifying the distribution of organisms (Kelly & Goulden, 2008). Thus, it is important to link functional traitbased responses to genetic diversity and environmental changes to improve our understanding of mechanisms that promote changes in the species distribution and adaptation (Nicotra et al., 2010). For example, a better understanding of the role of phenotypic plasticity in response to environmental variation will improve our ability to manage weed invasions and improve conservation efforts. In fact, it has been hypothesized that global change and biological invasions are synergistic processes that interact and impart negative ecological and economic impacts (Vilà, Corbin, Ducks, Pino, & Smith, 2007). Thus, studies that integrate biological invasions, climate change, and phenotypic plasticity are needed (Engel, Tollrian, & Jeschke, 2011). The invasion of the heptaploid austral cordgrass, Spartina densiflora Brongn. (syn. Sporobolus densiflorus following Peterson, Romaschenko, Arrieta, & Saarela, 2014), across a 12° wide latitudinal gradient along the Pacific coast of North America from San Francisco Bay (ca. 37°56′N, USA) to Vancouver Island (ca. 49°20′N, Canada) provides a natural model system for an integrated study of the potential mechanisms underlying the response of species traits to substantial variation in climate and other environmental variables and an opportunity to investigate how species establish and succeed in novel and contrasted environments. Previous studies in the field and in commongarden experiments on invasive S. densiflora populations sampled along the Pacific coast of North America revealed high levels of phenotypic plasticity for some foliar and tussock traits (Castillo, Grewell, Pickart, Figueroa, & Sytsma, 2016; Castillo et al., 2014; Grewell, Castillo, Skaer- Thomason, & Drenovsky, 2016). However, the levels of genetic diversity in these S. densiflora populations along the Pacific coast of North America experiencing high variability in sedimentary and climatic conditions have not been studied. In this study, we performed genetic analyses of S. densiflora using chloroplast (cp) and nuclear DNA sequences, involving simple sequence repeats (SSRs). Cp- DNA sequences (useful for identification of maternal lineages in a hybridization context and for phylogeographic analyses) and potentially highly polymorphic nuclear SSRs were employed to explore genetic diversity among populations established along the Pacific coast of North America. Phenotypic variation within and among populations of S. densiflora was also evaluated through analysis of 33 morphological, physiological, and biochemical traits, on plants sampled in the field that were compared to plants maintained in a commongarden. Five invasive populations of S. densiflora along the Pacific coast of North America from San Francisco Bay to Vancouver Island were analyzed to examine the mechanisms allowing this species to deal with changing environmental conditions and the relationships between these mechanisms. Analyzing the level of intrapopulation phenotypic variability and the phenotypic variation of many plant traits between populations of the same species invading contrasted environments and in a commongarden experiment, together with knowledge of genetic diversity, would allow us to know the nature of the abovementioned phenotypic variation 4994 | CASTILLO eT AL. and its level of phenotypic plasticity (Zhao et al., 2016). Combining these different levels of analysis (phenotypic and genetic, contrasted and common environments, and intra- and interpopulation) is particularly relevant and should shed light on the ways that a polyploid species with hybrid origin (intrinsic characteristics that can promote phenotypic variability at different levels; Grover et al., 2012) can respond to spatial and temporal environmental changes. The quantification of plasticity for several traits in several populations of the same species is a difficult task, and thus, available data in this regard are still very scarce (Molina- Montenegro & Naya, 2012). Integrated biogeographic approaches involving the study of natural populations along latitudinal gradients paired with commongarden experiments are a way to increase our understanding of plant species responses to environmental change, yet the spatial scale of most studies has been limited (Elsey- Quirk, Seliskar, & Gallagher, 2011; Hierro, Maron, & Callaway, 2005). Our overall objective was to assess whether trait responses were dependent on genetic differentiation or phenotypic plasticity or a combination of both factors. We hypothesized that (i) S. densiflora would show low genetic diversity and high phenotypic plasticity in response to contrasted environments and that they would be independent from each other, based on previous observations (Castillo et al., 2014, 2016; Grewell et al., 2016); (ii) due to the polyploid and hybrid nature of S. densiflora, at least certain foliar traits would show both high plasticity and high intrapopulation variability, as leaf (the organ of photosynthesis and transpiration) is highly sensitive to environmental conditions (Stephenson, Oliver, Burgos, & Gbur, 2006); and (iii) genetically based and environmentalbased interpopulation differences in response to a changing environment would increase with geographic distance. 2 | METHODS 2.1 | Studied sites Spartina densiflora populations were sampled in summer of 2010 at five invaded salt marshes along the Pacific coast of North America, covering its entire latitudinal distribution in North America (along 12° of latitude). See Castillo et al. (2014, 2016) for a full description of sampling sites. The southernmost population analyzed was at Corte Madera Creek in San Francisco Bay Estuary (SF; 37°56′33″N, 122°30′55″W; California, USA) in a middleelevation salt marsh dominated by S. densiflora, which was intentionally introduced from Humboldt Bay to stabilize tidal creek banks in a wetland restoration project in the late 1970s (Faber, 2000). The next population to the north was located at Vance Marsh, Humboldt Bay (HB; 40°49′53″N, 124°10′17″W; California, USA) in a middle intertidal elevation marsh where the species has been probably introduced accidentally by merchant vessels from Chile during the late nineteenth century (Spicher & Josselyn, 1985). The third population was near the mouth of the Mad River Estuary (MR; 40°56′10″N, 124°7′48″W; California, USA) fringing the narrow intertidal zone at a river bank, just north of Humboldt Bay. We also included a population of a middle marsh elevation from Grays Harbor Estuary (GH, 46°57′N, 124°8′08″W; Washington, USA), where S. densiflora invasion was first recorded in 2001 (Pfauth, Sytsma, & Isaacson, 2007). Finally, the fifth and northernmost population was in a fringing salt marsh along the narrow Baynes Sound channel on the east coast of Vancouver Island (VI; 49°33″N, 124°52′09″W; British Columbia, Canada), region where S. densiflora was first discovered in 2005 (Morgan & Sytsma, 2010). Spartina densiflora populations were visited during low tides. Ten adult tussocks with between 10 and 80 live shoots and a mean tussock area between 15 and 41 cm2 (ca. 2–3 years of age) were randomly selected, separated by a minimum of 1 m in attempt to sample potentially different genotypes, and individually marked in each of the five populations. Due to the low genetic diversity of S. densiflora populations recorded previously in San Francisco Bay and Humboldt Bay (Ayres et al., 2008), we assumed that n = 10 individuals would be representative of each population, a number of plants per population similar to that sampled in previous studies on genetic diversity of invasive species (Andreakis, Kooistra, & Procaccini, 2009; Guo et al., 2015; Moran, Reid, & Levine, 2017). The external limits of each S. densiflora tussock were clearly visible, so sampling tussocks separated by a minimum of 1 m ensured that plants coming from different establishment events were sampled. 2.2 | Chloroplast DNA sequencing DNA was extracted from dried leaves of the 10 marked individuals from each of the five populations by employing a NucleoSpin Plant Extraction Kit (Macherey- Nagel GmbH & Co, Düren, Germany). Chloroplast DNA sequences were analyzed and compared with corresponding sequences from GenBank, obtained on S. densiflora samples from various regions (see below). The trnT-trnF chloroplast DNA (Cp- DNA) region was chosen as it is well documented in Spartina sequence databases (and then useful for chlorotype comparison between samples from different species and geographic locations). This region was amplified by polymerase chain reaction (PCR) in a Mastercycle thermal cycler (Eppendorf AG, Hamburg, Germany) and underwent denaturation for 2 min at 94°C, followed by 35 cycles at 94°C for 30 s, 54°C for 30 s, and 72°C for 1 min 30 s, followed by an extension phase at 72°C for 7 min. PCR products were purified using a NucleoSpin Gel and PCR purification kit (Macherey- Nagel GmbH & Co) and sequenced directly (on both directions) by Sanger technique at Macrogen Europe sequencer (Amsterdam, The Netherlands). The amplification of the trnT-F region was performed in two sections using the primer pairs a and b to recover the trnT-trnL segment and primer pairs c and f to recover the trnL-trnF segment (Taberlet, Gielly, Pautou, & Bouvet, 1991). Sequence alignments and comparisons were performed with the Geneious 9.1.2 software (Kearse et al., 2012). Sequences were compared with previously published S. densiflora sequences (Ayres et al., 2008;, Baumel et al., 2001; Fortune et al., 2008). | 4995 CASTILLO eT AL. 2.3 | Simple sequence repeats Neutral molecular markers such as SSRs provide independent and complementary assessment of genetic diversity to phenotypic variation most likely submitted to selective constraints. They are thus commonly used for genotyping in studies of phenotypic variations and fitness (Agarwal, Shrivastava, & Padh, 2008; Chapman, Nakagawa, Coltman, Slate, & Sheldon, 2009). A total of 1435 microsatellite repeats (SSRs) of 4–6 bp were identified from Spartina maritima genomic data obtained by Illumina sequence assemblies (A. Salmon, unpublished), using the MISA (MIcroSAtellite Identification Tool) software (Thiel, Michalek, Varshney, & Graner, 2003). Primer pairs were designed using the default parameters (and optimized TM ranging between 55 and 60°C) of the PRIMER3 software (Koressaar & Remm, 2007; Untergasser et al., 2012), which resulted in a range of primer product sizes of 100–280 bp and SSR motifs of at least −4 to 6 nt long, repeated five to eight times. This resulted in 395 different microsatellite loci with successful primer design (A. Salmon, unpublished data). For the first test, we selected 20 primer pairs considering those with the highest number of expected amplicon size and number of motif repetitions. PCR was performed as described above with these 20 primer pairs in S. densiflora nuclear DNA and analyzed by electrophoresis in agarose gel (1%, 35 min, 100V). Eight primer pairs of which amplification products produced unique and clear bands were chosen for further analysis (Table 1). The 50 S. densiflora DNA samples previously extracted were diluted to a concentration of 10 ng/μl and amplified using the eight microsatellite primer pairs selected. The protocol used for microsatellite locus amplification by Baumel et al. (2016) was modified as follows: Amplifications were performed in 15 μl that contained 0.4 ng/ μl DNA, 1× GoTaq Flexi Buffer (Promega), 0.02 U/μl GoTaq (Flexi, Promega), 0.2 mmol/L dNTPs, 2 mmol/L MgCl2, and 0.04 μmol/L reverse and forward primers. FAM- labeled reverse primers were used. Amplifications were conducted in a Mastercycle Eppendorf thermal cycler, using the following program: 94°C 2 min., followed by 35 cycles 94°C 30 s, 50°C 30 s, and 72°C 1.30 min. The amplified products were diluted to 1/30 before subsequent analysis. An aliquot of 2 μl of these diluted PCR products was mixed with 10 μl of formamide solution (975 μl formamide + 25 μl of GeneScan- 500 LIZ size marker) and separated by electrophoresis in an ABI PRISM 16- capillary 3130xl Genetic Analyzer (Applied Biosystems Inc., Waltham, USA). Alleles were identified and binned using GeneMapper 4.1 software (Applied Biosystems Inc., Waltham, USA). 2.4 | Genetic diversity analyses SSR alleles were determined by comparison with the standard marker size (GeneScan- 500 LIZ size standard), and the different “genotypes” (harboring different allele combinations) were scored. S. densiflora is a heptaploid species (2n = 7x = 70), so up to seven alleles may be expected per locus. As allelic dosage could not be ascertained from the obtained chromatograms, SSRs were treated as “dominant markers” (i.e., alleles were scored as either present or absent) using the GenAlEx 6.502 software (Peakall & Smouse, 2006, 2012). Genetic diversity analyses are less affected by polyploidy using this approach (Obbard, Harris, & Pannell, 2006). The intrapopulation genetic diversity parameters, effective number of alleles, Shannon Information Diversity Index (Brown & Weir, 1983), expected heterozygosity (Hartl & Clark, 1997), and unbiased expected heterozygosity (Peakall & Smouse, 2012) were obtained for each population. A cluster analysis was performed using the phangorn R package (Schliep, 2011) with the matrix of pairwise genetic differences and the unweighted pairgroup method with arithmetic averages (UPGMA; Sneath & Sokal, 1973). The proportion of genetic diversity among and within populations was assessed by molecular variance analysis (AMOVA) using the GenAlEx 6.502 software for the five sampled populations of S. densiflora. Genetic differentiation among populations (ΦPT), analogous to Wright’s F’ST for dominant data, was obtained following standardization (Meirmans, 2006). 2.5 | Phenotypic variation analyses The following 33 tussock and foliar phenotypic traits were evaluated in the field for the genetically analyzed individuals of S. densiflora: shoot density in (1) alive, (2) spiked, and (3) dead stem density (shoot per m2), shoot (4) length (cm) and (5) diameter (mm), number of (6) live and (7) dead leaves, (8) total nonstructural carbohydrates in rhizomes (mg/g), (9) root biomass (g/m2), (10) rhizome biomass (g/ m2), (11) belowground biomass (BGB, g/m2), (12) dead shoot mass (g/ m2), (13) live shoot biomass (g/m2), (14) dead leaf mass (g/m2), (15) Locus Forward primer sequence Reverse primer sequence SSR MS2 ATATTCCGATCCCCTCCTTG TTCGATCGGTCATGTTTTGA (AAAG)n MS7 CAGAATCACCATCATCAGCG TTCCATTTTTCAGGGTGAGC (TGGCAG)n MS13 CTTGACCGCAACCAGTATGA CCCAGGGCAATGGTTATACA (TTCT)n MS14 TGAGTTTGAGTTCACGGTTCA ATGTGATGCCATTTCCACAA (AAAG)n MS15 TGCATTGCAGCAAGAGAATC CGCTAGCTGATCCTGGAAAC (GATG)n MS16 GGGACACGGGATAGGAAAGT CCGCCGTGCAATTATTTATC (GTGGA)n MS17 TTTGTTCAGCTTCAGCATGG TTCTTGCAGTCGTTCTGTGC (GAAA)n MS18 TCTTATGGACCCCTTGCAGT CATCCGATTGGCGTAAGATT (TGATA)n TABLE1 Nucleotide sequences of the selected primer pairs and their corresponding microsatellite repeated motif for the 8 loci analyzed in five invasive populations of Spartina densiflora along the Pacific coast of North America 4996 | CASTILLO eT AL. live leaf biomass (g/m2), (16) aboveground biomass (AGB, g/m2), (17) AGB: BGB ratio, (18) leaf area index (LAI) (m2/m2), (19) maximum leaf wide (cm), (20) leaf length (cm), (21) leaf area (cm2), (22) specific leaf area (SLA) (m2/g), (23) leaf longitudinal adaxial rolling (%), (24) leaf elongation rate (mm/day), (25) leaf water content (LWC) (%), content of (26) chlorophyll (Chl) a (mg/g DW), (27) Chl b (mg/g DW), (28) Chl a+b (mg/g DW), and (29) carotenoids (Car) (mg/g DW), (30) Chl:Car ratio, and (31) Chl a: Chl b ratio, and foliar content of (32) carbon and (33) nitrogen (mg/g). Leaves traits were measured in flag leaves to minimize differences due to the foliar ontogeny. These plant functional traits were chosen as they together describe appropriately the leaf characteristics and the tussock architecture of S. densiflora, being key factors in shaping the responses of this clonal species to its surrounding environment. Detailed methods and values of these plant traits for the studied populations of S. densiflora were reported previously by Castillo et al. (2014, 2016) and they are analyzed in this study using phenotypic plasticity indices and in relation to genetic analyses. A commongarden experiment was carried out to analyze whether the variation in the plant traits of the five studied populations of S. densiflora was the same at contrasted environmental conditions in the field than in a common environment. Thus, the commongarden experiment allowed us to study the nature of the changing responses of S. densiflora, analyzing whether its variations in different plant traits were fixed independently of the surrounding conditions or whether they depended on the environment. With this aim, every tussock of S. densiflora measured in the field was collected and rhizomes were separated and grown in perlite substrate for 27 months (plastic pots: 20 cm diameter × 18 cm height) in a commongarden experiment in the glasshouse facility of the University of Seville, Spain (37º21′42′′N, 5º59′15′′W). Pot bases were kept permanently flooded to a height of 2 cm and watered with 20% strength modified Hoagland’s nutrient solution (Epstein, 1972; Hoagland & Arnon, 1950). Solutions were changed once a week. After 27 months of growth, every tussock trait recorded in the field, except belowground biomass, rhizome TNC, and leaf carbon and nitrogen concentration, was recorded again for the same plants acclimated to common glasshouse conditions (n = 10 tussocks per population). Common glasshouse conditions included substrate redox potential of 234 ± 5 mV, pH 8.3 ± 0.2, electrical conductivity 0.5 ± 0.0 mS/cm, mean monthly air temperature 23 ± 2°C, and mean monthly air relative humidity 62 ± 2%. Light conditions of mean radiation of 700 μmol photon m−2 s−1 at canopy level and a daily photoperiod of 16 hr (that was extended with incandescent lights, Osram Vialox NAV- T (SON- T) 400 W, giving a continuum spectrum) were set up to imitate those recorded at higher latitude sites during summer in relation to photoperiod and light intensity. Spartina densiflora is a facultative halophyte that can germinate, establish, and develop in freshwater conditions (Castillo, Rubio- Casal et al., 2005; Nieva, Castellanos, & Figueroa, 2001). Freshwater was used to avoid salinity effects on tussock and foliar growth. Phenotypic distances between plants in the field were calculated as Gower similarity index ranging between 0 and 1 (Gower, 1971) using the package vegan (Oksanen, 2007) of R software to obtain a pairwise dissimilarity matrix. Correlations between genetic, phenotypic, and geographic distances in the field were performed using the Pearson’s correlation coefficient (p < .05). Comparisons of the dissimilarity matrices described above were tested by Mantel’s test for matrix correspondence (Mantel, 1967) using GenAlEx 6.502 software with 1000 permutations, following the method of Smouse and Long (1992); Smouse, Long, and Sokal (1986). A cluster analysis was performed using the package vegan (Oksanen, 2007) of the R software with the matrix of phenotypic pairwise dissimilarity. Interpopulation phenotypic variability index (or phenotypic plasticity index sensu Valladares, Sanchez- Gomez, & Zavala, 2006) (PIv) was calculated for each plant trait to evaluate differences between populations in phenotypic trait expression in the field and in the commongarden experiment. PIv for each trait was quantified as the percentage of change in absolute high and low mean trait values (Molina- Montenegro & Naya, 2012): where Xmax is the maximum mean value for a given trait and Xmin is the minimum mean value for a given trait at population level (Valladares et al., 2006). This index is adequate for comparison of plasticity sensu lato, which is more focused on the responses of species to variable environments and allows traits with divergent ranges in variation and/or that are expressed in different units to be compared (Valladares, Dobarro, Sanchez- Gomez, & Pearcy, 2005; Valladares et al., 2006). The PIv recorded in the commongarden experiment was not the product of phenotypic plasticity in response to different environments, but it was calculated as a way to measure to what extent the PIv recorded in the field was related to phenotypic plasticity or with interpopulation differences in intrapopulation phenotypic variability. Intrapopulation phenotypic variability index (IVI) was calculated also for each plant trait and population in the field and in the commongarden experiment as the percentage of change in absolute high and low trait values: where xmax is the highest value for a given trait and population and xmin is the lowest value for a given trait and population. To estimate the overall intrapopulation variability for each trait in the field and in the commongarden experiment, we calculated the arithmetic average of the intrapopulation variability observed for all five populations. To estimate the intrapopulation variability for each population, we calculated the arithmetic average of the values recorded for all 26 or 33 traits (in the commongarden experiment and in field, respectively). The overall mean intrapopulation variability in field and in the commongarden experiment was calculated as the arithmetic average of the variability for every plant trait. Intrapopulation variability was compared among traits and among populations using oneway analysis of variance (ANOVA) PIv=[(Xmax −Xmin)∕Xmax]×100 IVI =[(xmax −xmin)∕xmax]×100 | 4997 CASTILLO eT AL. (F- t e s t , p < .05) using the plant trait or the geographic location as the grouping factor, respectively. Dependent variables were tested for homogeneity of variance using Levene’s test and for normality using Kolmogorov–Smirnoff test. The relationships between PIv and intrapopulation variability index in the field and in the commongarden experiment were analyzed using Pearson’s correlation coefficient (r) and linear regressions. Total intrapopulation variability and total phenotypic plasticity (calculated as arithmetic means for all traits and populations together) were compared using pairedsamples ttest. Deviations were calculated as standard error of the mean (SE). The αlevel of significance was p < .05 for all tests. These analyses were carried out using SPSS 12.0 (SPSS Inc., Chicago, USA). 3 | RESULTS 3.1 | Genetic diversity 3.1.1 | DNA sequencing The chloroplast DNA sequences analyzed for the trnT-L (GenBank accessions MG869717 to MG869721) and trnL-F (GB accessions MG869712 to MG869716) regions were identical in the five populations examined in this study. No difference was encountered between these sequences and those previously reported for S. densiflora samples from California (Humboldt Bay, Baumel et al., 2001; San Francisco, Ayres et al., 2008) and Argentina–Chile (Fortune et al., 2008). 3.1.2 | Simple sequence repeats The eight microsatellites primer pairs yielded a total of 19 different alleles. Each primer pair amplified between one and five different alleles per locus, with an average of two alleles per locus. The eight loci examined exhibited different levels of genotype diversity: One single genotype was detected for all populations with the markers MS2, MS7, MS16, and MS17, two different genotypes at the MS13, MS14, and MS18 loci and three genotypes were encountered for MS15 (Table 2). We found seven different multilocus genotypes when all loci recorded were considered, with the most common multilocus genotype (I) present in 40% of the examined individuals (Table 3). The population from Grays Harbor presented two specific multilocus genotypes (II and VI) that were not present in any other population (Table 3). The population from Mad River had the lowest genetic diversity (Table 4), showing only three multilocus genotypes (Table 3) with seven indistinguishable plants following our molecular markers, nine of 10 plants grouped together in the same cluster and with individuals present only in two clusters in the UPGMA analysis (Figure 1). The other four populations (SF, HB, GH, and VI) presented higher intrapopulation genetic diversity (Table 4), all of them showing four multilocus genotypes, except HB with three genotypes (Table 3) and their individual plants scattered in many different clusters (Figure 1). The population from Vancouver Island presented one group of four plants and another group of three plants sharing the same genetic markers. The population from Grays Harbor showed two groups with four and two plants with the same markers, the population from Humboldt Bay had also two groups with three and two plants with identical markers, and the six of the 10 plants from San Francisco Bay were grouped in three pairs sharing the same genetic profile. Some genetic differentiation was observed among populations (ΦPT = 0.159), but the AMOVA showed that interpopulation variance contributed to only 16% of the total genetic variation, whereas withinpopulation variation represented 84% of the total variation (Table 5). In fact, the five populations of S. densiflora were genetically very similar as showed by Nei’s standard genetic distance (from 0 between SF and HB, SF and VI, and HB and VI to 0.15 between SF and GH; Table 4). The UPGMA analysis based on genetic distances revealed that each cluster was composed of individuals from different populations, which agreed with low interpopulation genetic differentiation (Figure 1). The Mantel test found a significant relationship between genetic distance and geographic distance (p = .030, n = 50; Figure 2a). 3.2 | Phenotypic variation Contrasting to genetic variation, consistent phenotypic differentiation among populations was encountered in the field. The dendrogram using phenotypic distances clustered together every individual of the same population, except for the population from San Francisco Bay of which individuals were grouped in different two clusters with plants from Mad River and Vancouver Island (Figure 1b). The Mantel TABLE2 Genotype frequencies at eight microsatellite loci in five populations of Spartina densiflora along the Pacific coast of North America Genotypes (fragment sizes in bp) VI GH MR HB SF MS2 186 1 1 1 1 1 MS7 248/261 1 1 1 1 1 MS13 114/245 0.5 0.4 0.3 00.5 245 0.5 0.6 0.7 10.5 MS14 263/270 0.1 00.1 00.4 270 0.9 10.9 10.6 MS15 137/140/223/230/253 0.1 0 0 0.14 0.11 223/230/253 0.8 10.89 0.72 0.78 223/230/239/253 0.1 00.11 0.14 0.11 MS16 240/249 1 1 1 1 1 MS17 275 1 1 1 1 1 MS18 218/254 10.25 1 1 1 218/254/275 00.75 0 0 0 Populations: VI, Vancouver Island; GH, Grays Harbor Estuary; MR, Mad River Estuary; HB, Humboldt Bay Estuary; SF, San Francisco Bay Estuary 4998 | CASTILLO eT AL. Genotype Multilocus genotype VI GH HB MR SF Total I0101001101111110111 0.33 0.13 0.43 0.78 0.29 0.400 II 0101001101111111111 0.00 0.50 0.00 0.00 0.00 0.100 III 0101001111111110111 0.00 0.00 0.14 0.11 0.00 0.050 IV 0101111101111110111 0.11 0.00 0.00 0.00 0.14 0.075 V1101001101111110111 0.45 0.12 0.43 0.00 0.28 0.225 VI 1101001101111111111 0.00 0.25 0.00 0.00 0.00 0.050 VII 1111001101111110111 0.11 0.00 0.00 0.11 0.29 0.100 Alleles were scored as present = 1 or absent = 0 for the eight microsatellite loci. Populations: VI, Vancouver Island; GH, Grays Harbor Estuary; MR, Mad River Estuary; HB, Humboldt Bay Estuary; SF, San Francisco Bay Estuary. TABLE3 Multilocus genotypes and their proportion for five populations of Spartina densiflora introduced along the Pacific coast of North America VI GH HB MR SF N 9.842 ± 0.086 9.368 ± 0.175 8.579 ± 0.279 9.684 ± 0.110 9.579 ± 0.116 Na 1.211 ± 0.123 0.895 ± 0.130 1.105 ± 0.130 1.000 ± 0.132 1.211 ± 0.123 Ne 1.060 ± 0.037 1.072 ± 0.055 1.053 ± 0.030 1.018 ± 0.010 1.085 ± 0.045 I0.074 ± 0.035 0.060 ± 0.042 0.070 ± 0.034 0.033 ± 0.018 0.095 ± 0.043 He 0.042 ± 0.023 0.041 ± 0.029 0.040 ± 0.021 0.016 ± 0.009 0.057 ± 0.028 uHe 0.045 ± 0.024 0.043 ± 0.031 0.043 ± 0.022 0.017 ± 0.009 0.060 ± 0.029 VI – GH 0.013 – HB 0.000 0.012 – MR 0.002 0.013 0.001 – SF 0.000 0.015 0.000 0.003 – Populations: VI, Vancouver Island; GH, Grays Harbor Estuary; HB, Humboldt Bay Estuary; MR, Mad River Estuary; SF, San Francisco Bay Estuary. TABLE4 Number of different alleles (N), number of effective alleles (Na), Shannon Information Diversity Index (I), expected heterozygosity (He), unbiased expected heterozygosity (uHe) (values are mean ± SE), and pairwise matrix of Nei unbiased genetic distances between five populations of Spartina densiflora introduced along the Pacific coast of North America FIGURE1 Cluster analysis (UPGM) of Spartina densiflora individuals based on (a) Huff genetic distance calculated on eight SSRs loci and (b) on Gower dissimilarity index based on 32 phenotypic traits for 50 individuals from five populations recorded in the field along the Pacific coast of North America. Populations: VI, Vancouver Island; GH, Grays Harbor Estuary; MR, Mad River Estuary; HB, Humboldt Bay Estuary; SF, San Francisco Bay Estuary 00.1/64.000.001/92.1 (a) (b) | 4999 CASTILLO eT AL. test showed that phenotypic distance increased with geographic distance (p = .010, n = 50; Figure 2b), whereas no significant relationship was found between genetic distance and phenotypic distance in the field (p = .403, n = 50; Figure 2c). Mean IVI for the 33 recorded plant traits ranged from 7 ± 2% for foliar carbon content in the field to 96 ± 3% for dead shoot mass in the field (Figure 3a). PIv in the field was from 9% for foliar carbon content to 86% for the leaf area index. In general, the traits with higher PIv were those related to tussock aboveground architecture (shoot densities, above- and belowground biomass, and the leaf area index), and leaf adaxial rolling and elongation. In contrast, those traits showing lower PIv in the field were related to shoot and leaf morphologies, leaf number per shoot, leaf water content, Chl:Car ratio, and foliar carbon content (Figure 3b). Every population showed similar mean IVI from trait expression in the field (F = 1.644, df = 4, p = .166) and in the commongarden experiment (F = 0.711, df = 4, p = .586), with similar mean values for the five studied populations in the field (56 ± 3%) and in the commongarden experiment (60 ± 2%) (tt e s t , p > .05) (Figure 4a,b). Mean IVI for plant traits was positively correlated with PIv in the field (r = .878, n = 33, p < .0001) (Figure 4c) and in the commongarden experiment (r = .804, n = 26, p < .0001) (Figure 4d). The IVI showed similar values than PIv for every plant trait in the field (Figure 4c). On the other hand, those plant traits that had higher IVI and higher PIv in the commongarden experiment also had higher values for both variables in the field (r = .923, p < .0001, n = 26; r = .688, n = 26, p < .0001, respectively; Figure 5a,b). In addition, the IVI had similar values for every plant trait in the field and in the commongarden experiment (Figure 5a), whereas PIv tended to be lower in the commongarden experiment than in the field (Figure 5b). Thus, total intrapopulation variability and total phenotypic plasticity indices for every plant trait and every population were similar in the field (56 ± 4% and 54 ± 4%, respectively; tt e s t , p > .05). In contrast, total intrapopulation variability was higher than total interpopulation variability in the commongarden experiment (60 ± 2% and 37 ± 4%, respectively; tt e s t , t = 8.036, df = 25, p < .0001). Those individuals with higher values in the field also had higher values in the commongarden experiment for contents of Chl b and Car and dead shoot mass (Chl b: r = .493, n = 20, p < .05; Car: r = .435, n = 20, p = .055; dead shoot mass: r = .328, n = 50, p < .05) (Figure 6). Moreover, those individuals with higher contents of Chl b also had higher contents of Car both in the field and in the glasshouse (field: r = .876, n = 25, p < .0001; glasshouse: r = .910, n = 25, p < .0001; Figure 6a,b). 4 | DISCUSSION Our results show that the populations of Spartina densiflora gradually invading the Pacific coast of North America along 12° of latitude from the late nineteenth century when it was first introduced to Humboldt Bay, until recently with the discovery of a new invasion in 2005 at Vancouver Island, are genetically very similar according to the chloroplast DNA sequences and the nuclear microsatellite loci investigated. This low genetic variation may be accentuated through the founder effect (Brzyski, Taylor, & McLetchie, 2014) for recent invasive populations that are going through bottlenecks due to genetic drift (Bilton, Paula, & Bishop, 2002; Dlugosch & Parker, 2008), which reduces diversity even in populations with high sexual reproduction rates (Kittelson & Boyd, 1997) promoting genetic recombination (Hartfield, 2016). Our results agree with those of Ayres et al. (2008) who recorded low genetic diversity among S. densiflora plants in San Francisco Bay and Humboldt Bay, which is consistent with the origin of San Francisco Bay population as transplants introduced from Humboldt Bay (Faber, 2000). Low genetic diversity is also encountered in S. versicolor populations from Europe which appear to be introduced plants from the Native American range of S. patens (Baumel et al., 2016). Our value of genetic differentiation among populations (ΦPT = 0.159) is in the range of those reported previously for other invasive species such as S. alterniflora in China (FST = 0.234, Guo et al., 2015; FST = 0.100, Bernik, Li, & Blum, 2016), the polyploid red algae Asparagopsis taxiformis (Delile) Trevis. in the Mediterranean (ΦPT = 0.249, being lower than in its native range, Andreakis et al., 2009), Ambrosia artemisiifolia L. in France (FST = 0.054, Genton, Shykoff, & Giraud, 2005), Solidago canadensis L. in Europe (GST (analogous to FST) = 0.109, Moran et al., 2017), or Mikania micrantha Kunth in China (FST = 0.044, Geng, Chen, Cai, Cao, & Ou- Yang, 2017). These genetically similar S. densiflora populations showed high levels of phenotypic variability when growing in dissimilar environments in the field. Most of the interpopulation differences in plant traits recorded in the field disappeared when growing in commongarden conditions, indicating that the observed variation in natural populations was the result of phenotypic plasticity (Castillo et al., 2014, 2016; Grewell et al., 2016). Nevertheless, the mean interpopulation phenotypic variability recorded for every plant trait in the commongarden experiment was 37% (in comparison with the 56% recorded in the field), which pointed to another process, beyond phenotypic plasticity, that was also responsible for a significant amount of the recorded interpopulation phenotypic variation. In this sense, our study revealed that invasive S. densiflora populations had a similar, TABLE5 Analysis of molecular variance (AMOVA) illustrating the proportion of genetic variation distributed within and among populations of Spartina densiflora along the Pacific coast of North America Source df SS MS Est. Var. Percentage Among populations 46.943 1.736 0.113 16 Within populations 45 27.086 0.602s 0.602 84 Total 49 34.029 0.715 100 5000 | CASTILLO eT AL. fixed, and high intrapopulation phenotypic variability for many plant traits in the field and in the commongarden experiment (where it was higher than phenotypic plasticity for most of the recorded traits). Moreover, this intrapopulation phenotypic variability was especially high, irrespective of environmental conditions, for those traits showing also high phenotypic plasticity, suggesting that intrapopulation phenotypic variability in S. densiflora would be regulated through the same mechanisms as those regulating phenotypic plasticity. Ensuring high levels of intrapopulation phenotypic variability for key plant functional traits would allow S. densiflora to respond adequately to the changing environment (Zhao et al., 2016) of salt marshes, where local temporal and spatial changes on salinity, soil moisture, oxygenation, and flooding period are frequent (Adam, 1990; Castillo et al., 2000). The presence of intrapopulation phenotypic variability within species can result in faster range expansions than if only a single phenotype is present in the landscape (Elliott & Cornell, 2012). Many of the studied plant traits had PIv values higher than 50% in the field, and the mean PIv for every trait was 54 ± 4%, within the range recorded for other invasive plant species (Funk, 2008) and native species along wide latitudinal ranges (Chun- Can et al., 2014; Molina- Montenegro & Naya, 2012; Zhao et al., 2012). Our results are also consistent with previous studies on other polyploid Spartina species (Castillo, Redondo, Wharmby, Luque, & Figueroa, 2005; Elsey- Quirk et al., 2011; Thompson et al., 1991; Trnka & Zedler, 2000). S. densiflora is a heptaploid species of hybrid origin (Fortune et al., 2008), which most likely has played a key role in FIGURE2 Mantel test for similarity between (a) genetic and geographic distances (p = .030), (b) phenotypic and geographic distances (p = .010), and (c) phenotypic and genetic distances (p = .403) for 50 individuals from five populations of Spartina densiflora recorded in the field along the Pacific coast of North America (n = 50) | 5007 CASTILLO eT AL. Peakall, R., & Smouse, P. E. (2012). GenAlEx 6.5: Genetic analysis in Excel. Population genetic software for teaching and researchan update. Bioinformatics, 28, 2537–2539. https://doi.org/10.1093/ bioinformatics/bts460 Peterson, P. M., Romaschenko, K., Arrieta, Y. H., & Saarela, J. M. (2014). A molecular phylogeny and new subgeneric classification of Sporobolus (Poaceae: Chloridoideae: Sporobolinae). Taxon, 63, 1212–1243. https://doi.org/10.12705/636.19 Pfauth, M., Sytsma, M., & Isaacson, D. (2007). Oregon Spartina response plan (p. 70). Portland, OR: Oregon Department of Agriculture. Quilot-Turion, B., Leppala, J., Leinonen, P. H., Waldmann, P., Savolainen, O., & Kuittinen, H. (2013). Genetic changes in flowering and morphology in response to adaptation to a highlatitude environment in Arabidopsis lyrata. Annals of Botany, 111, 957–968. https://doi. org/10.1093/aob/mct055 Richards, C. L. (2012). Epigenetics: Linking genotype and phenotype in development and evolution. Integrative and Comparative Biology, 52, 547–549. https://doi.org/10.1093/icb/ics108 Richards, C. L., Bossdorf, O., Muth, N. Z., Gurevitch, J., & Pigliucci, M. (2006). Jack of all trades, master of some? On the role of phenotypic plasticity in plant invasions. Ecological Letters, 9, 981–993. https:// doi.org/10.1111/j.1461-0248.2006.00950.x Salmon, A., Ainouche, M. L., & Wendel, J. F. (2005). Genetic and epigenetic consequences of recent hybridization and polyploidy in Spartina (Poaceae). Molecular Ecology, 14, 1163–1175. https://doi. org/10.1111/j.1365-294X.2005.02488.x Schierenbeck, M., Fleitas, M. C., Miralles, J. D., & Simon, M. R. (2016). Does radiation interception or radiation use efficiency limit the growth of wheat inoculated with tan spot or leaf rust? Field Crops Research, 199, 65–76. https://doi.org/10.1016/j.fcr.2016.09.017 Schliep, K. P. (2011). Phangorn: Phylogenetic analysis in R. Bioinformatics, 27, 592–593. https://doi.org/10.1093/bioinformatics/btq706 Smouse, P. E., & Long, J. C. (1992). Matrix correlation analysis in anthropology and genetics. American Journal of Physical Anthropology, 35, 187–213. https://doi.org/10.1002/(ISSN)1096-8644 Smouse, P. E., Long, J. C., & Sokal, R. R. (1986). Multiple regression and correlation extensions of the Mantel test of matrix correspondence. Systematic Zoology, 35, 627–632. https://doi.org/10.2307/2413122 Sneath, P. H. A., & Sokal, R. R. (1973). Numerical taxonomy (p. 549). San Francisco, CA: W. H Freeman. Spicher, D., & Josselyn, M. (1985). Spartina (Gramineae) in Northern California: Distribution and taxonomic notes. Madroño, 32, 158–167. Stephenson, D. O., Oliver, L. R., Burgos, N. R., & Gbur, E. E. (2006). Identification and characterization of pitted morning glory (Ipomoea lacunose) ecotypes. Weed Science, 54, 78–86. https://doi. org/10.1614/WS-04-205R.1 Stomp, M., van Dijk, M. A., van Overzee, M. J., Wortel, M. T., Sigon, C. A. M., Egas, M., … Huisman, J. (2008). The timescale of phenotypic plasticity and its impact on competition in fluctuating environments. The American Naturalist, 172, 169–185. https://doi. org/10.1086/591680 Strong, D. R., & Ayres, D. R. (2013). Ecological and evolutionary misadventures of Spartina. Annual Review of Ecology, Evolution and Systematics, 44, 389–410. https://doi.org/10.1146/ annurev-ecolsys-110512-135803 Taberlet, P., Gielly, L., Pautou, G., & Bouvet, J. (1991). Universal primers for amplification of three noncoding regions of chloroplast DNA. Plant Molecular Biology, 17, 1105–1109. https://doi.org/10.1007/ BF00037152 te Beest, M., Le Roux, J. J., Richardson, D. M., Brysting, A. K., Suda, J., Kubešová, M., & Pyšek, P. (2011). The more the better? The role of polyploidy in facilitating plant invasions. Annals of Botany, 109(1), 19–45. Thiel, T., Michalek, W., Varshney, R., & Graner, A. (2003). Exploiting EST databases for the development and characterization of genederived SSR- markers in barley (Hordeum vulgare L.). Theoretical and Applied Genetics, 106, 411–422. https://doi.org/10.1007/s00122-002-1031-0 Thompson, J. D., McNeilly, T., & Gay, A. J. (1991). Population variation in Spartina anglica C. E. Hubbard. New Phytologist, 117, 115–128. https://doi.org/10.1111/j.1469-8137.1991.tb00951.x Trnka, S., & Zedler, J. B. (2000). Site conditions, not parental phenotype, determine the height of Spartina foliosa. Estuaries, 23, 572–582. https://doi.org/10.2307/1353147 Untergasser, A., Cutcutache, I., Koressaar, T., Ye, J., Faircloth, B. C., Remm, M., & Rozen, S. G. (2012). Primer3 - new capabilities and interfaces. Nucleic Acids Research, 40(15), e115. https://doi.org/10.1093/nar/gks596 Valladares, F., Dobarro, I., Sanchez-Gomez, D., & Pearcy, R. W. (2005). Photoinhibition and drought in Mediterranean woody saplings: Scaling effects and interactions in sun and shade phenotypes. Journal of Experimental Botany, 56, 483–494. Valladares, F., Sanchez-Gomez, D., & Zavala, M. A. (2006). Quantitative estimation of phenotypic plasticity: Bridging the gap between the evolutionary concept and its ecological applications. Journal of Ecology, 94, 1103–1116. https://doi.org/10.1111/j.1365-2745.2006.01176.x Vilà, M., Corbin, J. D., Ducks, J. S., Pino, J., & Smith, S. D. (2007). Linking plant invasions to global change. In J. Canadell, D. Pataki, & L. Pitelka (Eds.), Terrestrial ecosystems in a changing world (pp. 93–102). New York, NY: Springer. https://doi.org/10.1007/978-3-540-32730-1 Weiss-Schneeweiss, H., Emadzade, K., Jang, T. S., & Schneeweiss, G. M. (2013). Evolutionary Consequences, Constraints and Potential of Polyploidy in Plants. Cytogenetic and Genome Research, 140, 137–150. https://doi.org/10.1159/000351727 Zhao, Y., Yang, X., Xi, X., Gao, X., & Sun, S. (2012). Phenotypic Plasticity in the Invasion of Crofton Weed (Eupatorium adenophorum) in China. Weed Science, 60, 431–439. https://doi.org/10.1614/ WS-D-11-00198.1 Zhao, N. X., Zhang, L. H., Zhao, T. T., Mo, L. D., Zhang, J. L., Gao, Y. B., & Wang, J. L. (2016). Trait differentiation among Stipa krylovii populations in the Inner Mongolia Steppe region. Flora, 223, 90–98. https:// doi.org/10.1016/j.flora.2016.05.004 How to cite this article: Castillo JM, Gallego-Tévar B, Figueroa E, et al. Low genetic diversity contrasts with high phenotypic variability in heptaploid Spartina densiflora populations invading the Pacific coast of North America. Ecol Evol. 2018;8:4992–5007. https://doi.org/10.1002/ece3.4063