Genomics of end-Pleistocene population replacement in a small mammal

Genetic admixture drives climate adaptation in the bank vole

Genetic admixture introduces new variants at relatively high frequencies, potentially aiding rapid responses to environmental changes. Here, we evaluate its role in adaptive variation related to climatic conditions in bank voles (Clethrionomys glareolus) in Britain, using whole-genome data. Our results reveal loci showing excess ancestry from one of the two postglacial colonist populations inconsistent with overall admixture patterns. Notably, loci associated with climate adaptation exhibit disproportionate amounts of excess ancestry, highlighting the impact of admixture between colonist populations on local adaptation. The results suggest strong and localized selection on climate-adaptive loci, as indicated by steep clines and/or shifted cline centres, during population replacement. A subset, including a haemoglobin gene, is associated with oxidative stress responses, underscoring a role of oxidative stress in local adaptation. Our study highlights the important contribution of admixture during secondary contact between populations from distinct climatic refugia enriching adaptive diversity. Understanding these dynamics is crucial for predicting future adaptive capacity to anthropogenic climate change.

Local adaptation and future climate vulnerability in a wild rodent

As climate change continues, species pushed outside their physiological tolerance limits must adapt or face extinction. When change is rapid, adaptation will largely harness ancestral variation, making the availability and characteristics of that variation of critical importance. Here, we used whole-genome sequencing and genetic-environment association analyses to identify adaptive variation and its significance in the context of future climates in a small Palearctic mammal, the bank vole (Clethrionomys glareolus). We found that peripheral populations of bank vole in Britain are already at the extreme bounds of potential genetic adaptation and may require an influx of adaptive variation in order to respond. Analyses of adaptive loci suggest regional differences in climate variables select for variants that influence patterns of population adaptive resilience, including genes associated with antioxidant defense, and support a pattern of thermal/hypoxic cross-adaptation. Our findings indicate that understanding potential shifts in genomic composition in response to climate change may be key to predicting species' fate under future climates.

Evolutionary history of the extinct wolf population from France in the context of global phylogeographic changes throughout the Holocene

Phylogeographic patterns in large mammals result from natural environmental factors and anthropogenic effects, which in some cases include domestication. The grey wolf was once widely distributed across the Holarctic, but experienced phylogeographic shifts and demographic declines during the Holocene. In the 19th-20th centuries, the species became extirpated from large parts of Europe due to direct extermination and habitat loss. We reconstructed the evolutionary history of the extinct Western European wolves based on the mitogenomic composition of 78 samples from France (Neolithic-20th century) in the context of other populations of wolves and dogs worldwide. We found a close genetic similarity of French wolves from ancient, medieval and recent populations, which suggests the long-term continuity of maternal lineages. MtDNA haplotypes of the French wolves showed large diversity and fell into two main haplogroups of modern Holarctic wolves. Our worldwide phylogeographic analysis indicated that haplogroup W1, which includes wolves from Eurasia and North America, originated in Northern Siberia. Haplogroup W2, which includes only European wolves, originated in Europe ~35 kya and its frequency was reduced during the Holocene due to an expansion of haplogroup W1 from the east. Moreover, we found that dog haplogroup D, currently restricted to Europe and the Middle East, was nested within the wolf haplogroup W2. This suggests European origin of haplogroup D, probably as a result of an ancient introgression from European wolves. Our results highlight the dynamic evolutionary history of European wolves during the Holocene, with a partial lineage replacement and introgressive hybridization with local dog populations.

Different waves of postglacial recolonisation and genomic structure of bank vole populations in NE Poland

Previous studies indicated that in some species phylogeographic patterns obtained in the analysis of nuclear and mitochondrial DNA (mtDNA) markers can be different. Such mitonuclear discordance can have important evolutionary and ecological consequences. In the present study, we aimed to check whether there was any discordance between mtDNA and nuclear DNA in the bank vole population in the contact zone of its two mtDNA lineages. We analysed the population genetic structure of bank voles using genome-wide genetic data (SNPs) and diversity of sequenced heart transcriptomes obtained from selected individuals from three populations inhabiting areas outside the contact zone. The SNP genetic structure of the populations confirmed the presence of at least two genetic clusters, and such division was concordant with the patterns obtained in the analysis of other genetic markers and functional genes. However, genome-wide SNP analyses revealed the more detailed structure of the studied population, consistent with more than two bank vole recolonisation waves, as recognised previously in the study area. We did not find any significant differences between individuals representing two separate mtDNA lineages of the species in functional genes coding for protein-forming complexes, which are involved in the process of cell respiration in mitochondria. We concluded that the contemporary genetic structure of the populations and the width of the contact zone were shaped by climatic and environmental factors rather than by genetic barriers. The studied populations were likely isolated in separate Last Glacial Maximum refugia for insufficient amount of time to develop significant genetic differentiation.

