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.     

Length polymorphisms at two candidate genes explain variation of migratory behaviors in blackpoll warblers (Setophaga striata)

Migratory behaviors such as the timing and duration of migration are genetically inherited and can be under strong natural selection, yet we still know very little about the specific genes or molecular pathways that control these behaviors. Studies in candidate genes Clock and Adcyap1 have revealed that both of these loci can be significantly correlated with migratory behaviors in birds, though observed relationships appear to vary across species. We investigated geographic genetic structure of Clock and Adcyap1 in four populations of blackpoll warblers (Setophaga striata), a Neotropical–Nearctic migrant that exhibits geographic variation in migratory timing and duration across its boreal breeding distribution. Further, we used data on migratory timing and duration, obtained from light‐level geolocator trackers to investigate candidate genotype–phenotype relationships at the individual level. While we found no geographic structure in either candidate gene, we did find evidence that candidate gene lengths are correlated with five of the six migratory traits. Maximum Clock allele length was significantly and negatively associated with spring arrival date. Minimum Adcyap1 allele length was significantly and negatively associated with spring departure date and positively associated with fall arrival date at the wintering grounds. Additionally, we found a significant interaction between Clock and Adcyap1 allele lengths on both spring and fall migratory duration. Adcyap1 heterozygotes also had significantly shorter migration duration in both spring and fall compared to homozygotes. Our results support the growing body of evidence that Clock and Adcyap1 allele lengths are correlated with migratory behaviors in birds.     

Molecular markers provide evidence for a broad‐fronted recolonisation of the widespread calcareous grassland species Sanguisorba minor from southern and cryptic northern refugia

• Calcareous grasslands belong to the most species‐rich and endangered habitats in Europe. However, little is known about the origin of the species typically occurring in these grasslands. In this study we analysed the glacial and post‐glacial history of Sanguisorba minor, a typical plant species frequently occurring in calcareous grasslands.
• The study comprised 38 populations throughout the whole distribution range of the species across Europe. We used molecular markers (AFLP) and applied Bayesian cluster analysis as well as spatial principal components analysis (sPCA) to identify glacial refugia and post‐glacial migration routes to Central Europe.
• Our study revealed significant differences in the level of genetic variation and the occurrence of rare fragments within populations of S. minor and a distinct separation of eastern and western lineages. The analyses uncovered traditional southern but also cryptic northern refugia and point towards a broad fronted post‐glacial recolonisation.
• Based on these results we postulate that incomplete lineage sorting may have contributed to the detected pattern of genetic variation and that S. minor recolonised Central Europe post‐glacially from Iberia and northern glacial refugia in France, Belgium or Germany. Our results highlight the importance of refugial areas for the conservation of intraspecific variation in calcareous grassland species.

     

Climate variables explain neutral and adaptive variation within salmonid metapopulations: the importance of replication in landscape genetics

Understanding how environmental variation influences population genetic structure is important for conservation management because it can reveal how human stressors influence population connectivity, genetic diversity and persistence. We used riverscape genetics modelling to assess whether climatic and habitat variables were related to neutral and adaptive patterns of genetic differentiation (population‐specific and pairwise F_ST) within five metapopulations (79 populations, 4583 individuals) of steelhead trout (Oncorhynchus mykiss) in the Columbia River Basin, USA. Using 151 putatively neutral and 29 candidate adaptive SNP loci, we found that climate‐related variables (winter precipitation, summer maximum temperature, winter highest 5% flow events and summer mean flow) best explained neutral and adaptive patterns of genetic differentiation within metapopulations, suggesting that climatic variation likely influences both demography (neutral variation) and local adaptation (adaptive variation). However, we did not observe consistent relationships between climate variables and F_ST across all metapopulations, underscoring the need for replication when extrapolating results from one scale to another (e.g. basin‐wide to the metapopulation scale). Sensitivity analysis (leave‐one‐population‐out) revealed consistent relationships between climate variables and F_ST within three metapopulations; however, these patterns were not consistent in two metapopulations likely due to small sample sizes (N = 10). These results provide correlative evidence that climatic variation has shaped the genetic structure of steelhead populations and highlight the need for replication and sensitivity analyses in land and riverscape genetics.     

