Seasonal Adaptation: Geographic Photoperiod-Temperature Patterns Explain Genetic Variation in the Common Vole Tsh Receptor

The vertebrate photoperiodic neuroendocrine system uses the photoperiod as a proxy to time the annual rhythms in reproduction. The thyrotropin receptor (TSHR) is a key protein in the mammalian seasonal reproduction pathway. Its abundance and function can tune sensitivity to the photoperiod. To investigate seasonal adaptation in mammals, the hinge region and the first part of the transmembrane domain of the Tshr gene were sequenced for 278 common vole (Microtus arvalis) specimens from 15 localities in Western Europe and 28 localities in Eastern Europe. Forty-nine single nucleotide polymorphisms (SNPs; twenty-two intronic and twenty-seven exonic) were found, with a weak or lack of correlation with pairwise geographical distance, latitude, longitude, and altitude. By applying a temperature threshold to the local photoperiod-temperature ellipsoid, we obtained a predicted critical photoperiod (pCPP) as a proxy for the spring onset of local primary food production (grass). The obtained pCPP explains the distribution of the genetic variation in Tshr in Western Europe through highly significant correlations with five intronic and seven exonic SNPs. The relationship between pCPP and SNPs was lacking in Eastern Europe. Thus, Tshr, which plays a pivotal role in the sensitivity of the mammalian photoperiodic neuroendocrine system, was targeted by natural selection in Western European vole populations, resulting in the optimized timing of seasonal reproduction.

Pan‐European phylogeography of the European roe deer ( Capreolus capreolus )

To provide the most comprehensive picture of species phylogeny and phylogeography of European roe deer (Capreolus capreolus), we analyzed mtDNA control region (610 bp) of 1469 samples of roe deer from Central and Eastern Europe and included into the analyses additional 1541 mtDNA sequences from GenBank from other regions of the continent. We detected two mtDNA lineages of the species: European and Siberian (an introgression of C. pygargus mtDNA into C. capreolus). The Siberian lineage was most frequent in the eastern part of the continent and declined toward Central Europe. The European lineage contained three clades (Central, Eastern, and Western) composed of several haplogroups, many of which were separated in space. The Western clade appeared to have a discontinuous range from Portugal to Russia. Most of the haplogroups in the Central and the Eastern clades were under expansion during the Weichselian glacial period before the Last Glacial Maximum (LGM), while the expansion time of the Western clade overlapped with the Eemian interglacial. The high genetic diversity of extant roe deer is the result of their survival during the LGM probably in a large, contiguous range spanning from the Iberian Peninsula to the Caucasus Mts and in two northern refugia.

A dynamic history of admixture from Mediterranean and Carpathian glacial refugia drives genomic diversity in the bank vole

Understanding the historical contributions of differing glacial refugia is key to evaluating the roles of microevolutionary forces, such as isolation, introgression, and selection in shaping genomic diversity in present-day populations. In Europe, where both Mediterranean and extra-Mediterranean (e.g., Carpathian) refugia of the bank vole (Clethrionomys glareolus) have been identified, mtDNA indicates that extra-Mediterranean refugia were the main source of colonization across the species range, while Mediterranean peninsulas harbor isolated, endemic lineages. Here, we critically evaluate this hypothesis using previously generated genomic data (>6,000 SNPs) for over 800 voles, focusing on genomic contributions to bank voles in central Europe, a key geographic area in considering range-wide colonization. The results provide clear evidence that both extra-Mediterranean (Carpathian) and Mediterranean (Spanish, Calabrian, and Balkan) refugia contributed to the ancestry and genomic diversity of bank vole populations across Europe. Few strong barriers to dispersal and frequent admixture events in central Europe have led to a prominent mid-latitude peak in genomic diversity. Although the genomic contribution of the centrally located Carpathian refugium predominates, populations in different parts of Europe have admixed origins from Mediterranean (28%-47%) and the Carpathian (53%-72%) sources. We suggest that the admixture from Mediterranean refugia may have provisioned adaptive southern alleles to more northern populations, facilitating the end-glacial spread of the admixed populations and contributing to increased bank vole diversity in central Europe. This study adds critical details to the complex end-glacial colonization history of this well-studied organism and underscores the importance of genomic data in phylogeographic interpretation.

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.

Phylogeographical structure of the pygmy shrew: revisiting the roles of southern and northern refugia in Europe

