Post‐glacial colonization routes coincide with a life‐history breakpoint along a latitudinal gradient

Although adaptive divergence along environmental gradients has repeatedly been demonstrated, the role of post‐glacial colonization routes in determining phenotypic variation along gradients has received little attention. Here, we used a hierarchical Q_ST–F_ST approach to separate the roles of adaptive and neutral processes in shaping phenotypic variation in moor frog (Rana arvalis) larval life histories along a 1,700 km latitudinal gradient across northern Europe. This species has colonized Scandinavia via two routes with a contact zone in northern Sweden. By using neutral SNP and common garden phenotypic data from 13 populations at two temperatures, we showed that most of the variation along the gradient occurred between the two colonizing lineages. We found little phenotypic divergence within the lineages; however, all phenotypic traits were strongly diverged between the southern and northern colonization routes, with higher growth and development rates and larger body size in the north. The Q_ST estimates between the colonization routes were four times higher than F_ST, indicating a prominent role for natural selection. Q_ST within the colonization routes did not generally differ from F_ST, but we found temperature‐dependent adaptive divergence close to the contact zone. These results indicate that lineage‐specific variation can account for much of the adaptive divergence along a latitudinal gradient.     

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.     

Cold adaptation in the Asian tiger mosquito's native range precedes its invasion success in temperate regions

Adaptation to environmental conditions within the native range of exotic species can condition the invasion success of these species outside their range. The striking success of the Asian tiger mosquito, Aedes albopictus, to invade temperate regions has been attributed to the winter survival of diapause eggs in cold environments. In this study, we evaluate genetic polymorphisms (SNPs) and wing morphometric variation among three biogeographical regions of the native range of A. albopictus. Reconstructed demographic histories of populations show an initial expansion in Southeast Asia and suggest that marine regression during late Pleistocene and climate warming after the last glacial period favored expansion of populations in southern and northern regions, respectively. Searching for genomic signatures of selection, we identified significantly differentiated SNPs among which several are located in or within 20 kb distance from candidate genes for cold adaptation. These genes involve cellular and metabolic processes and several of them have been shown to be differentially expressed under diapausing conditions. The three biogeographical regions also differ for wing size and shape, and wing size increases with latitude supporting Bergmann's rule. Adaptive genetic and morphometric variation observed along the climatic gradient of A. albopictus native range suggests that colonization of northern latitudes promoted adaptation to cold environments prior to its worldwide invasion.     

Neutral and functionally important genes shed light on phylogeography and the history of high‐altitude colonization in a widespread New World duck

Phylogeographic studies often infer historical demographic processes underlying species distributions based on patterns of neutral genetic variation, but spatial variation in functionally important genes can provide additional insights about biogeographic history allowing for inferences about the potential role of adaptation in geographic range evolution. Integrating data from neutral markers and genes involved in oxygen (O₂)‐transport physiology, we test historical hypotheses about colonization and gene flow across low‐ and high‐altitude regions in the Ruddy Duck (Oxyura jamaicensis), a widely distributed species in the New World. Using multilocus analyses that for the first time include populations from the Colombian Andes, we also examined the hypothesis that Ruddy Duck populations from northern South America are of hybrid origin. We found that neutral and functional genes appear to have moved into the Colombian Andes from both North America and southern South America, and that high‐altitude Colombian populations do not exhibit evidence of adaptation to hypoxia in hemoglobin genes. Therefore, the biogeographic history of Ruddy Ducks is likely more complex than previously inferred. Our new data raise questions about the hypothesis that adaptation via natural selection to high‐altitude conditions through amino acid replacements in the hemoglobin protein allowed Ruddy Ducks to disperse south along the high Andes into southern South America. The existence of shared genetic variation with populations from both North America and southern South America as well as private alleles suggests that the Colombian population of Ruddy Ducks may be of old hybrid origin. This study illustrates the breadth of inferences one can make by combining data from nuclear and functionally important loci in phylogeography, and underscores the importance of complete range‐wide sampling to study species history in complex landscapes.     

The east‐west‐north colonization history of the Mediterranean and Europe by the coastal plant Carex extensa (Cyperaceae)