Clock gene polymorphism,migratory behaviour and geographic distribution: a comparative study of trans‐Saharan migratory birds

Migratory behaviour is controlled by endogenous circannual rhythms that are synchronized by external cues, such as photoperiod. Investigations on the genetic basis of circannual rhythmicity in vertebrates have highlighted that variation at candidate ‘circadian clock’ genes may play a major role in regulating photoperiodic responses and timing of life cycle events, such as reproduction and migration. In this comparative study of 23 trans‐Saharan migratory bird species, we investigated the relationships between species‐level genetic variation at two candidate genes, Clock and Adcyap1, and species’ traits related to migration and geographic distribution, including timing of spring migration across the Mediterranean Sea, migration distance and breeding latitude. Consistently with previous evidence showing latitudinal clines in ‘circadian clock’ genotype frequencies, Clock allele size increased with breeding latitude across species. However, early‐ and late‐migrating species had similar Clock allele size. Species migrating over longer distances, showing delayed spring migration and smaller phenotypic variance in spring migration timing, had significantly reduced Clock (but not Adcyap1) gene diversity. Phylogenetic confirmatory path analysis suggested that migration date and distance were the most important variables directly affecting Clock gene diversity. Hence, our study supports the hypothesis that Clock allele size increases poleward as a consequence of adaptation to the photoperiodic regime of the breeding areas. Moreover, we show that long‐distance migration is associated with lower Clock diversity, coherently with strong stabilizing selection acting on timing of life cycle events in long‐distance migratory species, likely resulting from the time constraints imposed by late spring migration.     

Polymorphism at the Clock gene predicts phenology of long‐distance migration in birds

Dissecting phenotypic variance in life history traits into its genetic and environmental components is at the focus of evolutionary studies and of pivotal importance to identify the mechanisms and predict the consequences of human‐driven environmental change. The timing of recurrent life history events (phenology) is under strong selection, but the study of the genes that control potential environmental canalization in phenological traits is at its infancy. Candidate genes for circadian behaviour entrained by photoperiod have been screened as potential controllers of phenological variation of breeding and moult in birds, with inconsistent results. Despite photoperiodic control of migration is well established, no study has reported on migration phenology in relation to polymorphism at candidate genes in birds. We analysed variation in spring migration dates within four trans‐Saharan migratory species (Luscinia megarhynchos; Ficedula hypoleuca; Anthus trivialis; Saxicola rubetra) at a Mediterranean island in relation to Clock and Adcyap1 polymorphism. Individuals with larger number of glutamine residues in the poly‐Q region of Clock gene migrated significantly later in one or, respectively, two species depending on sex and whether the within‐individual mean length or the length of the longer Clock allele was considered. The results hinted at dominance of the longer Clock allele. No significant evidence for migration date to covary with Adcyap1 polymorphism emerged. This is the first evidence that migration phenology is associated with Clock in birds. This finding is important for evolutionary studies of migration and sheds light on the mechanisms that drive bird phenological changes and population trends in response to climate change.     

Conservation genetics of the pond bat (Myotis dasycneme) with special focus on the populations in northwestern Germany and in Jutland,Denmark

Conservation genetics is important in the management of endangered species, helping to understand their connectivity and long‐term viability, thus identifying populations of importance for conservation. The pond bat (Myotis dasycneme) is a rare species classified as “Near Threatened” with a wide but patchy Palearctic distribution. A total of 277 samples representing populations in Denmark, Germany, Latvia, Hungary, and Russia were used in the genetic analyses; 224 samples representing Denmark, Germany, and Russia were analyzed at 10 microsatellite loci; 241 samples representing all areas were analyzed using mitochondrial D‐loop and cytochrome B sequences. A Bayesian clustering approach revealed two poorly resolved clusters, one representing the Danish and German groups and the other the Russian group. However, significantly different pairwise F_ST and D_EST estimates were observed between the Danish and German groups and between the Danish and Russian groups suggesting a recent population structure. These conflicting results might be attributed to the effect of migration or low resolution due to the number of microsatellite markers used. After concatenating the two mitochondrial sequences, analysis detected significant genetic differentiation between all populations, probably due to genetic drift combined with a founder event. The phylogenetic tree suggested a closer relationship between the Russian and Northern European populations compared to the Hungarian population, implying that the latter belongs to an older ancestral population. This was supported by the observed haplotype network and higher nucleotide diversity in this population. The genetic structuring observed in the Danish/German pond bat stresses the need for a cross‐border management between the two countries. Further, the pronounced mtDNA structuring, together with the indicated migration between nearby populations suggest philopatric female behavior but male migration, emphasizes the importance of protecting suitable habitat mosaics to maintain a continuum of patches with dense pond bat populations across the species' distribution range.     

