The Celtic fringe of Britain: insights from small mammal phylogeography

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.

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.

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.

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.

Unexpected post-glacial colonisation route explains the white colour of barn owls (Tyto alba) from the British Isles

The climate fluctuations of the Quaternary shaped the movement of species in and out of glacial refugia. In Europe, the majority of species followed one of the described traditional postglacial recolonization routes from the southern peninsulas towards the north. Like most organisms, barn owls are assumed to have colonized the British Isles by crossing over Doggerland, a land bridge that connected Britain to northern Europe. However, while they are dark rufous in northern Europe, barn owls in the British Isles are conspicuously white, a contrast that could suggest selective forces are at play on the islands. Yet, our analysis of known candidate genes involved in coloration found no signature of selection. Instead, using whole genome sequences and species distribution modelling, we found that owls colonised the British Isles soon after the last glaciation, directly from a white coloured refugium in the Iberian Peninsula, before colonising northern Europe. They would have followed a hitherto unknown post‐glacial colonization route to the Isles over a westwards path of suitable habitat in now submerged land in the Bay of Biscay, thus not crossing Doggerland. As such, they inherited the white colour of their Iberian founders and maintained it through low gene flow with the mainland that prevents the import of rufous alleles. Thus, we contend that neutral processes probably explain this contrasting white colour compared to continental owls. With the barn owl being a top predator, we expect future research will show this unanticipated route was used by other species from its paleo community.

Evidence for multiple colonisations and Wolbachia infections shaping the genetic structure of the widespread butterfly Polyommatus icarus in the British Isles

The paradigm of isolation in southern refugia during glacial periods followed by expansions during interglacials, producing limited genetic differentiation in northern areas, dominates European phylogeography. However, the existence of complex structured populations in formerly glaciated areas, and islands connected to mainland areas during glacial maxima, call for alternative explanations. We reconstructed the mtDNA phylogeography of the widespread Polyommatus Icarus butterfly with an emphasis on the formerly glaciated and connected British Isles. We found distinct geographical structuring of CO1 haplogroups, with an ancient lineage restricted to the marginal European areas, including Northern Scotland and Outer Hebrides. Population genomic analyses, using ddRADSeq genomic markers, also reveal substantial genetic structuring within Britain. However, there is negligble mito-nuclear concordance consistent with independent demographic histories of mitochondrial versus nuclear DNA. While mtDNA-Wolbachia associations in northern Britain could account for the geographic structuring of mtDNA across most of the British Isles, for nuclear DNA markers (derived from ddRADseq data) butterflies from France cluster between northern and southern British populations - an observation consistent with a scenario of multiple recolonisation. Taken together our results suggest that contemporary mtDNA structuring in the British Isles (and potentially elsewhere in Europe) largely results from Wolbachia infections, however, nuclear genomic structuring suggests a history of at least two distinct colonisations. This two-stage colonisation scenario has previously been put forth to explain genetic diversity and structuring in other British flora and fauna. Additionally, we also present preliminary evidence for potential Wolbachia-induced feminization in the Outer Hebrides.

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.

Adaptation, ancestral variation and gene flow in a 'Sky Island' Drosophila species

Over time, populations of species can expand, contract, fragment and become isolated, creating subpopulations that must adapt to local conditions. Understanding how species maintain variation after divergence as well as adapt to these changes in the face of gene flow is of great interest, especially as the current climate crisis has caused range shifts and frequent migrations for many species. Here, we characterize how a mycophageous fly species, Drosophila innubila, came to inhabit and adapt to its current range which includes mountain forests in south-western USA separated by large expanses of desert. Using population genomic data from more than 300 wild-caught individuals, we examine four populations to determine their population history in these mountain forests, looking for signatures of local adaptation. In this first extensive study, establishing D. innubila as a key genomic "Sky Island" model, we find D. innubila spread northwards during the previous glaciation period (30-100 KYA) and have recently expanded even further (0.2-2 KYA). D. innubila shows little evidence of population structure, consistent with a recent establishment and genetic variation maintained since before geographic stratification. We also find some signatures of recent selective sweeps in chorion proteins and population differentiation in antifungal immune genes suggesting differences in the environments to which flies are adapting. However, we find little support for long-term recurrent selection in these genes. In contrast, we find evidence of long-term recurrent positive selection in immune pathways such as the Toll signalling system and the Toll-regulated antimicrobial peptides.

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.

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.

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.