Coastal plants are ideal models for studying the colonization routes of species because of the simple linear distributions of these species. Carex extensa occurs mainly in salt marshes along the Mediterranean and European coasts. Variation in cpDNA sequences, amplified fragment length polymorphisms (AFLPs) and simple sequence repeats (SSRs) of 24 populations were analysed to reconstruct its colonization history. Phylogenetic relationships indicate that C. extensa together with the South American Carex vixdentata and the southern African Carex ecklonii form a monophyletic group of halophilic species. Analyses of divergence times suggest that early lineage diversification may have occurred between the late Miocene and the late Pliocene (Messinian crisis). Phylogenetic and network analyses of cpDNA variation revealed the monophyly of the species and an ancestral haplotype contained in populations of the eastern Mediterranean. The AFLP and SSR analyses support a pattern of variation compatible with these two lineages. These analyses also show higher levels of genetic diversity and differentiation in the eastern population group, which underwent an east‐to‐west Mediterranean colonization. Quaternary climatic oscillations appear to have been responsible for the split between these two lineages. Secondary contacts may have taken place in areas near the Ligurian Sea in agreement with the gene flow detected in Corsican populations. The AFLP and SSR data accord with the ‘tabula rasa’ hypothesis in which a recent and rapid colonization of northern Europe took place from the western Mediterranean after the Last Glacial Maximum. The unbalanced west‐east vs. west‐north colonization may be as a result of ‘high density blocking’ effect.     

Caucasus Mountains divide postulated postglacial colonization routes in the white‐breasted hedgehog,Erinaceus concolor

For many European species, the mountains of the Alps and the Pyrenees have acted as significant barriers to northwards colonization from southern glacial refugia. To the east, the Caucasus Mountains would seem to have been a similar barrier to the white‐breasted hedgehog (Erinaceus concolor). A deep divergence among hedgehog mitochondrial sequences to the north and south of the Caucasus Mountains suggests two colonization routes, originating from separate refugial regions and divided by this mountain barrier. From a Balkan refugium, hedgehogs have colonized northwards into Russia and to the northern foothills of the Caucasus Mountains. The origins of hedgehogs colonizing the southern parts of the Caucasus are not entirely clear, although fossil and climatic data suggest a glacial refugium on the southern shores of the Black Sea. Divergence within the southern group indicates a long‐standing fragmentation within such a refugium or the presence of further cryptic refugia in Turkey and the Near East. The Caucasus barrier would seem to have been an important factor in structuring the late Pleistocene distribution of species.     

High genomic diversity in the bank vole at the northern apex of a range expansion: The role of multiple colonizations and end‐glacial refugia

The history of repeated northern glacial cycling and southern climatic stability has long dominated explanations for how genetic diversity is distributed within temperate species in Eurasia and North America. However, growing evidence indicates the importance of cryptic refugia for northern colonization dynamics. An important geographic region to assess this is Fennoscandia, where recolonization at the end of the last glaciation was restricted to specific routes and temporal windows. We used genomic data to analyse genetic diversity and colonization history of the bank vole (Myodes glareolus) throughout Europe (>800 samples) with Fennoscandia as the northern apex. We inferred that bank voles colonized Fennoscandia multiple times by two different routes; with three separate colonizations via a southern land‐bridge route deriving from a “Carpathian” glacial refugium and one via a north‐eastern route from an “Eastern” glacial refugium near the Ural Mountains. Clustering of genome‐wide SNPs revealed high diversity in Fennoscandia, with eight genomic clusters: three of Carpathian origin and five Eastern. Time estimates revealed that the first of the Carpathian colonizations occurred before the Younger Dryas (YD), meaning that the first colonists survived the YD in Fennoscandia. Results also indicated that introgression between bank and northern red‐backed voles (Myodes rutilus) took place in Fennoscandia just after end‐glacial colonization. Therefore, multiple colonizations from the same and different cryptic refugia, temporal and spatial separations and interspecific introgression have shaped bank vole genetic variability in Fennoscandia. Together, these processes drive high genetic diversity at the apex of the northern expansion in this emerging model species.     

Colonization routes of Pinus sylvestris inferred from distribution of mitochondrial DNA variation