Endemic or introduced? Phylogeography of Asparagopsis (Florideophyceae) in Australia reveals multiple introductions and a new mitochondrial lineage

The red seaweed Asparagopsis taxiformis embodies five cryptic mitochondrial lineages (lineage 1–5) introduced worldwide as a consequence of human mediated transport and climate change. We compared globally collected mitochondrial cox2‐3 intergenic spacer sequences with sequences produced from multiple Australian locations and South Korea to identify Asparagopsis lineages and to reveal cryptic introductions. We report A. taxiformis lineage 4 from Cocos (Keeling) Islands, Australia, and the highly invasive Indo‐Pacific Mediterranean lineage 2 from South Korea and Lord Howe Island, Australia. Phylogeographic analysis showed a clear haplotype and geographic separation between western Australian and Great Barrier Reef (GBR) isolates belonging to the recently described lineage 5. The same lineage, however, was characterized by a substantial genetic and geographic break between the majority of Australian specimens and Asparagopsis collections from South Solitary Island, Southern GBR, Lord Howe Island, Kermadec Islands, Norfolk Island, New Caledonia and French Polynesia. The disjunct geographic distribution and sequence divergence between these two groups supports the recognition of a sixth cryptic A. taxiformis mitochondrial lineage. As climatic changes accelerate the relocation of biota and offer novel niches for colonization, periodic surveys for early detection of cryptic invasive seaweeds will be critical in determining whether eradication or effective containment of the aliens are feasible.     

Adcyap1 polymorphism covaries with breeding latitude in a Nearctic migratory songbird,the Wilson's warbler (Cardellina pusilla)

Understanding the genetic background of complex behavioral traits, showing multigenic control and extensive environmental effects, is a challenging task. Among such traits, migration is known to show a large additive genetic component. Yet, the identification of specific genes or gene regions explaining phenotypic variance in migratory behavior has received less attention. Migration ultimately depends on seasonal cycles, and polymorphism at phenological candidate genes may underlie variation in timing of migration or other aspects of migratory behavior. In this study of a Nearctic–Neotropical migratory songbird, the Wilson's warbler (Cardellina pusilla), we investigated the association between polymorphism at two phenological candidate genes, Clock and Adcyap1, and two aspects of the migratory phenotype, timing of spring migration through a stopover site and inferred latitude of the breeding destination. The breeding destination of migrating individuals was identified using feather deuterium ratio (δ²H), which reliably reflects breeding latitude throughout the species' western breeding range. Ninety‐eight percent of the individuals were homozygous at Clock, and the rare heterozygotes did not deviate from homozygous migration phenology. Adcyap1 was highly polymorphic, and allele size was not significantly associated with migration date. However, Adcyap1 allele size significantly positively predicted the inferred breeding latitude of males but not of females. Moreover, we found a strong positive association between inferred breeding latitude and Adcyap1 allele size in long‐distance migrating birds from the northern sector of the breeding range (western Canada), while this was not the case in short‐distance migrating birds from the southern sector of the breeding range (coastal California). Our findings support previous evidence for a role of Adcyap1 in shaping the avian migratory phenotype, while highlighting that patterns of phenological candidate gene–phenotype associations may be complex, significantly varying between geographically distinct populations and even between the sexes.     

Local adaptation through genetic differentiation in highly fragmented Tilia cordata populations