Playing Hide-and-Seek in Beta-Globin Genes: Gene Conversion Transferring a Beneficial Mutation between Differentially Expressed Gene Duplicates

Increasing evidence suggests that adaptation to diverse environments often involves selection on existing variation rather than new mutations. A previous study identified a nonsynonymous single nucleotide polymorphism (SNP) in exon 2 of two paralogous β-globin genes of the bank vole (Clethrionomys glareolus) in Britain in which the ancestral serine (Ser) and the derived cysteine (Cys) allele represent geographically partitioned functional variation affecting the erythrocyte antioxidative capacity. Here we studied the geographical pattern of the two-locus Ser/Cys polymorphism throughout Europe and tested for the geographic correlation between environmental variables and allele frequency, expected if the polymorphism was under spatially heterogeneous environment-related selection. Although bank vole population history clearly is important in shaping the dispersal of the oxidative stress protective Cys allele, analyses correcting for population structure suggest the Europe-wide pattern is affected by geographical variation in environmental conditions. The β-globin phenotype is encoded by the major paralog HBB-T1 but we found evidence of bidirectional gene conversion of exon 2 with the low-expression paralog HBB-T2. Our data support the model where gene conversion reshuffling genotypes between high- and low- expressed paralogs enables tuning of erythrocyte thiol levels, which may help maintain intracellular redox balance under fluctuating environmental conditions. Therefore, our study suggests a possible role for gene conversion between differentially expressed gene duplicates as a mechanism of physiological adaptation of populations to new or changing environments.

The colonization history of British water vole (Arvicola amphibius (Linnaeus, 1758)): origins and development of the Celtic fringe

The terminal Pleistocene and Early Holocene, a period from 15 000 to 18 000 Before Present (BP), was critical in establishing the current Holarctic fauna, with temperate-climate species largely replacing cold-adapted ones at mid-latitudes. However, the timing and nature of this process remain unclear for many taxa, a point that impacts on current and future management strategies. Here, we use an ancient DNA dataset to test more directly postglacial histories of the water vole (Arvicola amphibius, formerly A terrestris), a species that is both a conservation priority and a pest in different parts of its range. We specifically examine colonization of Britain, where a complex genetic structure can be observed today. Although we focus on population history at the limits of the species' range, the inclusion of additional European samples allows insights into European postglacial colonization events and provides a molecular perspective on water vole taxonomy.

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.

Frequent, geographically structured heteroplasmy in the mitochondria of a flowering plant, ribwort plantain (Plantago lanceolata)

Recent research has convincingly documented cases of mitochondrial heteroplasmy in a small set of wild and cultivated plant species. Heteroplasmy is suspected to be common in flowering plants and investigations of additional taxa may help understand the mechanisms generating heteroplasmy as well as its effects on plant phenotypes. The role of mitochondrial heteroplasmy is of particular interest in plants as cytoplasmic male sterility is controlled by mitochondrial genotypes, sometimes leading to co-occurring female and hermaphroditic individuals (gynodioecy). Paternal leakage may be important in the evolution of mating systems in such populations. We conducted a genetic survey of the gynodioecious plant Plantago lanceolata, in which heteroplasmy has not previously been reported, and estimated the frequencies of mitochondrial genotypes and heteroplasmy. Sanger sequence genotyping of 179 individuals from 15 European populations for two polymorphic mitochondrial loci, atp6 and rps12, identified 15 heteroplasmic individuals. These were distributed among 6 of the 10 populations that had polymorphisms in the target loci and represented 8% of all sampled individuals and 15% of the individuals in those 6 populations. The incidence was highest in Northern England and Scotland. Our results are consistent with geographic differences in the incidence of paternal leakage and/or the rates of nuclear restoration of male fertility.

New insights on postglacial colonization in western Europe: the phylogeography of the Leisler's bat (Nyctalus leisleri)

