Conservation genomics of anadromous Atlantic salmon across its North American range: outlier loci identify the same patterns of population structure as neutral loci

Anadromous Atlantic salmon (Salmo salar) is a species of major conservation and management concern in North America, where population abundance has been declining over the past 30 years. Effective conservation actions require the delineation of conservation units to appropriately reflect the spatial scale of intraspecific variation and local adaptation. Towards this goal, we used the most comprehensive genetic and genomic database for Atlantic salmon to date, covering the entire North American range of the species. The database included microsatellite data from 9142 individuals from 149 sampling locations and data from a medium‐density SNP array providing genotypes for >3000 SNPs for 50 sampling locations. We used neutral and putatively selected loci to integrate adaptive information in the definition of conservation units. Bayesian clustering with the microsatellite data set and with neutral SNPs identified regional groupings largely consistent with previously published regional assessments. The use of outlier SNPs did not result in major differences in the regional groupings, suggesting that neutral markers can reflect the geographic scale of local adaptation despite not being under selection. We also performed assignment tests to compare power obtained from microsatellites, neutral SNPs and outlier SNPs. Using SNP data substantially improved power compared to microsatellites, and an assignment success of 97% to the population of origin and of 100% to the region of origin was achieved when all SNP loci were used. Using outlier SNPs only resulted in minor improvements to assignment success to the population of origin but improved regional assignment. We discuss the implications of these new genetic resources for the conservation and management of Atlantic salmon in North America.     

Conserving marine biodiversity: insights from life-history trait candidate genes in Atlantic cod (Gadus morhua)

Recent technological developments have facilitated an increased focus on identifying genomic regions underlying adaptive trait variation in natural populations, and it has been advocated that this information should be important for designating population units for conservation. In marine fishes, phenotypic studies have suggested adaptation through divergence of life-history traits among natural populations, but the distribution of adaptive genetic variation in these species is still relatively poorly known. In this study, we extract information about the geographical distribution of genetic variation for 33 single nucleotide polymorphisms (SNPs) associated with life-history trait candidate genes, and compare this to variation in 70 putatively neutral SNPs in Atlantic cod (Gadus morhua). We analyse samples covering the major population complexes in the eastern Atlantic and find strong evidence for non-neutral levels and patterns of population structuring for several of the candidate gene-associated markers, including two SNPs in the growth hormone 1 gene. Thus, this study aligns with findings from phenotypic studies, providing molecular data strongly suggesting that these or closely linked genes are under selection in natural populations of Atlantic cod. Furthermore, we find that patterns of variation in outlier markers do not align with those observed at selectively neutral markers, and that outlier markers identify conservation units on finer geographical scales than those revealed when analysing only neutral markers. Accordingly, results also suggest that information about adaptive genetic variation will be useful for targeted conservation and management in this and other marine species.     

Population differentiation determined from putative neutral and divergent adaptive genetic markers in Eulachon (Thaleichthys pacificus,Osmeridae), an anadromous Pacific smelt

Twelve eulachon (Thaleichthys pacificus, Osmeridae) populations ranging from Cook Inlet, Alaska and along the west coast of North America to the Columbia River were examined by restriction‐site‐associated DNA (RAD) sequencing to elucidate patterns of neutral and adaptive variation in this high geneflow species. A total of 4104 single‐nucleotide polymorphisms (SNPs) were discovered across the genome, with 193 putatively adaptive SNPs as determined by F_ST outlier tests. Estimates of population structure in eulachon with the putatively adaptive SNPs were similar, but provided greater resolution of stocks compared with a putatively neutral panel of 3911 SNPs or previous estimates with 14 microsatellites. A cline of increasing measures of genetic diversity from south to north was found in the adaptive panel, but not in the neutral markers (SNPs or microsatellites). This may indicate divergent selective pressures in differing freshwater and marine environments between regional eulachon populations and that these adaptive diversity patterns not seen with neutral markers could be a consideration when determining genetic boundaries for conservation purposes. Estimates of effective population size (N_e) were similar with the neutral SNP panel and microsatellites and may be utilized to monitor population status for eulachon where census sizes are difficult to obtain. Greater differentiation with the panel of putatively adaptive SNPs provided higher individual assignment accuracy compared to the neutral panel or microsatellites for stock identification purposes. This study presents the first SNPs that have been developed for eulachon, and analyses with these markers highlighted the importance of integrating genome‐wide neutral and adaptive genetic variation for the applications of conservation and management.     