Diverse susceptibilities and responses of human and rodent cells to orthohantavirus infection reveal different levels of cellular restriction

Orthohantaviruses are rodent-borne emerging viruses that may cause severe diseases in humans but no apparent pathology in their small mammal reservoirs. However, the mechanisms leading to tolerance or pathogenicity in humans and persistence in rodent reservoirs are poorly understood, as is the manner in which they spread within and between organisms. Here, we used a range of cellular and molecular approaches to investigate the interactions of three different orthohantaviruses-Puumala virus (PUUV), responsible for a mild to moderate form of hemorrhagic fever with renal syndrome in humans, Tula virus (TULV) with low pathogenicity, and non-pathogenic Prospect Hill virus (PHV)-with human and rodent host cell lines. Besides the fact that cell susceptibility to virus infection was shown to depend on the cell type and virus strain, the three orthohantaviruses were able to infect Vero E6 and HuH7 human cells, but only the former secreted infectious particles. In cells derived from PUUV reservoir, the bank vole (Myodes glareolus), PUUV achieved a complete viral cycle, while TULV did not enter the cells and PHV infected them but did not produce infectious particles, reflecting differences in host specificity. A search for mature virions by electron microscopy (EM) revealed that TULV assembly occurred in part at the plasma membrane, whereas PHV particles were trapped in autophagic vacuoles in cells of the heterologous rodent host. We described differential interactions of orthohantaviruses with cellular factors, as supported by the cellular distribution of viral nucleocapsid protein with cell compartments, and proteomics identification of cellular partners. Our results also showed that interferon (IFN) dependent gene expression was regulated in a cell and virus species dependent manner. Overall, our study highlighted the complexity of the host-virus relationship and demonstrated that orthohantaviruses are restricted at different levels of the viral cycle. In addition, the study opens new avenues to further investigate how these viruses differ in their interactions with cells to evade innate immunity and how it depends on tissue type and host species.

Genic distribution modelling predicts adaptation of the bank vole to climate change

The most likely pathway for many species to survive future climate change is by pre-existing trait variation providing a fitness advantage under the new climate. Here we evaluate the potential role of haemoglobin (Hb) variation in bank voles under future climate change. We model gene-climate relationships for two functionally distinct Hb types, HbS and HbF, which have a north-south distribution in Britain presenting an unusually tractable system linking genetic variation in physiology to geographical and temporal variation in climate. Projections to future climatic conditions suggest a change in relative climatic suitability that would result in HbS being displaced by HbF in northern Britain. This would facilitate local adaptation to future climate—without Hb displacement, populations in northern Britain would likely be suboptimally adapted because their Hb would not match local climatic conditions. Our study shows how pre-existing physiological differences can influence the adaptive capacity of species to climate change.

Haemoglobin variation in British bank voles combined with climate models predict future regional allelic replacement reflecting capacity for adaptation to climate change.

Deep structure, long-distance migration and admixture in the colour polymorphic land snail Cepaea nemoralis

Although snails of the genus Cepaea have historically been important in studying colour polymorphism, an ongoing issue is that there is a lack of knowledge of the underlying genetics of the polymorphism, as well as an absence of genomic data to put findings in context. We, therefore, used phylogenomic methods to begin to investigate the post‐glacial history of Cepaea nemoralis, with a long‐term aim to understand the roles that selection and drift have in determining both European‐wide and local patterns of colour polymorphism. By combining prior and new mitochondrial DNA data from over 1500 individuals with ddRAD genomic data from representative individuals across Europe, we show that patterns of differentiation are primarily due to multiple deeply diverged populations of snails. Minimally, there is a widespread Central European population and additional diverged groups in Northern Spain, the Pyrenees, as well as likely Italy and South Eastern Europe. The genomic analysis showed that the present‐day snails in Ireland and possibly some other locations are likely descendants of admixture between snails from the Pyrenees and the Central European group, an observation that is consistent with prior inferences from mitochondrial DNA alone. The interpretation is that C. nemoralis may have arrived in Ireland via long‐distance migration from the Pyrenean region, subsequently admixing with arrivals from elsewhere. This work, therefore, provides a baseline expectation for future studies on the genetics of the colour polymorphism, as well as providing a comparator for similar species.

The Bank Vole (Clethrionomys glareolus) as a Model System for Adaptive Phylogeography in the European Theater