Understanding the present-day distribution of molecular variation requires knowledge about the history of the species. Past colonization routes and locations of refugia of Scots pine (Pinus sylvestris) were inferred from variation in mitochondrial DNA in material collected from 37 populations located in countries within, and immediately adjacent to the continent of Europe. Two mitochondrial regions, nad1 intron (exon B/C) and nad7 intron 1, were included in the study. Differentiation in maternally inherited mitochondria was high (G _ST′ = 0.824). Two new haplotypes were found at the nad7 intron 1. The occurrence of a 5-bp indel variant was restricted to the Turkish Kalabak population and a 32 bp only found in Central, Eastern, and Northern Europe. The complete absence of the 32-bp indel from the Mediterranean peninsulas supports the view that coniferous forests existed outside these areas during the last glacial maximum, and these populations contributed to the subsequent colonization of the northern parts of Europe. P. sylvestris shares features of its glacial and postglacial history with two other northern, cold-tolerant tree species, Picea abies and Betula sp. These three species differ from many other European trees for which pollen core and molecular evidence indicate colonization from southern refugia after the last glacial period.     

A highly pleiotropic amino acid polymorphism in the Drosophila insulin receptor contributes to life‐history adaptation

Finding the specific nucleotides that underlie adaptive variation is a major goal in evolutionary biology, but polygenic traits pose a challenge because the complex genotype–phenotype relationship can obscure the effects of individual alleles. However, natural selection working in large wild populations can shift allele frequencies and indicate functional regions of the genome. Previously, we showed that the two most common alleles of a complex amino acid insertion–deletion polymorphism in the Drosophila insulin receptor show independent, parallel clines in frequency across the North American and Australian continents. Here, we report that the cline is stable over at least a five‐year period and that the polymorphism also demonstrates temporal shifts in allele frequency concurrent with seasonal change. We tested the alleles for effects on levels of insulin signaling, fecundity, development time, body size, stress tolerance, and life span. We find that the alleles are associated with predictable differences in these traits, consistent with patterns of Drosophila life‐history variation across geography that likely reflect adaptation to the heterogeneous climatic environment. These results implicate insulin signaling as a major mediator of life‐history adaptation in Drosophila, and suggest that life‐history trade‐offs can be explained by extensive pleiotropy at a single locus.     

Detecting short spatial scale local adaptation and epistatic selection in climate‐related candidate genes in European beech (Fagus sylvatica) populations

Detecting signatures of selection in tree populations threatened by climate change is currently a major research priority. Here, we investigated the signature of local adaptation over a short spatial scale using 96 European beech (Fagus sylvatica L.) individuals originating from two pairs of populations on the northern and southern slopes of Mont Ventoux (south‐eastern France). We performed both single and multilocus analysis of selection based on 53 climate‐related candidate genes containing 546 SNPs. F_ST outlier methods at the SNP level revealed a weak signal of selection, with three marginally significant outliers in the northern populations. At the gene level, considering haplotypes as alleles, two additional marginally significant outliers were detected, one on each slope. To account for the uncertainty of haplotype inference, we averaged the Bayes factors over many possible phase reconstructions. Epistatic selection offers a realistic multilocus model of selection in natural populations. Here, we used a test suggested by Ohta based on the decomposition of the variance of linkage disequilibrium. Overall populations, 0.23% of the SNP pairs (haplotypes) showed evidence of epistatic selection, with nearly 80% of them being within genes. One of the between gene epistatic selection signals arose between an F_ST outlier and a nonsynonymous mutation in a drought response gene. Additionally, we identified haplotypes containing selectively advantageous allele combinations which were unique to high or low elevations and northern or southern populations. Several haplotypes contained nonsynonymous mutations situated in genes with known functional importance for adaptation to climatic factors.     

Tracing geographical patterns of population differentiation in a widespread mangrove gastropod: genetic and geometric morphometrics surveys along the eastern African coast

In the present study, we assessed the inter‐ and intrapopulation genetic and morphological variation of Cerithidea decollata along the eastern coast of Africa. The population structure of C. decollata along the latitudinal gradient was examined by sequencing 420 bp of the mitochondrial cytochrome c oxidase I (COI) gene in 172 snails from 29 sites, in a combined analysis with geometric morphometrics in 1799 snails from 32 sites. Analysis of molecular variance and spatial analysis of molecular variance showed a moderate spatial population differentiation from Kenya to the Republic of South Africa, suggesting genetic divergence between the northern, central, and southern regions. This structure appears to be the result of life‐history traits combined with oceanographic features. Haplotype network and mismatch analysis suggest a recent population expansion during the Holsteinian interglacial period in the northern region and several colonization events in the central and southern regions. The morphometric approach suggests that morphological variation in shell shape is somewhat independent of the genetic divergence, revealing an overlap of shape across the latitudinal gradient but significant differences among‐population at a local level. This may indicate that similar ecological pressures are acting along the coast, leading to the occurrence of similar morphological characters. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, 107 , 647–663.     