Latitudinal divergence in a widespread amphibian: Contrasting patterns of neutral and adaptive genomic variation

Stochastic effects from demographic processes and selection are expected to shape the distribution of genetic variation in spatially heterogeneous environments. As the amount of genetic variation is central for long‐term persistence of populations, understanding how these processes affect variation over large‐scale geographical gradients is pivotal. We investigated the distribution of neutral and putatively adaptive genetic variation, and reconstructed demographic history in the moor frog (Rana arvalis) using 136 individuals from 15 populations along a 1,700‐km latitudinal gradient from northern Germany to northern Sweden. Using double digest restriction‐site associated DNA sequencing we obtained 27,590 single nucleotide polymorphisms (SNPs), and identified differentiation outliers and SNPs associated with growing season length. The populations grouped into a southern and a northern cluster, representing two phylogeographical lineages from different post‐glacial colonization routes. Hybrid index estimation and demographic model selection showed strong support for a southern and northern lineage and evidence of gene flow between regions located on each side of a contact zone. However, patterns of past gene flow over the contact zone differed between neutral and putatively adaptive SNPs. While neutral nucleotide diversity was higher along the southern than the northern part of the gradient, nucleotide diversity in differentiation outliers showed the opposite pattern, suggesting differences in the relative strength of selection and drift along the gradient. Variation associated with growing season length decreased with latitude along the southern part of the gradient, but not along the northern part where variation was lower, suggesting stronger climate‐mediated selection in the north. Outlier SNPs included loci involved in immunity and developmental processes.     

Genomic footprints of adaptation in a cooperatively breeding tropical bird across a vegetation gradient

Identifying the genetic basis of phenotypic variation and its relationship with the environment is key to understanding how local adaptations evolve. Such patterns are especially interesting among populations distributed across habitat gradients, where genetic structure can be driven by isolation by distance (IBD) and/or isolation by environment (IBE). Here, we used variation in ~1,600 high‐quality SNPs derived from paired‐end sequencing of double‐digest restriction site‐associated DNA (ddRAD‐Seq) to test hypotheses related to IBD and IBE in the Yucatan jay (Cyanocorax yucatanicus), a tropical bird endemic to the Yucatán Peninsula. This peninsula is characterized by a precipitation and vegetation gradient—from dry to evergreen tropical forests—that is associated with morphological variation in this species. We found a moderate level of nucleotide diversity (π = .008) and little evidence for genetic differentiation among vegetation types. Analyses of neutral and putatively adaptive SNPs (identified by complementary genome‐scan approaches) indicate that IBD is the most reliable explanation to account for frequency distribution of the former, while IBE has to be invoked to explain those of the later. These results suggest that selective factors acting along a vegetation gradient can promote local adaptation in the presence of gene flow in a vagile, nonmigratory and geographically restricted species. The putative candidate SNPs identified here are located within or linked to a variety of genes that represent ideal targets for future genomic surveys.     

Riverscape genomics of a threatened fish across a hydroclimatically heterogeneous river basin

Understanding how natural selection generates and maintains adaptive genetic diversity in heterogeneous environments is key to predicting the evolutionary response of populations to rapid environmental change. Detecting selection in complex spatial environments remains challenging, especially for threatened species where the effects of strong genetic drift may overwhelm signatures of selection. We carried out a basinwide riverscape genomic analysis in the threatened southern pygmy perch (Nannoperca australis), an ecological specialist with low dispersal potential. High‐resolution environmental data and 5162 high‐quality filtered SNPs were used to clarify spatial population structure and to assess footprints of selection associated with a steep hydroclimatic gradient and with human disturbance across the naturally and anthropogenically fragmented Murray–Darling Basin (Australia). Our approach included F_ST outlier tests to define neutral loci, and a combination of spatially explicit genotype–environment association analyses to identify candidate adaptive loci while controlling for the effects of landscape structure and shared population history. We found low levels of genetic diversity and strong neutral population structure consistent with expectations based on spatial stream hierarchy and life history. In contrast, variables related to precipitation and temperature appeared as the most important environmental surrogates for putatively adaptive genetic variation at both regional and local scales. Human disturbance also influenced the variation in candidate loci for adaptation, but only at a local scale. Our study contributes to understanding of adaptive evolution along naturally and anthropogenically fragmented ecosystems. It also offers a tangible example of the potential contributions of landscape genomics for informing in situ and ex situ conservation management of biodiversity.     