We assessed the level of geographic differentiation of Tilia cordata in Denmark based on tests of 91 trees selected from 12 isolated populations. We used quantitative analysis of spring phenology and population genetic analysis based on SSR markers to infer the likely historical genetic processes within and among populations. High genetic variation within and among populations was observed in spring phenology, which correlated with spring temperatures at the origin of the tested T. cordata trees. The population genetic analysis revealed significant differentiation among the populations, but with no clear sign of isolation by distance. We infer the findings as indications of ongoing fine scale selection in favor of local growth conditions made possible by limited gene flow among the small and fragmented populations. This hypothesis fits well with reports of limited fruiting in the investigated Danish T. cordata populations, while the species is known for its ability to propagate vegetatively by root suckers. Our results suggest that both divergent selection and genetic drift may have played important roles in forming the genetic patterns of T. cordata at its northern distribution limit. However, we also speculate that epigenetic mechanism arising from the original population environment could have created similar patterns in regulating the spring phenology.     

Candidate gene polymorphisms for behavioural adaptations during urbanization in blackbirds

Successful urban colonization by formerly rural species represents an ideal situation in which to study adaptation to novel environments. We address this issue using candidate genes for behavioural traits that are expected to play a role in such colonization events. We identified and genotyped 16 polymorphisms in candidate genes for circadian rhythms, harm avoidance and migratory and exploratory behaviour in 12 paired urban and rural populations of the blackbird Turdus merula across the Western Palaearctic. An exonic microsatellite in the SERT gene, a candidate gene for harm avoidance behaviour, exhibited a highly significant association with habitat type in an analysis conducted across all populations. Genetic divergence at this locus was consistent in 10 of the 12 population pairs; this contrasts with previously reported stochastic genetic divergence between these populations at random markers. Our results indicate that behavioural traits related to harm avoidance and associated with the SERT polymorphism experience selection pressures during most blackbird urbanization events. These events thus appear to be influenced by homogeneous adaptive processes in addition to previously reported demographic founder events.     

Population genetics of Wolbachia‐infected,parthenogenetic and uninfected,sexual populations of Tetrastichus coeruleus (Hymenoptera: Eulophidae)

Wolbachia are endosymbiotic bacteria known to manipulate the reproduction of their hosts. These manipulations are expected to have consequences on the population genetics of the host, such as heterozygosity levels, genetic diversity and gene flow. The parasitoid wasp Tetrastichus coeruleus has populations that are infected with parthenogenesis‐inducing Wolbachia and populations that are not infected. We studied the population genetics of T. coeruleus between and within Wolbachia‐infected and uninfected populations, using nuclear microsatellites and mitochondrial DNA. We expected reduced genetic diversity in both DNA types in infected populations. However, migration and gene flow could introduce new DNA variants into populations. We therefore paid special attention to individuals with unexpected (genetic) characteristics. Based on nuclear and mitochondrial DNA, two genetic clusters were evident: a thelytokous cluster containing all Wolbachia‐infected, parthenogenetic populations and an arrhenotokous cluster containing all uninfected, sexual populations. Nuclear and mitochondrial DNA did not exhibit concordant patterns of variation, although there was reduced genetic diversity in infected populations for both DNA types. Within the thelytokous cluster, there was nuclear DNA variation, but no mitochondrial DNA variation. This nuclear DNA variation may be explained by occasional sex between infected females and males, by horizontal transmission of Wolbachia, and/or by novel mutations. Several females from thelytokous populations were uninfected and/or heterozygous for microsatellite loci. These unexpected characteristics may be explained by migration, by inefficient transmission of Wolbachia, by horizontal transmission of Wolbachia, and/or by novel mutations. However, migration has not prevented the build‐up of considerable genetic differentiation between thelytokous and arrhenotokous populations.     

Microsatellite DNA Population Structure and Stock Identification of Steelhead Trout (<Emphasis Type="BoldItalic">Oncorhynchus mykiss</Emphasis>) in the Nass and Skeena Rivers in Northern British Columbia