Following the rivers: historical reconstruction of California voles Microtus californicus (Rodentia: Cricetidae) in the deserts of eastern California

The California vole, Microtus californicus, restricted to habitat patches where water is available nearly year‐round, is a remnant of the mesic history of the southern Great Basin and Mojave deserts of eastern California. The history of voles in this region is a model for species‐edge population dynamics through periods of climatic change. We sampled voles from the eastern deserts of California and examined variation in the mitochondrial cytb gene, three nuclear intron regions, and across 12 nuclear microsatellite markers. Samples are allocated to two mitochondrial clades: one associated with southern California and the other with central and northern California. The limited mtDNA structure largely recovers the geographical distribution, replicated by both nuclear introns and microsatellites. The most remote population, Microtus californicus scirpensis at Tecopa near Death Valley, was the most distinct. This population shares microsatellite alleles with both mtDNA clades, and both its northern clade nuclear introns and southern clade mtDNA sequences support a hybrid origin for this endangered population. The overall patterns support two major invasions into the desert through an ancient system of riparian corridors along streams and lake margins during the latter part of the Pleistocene followed by local in situ divergence subsequent to late Pleistocene and Holocene drying events. Changes in current water resource use could easily remove California voles from parts of the desert landscape.     

Genetic and morphologic diversity of European fan palm (Chamaerops humilis L.) populations from different environments from Sicily

Chamaerops humilis is decreasing in abundance in Mediterranean Europe, which has induced the European community to call for its protection in Special Areas of Conservation. However, information about its genetic and morphological variability, which is crucial to the development of any conservation strategies, is insufficient. The present study aimed to investigate the genetic and morphological variability of C. humilis in Sicily, which was selected as a model because of the high number of dense populations. The relationships between morphological traits and climatic variables were studied to highlight patterns of adaptation to the environment, along with the genetic similarity among the populations. Ten natural populations were sampled, analyzed using 28 specifically designed SSR primers, and evaluated based on 29 morphological traits. The populations were clustered similarly based on genetic and morphological traits. Heterozygosity was high and inbreeding coefficients were low. These results, along with higher intra‐ than inter‐population differentiation, suggest that C. humilis populations in Sicily differentiated from a common ancestor and that inter‐population variation arose from secondary evolution processes induced by ecological adaptation. The correlations between climatic variables and morphological traits suggest that the morphological adaptation to arid environments depends more on summer temperatures than on evapotranspiration or rainfall and that autumn and winter temperatures are determinants of the species establishment at new sites. Considering the response of C. humilis to seasonal temperatures, the present results indicate this species as a candidate for tracking climatic changes in Europe. Further studies are needed to highlight the adaptation of C. humilis to cold environments. Palaeo‐climatological and ‐ecological studies could help clarify its strategies for the colonization of new sites. © 2014 The Linnean Society of London, Botanical Journal of the Linnean Society, 2014, 176 , 66–81.     

Complex patterns of local adaptation in heat tolerance in Drosophila simulans from eastern Australia

Latitudinal clines are considered a powerful means of investigating evolutionary responses to climatic selection in nature. However, most clinal studies of climatic adaptation in Drosophila have involved species that contain cosmopolitan inversion polymorphisms that show clinal patterns themselves, making it difficult to determine whether the traits or inversions are under selection. Further, although climatic selection is unlikely to act on only one life stage in metamorphic organisms, a few studies have examined clinal patterns across life stages. Finally, clinal patterns of heat tolerance may also depend on the assay used. To unravel these potentially confounding effects on clinal patterns of thermal tolerance, we examined adult and larval heat tolerance traits in populations of Drosophila simulans from eastern Australia using static and dynamic (ramping 0.06 °C min^?1) assays. We also used microsatellites markers to clarify whether demographic factors or selection are responsible for population differentiation along clines. Significant cubic clinal patterns were observed for adult static basal, hardened and dynamic heat knockdown time and static basal heat survival in larvae. In contrast, static, hardened larval heat survival increased linearly with latitude whereas no clinal association was found for larval ramping survival. Significant associations between adult and larval traits and climatic variables, and low population differentiation at microsatellite loci, suggest a role for climatic selection, rather than demographic processes, in generating these clinal patterns. Our results suggest that adaptation to thermal stress may be species and life‐stage specific, complicating our efforts to understand the evolutionary responses to selection for increasing thermotolerance.     