Seascape genomics reveals adaptive divergence in a connected and commercially important mollusc,the greenlip abalone (Haliotis laevigata), along a longitudinal environmental gradient

Populations of broadcast spawning marine organisms often have large sizes and are exposed to reduced genetic drift. Under such scenarios, strong selection associated with spatial environmental heterogeneity is expected to drive localized adaptive divergence, even in the face of connectivity. We tested this hypothesis using a seascape genomics approach in the commercially important greenlip abalone (Haliotis laevigata). We assessed how its population structure has been influenced by environmental heterogeneity along a zonal coastal boundary in southern Australia linked by strong oceanographic connectivity. Our data sets include 9,109 filtered SNPs for 371 abalones from 13 localities and environmental mapping across ~800 km. Genotype–environment association analyses and outlier tests defined 8,786 putatively neutral and 323 candidate adaptive loci. From a neutral perspective, the species is better represented by a metapopulation with very low differentiation (global F_ST = 0.0081) and weak isolation by distance following a stepping‐stone model. For the candidate adaptive loci, however, model‐based and model‐free approaches indicated five divergent population clusters. After controlling for spatial distance, the distribution of putatively adaptive variation was strongly correlated to selection linked to minimum sea surface temperature and oxygen concentration. Around 80 candidates were annotated to genes with functions related to high temperature and/or low oxygen tolerance, including genes that influence the resilience of abalone species found in other biogeographic regions. Our study includes a documented example about the uptake of genomic information in fisheries management and supports the hypothesis of adaptive divergence due to coastal environmental heterogeneity in a connected metapopulation of a broadcast spawner.     

Selection and sex‐biased dispersal in a coastal shark: the influence of philopatry on adaptive variation

Sex‐biased dispersal is expected to homogenize nuclear genetic variation relative to variation in genetic material inherited through the philopatric sex. When site fidelity occurs across a heterogeneous environment, local selective regimes may alter this pattern. We assessed spatial patterns of variation in nuclear‐encoded, single nucleotide polymorphisms (SNPs) and sequences of the mitochondrial control region in bonnethead sharks (Sphyrna tiburo), a species thought to exhibit female philopatry, collected from summer habitats used for gestation. Geographic patterns of mtDNA haplotypes and putatively neutral SNPs confirmed female philopatry and male‐mediated gene flow along the northeastern coast of the Gulf of Mexico. A total of 30 outlier SNP loci were identified; alleles at over half of these loci exhibited signatures of latitude‐associated selection. Our results indicate that in species with sex‐biased dispersal, philopatry can facilitate sorting of locally adaptive variation, with the dispersing sex facilitating movement of potentially adaptive variation among locations and environments.     

Landscape heterogeneity and local adaptation define the spatial genetic structure of Pacific salmon in a pristine environment

Identifying the spatial distribution of genetic variation across the landscape is an essential step in informing species conservation. Comparison of closely related and geographically overlapping species can be particularly useful in cases where landscape may similarly influence genetic structure. Congruent patterns among species highlight the importance that landscape heterogeneity plays in determining genetic structure whereas contrasting patterns emphasize differences in species-specific ecology and life-history or the importance of species-specific adaptation to local environments. We examined the interacting roles of demography and adaptation in determining spatial genetic structure in two closely related and geographically overlapping species in a pristine environment. Using single nucleotide polymorphism (SNP) loci exhibiting both neutral and putative adaptive variation, we evaluated the genetic structure of sockeye salmon in the Copper River, Alaska; these data were compared to existing data for Chinook salmon from the same region. Overall, both species exhibited patterns of isolation by distance; the spatial distribution of populations largely determined the distribution of genetic variation across the landscape. Further, both species exhibited largely congruent patterns of within- and among-population genetic diversity, highlighting the role that landscape heterogeneity and historical processes play in determining spatial genetic structure. Potential adaptive differences among geographically proximate sockeye salmon populations were observed when high F_ST outlier SNPs were evaluated in a landscape genetics context. Results were evaluated in the context of conservation efforts with an emphasis on reproductive isolation, historical processes, and local adaptation.     

Conserving adaptive genetic diversity in dynamic landscapes