Despite recent advances in the understanding of the interplay between a dynamic physical environment and phylogeography in Europe, the origins of contemporary Irish biota remain uncertain. Current thinking is that Ireland was colonized post-glacially from southern European refugia, following the end of the last glacial maximum (LGM), some 20 000 years BP. The Leisler's bat (Nyctalus leisleri), one of the few native Irish mammal species, is widely distributed throughout Europe but, with the exception of Ireland, is generally rare and considered vulnerable. We investigate the origins and phylogeographic relationships of Irish populations in relation to those across Europe, including the closely related species N. azoreum. We use a combination of approaches, including mitochondrial and nuclear DNA markers, in addition to approximate Bayesian computation and palaeo-climatic species distribution modelling. Molecular analyses revealed two distinct and diverse European mitochondrial DNA lineages, which probably diverged in separate glacial refugia. A western lineage, restricted to Ireland, Britain and the Azores, comprises Irish and British N. leisleri and N. azoreum specimens; an eastern lineage is distributed throughout mainland Europe. Palaeo-climatic projections indicate suitable habitats during the LGM, including known glacial refugia, in addition to potential novel cryptic refugia along the western fringe of Europe. These results may be applicable to populations of many species.

Adaptive phylogeography: functional divergence between haemoglobins derived from different glacial refugia in the bank vole

Over the years, researchers have used presumptively neutral molecular variation to infer the origins of current species' distributions in northern latitudes (especially Europe). However, several reported examples of genic and chromosomal replacements suggest that end-glacial colonizations of particular northern areas may have involved genetic input from different source populations at different times, coupled with competition and selection. We investigate the functional consequences of differences between two bank vole (Clethrionomys glareolus) haemoglobins deriving from different glacial refugia, one of which partially replaced the other in Britain during end-glacial climate warming. This allows us to examine their adaptive divergence and hence a possible role of selection in the replacement. We determine the amino acid substitution Ser52Cys in the major expressed β-globin gene as the allelic difference. We use structural modelling to reveal that the protein environment renders the 52Cys thiol a highly reactive functional group and we show its reactivity in vitro. We demonstrate that possessing the reactive thiol in haemoglobin increases the resistance of bank vole erythrocytes to oxidative stress. Our study thus provides striking evidence for physiological differences between products of genic variants that spread at the expense of one another during colonization of an area from different glacial refugia.

Mitogenomic phylogenetics of the bank vole Clethrionomys glareolus, a model system for studying end-glacial colonization of Europe

We have revisited the mtDNA phylogeny of the bank vole Clethrionomys glareolus based on Sanger and next-generation Illumina sequencing of 32 complete mitochondrial genomes. The bank vole is a key study species for understanding the response of European fauna to the climate change following the Last Glacial Maximum (LGM) and one of the most convincing examples of a woodland mammal surviving in cryptic northern glacial refugia in Europe. The genomes sequenced included multiple representatives of each of the eight bank vole clades previously described based on cytochrome b (cob) sequences. All clades with the exception of the Basque - likely a misidentified pseudogene clade - were highly supported in all phylogenetic analyses and the relationships between the clades were resolved with high confidence. Our data extend the distribution of the Carpathian clade, the marker of a northern glacial refugium in the Carpathian Mountains, to include Britain and Fennoscandia (but not adjacent areas of continental Europe). The Carpathian sub-clade that colonized Britain and Fennoscandia had a somewhat different history from the sub-clade currently found in or close to the Carpathians and may have derived from a more north-westerly refugial area. The two bank vole populations that colonized Britain at the end of the last glaciation are for the first time linked with particular continental clades, the first colonists with the Carpathian clade and the second colonists with the western clade originating in a more southerly refugium in the vicinity of the Alps. We however found no evidence that a functional divergence of proteins encoded in the mitochondrial genome promoted the partial genetic replacement of the first colonists by the second colonists detected previously in southern Britain. We did identify one codon site that changed more often and more radically in the tree than expected and where the observed amino acid change may affect the reductase activity of the cytochrome bc1 complex, but the change was not specific to a particular clade. We also found an excess of radical changes to the primary protein structure for geographically restricted clades from southern Italy and Norway, respectively, possibly related to stronger selective pressure at the latitudinal extremes of the bank vole distribution. However, overall, we find little evidence of pervasive effects of deviation from neutrality on bank vole mtDNA phylogeography.