Phylogeography of Pulsatilla vernalis (L.) Mill. (Ranunculaceae): chloroplast DNA reveals two evolutionary lineages across central Europe and Scandinavia

Aim The aim of this study was to test hypotheses regarding some of the main phylogeographical patterns proposed for European plants, in particular the locations of glacial refugia, the post‐glacial colonization routes, and genetic affinities between southern (alpine) and northern (boreal) populations. Location The mountains of Europe (Alps, Balkans, Carpathians, Central Massif, Pyrenees, Scandinavian chain, Sudetes), and central European/southern Scandinavian lowlands. Methods As our model system we used Pulsatilla vernalis, a widely distributed European herbaceous plant occurring both in the high‐mountain environments of the Alps and other European ranges and in lowlands north of these ranges up to Scandinavia. Based on a distribution‐wide sampling of 61 populations, we estimated chloroplast DNA (cpDNA) variation along six regions using polymerase chain reaction–restriction fragment‐length polymorphisms (PCR–RFLPs) (trnH–trnK, trnK–trnK, trnC–trnD, psbC–trnS, psaA–trnS, trnL–trnF) and further sequencing of trnL–trnF and trnH–psbA. In addition, 11 samples of other European species of Pulsatilla were sequenced to survey the genus‐scale cpDNA variation. Results Eleven PCR–RFLP polymorphisms were detected in P. vernalis, revealing seven haplotypes. They formed two distinct genetic groups. Three haplotypes representing both groups dominated and were widely distributed across Europe, whereas the others were restricted to localized regions (central Alps, Tatras/Sudetes mountains) or single populations. Sequencing analysis confirmed the reliability of PCR–RFLPs and homology of haplotypes across their distribution. The chloroplast DNA variation across the section Pulsatilla was low, but P. vernalis did not share haplotypes with other species. Main conclusions The genetic distinctiveness of P. vernalis populations from the south‐western Alps with respect to other Alpine populations, as well as the affinities between the former populations and those from the eastern Pyrenees, is demonstrated, thus providing support for the conclusions of previous studies. Glacial refugia in the Dolomites are also suggested. Isolation is inferred for the high‐mountain populations from the Tatras and Sudetes; this is in contrast to the case for the Balkans, which harboured the common haplotype. Specific microsatellite variation indicates the occurrence of periglacial lowland refugia north of the Alps, acting as a source for the post‐glacial colonization of Scandinavia. The presence of different fixed haplotypes in eastern and western Scandinavia, however, suggests independent post‐glacial colonization of these two areas, with possible founder effects.     

Genetic consequences of glacial survival and postglacial colonization in Norway spruce: combined analysis of mitochondrial DNA and fossil pollen

Norway spruce (Picea abies [L.] Karst.) is a broadly distributed European conifer tree whose history has been intensively studied by means of fossil records to infer the location of full‐glacial refugia and the main routes of postglacial colonization. Here we use recently compiled fossil pollen data as a template to examine how past demographic events have influenced the species’ modern genetic diversity. Variation was assessed in the mitochondrial nad1 gene containing two minisatellite regions. Among the 369 populations (4876 trees) assayed, 28 mitochondrial variants were identified. The patterns of population subdivision superimposed on interpolated fossil pollen distributions indicate that survival in separate refugia and postglacial colonization has led to significant structuring of genetic variation in the southern range of the species. The populations in the northern range, on the other hand, showed a shallow genetic structure consistent with the fossil pollen data, suggesting that the vast northern range was colonized from a single refugium. Although the genetic diversity decreased away from the putative refugia, there were large differences between different colonization routes. In the Alps, the diversity decreased over short distances, probably as a result of population bottlenecks caused by the presence of competing tree species. In northern Europe, the diversity was maintained across large areas, corroborating fossil pollen data in suggesting that colonization took place at high population densities. The genetic diversity increased north of the Carpathians, probably as a result of admixture of expanding populations from two separate refugia.     