Genetic variation supplies the raw material for adaptation, evolution and survival of populations and has therefore been a key focus of conservation biology ever since its foundation (Soulé 1985). In previous decades, the neutral component of genetic diversity (generated by mutation and shaped by drift) has been the subject of intense scientific research, fuelled by the increasing availability of molecular markers. On the other hand, the adaptive component of genetic diversity, which is shaped by the action of natural selection, has long remained elusive and difficult to assess, especially at small spatial or temporal scales (Ouborg et al. 2010). Fortunately, new technological and methodological developments now make it possible to identify loci in the genome that are influenced by selection, and thus to get a more complete view of genetic diversity. One article featured in this issue of Molecular Ecology is a good example of this recent breakthrough. Richter-Boix et al. (2011) examined a network of moor frog populations breeding in contrasting habitats in order to understand how landscape features influence patterns of genetic variation. They combined information from both neutral markers and loci putatively under selection to quantify the relative roles of selection and isolation in the evolution of fine-scale local adaptations in these populations. This study nicely illustrates how data on polymorphisms of neutral and adaptive loci can now be judiciously synthesized to help identify the best strategies for preserving adaptive variation, and more generally to enlighten conservation and population-management plans.     

Adaptive and neutral markers both show continent‐wide population structure of mountain pine beetle (Dendroctonus ponderosae)

Assessments of population genetic structure and demographic history have traditionally been based on neutral markers while explicitly excluding adaptive markers. In this study, we compared the utility of putatively adaptive and neutral single‐nucleotide polymorphisms (SNPs) for inferring mountain pine beetle population structure across its geographic range. Both adaptive and neutral SNPs, and their combination, allowed range‐wide structure to be distinguished and delimited a population that has recently undergone range expansion across northern British Columbia and Alberta. Using an equal number of both adaptive and neutral SNPs revealed that adaptive SNPs resulted in a stronger correlation between sampled populations and inferred clustering. Our results suggest that adaptive SNPs should not be excluded prior to analysis from neutral SNPs as a combination of both marker sets resulted in better resolution of genetic differentiation between populations than either marker set alone. These results demonstrate the utility of adaptive loci for resolving population genetic structure in a nonmodel organism.     

Broad Concordance in the Spatial Distribution of Adaptive and Neutral Genetic Variation across an Elevational Gradient in Deer Mice

When species are continuously distributed across environmental gradients, the relative strength of selection and gene flow shape spatial patterns of genetic variation, potentially leading to variable levels of differentiation across loci. Determining whether adaptive genetic variation tends to be structured differently than neutral variation along environmental gradients is an open and important question in evolutionary genetics. We performed exome-wide population genomic analysis on deer mice sampled along an elevational gradient of nearly 4,000 m of vertical relief. Using a combination of selection scans, genotype−environment associations, and geographic cline analyses, we found that a large proportion of the exome has experienced a history of altitude-related selection. Elevational clines for nearly 30% of these putatively adaptive loci were shifted significantly up- or downslope of clines for loci that did not bear similar signatures of selection. Many of these selection targets can be plausibly linked to known phenotypic differences between highland and lowland deer mice, although the vast majority of these candidates have not been reported in other studies of highland taxa. Together, these results suggest new hypotheses about the genetic basis of physiological adaptation to high altitude, and the spatial distribution of adaptive genetic variation along environmental gradients.     

Landscape features along migratory routes influence adaptive genomic variation in anadromous steelhead (Oncorhynchus mykiss)

Organisms typically show evidence of adaptation to features within their local environment. However, many species undergo long‐distance dispersal or migration across larger geographic regions that consist of highly heterogeneous habitats. Therefore, selection may influence adaptive genetic variation associated with landscape features at residing sites and along migration routes in migratory species. We tested for genomic adaptation to landscape features at natal spawning sites and along migration paths to the ocean of anadromous steelhead trout (Oncorhynchus mykiss) in the Columbia River Basin. Results from multivariate ordination, gene–environment association and outlier analyses using 24,526 single nucleotide polymorphisms (SNPs) provided evidence that adaptive allele frequencies were more commonly associated with landscape features along migration paths than features at natal sites (91.8% vs. 8.2% of adaptive loci, respectively). Among the 45 landscape variables tested, migration distance to the ocean and mean annual precipitation along migration paths were significantly associated with adaptive genetic variation in three distinct genetic groups. Additionally, variables such as minimum migration water temperature and mean migration slope were significant only in inland stocks of steelhead that migrate up to 1,200 km farther than those near the coast, indicating regional differences in migratory selective pressures. This study provides novel approaches for investigating migratory corridors and some of the first evidence that environment along migration paths can lead to substantial divergent selection. Consequently, our approach to understand genetic adaptation to migration conditions can be applied to other migratory species when migration or dispersal paths are generally known.     

Genome-Wide Patterns of Adaptation to Temperate Environments Associated with Transposable Elements in Drosophila