The Heterogeneity, Distribution, and Environmental Associations of Borrelia burgdorferi Sensu Lato, the Agent of Lyme Borreliosis, in Scotland

Lyme borreliosis is an emerging infectious human disease caused by the Borrelia burgdorferi sensu lato complex of bacteria with reported cases increasing in many areas of Europe and North America. To understand the drivers of disease risk and the distribution of symptoms, which may improve mitigation and diagnostics, here we characterize the genetics, distribution, and environmental associations of B. burgdorferi s.l. genospecies across Scotland. In Scotland, reported Lyme borreliosis cases have increased almost 10-fold since 2000 but the distribution of B. burgdorferi s.l. is so far unstudied. Using a large survey of over 2200 Ixodes ricinus tick samples collected from birds, mammals, and vegetation across 25 sites we identified four genospecies: Borrelia afzelii (48%), Borrelia garinii (36%), Borrelia valaisiana (8%), and B. burgdorferi sensu stricto (7%), and one mixed genospecies infection. Surprisingly, 90% of the sequence types were novel and, importantly, up to 14% of samples were mixed intra-genospecies co-infections, suggesting tick co-feeding, feeding on multiple hosts, or multiple infections in hosts. B. garinii (hosted by birds) was considerably more genetically diverse than B. afzelii (hosted by small mammals), as predicted since there are more species of birds than small mammals and birds can import strains from mainland Europe. Higher proportions of samples contained B. garinii and B. valaisiana in the west, while B. afzelii and B. garinii were significantly more associated with mixed/deciduous than with coniferous woodlands. This may relate to the abundance of transmission hosts in different regions and habitats. These data on the genetic heterogeneity within and between Borrelia genospecies are a first step to understand pathogen spread and could help explain the distribution of patient symptoms, which may aid local diagnosis. Understanding the environmental associations of the pathogens is critical for rational policy making for disease risk mitigation and land management.

Land-bridge calibration of molecular clocks and the post-glacial Colonization of Scandinavia by the Eurasian field vole Microtus agrestis

Phylogeography interprets molecular genetic variation in a spatial and temporal context. Molecular clocks are frequently used to calibrate phylogeographic analyses, however there is mounting evidence that molecular rates decay over the relevant timescales. It is therefore essential that an appropriate rate is determined, consistent with the temporal scale of the specific analysis. This can be achieved by using temporally spaced data such as ancient DNA or by relating the divergence of lineages directly to contemporaneous external events of known time. Here we calibrate a Eurasian field vole (Microtus agrestis) mitochondrial genealogy from the well-established series of post-glacial geophysical changes that led to the formation of the Baltic Sea and the separation of the Scandinavian peninsula from the central European mainland. The field vole exhibits the common phylogeographic pattern of Scandinavian colonization from both the north and the south, however the southernmost of the two relevant lineages appears to have originated in situ on the Scandinavian peninsula, or possibly in the adjacent island of Zealand, around the close of the Younger Dryas. The mitochondrial substitution rate and the timescale for the genealogy are closely consistent with those obtained with a previous calibration, based on the separation of the British Isles from mainland Europe. However the result here is arguably more certain, given the level of confidence that can be placed in one of the central assumptions of the calibration, that field voles could not survive the last glaciation of the southern part of the Scandinavian peninsula. Furthermore, the similarity between the molecular clock rate estimated here and those obtained by sampling heterochronous (ancient) DNA (including that of a congeneric species) suggest that there is little disparity between the measured genetic divergence and the population divergence that is implicit in our land-bridge calibration.

The complete mitochondrial genome of the bank vole Clethrionomys glareolus (Rodentia: Arvicolinae)

We present the first complete sequence of the bank vole (Clethrionomys glareolus) mitochondrial genome (GenBank accession no. KF918859). The bank vole mitogenome is 16,353 base pairs long and shows the gene content, genome architecture and gene strand asymmetry typical for mammals. The sequence provides an important new genomic resource for the bank vole, which is a popular study species in ecological and evolutionary research.

Multilocus sequence typing using mitochondrial genes (mtMLST) reveals geographic population structure of Ixodes ricinus ticks