The legacy of climatic changes during the Pleistocene glaciations allows inferences to be made about the patterns and processes associated with range expansion/colonization, including evolutionary adaptation. With the increasing availability of population genomic data, we have the opportunity to examine these questions in detail and in a variety of non-traditional model species. As an exemplar, here we review more than two decades of work by our group and others that illustrate the potential of a single “non-model model” mammal species - the bank vole (Clethrionomys glareolus), which is particularly well suited to illustrate the complexities that may be associated with range expansion and the power of genomics (and other datasets) to uncover them. We first summarize early phylogeographic work using mitochondrial DNA and then describe new phylogeographic insights gained from population genomic analysis of genome-wide SNP data to highlight the bank vole as one of the most compelling examples of a forest mammal, that survived in cryptic extra-Mediterranean (“northern”) glacial refugia in Europe, and as one of the species in which substantial replacement and mixing of lineages originating from different refugia occurred during end-glacial colonization. Our studies of bank vole hemoglobin structure and function, as well as our recent ecological niche modeling study examining differences among bank vole lineages, led us to develop the idea of “adaptive phylogeography.” This is what we call the study of the role of adaptive differences among populations in shaping phylogeographic patterns. Adaptive phylogeography provides a link between past population history and adaptation that can ultimately help predict the potential of future species responses to climate change. Because the bank vole is part of a community of organisms whose range has repeatedly contracted and then expanded in the past, what we learn from the bank vole will be useful for our understanding of a broad range of species.

Genomic insights into the historical and contemporary demographics of the grey reef shark

Analyses of genetic diversity can shed light on both the origins of biodiversity hotspots, as well as the conservation status of species that are impacted by human activities. With these objectives, we assembled a genomic dataset of 14,935 single nucleotide polymorphisms from 513 grey reef sharks (Carcharhinus amblyrhynchos) sampled across 17 locations in the tropical Indo-Pacific. We analysed geographic variation in genetic diversity, estimated ancient and contemporary effective population size (N_e) across sampling locations (using coalescent and linkage disequilibrium methods) and modelled the history of gene flow between the Coral Triangle and the Coral Sea. Genetic diversity decreased with distance away from the Coral Triangle and north-western Australia, implying that C. amblyrhynchos may have originated in this region. Increases in N_e were detected across almost all sampling locations 40,000-90,000 generations ago (approximately 0.6-1.5 mya, given an estimated generation time of 16.4 years), suggesting a range expansion around this time. More recent, secondary increases in N_e were inferred for the Misool and North Great Barrier Reef sampling locations, but joint modelling did not clarify whether these were due to population growth, migration, or both. Despite the greater genetic diversity and ancient N_e observed at sites around Australia and the Coral Triangle, remote reefs around north-western New Caledonia had the highest contemporary N_e, demonstrating the importance of using multiple population size assessment methods. This study provides insight into both the past and present demographics of C. amblyrhynchos and contributes to our understanding of evolution in marine biodiversity hotspots.

Next-generation phylogeography resolves post-glacial colonization patterns in a widespread carnivore, the red fox (Vulpes vulpes), in Europe

Carnivores tend to exhibit a lack of (or less pronounced) genetic structure at continental scales in both a geographic and temporal sense and this can confound the identification of post-glacial colonization patterns in this group. In this study we used genome-wide data (using genotyping by sequencing [GBS]) to reconstruct the phylogeographic history of a widespread carnivore, the red fox (Vulpes vulpes), by investigating broad-scale patterns of genomic variation, differentiation and admixture amongst contemporary populations in Europe. Using 15,003 single nucleotide polymorphisms (SNPs) from 524 individuals allowed us to identify the importance of refugial regions for the red fox in terms of endemism (e.g., Iberia). In addition, we tested multiple post-glacial recolonization scenarios of previously glaciated regions during the Last Glacial Maximum using an Approximate Bayesian Computation (ABC) approach that were unresolved from previous studies. This allowed us to identify the role of admixture from multiple source population post-Younger Dryas in the case of Scandinavia and ancient land-bridges in the colonization of the British Isles. A natural colonization of Ireland was deemed more likely than an ancient human-mediated introduction as has previously been proposed and potentially points to a larger mammalian community on the island in the early post-glacial period. Using genome-wide data has allowed us to tease apart broad-scale patterns of structure and diversity in a widespread carnivore in Europe that was not evident from using more limited marker sets and provides a foundation for next-generation phylogeographic studies in other non-model species.

Niche differentiation in a postglacial colonizer, the bank vole Clethrionomys glareolus