Investigating spatial patterns of loci under selection can give insight into how populations evolved in response to selective pressures and can provide monitoring tools for detecting the impact of environmental changes on populations. Drosophila is a particularly good model to study adaptation to environmental heterogeneity since it is a tropical species that originated in sub-Saharan Africa and has only recently colonized the rest of the world. There is strong evidence for the adaptive role of Transposable Elements (TEs) in the evolution of Drosophila, and TEs might play an important role specifically in adaptation to temperate climates. In this work, we analyzed the frequency of a set of putatively adaptive and putatively neutral TEs in populations with contrasting climates that were collected near the endpoints of two known latitudinal clines in Australia and North America. The contrasting results obtained for putatively adaptive and putatively neutral TEs and the consistency of the patterns between continents strongly suggest that putatively adaptive TEs are involved in adaptation to temperate climates. We integrated information on population behavior, possible environmental selective agents, and both molecular and functional information of the TEs and their nearby genes to infer the plausible phenotypic consequences of these insertions. We conclude that adaptation to temperate environments is widespread in Drosophila and that TEs play a significant role in this adaptation. It is remarkable that such a diverse set of TEs located next to a diverse set of genes are consistently adaptive to temperate climate-related factors. We argue that reverse population genomic analyses, as the one described in this work, are necessary to arrive at a comprehensive picture of adaptation.     

Genetic analysis of differentiation among breeding ponds reveals a candidate gene for local adaptation in Rana arvalis

One of the main questions in evolutionary and conservation biology is how geographical and environmental features of the landscape shape neutral and adaptive genetic variation in natural populations. The identification of genomic polymorphisms that account for adaptive variation can aid in finding candidate loci for local adaptation. Consequently, a comparison of spatial patterns in neutral markers and loci under selection may help disentangle the effects of gene flow, genetic drift and selection at the landscape scale. Many amphibians breed in wetlands, which differ in environmental conditions and in the degree of isolation, enhancing the potential for local adaptation. We used microsatellite markers to measure genetic differentiation among 17 local populations of Rana arvalis breeding in a network of wetlands. We found that locus RC08604 deviated from neutral expectations, suggesting that it is a good candidate for directional selection. We used a genetic network analysis to show that the allele distribution in this locus is correlated with habitat characteristics, whereas this was not the case at neutral markers that displayed a different allele distribution and population network in the study area. The graph approach illustrated the genomic heterogeneity (neutral loci vs. the candidate locus for directional selection) of gene exchange and genetic divergence among populations under directional selection. Limited gene flow between wetlands was only observed at the candidate genomic region under directional selection. RC08604 is partially located inside an up‐regulated thyroid‐hormone receptor (TRβ) gene coordinating the expression of other genes during metamorphosis and appears to be linked with variation in larval life‐history traits found among R. arvalis populations. We suggest that directional selection on genes coding larval life‐history traits is strong enough to maintain the divergence in these genomic regions, reducing the effective recombination of locally adapted alleles but not in other regions of the genome. Integrating this knowledge into conservation plans at the landscape scale will improve the design of management strategies to preserve adaptive genetic diversity in wetland networks.     

A conceptual framework for the spatial analysis of landscape genetic data

Understanding how landscape heterogeneity constrains gene flow and the spread of adaptive genetic variation is important for biological conservation given current global change. However, the integration of population genetics, landscape ecology and spatial statistics remains an interdisciplinary challenge at the levels of concepts and methods. We present a conceptual framework to relate the spatial distribution of genetic variation to the processes of gene flow and adaptation as regulated by spatial heterogeneity of the environment, while explicitly considering the spatial and temporal dynamics of landscapes, organisms and their genes. When selecting the appropriate analytical methods, it is necessary to consider the effects of multiple processes and the nature of population genetic data. Our framework relates key landscape genetics questions to four levels of analysis: (i) node-based methods, which model the spatial distribution of alleles at sampling locations (nodes) from local site characteristics; these methods are suitable for modeling adaptive genetic variation while accounting for the presence of spatial autocorrelation. (ii) Link-based methods, which model the probability of gene flow between two patches (link) and relate neutral molecular marker data to landscape heterogeneity; these methods are suitable for modeling neutral genetic variation but are subject to inferential problems, which may be alleviated by reducing links based on a network model of the population. (iii) Neighborhood-based methods, which model the connectivity of a focal patch with all other patches in its local neighborhood; these methods provide a link to metapopulation theory and landscape connectivity modeling and may allow the integration of node- and link-based information, but applications in landscape genetics are still limited. (iv) Boundary-based methods, which delineate genetically homogeneous populations and infer the location of genetic boundaries; these methods are suitable for testing for barrier effects of landscape features in a hypothesis-testing framework. We conclude that the power to detect the effect of landscape heterogeneity on the spatial distribution of genetic variation can be increased by explicit consideration of underlying assumptions and choice of an appropriate analytical approach depending on the research question.     