Population structure and the application to genetic stock identification for steelhead (Oncorhynchus mykiss) in the Nass and Skeena Rivers in northern British Columbia was examined using microsatellite markers. Variation at 8 microsatellite loci (Oki200, Omy77, Ots1, Ots3, Ssa85, Ots100, Ots103, and Ots108) was surveyed for approximately 930 steelhead from 7 populations in the Skeena River drainage and 850 steelhead from 10 populations in the Nass River drainage, as well as 1550 steelhead from test fisheries near the mouth of each river. Differentiation among populations within rivers accounted for about 1.9 times the variation observed among years within populations, with differences between drainages less than variation among populations within drainages. In the Nass River, winter-run populations formed a distinct group from the summer-run populations. Winter-run populations were not assessed in the Skeena River watershed. Simulated mixed-stock samples suggested that variation at the 8 microsatellite loci surveyed should provide relatively accurate and precise estimates of stock composition for fishery management applications within drainages. In the Skeena River drainage in 1998, Babine River (27%) and Bulkley drainage populations (31%) comprised the main components of the returns. For the Nass River in 1998 steelhead returning to Bell-Irving River were estimated to have comprised 39% of the fish sampled in the test fishery, with another 27% of the returns estimated to be derived from Cranberry River. The survey of microsatellite variation did not reveal enough differentiation between Nass River and Skeena River populations to be applied confidently in estimation of stock composition in marine fisheries at this time. Received January 14, 2000; accepted July 13, 2000.     

A clever solution to a vexing problem

F_ST (as well as related measures such as G_ST) has long been used both as a measure of the relative amount of genetic variation between populations and as an indicator of the amount of gene flow among populations. Unfortunately, F_ST and its clones are also sensitive to mutation, particularly when the mutation rate per locus is greater than the migration rate among populations. Relatively high mutation rates cause estimates of F_ST and G_ST to be much lower than researchers sometimes expect, when migration rates are low in the studied species. Several recent suggestions for dealing with this problem have been unsatisfactory for one reason or another, and no general solution exists (if we are not to abandon these otherwise useful measures of differentiation). In an important article in this issue, Jinliang Wang (2015) shows that it is possible to identify whether the genetic markers in a given study are likely to give estimates of F_ST that are strongly affected by mutation. The proposed test is simple and elegant, and with it, molecular ecologists can determine whether the F_ST from their makers can be depended on for further inference about their species’ genome and the demographic forces which shaped its patterns.     

Peripheral genetic structure of Helicoverpa zea indicates asymmetrical panmixia

Seasonal climatic shifts create peripheral habitats that alternate between habitable and uninhabitable for migratory species. Such dynamic peripheral habitats are potential sites where migratory species could evolve high genetic diversity resulting from convergence of immigrants from multiple regionally distant areas. Migrant populations of Helicoverpa zea (Boddie) captured during two different seasons were assessed for genetic structure using microsatellite markers and for host plant type using stable carbon isotope analysis. Individuals (N = 568) were genotyped and divided into 13 putative populations based on collection site and time. Fixation indices (F‐statistics), analysis of molecular variance (AMOVA), and discriminant analysis of principal components (DAPC) were used to examine within and among population genetic variation. Mean number of alleles per locus was 10.25 (± 3.2 SD), and allelic richness ranged from 2.38 to 5.13 (± 3.2 SD). The observed and expected heterozygosity ranged from 0.07 to 0.48 and 0.08 to 0.62, respectively. Low F_ST (0.01 to 0.02) and high F_IS (0.08 to 0.33) values suggest captured migrants originated from breeding populations with different allele frequencies. We postulate that high genetic diversity within migrant populations and low genetic differentiation among migrant populations of H. zea are the result of asymmetrical immigration due to the high dispersal and reproductive behavior of H. zea, which may hinder the adaptation and establishment of H. zea to peripheral habitat. These findings highlight the importance of assessing peripheral population structure in relation to ecological and evolutionary dynamics of this and other highly reproductive and dispersive species.     

Using soil seed banks to assess temporal patterns of genetic variation in invasive plant populations

Most research on the genetics of invasive plant species has focused on analyzing spatial differences among existing populations. Using a long‐established Gunnera tinctoria population from Ireland, we evaluated the potential of using plants derived from seeds associated with different soil layers to track genetic variation through time. This species and site were chosen because (1) G. tinctoria produces a large and persistent seed bank; (2) it has been present in this locality, Sraheens, for ~90 years; (3) the soil is largely undisturbed; and (4) the soil's age can be reliably determined radiometrically at different depths. Amplified fragment length polymorphic markers (AFLPs) were used to assess differences in the genetic structure of 75 individuals sampled from both the standing population and from four soil layers, which spanned 18 cm (estimated at ~90 years based on ²¹⁰Pb and ¹³⁷Cs dating). While there are difficulties in interpreting such data, including accounting for the effects of selection, seed loss, and seed migration, a clear pattern of lower total allele counts, percentage polymorphic loci, and genetic diversity was observed in deeper soils. The greatest percentage increase in the measured genetic variables occurred prior to the shift from the lag to the exponential range expansion phases and may be of adaptive significance. These findings highlight that seed banks in areas with long‐established invasive populations can contain valuable genetic information relating to invasion processes and as such, should not be overlooked.     