Southern and northern glacial refugia are considered paradigms that explain the complex phylogeographical patterns and processes of European biota. Here, we provide a revisited statistical phylogeographical analysis of the pygmy shrew Sorex minutus Linnaeus, 1766 (Eulipotyphla, Soricidae), examining its genetic diversity, genetic differentiation and demographic history in the Mediterranean peninsulas and in Western and Central Europe. The results showed support for genetically distinct and diverse phylogeographical groups consistent with southern and northern glacial refugia, as expected from previous studies. We also identified geographical barriers concordant with glaciated mountain ranges during the Last Glacial Maximum (LGM), early diversification events dated between the Late Pleistocene and Early Holocene for the main phylogeographical groups, and recent (post-LGM) patterns of demographic expansions. This study is the most comprehensive investigation of this species to date, and the results have implications for the conservation of intraspecific diversity and the preservation of the evolutionary potential of S. minutus.

New record of Microtusmystacinus in eastern Kazakhstan: phylogeographical considerations

The Eastern European vole (Microtusmystacinus) is an arvicoline rodent distributed across northern and eastern Europe, the Balkans, Turkey, Armenia, NW and N Iran, Russia as far east as the Tobol River in W Siberia, and W and N Kazakhstan. We present a novel records from eastern Kazakhstan (the village of Dzhambul – 49°14'21.3"N, 86°18'29.9"E and the village of Sekisovka – 50°21'9.18"N, 82°35'46.5"E) based on mtDNA and we discuss implications of this findings on biogeography of eastern Kazakhstan populations. Marine Isotope Stage 11 is considered an important period for the diversification of the arvalis species group. In the context of our study, it is important to analyse genetically discontinuous Siberian populations, and the current distribution of Microtusmystacinus in new localities in eastern Kazakhstan.

Historical changes in grassland area determined the demography of semi-natural grassland butterflies in Japan

Semi-natural grassland areas expanded worldwide several thousand years ago following an increase in anthropogenic activities. However, semi-natural grassland habitat areas have been declining in recent decades due to changes in landuse, which have caused a loss of grassland biodiversity. Reconstructing historical and recent demographic changes in semi-natural grassland species will help clarify the factors affecting their population decline. Here we quantified past and recent demographic histories of Melitaea ambigua (Lepidoptera; Nymphalidae), an endangered grassland butterfly species in Japan. We examined changes in demography over the past 10,000 years based on 1378 bp of mitochondrial COI gene. We then examined changes in its genetic diversity and structure during the last 30 years using nine microsatellite DNA markers. The effective population size of M. ambigua increased about 3000-6000 years ago. In contrast, the genetic diversity and effective population sizes of many populations significantly declined from the 1980s to 2010s, which is consistent with a recent decline in the species population size. Our data suggest that the M. ambigua demography can be traced to changes in area covered by semi-natural grasslands throughout the Holocene.

Between the Balkans and the Baltic: Phylogeography of a Common Vole Mitochondrial DNA Lineage Limited to Central Europe

The common vole (Microtus arvalis) has been a model species of small mammal for studying end-glacial colonization history. In the present study we expanded the sampling from central and eastern Europe, analyzing contemporary genetic structure to identify the role of a potential 'northern glacial refugium', i.e. a refugium at a higher latitude than the traditional Mediterranean refugia. Altogether we analyzed 786 cytochrome b (cytb) sequences (representing mitochondrial DNA; mtDNA) from the whole of Europe, adding 177 new sequences from central and eastern Europe, and we conducted analyses on eight microsatellite loci for 499 individuals (representing nuclear DNA) from central and eastern Europe, adding data on 311 new specimens. Our new data fill gaps in the vicinity of the Carpathian Mountains, the potential northern refugium, such that there is now dense sampling from the Balkans to the Baltic Sea. Here we present evidence that the Eastern mtDNA lineage of the common vole was present in the vicinity of this Carpathian refugium during the Last Glacial Maximum and the Younger Dryas. The Eastern lineage expanded from this refugium to the Baltic and shows low cytb nucleotide diversity in those most northerly parts of the distribution. Analyses of microsatellites revealed a similar pattern but also showed little differentiation between all of the populations sampled in central and eastern Europe.

Retreat and extinction of the Late Pleistocene cave bear (Ursus spelaeus sensu lato)

The cave bear (Ursus spelaeus sensu lato) is a typical representative of Pleistocene megafauna which became extinct at the end of the Last Glacial. Detailed knowledge of cave bear extinction could explain this spectacular ecological transformation. The paper provides a report on the youngest remains of the cave bear dated to 20,930 ± 140 ¹⁴C years before present (BP). Ancient DNA analyses proved its affiliation to the Ursus ingressus haplotype. Using this record and 205 other dates, we determined, following eight approaches, the extinction time of this mammal at 26,100–24,300 cal. years BP. The time is only slightly earlier, i.e. 27,000–26,100 cal. years BP, when young dates without associated collagen data are excluded. The demise of cave bear falls within the coldest phase of the last glacial period, Greenland Stadial 3. This finding and the significant decrease in the cave bear records with cooling indicate that the drastic climatic changes were responsible for its extinction. Climate deterioration lowered vegetation productivity, on which the cave bear strongly depended as a strict herbivore. The distribution of the last cave bear records in Europe suggests that this animal was vanishing by fragmentation into subpopulations occupying small habitats. One of them was the Kraków-Częstochowa Upland in Poland, where we discovered the latest record of the cave bear and also two other, younger than 25,000 ¹⁴C years BP. The relatively long survival of this bear in karst regions may result from suitable microclimate and continuous access to water provided by deep aquifers, indicating a refugial role of such regions in the Pleistocene for many species.

Purifying Selection, Density Blocking and Unnoticed Mitochondrial DNA Diversity in the Red Deer, Cervus elaphus

The trajectories of postglacial range expansions, the occurrence of lineage patches and the formation and maintenance of secondary contact between lineages may mostly reflect neutral demographic processes, including density blocking, that may leave long-lasting genetic signatures. However, a few studies have recently shown that climate may also play a role. We used red deer, a large, mobile herbivore that is assumed to be sensitive to climate change, to test hypotheses of possible selection on the mitochondrial DNA cytochrome b gene (mtDNA cytb) and competitive and/or density-blocking (using mtDNA control region). We searched for a possible link between the phylogeographic structure and abiotic climatic variables. Finally, we tested for isolation by distance and isolation by environment and assessed the impact of human-mediated translocations on the genetic structure of red deer. Our analysis of 30 red deer populations in Poland using the mtDNA control region (N = 357) and cytochrome b (N = 50) markers not only confirmed the presence of the Western and South-Eastern lineages of the species but also indicated the presence of a previously unnoticed, rare relic haplotype that grouped together C. e. italicus from Italy (the Mesola deer). No significant signs of positive selection were detected for the mtDNA cytb gene in the studied red deer. However, a significant signal for purifying selection was found in our study that may explain the narrowness of the contact zone because gene flow between the Western and South-Eastern lineages should drive relatively strong mito-nuclear incompatibilities. MtDNA control region differentiation among red deer populations in Poland correlated with different abiotic climatic variables. Strikingly, the southernmost ice sheet limits during the Elsterian was the most important factor, and it explained the largest amount of variation. However, neither isolation by distance (IBD) nor isolation by environment (IBE) were recorded, and a very limited impact of human translocations was evident. The above-mentioned results suggest that in contemporary red deer populations in Poland, the phylogeographic pattern is well preserved, and long-term processes (density and/or competitive blocking) still play a major role.

Contrasting evolutionary histories of the legless lizards slow worms (Anguis) shaped by the topography of the Balkan Peninsula

Background

Genetic architecture of a species is a result of historical changes in population size and extent of distribution related to climatic and environmental factors and contemporary processes of dispersal and gene flow. Population-size and range contractions, expansions and shifts have a substantial effect on genetic diversity and intraspecific divergence, which is further shaped by gene-flow limiting barriers. The Balkans, as one of the most important sources of European biodiversity, is a region where many temperate species persisted during the Pleistocene glaciations and where high topographic heterogeneity offers suitable conditions for local adaptations of populations. In this study, we investigated the phylogeographical patterns and demographic histories of four species of semifossorial slow-worm lizards (genus Anguis) present in the Balkan Peninsula, and tested the relationship between genetic diversity and topographic heterogeneity of the inhabited ranges.

Results

We inferred phylogenetic relationships, compared genetic structure and historical demography of slow worms using nucleotide sequence variation of mitochondrial DNA. Four Anguis species with mostly parapatric distributions occur in the Balkan Peninsula. They show different levels of genetic diversity. A signature of population growth was detected in all four species but with various courses in particular populations. We found a strong correlation between genetic diversity of slow-worm populations and topographic ruggedness of the ranges (mountain systems) they inhabit. Areas with more rugged terrain harbour higher genetic diversity.

Conclusions

Phylogeographical pattern of the genus Anguis in the Balkans is concordant with the refugia-within-refugia model previously proposed for both several other taxa in the region and other main European Peninsulas. While slow-worm populations from the southern refugia mostly have restricted distributions and have not dispersed much from their refugial areas, populations from the extra-Mediterranean refugia in northern parts of the Balkans have colonized vast areas of eastern, central, and western Europe. Besides climatic historical events, the heterogeneous topography of the Balkans has also played an important role in shaping genetic diversity of slow worms.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0669-1) contains supplementary material, which is available to authorized users.