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

Genomic signatures of climate adaptation in bank voles

Evidence for divergent selection and adaptive variation across the landscape can provide insight into a species' ability to adapt to different environments. However, despite recent advances in genomics, it remains difficult to detect the footprints of climate-mediated selection in natural populations. Here, we analysed ddRAD sequencing data (21,892 SNPs) in conjunction with geographic climate variation to search for signatures of adaptive differentiation in twelve populations of the bank vole (Clethrionomys glareolus) distributed across Europe. To identify the loci subject to selection associated with climate variation, we applied multiple genotype-environment association methods, two univariate and one multivariate, and controlled for the effect of population structure. In total, we identified 213 candidate loci for adaptation, 74 of which were located within genes. In particular, we identified signatures of selection in candidate genes with functions related to lipid metabolism and the immune system. Using the results of redundancy analysis, we demonstrated that population history and climate have joint effects on the genetic variation in the pan-European metapopulation. Furthermore, by examining only candidate loci, we found that annual mean temperature is an important factor shaping adaptive genetic variation in the bank vole. By combining landscape genomic approaches, our study sheds light on genome-wide adaptive differentiation and the spatial distribution of variants underlying adaptive variation influenced by local climate in bank voles.

Shrew Communities in Mediterranean Agro-Ecosystems of Central Greece: Associations with Crop Types, Land Uses, and Soil Parameters

Shrew communities play a crucial role in a diverse range of natural, urban, and agricultural ecosystems. We used Barn owl diet analysis as the ideal proxy to assess small-mammal distribution patterns on large spatial scales. More than 10,000 pellets were analyzed from Thessaly, the largest agricultural prefecture located in central Greece. A total of more than 29,000 prey items were identified, one of the largest datasets used in similar analyses in Europe. Three discrete shrew species were present in Thessaly agricultural plains, central Greece (Güldenstädt's shrew Crocidura gueldenstaedtii, Bi-coloured shrew Crocidura leucodon, and Pygmy white-toothed shrew Suncus etruscus), which comprised a total of 7452 shrews, representing 25.64% of the total small-mammals' dataset. C. gueldenstaedtii and S. etruscus demonstrated strong associations with heavy argillaceous-clay soils and Vertisol soil types, whereas S. etruscus was also associated with non-irrigated land and non-intensive cultivated plots. C. leucodon demonstrated no significant associations to any environmental gradient and demonstrated habitat plasticity, most possibly shaped by existing resources and competition. Our study highlights the important insights gained from Barn owl diet analysis in respect of small-mammal assemblages on broad geographical scales, and the inclusion of soil parameters as drivers of habitat suitability and distribution patterns for small-mammal responses.

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.

Population genetic analysis of the microsporidium Ordospora colligata reveals the role of natural selection and phylogeography on its extremely compact and reduced genome

The determinants of variation in a species' genome-wide nucleotide diversity include historical, environmental, and stochastic aspects. This diversity can inform us about the species' past and present evolutionary dynamics. In parasites, the mode of transmission and the interactions with the host might supersede the effects of these aspects in shaping parasite genomic diversity. We used genomic samples from 10 populations of the microsporidian parasite Ordospora colligata to investigate present genomic diversity and how it was shaped by evolutionary processes, specifically, the role of phylogeography, co-phylogeography (with the host), natural selection, and transmission mode. Although very closely related microsporidia cause diseases in humans, O. colligata is specific to the freshwater crustacean Daphnia magna and has one of the smallest known eukaryotic genomes. We found an overlapping phylogeography between O. colligata and its host highlighting the long-term, intimate relationship between them. The observed geographic distribution reflects previous findings that O. colligata exhibits adaptations to colder habitats, which differentiates it from other microsporidian gut parasites of D. magna predominantly found in warmer areas. The co-phylogeography allowed us to calibrate the O. colligata phylogeny and thus estimate its mutation rate. We identified several genetic regions under potential selection. Our whole-genome study provides insights into the evolution of one of the most reduced eukaryotic genomes and shows how different processes shape genomic diversity of an obligate parasite.

Distribution of Four Vole Species through the Barn Owl Tyto alba Diet Spectrum: Pattern Responses to Environmental Gradients in Intensive Agroecosystems of Central Greece

Voles are the most common vertebrate pests in European agriculture. Identifying their distribution and abundance patterns provides valuable information for future management. Barn Owl diet analysis is one of the optimum methods used to record small mammal distribution patterns on large spatial scales. From 2003 to 2005, a total of 10,065 Barn Owl pellets were collected and analyzed from 31 breeding sites in the largest agroecosystem in Greece, the Thessaly plains. A total of 29,061 prey items were identified, offering deep insight into small mammal distribution, specifically voles. Four discrete vole species (Harting's vole Microtus hartingi, East European vole Microtus levis, Thomas's pine vole Microtus thomasi, and Grey dwarf hamster Cricetulus migratorius) comprised 40.5% (11,770 vole prey items) of the total Barn Owl prey intake. The presence and abundance of the voles varied according to underlying environmental gradients, with soil texture and type playing a major role. M. levis showed no significant attachments to gradients, other than a mild increase in Mollisol soils. It was syntopic in all sites with M. hartingi, which was the dominant and most abundant small mammal species, preferring non-arable cultivated land, natural grasslands, set-aside fields, and fallow land. M. thomasi was strictly present in western Thessaly and strongly associated with a sandy-clay soil texture and Alfisol soils. C. migratorius was the least represented vole (162 items), exclusively present in eastern Thessaly and demonstrating a stronger association with cereals, Mollisol soils, and an argillaceous-clay soil texture. This is the first study in Greece at such a large spatial scale, offering insights for pest rodents' distribution in intensive agroecosystems and their response to environmental gradients including soil parameters.

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.

Complete mitogenomes reveal limited genetic variability in the garden dormouse Eliomys quercinus of the Iberian Peninsula

The garden dormouse Eliomys quercinus is a poorly known Western Palearctic species experiencing a global decline. Even though the availability of genetic information is key to assess the driversunderlying demographic changes in wild populations and plan adequate management, data on E. quercinus are still scant. In this study, we reconstructed the complete mitogenomes of four E. quercinus individuals from southern Spain using in–solution enriched libraries, and found evidence of limited genetic variability. We then compared their cytochrome b sequences to those of conspecifics from other countries and supported the divergent but genetically depauperate position of this evolutionarily significant unit (ESU). The information produced will assist future conservation studies on this little–studied rodent.

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.