The hard tick Ixodes ricinus is the principal vector of Lyme borreliosis (LB) group spirochaetes in Europe, but it also transmits a large number of other microbial pathogens that are of importance to animal and human health. Here, we characterise geographically distinct populations of this important ectoparasite based on multilocus sequence typing (MLST) of multiple mitochondrial (mt) genes (mtMLST). Internal fragments of approximately 500 bp were amplified and sequenced for 6 protein-encoding and ribosomal genes (atp6, coi, coii, coiii, cytB, and 12s). The samples analysed consisted of 506 questing nymphs collected in Britain and Latvia in 2006-2008 and in Latvia in 2002. Although little genetic structure has previously been observed in I. ricinus ticks among Europe, our data could clearly differentiate these 2 populations. Here, we argue that this novel scheme provides additional phylogenetic resolution which is important for understanding the genetic and geographic structure of I. ricinus populations. This in turn will benefit monitoring and management of tick-borne diseases.

Irish Cepaea nemoralis Land Snails Have a Cryptic Franco-Iberian Origin That Is Most Easily Explained by the Movements of Mesolithic Humans

The origins of flora and fauna that are only found in Ireland and Iberia, but which are absent from intervening countries, is one of the enduring questions of biogeography. As Southern French, Iberian and Irish populations of the land snail Cepaea nemoralis sometimes have a similar shell character, we used mitochondrial phylogenies to begin to understand if there is a shared “Lusitanian” history. Although much of Europe contains snails with A and D lineages, by far the majority of Irish individuals have a lineage, C, that in mainland Europe was only found in a restricted region of the Eastern Pyrenees. A past extinction of lineage C in the rest of Europe cannot be ruled out, but as there is a more than 8000 year continuous record of Cepaea fossils in Ireland, the species has long been a food source in the Pyrenees, and the Garonne river that flanks the Pyrenees is an ancient human route to the Atlantic, then we suggest that the unusual distribution of the C lineage is most easily explained by the movements of Mesolithic humans. If other Irish species have a similarly cryptic Lusitanian element, then this raises the possibility of a more widespread and significant pattern.

Population genetics, taxonomy, phylogeny and evolution of Borrelia burgdorferi sensu lato

In order to understand the population structure and dynamics of bacterial microorganisms, typing systems that accurately reflect the phylogenetic and evolutionary relationship of the agents are required. Over the past 15 years multilocus sequence typing schemes have replaced single locus approaches, giving novel insights into phylogenetic and evolutionary relationships of many bacterial species and facilitating taxonomy. Since 2004, several schemes using multiple loci have been developed to better understand the taxonomy, phylogeny and evolution of Lyme borreliosis spirochetes and in this paper we have reviewed and summarized the progress that has been made for this important group of vector-borne zoonotic bacteria.

Human-aided dispersal has altered but not erased the phylogeography of the tench

Human-aided dispersal can result in phylogeographic patterns that do not reflect natural historical processes, particularly in species prone to intentional translocations by humans. Here, we use a multiple-gene sequencing approach to assess the effects of human-aided dispersal on phylogeography of the tench Tinca tinca, a widespread Eurasian freshwater fish with a long history in aquaculture. Spatial genetic analysis applied to sequence data from four unlinked loci and 67 geographic localities (38-382 gene copies per locus) defined two groups of populations that were little structured geographically but were significantly differentiated from each other, and it identified locations of major genetic breaks, which were concordant across genes and were driven by distributions of two phylogroups. This pattern most reasonably reflects isolation in two major glacial refugia and subsequent range expansions, with the Eastern and Western phylogroups remaining largely allopatric throughout the tench range. However, this phylogeographic variation was also present in all 17 cultured breeds studied, and some populations at the western edge of the native range contained the Eastern phylogroup. Thus, natural processes have played an important role in structuring tench populations, but human-aided dispersal has also contributed significantly, with the admixed genetic composition of cultured breeds most likely contributing to the introgression.

Colonization of Ireland: revisiting 'the pygmy shrew syndrome' using mitochondrial, Y chromosomal and microsatellite markers