Insights into the neutral and adaptive processes shaping the spatial distribution of genomic variation in the economically important Moroccan locust (Dociostaurus maroccanus)

Understanding the processes that shape neutral and adaptive genomic variation is a fundamental step to determine the demographic and evolutionary dynamics of pest species. Here, we use genomic data obtained via restriction site‐associated DNA sequencing to investigate the genetic structure of Moroccan locust (Dociostaurus maroccanus) populations from the westernmost portion of the species distribution (Iberian Peninsula and Canary Islands), infer demographic trends, and determine the role of neutral versus selective processes in shaping spatial patterns of genomic variation in this pest species of great economic importance. Our analyses showed that Iberian populations are characterized by high gene flow, whereas the highly isolated Canarian populations have experienced strong genetic drift and loss of genetic diversity. Historical demographic reconstructions revealed that all populations have passed through a substantial genetic bottleneck around the last glacial maximum (~21 ka BP) followed by a sharp demographic expansion at the onset of the Holocene, indicating increased effective population sizes during warm periods as expected from the thermophilic nature of the species. Genome scans and environmental association analyses identified several loci putatively under selection, suggesting that local adaptation processes in certain populations might not be impeded by widespread gene flow. Finally, all analyses showed few differences between outbreak and nonoutbreak populations. Integrated pest management practices should consider high population connectivity and the potential importance of local adaptation processes on population persistence.     

Comparison of F(ST) outlier tests for SNP loci under selection

Genome scans with many genetic markers provide the opportunity to investigate local adaptation in natural populations and identify candidate genes under selection. In particular, SNPs are dense throughout the genome of most organisms and are commonly observed in functional genes making them ideal markers to study adaptive molecular variation. This approach has become commonly employed in ecological and population genetics studies to detect outlier loci that are putatively under selection. However, there are several challenges to address with outlier approaches including genotyping errors, underlying population structure and false positives, variation in mutation rate and limited sensitivity (false negatives). In this study, we evaluated multiple outlier tests and their type I (false positive) and type II (false negative) error rates in a series of simulated data sets. Comparisons included simulation procedures (FDIST2, ARLEQUIN v.3.5 and BAYESCAN) as well as more conventional tools such as global F(ST) histograms. Of the three simulation methods, FDIST2 and BAYESCAN typically had the lowest type II error, BAYESCAN had the least type I error and Arlequin had highest type I and II error. High error rates in Arlequin with a hierarchical approach were partially because of confounding scenarios where patterns of adaptive variation were contrary to neutral structure; however, Arlequin consistently had highest type I and type II error in all four simulation scenarios tested in this study. Given the results provided here, it is important that outlier loci are interpreted cautiously and error rates of various methods are taken into consideration in studies of adaptive molecular variation, especially when hierarchical structure is included.     

Longitudinal trends in climate drive flowering time clines in North American Arabidopsis thaliana

Introduced species frequently show geographic differentiation, and when differentiation mirrors the ancestral range, it is often taken as evidence of adaptive evolution. The mouse-ear cress (Arabidopsis thaliana) was introduced to North America from Eurasia 150-200 years ago, providing an opportunity to study parallel adaptation in a genetic model organism. Here, we test for clinal variation in flowering time using 199 North American (NA) accessions of A. thaliana, and evaluate the contributions of major flowering time genes FRI, FLC, and PHYC as well as potential ecological mechanisms underlying differentiation. We find evidence for substantial within population genetic variation in quantitative traits and flowering time, and putatively adaptive longitudinal differentiation, despite low levels of variation at FRI, FLC, and PHYC and genome-wide reductions in population structure relative to Eurasian (EA) samples. The observed longitudinal cline in flowering time in North America is parallel to an EA cline, robust to the effects of population structure, and associated with geographic variation in winter precipitation and temperature. We detected major effects of FRI on quantitative traits associated with reproductive fitness, although the haplotype associated with higher fitness remains rare in North America. Collectively, our results suggest the evolution of parallel flowering time clines through novel genetic mechanisms.