Analysis of population genetic structure and gene flow in an annual plant before and after a rapid evolutionary response to drought

Climate change can have widespread and devastating impacts, but little is known about the effects of this change on the genetics of natural populations. We examined the genetics of two populations of an annual plant (Brassica rapa) before and after an extended drought in California. We found limited genetic differences between pre-drought ancestors and post-drought descendants, indicating that the two populations did not become more genetically similar, despite a greater overlap in reproductive timing after the drought. The findings of this study contribute to an understanding of the genetic impacts of climate change, which is critical in optimizing conservation efforts.


Introduction
Climate change is a major environmental concern and one of the greatest threats to biodiversity Bellard et al. 2012). Climate can influence natural systems from a genetic to an ecosystem level, affecting individuals' fitness and species' survival, and is a major factor in determining community composition and ecosystem function (Bellard et al. 2012). The complex relationships between biotic and abiotic factors, combined with the intricate responses of species to a changing climate and other threats to biodiversity, make understanding the current and future impacts of global climate change difficult as well as of great importance (Parmesan and Yohe 2003;Pearson and Dawson 2003;Walther 2010).
Climate change can act as a strong selective pressure in natural populations (Jump and Peñ uelas 2005;Gienapp et al. 2008;Hoffmann and Sgró 2011). Several recent studies have illustrated the potential for rapid evolution in response to changes in species' surroundings (Huey et al. 2000;Stockwell et al. 2003;Franks et al. 2007). Adaptive responses in both plants and animals frequently involve the evolution of phenological traits (Fitter and Fitter 2002;Peñ uelas et al. 2002;Menzel et al. 2006;Hoffmann et al. 2010), which often depend on seasonal cues, such as temperature and/or precipitation conditions (Franks et al. 2007;Jarrad et al. 2008). Consequently, changes in climate patterns can have major impacts on phenology, and evolutionary changes in phenology may further affect adaptive ability by altering patterns of gene flow between populations (Garant et al. 2007;Franks and Weis 2009). For example, in populations that are differentially affected by changes in environmental conditions, phenology (i.e. timing of reproduction) may become more or less synchronized, impacting the opportunity for reproduction across populations.
Generally, gene flow is expected to reduce local adaptation by homogenizing populations found in differing environments, or by spreading detrimental alleles across populations (Ellstrand et al. 1989;Holt and Gomulkiewicz 1997). However, gene flow may also serve to introduce potentially adaptive alleles to populations, and increase genetic diversity, which natural selection can act upon to provide an evolutionary response (Ellstrand et al. 1989;Räsänen and Hendry 2008;North et al. 2011;Sexton et al. 2011). An improved understanding of both the effects of climate on evolutionary changes in phenology, and the consequences of these phenological changes for gene flow and local adaptation, will be valuable in predicting how species may respond to changes in climate (Garant et al. 2007).
Adaptation can provide a response to changes in climate in some natural populations (Franks et al. 2014). Franks et al. (2007) found that an extended drought in southern California led to an adaptive response in the flowering times of two populations of Brassica rapa (field mustard). The drought, lasting from 2000 to 2004, corresponded to a shift to earlier mean flowering time in both populations, causing an abbreviated growing season and resulting in more synchronous flowering times between the populations following the drought (Franks et al. 2007;Franks and Weis 2009).
Seeds were collected from both sites in May 1997 (prior to the drought) and June 2004 (immediately following the drought), at the end of the growing season for both populations in these years. A parent-offspring analysis from a greenhouse study demonstrated an additive genetic basis for this evolutionary change in flowering time (Franks et al. 2007). Additionally, further analysis has shown a decrease in phenological isolation of the two sites between 1997 (Isolation Index ¼ 33 %) and 2004 (Isolation Index ¼ 14 %) as a result of the change in flowering time (Franks and Weis 2009).
The current study focuses on the impacts of the change in flowering times on the genetics of these natural populations. We analysed the population genetics of this system before and after the drought to examine genetic structuring and gene flow across this phenological change. While previous studies have used population genetic data to infer historic levels of gene flow and genetic structure, this study is unique in using the 'resurrection approach' (Franks et al. 2008) of comparing population genetic patterns in ancestors and descendants using stored seeds.
Based on previous analyses of this system (Franks et al. 2007;Franks and Weis 2009), we expected to find an increase in gene flow between the populations following the drought, corresponding to decreased phenological, and potentially reproductive, isolation. Given the longdistance (.5 km) pollinator-mediated pollen dispersal found in related species (Devaux et al. 2005(Devaux et al. , 2007, we expected that our sites (3 km apart) would be connected by some degree of gene flow. Regardless, our approach allowed us to determine genetic changes that occurred both within and between populations over time. Changes in population structure would not only elucidate how an adaptive evolutionary change influenced genetic patterns, but may also have implications for the ability of these populations to further evolve in response to continued changes in climate. The results of this study can contribute to the growing understanding of natural responses to climate change events, which will be necessary in aiding in the conservation of species as the environment continues to change.

Study system
Brassica rapa is a weedy plant, native to Eurasia but currently distributed worldwide with cultivars and naturalized populations (Suwabe et al. 2002). The populations of interest in this study are found in southern California, where B. rapa was introduced 300 years ago and has become naturalized (Franks et al. 2007). In these populations, B. rapa is an annual plant, produces perfect flowers and is self-incompatible. Once an individual plant has begun to flower, many new flowers may open each day, each flower lasting 3 days, for the remainder of the growing season. In southern California, B. rapa germinates in the winter, flowers in the spring and continues to produce flowers until the end of the rainy season (late spring to summer). Flowers are pollinated by many species of vagile insects, such as bees and flies, allowing for the dispersal of pollen across these insects' foraging range (Davis et al. 1994;Warwick et al. 2003).
The populations were located in upper Newport Harbor Back Bay (BB) and near the San Joaquin Marsh Arboretum (Arb), in Irvine, CA, USA (Fig. 1). The two locations are 3 km apart and are characterized by differing soil conditions. Back Bay soil is sandy and dry in comparison to Arb, where the soil has a higher water table and clay content, aiding in moisture retention (Franke et al. 2006). These two populations also differed in their flowering time ranges, with BB flowering earlier than Arb prior to the drought (Franks et al. 2007). Both populations comprised hundreds to many thousands of individuals both before and after the drought.

Samples
Seeds collected from individuals in the Arb and BB populations in 1997 and 2004 were stored at 4 8C. Over 10 000 seeds were collected from individuals across both populations to obtain representative samples of the gene pools. Seeds from this stock were grown for tissue collection on light carts in the Franks laboratory at Fordham University. Plants were watered daily and fertilized weekly yielding 259 total samples (56 Arb 1997, 39 BB 1997, 89 Arb 2004and 75 BB 2004 (Table 1).

DNA extraction
Leaf tissue was collected from individuals following the development of true leaves. Approximately 50 mg of tissue per individual was removed, dried and stored in silica gel until extraction. Dried tissue was homogenized using Qiagen's TissueLyser system (Qiagen, Valencia, CA, USA). DNA was extracted using a Qiagen DNeasy plant kit following manufacturer's DNeasy 96 Protocol. This extraction process yielded 200 mL of eluted DNA, which was stored at 220 8C.

Microsatellite analysis
Ten microsatellite loci were chosen based on repeatable amplification success and high degree of polymorphism: BN12A (Szewc-McFadden et al. 1996), Na10-A08, Na10-D09, Na10-G10, Ni4-A03, Ol10-D08, Ra2-E04, Ra2-E12 (Lowe et al. 2002(Lowe et al. , 2004, BRMS-040 and BRMS-037 (Suwabe et al. 2002). PCR was carried out in 12.5 mL reactions using TC-5000 thermal cyclers (Techne, Bibby Scientific Ltd, Staffordshire, UK) to isolate and amplify microsatellite loci. These reactions included 0.5 mL of unlabelled forward and reverse primers, 0.5 mL of  A second PCR was run to incorporate a fluorescent label (WellRED Dye D2, D3 or D4) (Sigma-Aldrich Co., St. Louis, MO, USA) using 0.5 mL reverse primer and 0.5 mL fluorescently labelled primer (M13F (-20) ) corresponding to the tag sequence added to the 5 ′ end of the initial forward primer. This was included in a reaction with 0.5 mL of 1 : 100 diluted PCR product from the first PCR, along with RNase-free water and Am-pliTaq Gold 360 Master Mix. This second PCR followed the same cycling protocol as the first. PCR products were run on 1 % agarose gel to confirm the presence of target microsatellite fragments. Fragment length was then analysed with a Beckman Coulter CEQ 8800 using Genetic Analysis System v.9.0 (Beckman Coulter, Fullerton, CA, USA). Samples were run with 0.5 mL of Beckman Coulter 400 bp Size Standard and PCR product at volumes according to their fluorescent label: 2.5 mL of M13F(-20), 2.5 mL of M13F(-40) and 1.5 mL of M13R. Alleles were scored using GENEMARKER v.1.90 (SoftGenetics, State College, PA, USA).

Statistical analysis
GENALEX v.6.41 (Peakall and Smouse 2006) was used to determine measures of microsatellite diversity, including number of alleles (A), expected heterozygosity (H e ), observed heterozygosity (H o ) and inbreeding coefficient (F IS ) for each locus. ARLEQUIN v.3.5.1.3 (Excoffier and Lischer 2010) was then used to calculate Hardy -Weinberg Equilibrium (HWE) and to test for linkage disequilibrium (LD) across all loci. Departure from HWE was determined for each locus using the Markov chain method with a length of 1.0 × 10 5 steps (Guo and Thompson 1992). Linkage disequilibrium was determined using an Expectation Maximization algorithm with 1.0 × 10 4 permutations (Lewontin and Kojima 1960;Slatkin 1994;Slatkin and Excoffier 1996). The frequency of null alleles for a locus was estimated using an Expectation Maximization algorithm (Dempster et al. 1977) in GENEPOP v 4.0.10 (Raymond and Rousset 1995;Rousset 2008). These parameters describe the loci and can indicate the influence that each has in defining measures of structure in these populations.
GENALEX was also used to determine measures of allelic diversity within these populations including polymorphism (P) across all loci, the number of private alleles (PA), A, H e and H o . F IS was calculated using ARLE-QUIN for each population. F-statistics, such as F IS (inbreeding) and F ST (fixation), are used to describe the fixation of alleles for individuals in a population, and for populations in a metapopulation, respectively (Wright 1951;Weir and Cockerham 1984). Several analogues of F ST have been developed to account for the mutational characteristics of microsatellites (Slatkin 1995;Michalakis and Excoffier 1996;Hedrick 2005;Jost 2008). In addition to F ST , we also calculated f ST (Michalakis and Excoffier 1996), which takes into account the step-wise mutation model believed to be characteristic of many microsatellites. Because F ST and gene flow are inversely related (Wright 1951), we used F ST to test the hypothesis that gene flow increased following the evolutionary shift to greater flowering synchrony between the populations.
Pairwise comparisons of F ST and f ST , as well as AMOVA were performed in GENALEX. To determine whether levels of fixation and gene flow changed across this period, we compared F ST values across these years. Significance was established using 95 % confidence intervals (CI) as determined through 1000 bootstrapping replicates to estimate F ST using GDA v.1.1 (Lewis and Zaykin 2002).
Bayesian estimates of genetic clustering probabilistically assigns individuals to populations defined by allele frequencies at multiple loci, and were determined using STRUCTURE v.2.3.2.1 (Pritchard et al. 2000) for 5.0 × 10 4 burn-in repetitions and 1.0 × 10 6 MCMC simulations at four iterations. This allowed for an estimate of the number of genetic units, K, following Evanno et al. (2005) using STRUCTURE HARVESTER web v.0.6.92 (Earl and vonHoldt 2012). Outputs of four iterations at K ¼ 3 were run in  (Jakobsson and Rosenberg 2007) to align these clustering results, and were then visualized using DISTRUCT v.1.1 (Rosenberg 2004). Data were reanalysed using only loci in HWE in each population [see Supporting Information], as the analyses assume that these microsatellites are acting as neutral markers and a deviation from HWE may be an indication of an effect of selection. The same procedures were used to estimate F IS within each population, F ST and f ST between Arb and BB in 1997 and 2004 and to compare F ST between the 1997 and 2004 generations. Results for these analyses are noted along with the results of analyses run with all 10 loci.

Allelic characteristics
Across the populations, the loci analysed showed a mean polymorphism of 95 % per population (Table 1). Between two and nine alleles were found per locus (mean ¼ 5.8 alleles per locus) [see Supporting Information] and between zero and eight PA were found in each population (Table 1). Of the 38 polymorphic loci (for 10 loci across four populations), 26 exhibited departure from HWE (P , 0.05). Out of 180 total locus pairs, 22 pairs across all populations showed LD. However, these locus pairs were not linked in each population. The estimated frequency of null alleles ranged from ,1.0 × 10 24 to 0.3558 for all loci. The loci in HWE were found to have a lower frequency of null alleles. Therefore, supplementary analyses performed using only loci in HWE, also minimized the influence of allelic dropout.

Intra-population characteristics
We found high genetic diversity, and some changes in the number of PA, but no significant change in overall genetic diversity over time [see Supporting Information]. In the BB population, the number of PA increased from zero to eight between the sampling years ( All populations showed fewer than expected heterozygotes and high levels of inbreeding (F IS ) (Table 1). F IS was significantly greater than zero for all populations and generations ( Table 1). Levels of F IS did not change significantly over time [see Supporting Information]. Measures of F IS using only loci in HWE showed lower levels of inbreeding (between 0.065 for Arb97 and 0.085 for BB04) [see Supporting Information] than the analyses performed using all 10 loci (Table 1), which is expected as random mating is an assumption of HWE.

Genetic structure and gene flow
We estimated fixation between all pairs of populations (BB 97, BB 04, Arb 97 and Arb 04) as F ST and f ST . Fixation between populations was low but significant for all pairs except BB 97 and BB 04. F ST ranged from 0.027 (between Arb 97 and Arb 04) to 0.039 (BB 97 and Arb 04) ( Table 2).
In addition to pairwise comparisons, we used AMOVA to measure F ST and f ST in both 1997 and 2004. We found that estimates of F ST and f ST were low but significant and increased slightly from 1997 to 2004 (Table 3)  To compare fixation between the two populations before and after the drought, we permuted F ST (Fig. 2) using 1000 bootstrap replications. We estimated F ST to be 0.053 (0.019 -0.082; 95 % CI) between BB and Arb in 1997 and 0.076 (0.024 -0.136) in 2004. As these estimates of F ST do not overlap with zero, this suggests significant levels of differentiation between the populations in both years. However, the estimates of F ST were not significantly different between these 2 years (Fig. 2), indicating a lack of change in fixation or gene flow between the generations. Analysis using loci in HWE similarly shows no    Bayesian analysis of this data suggests the existence of three distinct clusters (K ¼ 3) in this system [see Supporting Information]. The greater proportion of individuals assigned to the genetic cluster represented by yellow in BB as compared with Arb indicates some genetic distinctiveness between these locations in both years. However, every individual sampled is estimated, in some proportion, to be assigned to each of these three clusters (yellow, orange and blue) (Fig. 3) suggesting low levels of structure, and consequently high levels of gene flow or genetic similarity in 1997 and 2004.

Discussion
The results of this study show relatively high levels of inbreeding within the Arb and BB populations, and low fixation between these locations in 1997 and 2004. Despite a previous finding of reduced phenological, and potentially reproductive isolation between these populations following a 5-year drought, we did not find a significant change in gene flow or in genetic structuring across this time.

Genetic diversity
Microsatellite data from these populations revealed high allelic diversity. Across all loci, we found 58 alleles, com-  (Selkoe and Toonen 2006), as no nearby populations are known to exist to introduce new alleles via gene flow. Mutation rate is dependent on the effective population size and can vary due to environmental stresses, allowing for different trends in allelic diversity at each site (Hoffmann and Hercus 2000). However, the increase in PA in BB may also be due, fully or in part, to the relatively small sample size of BB 97 compared with BB 04 and a  Table 3. Hierarchical AMOVA tables. The proportion of genetic variation is partitioned among populations, within populations and within individuals. Overall fixation, F ST , and corresponding Nm provided for both years (1997 in (a); 2004 in (b)  failure to detect low frequency alleles with our sampling effort. Thus, these populations appear to have high diversity at these loci, but there is little evidence to suggest changes in diversity over time.

Intra-population genetic structuring
We found evidence for high levels of inbreeding in our populations despite the genetic self-incompatibility of B. rapa (Sobotka et al. 2000). This high level of inbreeding could be due to local seed and pollen dispersal causing mating between relatives, or to assortative mating. Some level of assortative mating was likely in this system due to the range of flowering times, as individuals may only cross with other individuals contemporaneously in flower (Weis and Kossler 2004;Franks and Weis 2009), and individuals with similar flowering times may be more likely related to one another than expected by chance. Levels of inbreeding did not change significantly between these years despite changes in the flowering phenology of each population (Franks and Weis 2009). As the drought resulted in an abbreviated growing season in 2004, this could allow for a greater proportion of individuals to be in flower at the same time thereby reducing the likelihood of potential crosses between closely related individuals. However, genetic similarity among all individuals may have allowed these populations to maintain high F IS across these generations. Ultimately, this shift in phenology did not appear to influence the genetic structuring within each population.

Fixation and gene flow over time
We found low levels of fixation, and correspondingly high levels of gene flow or shared common ancestry, between the two populations both before and after the drought.
The lowest values of fixation were found between years within the same population, which is likely due to the 2004 generation maintaining much of the genetic character of its ancestral 1997 populations. However, we found very little deviation from these F ST values between all other population comparisons, indicating maintenance of, and no large change in, population connectivity over these generations. We calculated F ST and f ST in both of these generations to investigate whether a change in gene flow occurred with the decreased phenological isolation between the populations. Although these parameters have been critiqued as estimates of genetic differentiation, particularly for microsatellite data (Whitlock and McCauley 1999;Jost 2008), the focus of our analysis is to determine if the populations became more genetically similar over time, as would occur if gene flow increased. Thus, examining these parameters and their change over time is appropriate for our goals. Both F ST and f ST indicate similar patterns of fixation between the BB and Arb populations with lower, significant fixation occurring prior to the drought in 1997 and higher, significant fixation following the drought in 2004 (Table 3) [see Supporting Information]. These values of fixation are very low and represent minimal structure between the populations. The statistical significance of these measures suggests a deviation from no fixation, however the low estimates of F ST may not represent biologically meaningful levels of structure. The 95 % bootstrapped CI indicate no significant change in F ST or gene flow over time (Fig. 2). This result was not expected given the increased phenological opportunity for gene flow. However, in some systems, a decrease in phenological isolation may not cause changes in gene flow, particularly when minimal structure is present prior to the phenological change.
The data were also analysed using Bayesian statistics, which showed little genetic distinctiveness between these populations and generations. Three genetic clusters were recovered for this dataset, and all individuals in each population were assigned to these clusters to varying degrees (Fig. 3). This indicates a high degree of genetic relatedness of the BB and Arb sites across these years. Despite some variation across the two sites in terms of the relative proportions of each cluster, the assignment of all individuals sampled to each of the three clusters suggests that these populations had a recent shared ancestry. Furthermore, no major shift in proportion of these genetic units is evident across time (Fig. 3).
Overall, this assessment of the data supports the results of F ST analyses. The minimal observed change in structure and gene flow may have been due to an insignificant change in migration rate over time, or these results may have instead been impacted by the populations' genetic similarity or large effective population sizes (Neigel 2002).

Evolutionary, environmental and community considerations
The lack of an expected change in gene flow over time could be due to several factors including the evolutionary history of these populations, seed banks, a low initial degree of genetic structure between the populations or a direct effect of the drought on the movement of pollen between populations. While research has been done to examine the phenological traits of these populations, little is known about their evolutionary history preceding the initial Franks et al. (2007) study, and pollination dynamics are not well understood. The low levels of fixation maintained between the two populations across seven generations (Table 3) [see Supporting Information] could have resulted if they had been very genetically similar across the data collection period. For example, if both populations had been established by individuals from the same source population, one population had been established from the other, or they were effectively a single population in the recent past, they would exhibit low levels of fixation. In any of these cases, fixation could have increased over time, given the 3 km separating these populations, allowing for divergence through local effects of genetic drift and mutation at neutral loci. However, it is possible that neither of these processes had an effect on fixation in this time, or that high levels of gene flow diminished a trend toward fixation. Seed banks present in these populations might influence phenology and genetics. Both populations likely have persistent seed banks, as seed viability remains high over time, especially under the dry conditions that these populations experience outside of the growing season, and as persistent seed banks were found in the related species, Brassica napus (Simard et al. 2002). The presence of a seed bank would serve to hinder the evolutionary response to selection across this climatic shift, as genotypes adapted to previous conditions could emerge after conditions changed. However, evolutionary change in phenology was previously documented in these populations across these years (Franks et al. 2007) despite the likely presence of a seed bank. Nevertheless, a seed bank could have been a contributing factor to the lack of change in gene flow seen in these populations.
Additionally, given that little information exists regarding the pollinators of this system in these years, it is unknown whether they were able to easily traverse the 3 km distance separating the populations. Some sources define 3 km as the maximum distance that insects will forage (Rieger et al. 2002), while others believe that flowering plant populations within a range of 10 -20 km can still exhibit substantial levels of gene flow, depending on the type of insect pollinator (dePamphilis and Wyatt 1989;Ellstrand 1992;Cresswell et al. 2000). Therefore, understanding whether the distance between these populations is a reasonable foraging range for this system's pollinators is necessary to properly interpret these results. The observed levels of fixation suggest very high levels of gene flow but this may instead be a consequence of historic patterns of population connectivity (Heywood 1991;Slatkin 1995;Whitlock and McCauley 1999).
The landscape surrounding these two populations is highly urbanized and developed (Fig. 1), but contains managed conservation areas surrounding the populations, including the UCI Arboretum grounds and the Upper Newport Bay Nature Preserve and Ecological Reserve. Although weedy populations of other wild and landscaped flowering plants are present in the area and could potentially serve as stepping stones in pollinators moving between the two study populations, the extent to which this occurs, and the effects of urbanization in this area on pollinator movement, is not known.
Another potential explanation for the absence of a significant change in gene flow with decreased phenological isolation is the effect of the drought, or other unknown changes, on pollinators. As B. rapa is primarily pollinated by insects, gene flow would require the foraging activity of these pollinators. As climate change is negatively impacting many systems across the globe (Hannah et al. 2002;Parmesan and Yohe 2003), the effect of this drought may have altered the population sizes of insects in the area or caused some modification of their behaviour that impacted their ability to transfer pollen across this range (Ollerton 1996;Rafferty and Ives 2011). Therefore, if the potential for increased gene flow existed in these populations over this drought period, negative impacts of the drought on pollinators could have had effects throughout the community, which may have minimized the level of gene flow achieved through this change (Ollerton 1996;Williams et al. 2008;Walther 2010;Rafferty and Ives 2011). As little is known about the evolutionary history of this system and about the community at the time of this drought, any variation in gene flow or genetic structuring cannot be directly attributed to a change in climate. However, an effect of the drought on the pollinators in this area would likely have consequences for gene flow (Cresswell et al. 2000).

Conclusions
Despite the lack of a significant change in genetic structuring of these populations in response to a climate change event, this case highlights the importance of understanding the effects of environmental changes on gene flow between natural populations. Further research examining the genetic responses of natural populations to climate change across multiple systems will be useful in increasing our understanding of the evolutionary consequences of climate change (Franks and Hoffmann 2012) and establishing a baseline understanding with which to advise management decisions.

Sources of Funding
Funding was provided by a Sigma Xi Grant-in-Aidof-Research (G20101015154375) (USA) to R.S.W., a Faculty Research Grant from Fordham and a grant (DEB-1142784) from the National Science Foundation (USA) to S.J.F. and by the New York Botanical Garden.

Supporting Information
The following additional information is available in the online version of this article - Table S1. 'HWE estimates of F IS ' provides estimates of F IS , and corresponding 95 % CI, from 1000 bootstrap replicates using only loci found to be in HWE. Loci used are listed for each population and generation. Table S2. 'HWE estimates of F ST ' provides estimates of F ST , and corresponding 95 % CI (when available), from 1000 bootstrap replicates using only loci found to be in HWE. Loci used are listed for each generation. Table S3. 'Hierarchical AMOVA tables using loci in HWE' (F ST ) provides estimates of F ST and the distribution of molecular variance for both generations, using only loci in HWE. Table S4. 'Hierarchical AMOVA tables using loci in HWE' (f ST ) provides estimates of f ST and the distribution of molecular variance for both generations, using only loci in HWE. Table S5. 'Loci in HWE' provides a list of the loci estimated to be in HWE for each population and generation. Table S6. 'Microsatellite loci characteristics' lists allelic parameters calculated for each of the 10 loci used in this study. Table S7. 'Hierarchical AMOVA tables' provides estimates of f ST and the distribution of molecular variance for both generations, using all 10 loci. Table S8. 'Within population parameters' provides estimates of standard deviation for genetic parameters of each population and generation. Table S9. 'STRUCTURE likelihood outputs' provides likelihood scores used to determine the optimal number of clusters comprising these populations. Table S10. 'Within population parameters over time' provides the results of Student's t-test comparing within population genetic parameters across generations. Table S11. 'CI for within population fixation' provides the 95 % CI corresponding to estimates of F IS from 1000 bootstrap replicates, for each population and generation. Table S12. 'Estimates of F IS across populations' provides global estimates of F IS for each generation, and across generations.