Species-level environmental niche modeling has been crucial in efforts to understand how species respond to climate variation and change. However, species often exhibit local adaptation and intraspecific niche differences that may be important to consider in predicting responses to climate. Here, we explore whether phylogeographic lineages of the bank vole originating from different glacial refugia (Carpathian, Western, Eastern, and Southern) show niche differentiation, which would suggest a role for local adaptation in biogeography of this widespread Eurasian small mammal. We first model the environmental requirements for the bank vole using species-wide occurrences (210 filtered records) and then model each lineage separately to examine niche overlap and test for niche differentiation in geographic and environmental space. We then use the models to estimate past [Last Glacial Maximum (LGM) and mid-Holocene] habitat suitability to compare with previously hypothesized glacial refugia for this species. Environmental niches are statistically significantly different from each other for all pairs of lineages in geographic and environmental space, and these differences cannot be explained by habitat availability within their respective ranges. Together with the inability of most of the lineages to correctly predict the distributions of other lineages, these results support intraspecific ecological differentiation in the bank vole. Model projections of habitat suitability during the LGM support glacial survival of the bank vole in the Mediterranean region and in central and western Europe. Niche differences between lineages and the resulting spatial segregation of habitat suitability suggest ecological differentiation has played a role in determining the present phylogeographic patterns in the bank vole. Our study illustrates that models pooling lineages within a species may obscure the potential for different responses to climate change among populations.

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.

An amphibian species pushed out of Britain by a moving hybrid zone

Classical theory states that hybrid zones will be stable in troughs of low population density where dispersal is hampered. Yet, evidence for moving hybrid zones is mounting. One possible reason that moving zones have been underappreciated is that they may drive themselves into oblivion and with just the superseding species remaining, morphological and genetic signals of past species replacement may be difficult to appreciate. Using genetic data (32 diagnostic single nucleotide polymorphisms) from a clinal hybrid zone of the common toad (Bufo bufo) and the spined toad (Bufo spinosus) in France for comparison, alleles of the latter species were documented in common toads in the south of Great Britain, at frequencies in excess of 10%. Because long distance dispersal across the Channel is unlikely, the conclusion reached was that the continental toad hybrid zone which previously extended into Britain, moved southwards and extirpated B. spinosus. Species distribution models for the mid-Holocene and the present support that climate has locally changed in favour of B. bufo. The system bears resemblance with the demise of Homo neanderthalensis and the rise of Homo sapiens and provides an example that some paleoanthropologists demanded in support of a hominin "leaky replacement" scenario. The toad example is informative just because surviving pure B. spinosus and an extant slowly moving interspecific hybrid zone are available for comparison.

Field vole-associated Traemmersee hantavirus from Germany represents a novel hantavirus species

Vole-associated hantaviruses occur in the Old and New World. Tula orthohantavirus (TULV) is widely distributed throughout the European continent in its reservoir, the common vole (Microtus arvalis), but the virus was also frequently detected in field voles (Microtus agrestis) and other vole species. TULV and common voles are absent from Great Britain. However, field voles there harbor Tatenale and Kielder hantaviruses. Here we screened 126 field voles and 13 common voles from Brandenburg, Germany, for hantavirus infections. One common vole and four field voles were anti-TULV antibody and/or TULV RNA positive. In one additional, seropositive field vole a novel hantavirus sequence was detected. The partial S and L segment nucleotide sequences were only 61.1% and 75.6% identical to sympatrically occurring TULV sequences, but showed highest similarity of approximately 80% to British Tatenale and Kielder hantaviruses. Subsequent determination of the entire nucleocapsid (N), glycoprotein (GPC), and RNA-dependent RNA polymerase encoding sequences and determination of the pairwise evolutionary distance (PED) value for the concatenated N and GPC amino acid sequences confirmed a novel orthohantavirus species, tentatively named Traemmersee orthohantavirus. The identification of this novel hantavirus in a field vole from eastern Germany underlines the necessity of a large-scale, broad geographical hantavirus screening of voles to understand evolutionary processes of virus-host associations and host switches.

Genomics of end-Pleistocene population replacement in a small mammal

Current species distributions at high latitudes are the product of expansion from glacial refugia into previously uninhabitable areas at the end of the last glaciation. The traditional view of postglacial colonization is that southern populations expanded their ranges into unoccupied northern territories. Recent findings on mitochondrial DNA (mtDNA) of British small mammals have challenged this simple colonization scenario by demonstrating a more complex genetic turnover in Britain during the Pleistocene-Holocene transition where one mtDNA clade of each species was replaced by another mtDNA clade of the same species. Here, we provide evidence from one of those small mammals, the bank vole (Clethrionomys glareolus), that the replacement was genome-wide. Using more than 10 000 autosomal SNPs we found that similar to mtDNA, bank vole genomes in Britain form two (north and south) clusters which admix. Therefore, the genome of the original postglacial colonists (the northern cluster) was probably replaced by another wave of migration from a different continental European population (the southern cluster), and we gained support for this by modelling with approximate Bayesian computation. This finding emphasizes the importance of analysis of genome-wide diversity within species under changing climate in creating opportunities for sophisticated testing of population history scenarios.