There is great uncertainty about how Ireland attained its current fauna and flora. Long-distance human-mediated colonization from southwestern Europe has been seen as a possible way that Ireland obtained many of its species; however, Britain has (surprisingly) been neglected as a source area for Ireland. The pygmy shrew has long been considered an illustrative model species, such that the uncertainty of the Irish colonization process has been dubbed 'the pygmy shrew syndrome'. Here, we used new genetic data consisting of 218 cytochrome (cyt) b sequences, 153 control region sequences, 17 Y-intron sequences and 335 microsatellite multilocus genotypes to distinguish between four possible hypotheses for the colonization of the British Isles, formulated in the context of previously published data. Cyt b sequences from western Europe were basal to those found in Ireland, but also to those found in the periphery of Britain and several offshore islands. Although the central cyt b haplotype in Ireland was found in northern Spain, we argue that it most likely occurred in Britain also, from where the pygmy shrew colonized Ireland as a human introduction during the Holocene. Y-intron and microsatellite data are consistent with this hypothesis, and the biological traits and distributional data of pygmy shrews argue against long-distance colonization from Spain. The compact starburst of the Irish cyt b expansion and the low genetic diversity across all markers strongly suggests a recent colonization. This detailed molecular study of the pygmy shrew provides a new perspective on an old colonization question.

Post-glacial partitioning of mitochondrial genetic variation in the field vole

Genetic markers are often used to examine population history. There is considerable debate about the behaviour of molecular clock rates around the population-species transition. Nevertheless, appropriate calibration is critical to any inference regarding the absolute timing and scale of demographic changes. Here, we use a mitochondrial cytochrome b gene genealogy, based entirely on modern sequences and calibrated from recent geophysical events, to date the post-glacial expansion of the Eurasian field vole (Microtus agrestis), a widespread temperate mammal species. The phylogeographic structure reflects the subsequent expansion of populations that went through bottlenecks at the time of the Younger Dryas (ca 12,000 years BP) rather than the Last Glacial Maximum (LGM, ca 24,000 years BP), which is usually seen as the time when present-day patterns were determined. The nucleotide substitution rate that was estimated here, ca 4 × 10(-7) substitutions/site/year, remains extremely high throughout the relevant time frame. Calibration with similarly high population-based substitution rates, rather than long-term rates derived from species divergence times, will show that post-LGM climatic events generated current phylogeographic structure in many other organisms from temperate latitudes.

House mouse colonization patterns on the sub-Antarctic Kerguelen Archipelago suggest singular primary invasions and resilience against re-invasion

BACKGROUND

Starting from Western Europe, the house mouse (Mus musculus domesticus) has spread across the globe in historic times. However, most oceanic islands were colonized by mice only within the past 300 years. This makes them an excellent model for studying the evolutionary processes during early stages of new colonization. We have focused here on the Kerguelen Archipelago, located within the sub-Antarctic area and compare the patterns with samples from other Southern Ocean islands.

RESULTS

We have typed 18 autosomal and six Y-chromosomal microsatellite loci and obtained mitochondrial D-loop sequences for a total of 534 samples, mainly from the Kerguelen Archipelago, but also from the Falkland Islands, Marion Island, Amsterdam Island, Antipodes Island, Macquarie Island, Auckland Islands and one sample from South Georgia. We find that most of the mice on the Kerguelen Archipelago have the same mitochondrial haplotype and all share the same major Y-chromosomal haplotype. Two small islands (Cochons Island and Cimetière Island) within the archipelago show a different mitochondrial haplotype, are genetically distinct for autosomal loci, but share the major Y-chromosomal haplotype. In the mitochondrial D-loop sequences, we find several single step mutational derivatives of one of the major mitochondrial haplotypes, suggesting an unusually high mutation rate, or the occurrence of selective sweeps in mitochondria.

CONCLUSIONS

Although there was heavy ship traffic for over a hundred years to the Kerguelen Archipelago, it appears that the mice that have arrived first have colonized the main island (Grande Terre) and most of the associated small islands. The second invasion that we see in our data has occurred on islands that are detached from Grande Terre and were likely to have had no resident mice prior to their arrival. The genetic data suggest that the mice of both primary invasions originated from related source populations. Our data suggest that an area colonized by mice is refractory to further introgression, possibly due to fast adaptations of the resident mice to local conditions.