Climatic adaptation in an isolated and genetically impoverished amphibian population

The capacity of populations to respond adaptively to environmental change is essential for their persistence. Isolated populations often harbour reduced genetic variation, which is predicted to decrease adaptive potential, and can be detrimental under the current scenarios of global change. In this study, we examined climatic adaptation in larval life history traits in the pool frog Rana lessonae along a latitudinal gradient across the northern distribution area of the species, paying special attention to the isolated and genetically impoverished fringe populations in central Sweden. Larvae from eight populations within three geographic areas (Poland, Latvia and Sweden) were reared under three temperatures (19, 22 and 26°C) in a common garden laboratory experiment. We found clear evidence for latitudinal adaptation in R. lessonae populations, larvae from higher latitudes growing and developing faster than low‐latitude ones. Larvae from the Swedish populations were able to compensate for the effects of cooler temperatures and a shorter growth season with genetically higher growth and development rates (i.e. countergradient variation) in the two higher temperature treatments, but there was no difference among the populations at the lowest temperature treatment, which is likely to be close to the temperature limiting growth in R. lessonae. Our results demonstrate that isolated and genetically impoverished populations can be locally adapted, and identify the Swedish fringe populations as a significant conservation unit adapted to the northern environmental conditions.     

Let's mate here and now – seasonal constraints increase mating efficiency

1. Latitudinal climatic conditions shape the length of the mating season and could thus influence reproductive traits. Knowledge of how animals behave along latitudinal clines will increase understanding of the impact of climate on sexual selection and might help in the prediction of whether peripheral populations will spread or shrink in response to changes in climate. 2. This study investigated variation in the mating efficiency of a temperate insect, the emerald damselfly Lestes sponsa, under semi‐natural field conditions along a latitudinal gradient covering three regions of the species' distribution: south, central and north. 3. A comparison was done of the proportion of copulating males, the proportion of males that formed tandems but did not copulate (unsuccessful males), and the proportion of males that did not attempt to form a tandem (passive males) in these three regions. 4. It was found that the proportion of copulations was significantly higher at northern latitudes than in the southern and central regions. Southern latitudes had a higher proportion of successful copulations compared with central latitudes. The northern region had a significantly lower frequency of passive males. The southern region had an intermediate proportion of passive males, and the central region had the highest proportion. The proportion of unsuccessful males did not differ between regions. The population density across sites did not affect these results. 5. The study shows that damselflies inhabiting northern populations mate more intensively than individuals from southern and central populations. This suggests that more restrictive environmental conditions during a brief mating season select for higher mating efficiency.     

Selection for population‐specific adaptation shaped patterns of variation in the photoperiod pathway genes in Arabidopsis lyrata during post‐glacial colonization

Spatially varying selection can lead to population‐specific adaptation, which is often recognized at the phenotypic level; however, the genetic evidence is weaker in many groups of organisms. In plants, environmental shifts that occur due to colonization of a novel environment may require adaptive changes in the timing of growth and flowering, which are often governed by location‐specific environmental cues such as day length. We studied locally varying selection in 19 flowering time loci in nine populations of the perennial herb Arabidopsis lyrata, which has a wide but patchy distribution in temperate and boreal regions of the northern hemisphere. The populations differ in their recent population demographic and colonization histories and current environmental conditions, especially in the growing season length. We searched for population‐specific molecular signatures of directional selection by comparing a set of candidate flowering time loci with a genomic reference set within each population using multiple approaches and contrasted the patterns of different populations. The candidate loci possessed approximately 20% of the diversity of the reference loci. On average the flowering time loci had more rare alleles (a smaller Tajima's D) and an excess of highly differentiated sites relative to the reference, suggesting positive selection. The strongest signal of selection was detected in photoperiodic pathway loci in the colonizing populations of Northwestern Europe, whereas no evidence of positive selection was detected in the Central European populations. These findings emphasized the population‐specific nature of selection and suggested that photoperiodic adaptation was important during postglacial colonization of the species.