A comparative plastomics approach reveals available molecular markers for the phylogeographic study of Dendrobium huoshanense,an endangered orchid with extremely small populations

Comparative plastomics approaches have been used to identify available molecular markers for different taxonomic level studies of orchid species. However, the adoption of such methods has been largely limited in phylogeographic studies. Therefore, in this study, Dendrobium huoshanense, an endangered species with extremely small populations, was used as a model system to test whether the comparative plastomic approaches could screen available molecular markers for the phylogeographic study. We sequenced two more plastomes of D. huoshanense and compared them with our previously published one. A total of 27 mutational hotspot regions and six polymorphic cpSSRs have been screened for the phylogeographic studies of D. huoshanense. The cpDNA haplotype data revealed that the existence of haplotype distribution center was located in Dabieshan Mts. (Huoshan). The genetic diversity and phylogenetic analyses showed that the populations of D. huoshanense have been isolated and evolved independently for long period. On the contrary, based on cpSSR data, the genetic structure analysis revealed a mixed structure among the populations in Anhui and Jiangxi province, which suggested that the hybridization or introgression events have occurred among the populations of D. huoshanense. These results indicated that human activities have played key roles in shaping the genetic diversity and distributional patterns of D. huoshanense. According to our results, both two markers showed a high resolution for the phylogeographic studies of D. huoshanense. Therefore, we put forth that comparative plastomic approaches could revealed available molecular markers for phylogeographic study, especially for the species with extremely small populations.     

Variation in lipid synthesis,but genetic homogeneity,among Leptopilina parasitic wasp populations

Lipid synthesis can have a major effect on survival and reproduction, yet most insect parasitoids fail to synthesize lipids. For parasitic wasps in the genus Leptopilina, however, studies have suggested that there is intraspecific variation in the ability for lipid synthesis. These studies were performed on only few populations, and a large‐scale investigation of both lipogenic ability and population genetic structure is now needed. Here, we first examined lipogenic ability of nine Leptopilina heterotoma populations collected in 2013 and found that five of nine populations synthesized lipids. The 2013 populations could not be used to determine genetic structure; hence, we obtained another 20 populations in 2016 that were tested for lipogenic ability. Thirteen of 20 populations (all Leptopilina heterotoma) were then used to determine the level of genetic differentiation (i.e., haplotype and nucleotide diversity) by sequencing neutral mitochondrial (COI) and nuclear (ITS2) markers. None of the 2016 populations synthesized lipids, and no genetic differentiation was found. Our results did reveal a nearly twofold increase in mean wasp lipid content at emergence in populations obtained in 2016 compared to 2013. We propose that our results can be explained by plasticity in lipid synthesis, where lipogenic ability is determined by environmental factors, such as developmental temperature and/or the amount of lipids carried over from the host.     

Higher genetic diversity on mountain tops: the role of historical and contemporary processes in shaping genetic variation in the bank vole

Glacial phases during the Pleistocene caused remarkable changes in species range distributions, with inevitable genetic consequences. Specifically, during interglacial phases, when the ice melted and new habitats became suitable again, species could recolonize regions that were previously covered by ice, such as high latitudes and elevations. Based on theoretical models and empirical data, a decrease in genetic variation is predicted along recolonization routes as a result of the consecutive founder effects that characterize the recolonization process. In the present study, we assessed the relative importance of historical and contemporary processes in shaping genetic diversity and differentiation of bank vole (Myodes glareolus) populations at different elevations in the Swiss Alps. By contrast to expectations, we found that genetic variation increased with elevation. Estimates of recent migration rates and a contrasting pattern of genetic differentiation observed at the mitochondrial cytochrome b gene and nuclear microsatellites support the hypothesis that higher genetic diversity at high elevation results from contemporary gene flow. Although historical recolonization processes can have marked effects on the genetic structure of populations, the present study provides an example where contemporary processes along an environmental gradient can reverse predicted patterns of genetic variation.