DNA footprints of European hedgehogs, Erinaceus europaeus and E. concolor: Pleistocene refugia, postglacial expansion and colonization routes

Hematology, Biochemistry, and Protein Electrophoresis Reference Intervals of Western European Hedgehog (Erinaceus europaeus) from a Rehabilitation Center in Northern Portugal

The Western European hedgehog (Erinaceus europaeus) can work as a bioindicator of environmental pollution and be a host for multiple zoonotic agents, making it relevant in terms of One Health studies. It is essential to deepen the knowledge on this species and calculate reference intervals (RIs) for the usual hematological and biochemical parameters. For this retrospective study (2017–2022), the archives of the Clinical Pathology Laboratory (LPC) of University of Trás-os-Montes and Alto Douro (UTAD) Veterinary Teaching Hospital were analyzed. Data of hematology, clinical biochemistry, and protein electrophoresis from 37 healthy hedgehogs of the Wild Animal Rehabilitation Center at UTAD, Northern Portugal, were included. It was possible to calculate RIs for almost all of the variables in the study, using Reference Value Advisor V2.1. Moreover, sex and age effects were investigated: alkaline phosphatase (p = 0.012, higher in males); total proteins (p = 0.034, higher in adults); mean cell volume (p = 0.007) and mean corpuscular hemoglobin (p = 0.010) (both higher in juveniles); and red blood cell distribution width (p = 0.021, higher in adults). Our study allowed for the first time to define RIs for a population of hedgehogs in Portugal, having a potentially relevant impact on species conservation and in the human–animal health interface.

Contrasting levels of hybridization across the two contact zones between two hedgehog species revealed by genome-wide SNP data

Hybridization and introgression have played important roles in the history of various species, including lineage diversification and the evolution of adaptive traits. Hybridization can accelerate the development of reproductive isolation between diverging species, and thus valuable insight into the evolution of reproductive barrier formation may be gained by studying secondary contact zones. Hedgehogs of the genus Erinaceus, which are insectivores sensitive to changes in climate, are a pioneer model in Pleistocene phylogeography. The present study provides the first genome-wide SNP data regarding the Erinaceus hedgehogs species complex, offering a unique comparison of two secondary contact zones between Erinaceus europaeus and E. roumanicus. Results confirmed diversification of the genus during the Pleistocene period, and detected a new refugial lineage of E. roumanicus outside the Mediterranean basin, most likely in the Ponto-Caspian region. In the Central European zone, the level of hybridization was low, whereas in the Russian-Baltic zone, both species hybridise extensively. Asymmetrical gene flow from E. europaeus to E. roumanicus suggests that reproductive isolation varies according to the direction of the crosses in the hybrid zones. However, no loci with significantly different patterns of introgression were detected. Markedly different pre- and post-zygotic barriers, and thus diverse modes of species boundary maintenance in the two contact zones, likely exist. This pattern is probably a consequence of the different age and thus of the different stage of evolution of reproductive isolating mechanisms in each hybrid zone.

The European Pine Marten Martes martes (Linnaeus, 1758) Is Autochthonous in Sicily and Constitutes a Well-Characterised Major Phylogroup within the Species (Carnivora, Mustelidae)

Simple Summary

The faunal assemblages currently occurring on islands are often a mélange of native and non-native species, and the actual status of some of them is doubtful at present. Since different laws are enforced for native and non-native species, sound knowledge about their status is pivotal for both their management and our understanding of the natural history of the studied areas. In the frame of this work, the Sicilian population of the European pine marten is genetically characterized for the first time, based on mitochondrial DNA sequences. Our results prove that the European pine marten is native on the island, where it arrived during the Pleistocene, and is represented there by a well-differentiated and endemic evolutionary lineage.

Abstract

No molecular data are currently available for the Sicilian populations of the European pine marten Martes martes, thus preventing any sound inference about its native or non-native status on the island, as well as the local phylogeography of the species. In order to investigate these issues, we sequenced two mtDNA markers in road-killed specimens collected in Sicily. Both markers consistently demonstrated the existence of a well-characterised Sicilian clade of the species, which is endemic to the island and constitutes the sister group of a clade including the Mediterranean and Central–North European major phylogroups of the European pine marten. Such evidence supports the autochthony of Martes martes in Sicily and points to a natural Pleistocene colonisation of the island followed by isolation. The occurrence of a, to date undetected, major phylogroup of the species in Sicily calls for the dedicated monitoring of the Sicilian populations of the species in order to preserve this evolutionarily significant unit.

Genetic evidence for post-glacial expansion from a southern refugium in the eastern moa (Emeus crassus)

Cycles of glacial expansion and contraction throughout the Pleistocene drove increases and decreases, respectively, in the geographical range and population size of many animal species. Genetic data have revealed that during glacial maxima the distribution of many Eurasian animals was restricted to small refugial areas, from which species expanded to reoccupy parts of their former range as the climate warmed. It has been suggested that the extinct eastern moa (Emeus crassus)-a large, flightless bird from New Zealand-behaved analogously during glacial maxima, possibly surviving only in a restricted area of lowland habitat in the southern South Island of New Zealand during the Last Glacial Maximum (LGM). However, previous studies have lacked the power and geographical sampling to explicitly test this hypothesis using genetic data. Here we analyse 46 ancient mitochondrial genomes from Late Pleistocene and Holocene bones of the eastern moa from across their post-LGM distribution. Our results are consistent with a post-LGM increase in the population size and genetic diversity of eastern moa. We also demonstrate that genetic diversity was higher in eastern moa from the southern extent of their range, supporting the hypothesis that they expanded from a single glacial refugium following the LGM.

Genetic diversity and population structure of the Western European hedgehog, Erinaceus europaeus: conservation status of populations in the Iberian Peninsula

Anthropogenic habitat fragmentation and roadkill mortality are considered important threats to European hedgehogs. Habitat fragmentation isolates hedgehog populations and, as a consequence, reduces their genetic diversity and leads the populations to vulnerable situations. The hedgehog populations in the Iberian Peninsula represent the southern limit of the species. We used microsatellite markers to estimate the genetic diversity and population structure of Erinaceus europaeus on the Iberian Peninsula. The obtained results indicated the presence of two differentiated groups, north-western and north-eastern, which coincided with the distribution of the two phylogeographic mitochondrial lineages described in the Peninsula. Moreover, in the north-eastern group, three genetically different clusters (Girona, Central Catalonia and Zoo) were identified. The highest genetic diversity (Hs = 0.696) was detected in the north-western region. Significant genetic differentiation (FST range = 0.072–0.224) was found among the clusters, indicating that these groups are well differentiated and present low gene flow. We concluded that the north-western group is genetically stable, whereas in the north-eastern region, despite some contact among groups, some populations are isolated and vulnerable.

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.

Biogeography of Italy revisited: genetic lineages confirm major phylogeographic patterns and a pre-Pleistocene origin of its biota

Warm-adapted species survived the cold stages of the past glacial-interglacial cycles in southern European peninsulas and recolonized major parts of Central and Northern Europe in the wake of postglacial warming. However, many of the genetic lineages which differentiated within these refugia predate even the Pleistocene. One of these refugia is the Italian Peninsula with its surrounding islands. In this meta-analysis, we compiled phylogeographic patterns of multiple species across this region. We transformed existing geographic information on 78 animal and plant species (with a total of 471 genetic lineages) within 17 predefined areas into presence/absence matrices. We elaborated three matrices: (i) only old lineages exclusively distinguished by deep splits, (ii) only young lineages distinguished by shallow (i.e. recent) splits, and (iii) presence/absence of the respective species. To infer biogeographic relationships between the predefined areas, we performed bootstrapped neighbour joining cluster analyses on these three matrices. In addition, we reviewed the geological history of Italy to identify causes of the observed biogeographic patterns. We found Sardinia and Corsica to be biogeographically closely linked with each other, and that they diverge strongly from all other regions. Sicily also diverges strongly from all other regions, while the intra-island differentiation was comparatively low. On the Italian mainland, Calabria exhibited the most pronounced biogeographic differentiation, often with several lineages present, resulting from old vicariance events within the region. Furthermore, southern Apulia and the Po Plain with adjoining areas of northern peninsular Italy displayed considerable distinctiveness. Admixture prevailed in the areas between these three regions. The ancient isolation of Sicily, as well as Sardinia plus Corsica, resulted in endemic lineages with only moderate recent exchange with adjacent mainland regions. Pronounced diversification occurs within the Italian Peninsula. The complex tectonic activities, such as shifting (micro)plates, submergence of major parts of peninsular Italy with the genesis of numerous Pliocene islands, in combination with the climatic cycles during the Pleistocene have most likely generated the current biogeographic pattern of species. Immigrations from the Balkan Peninsula into northern Italy partly accounted for the distinctiveness of this region.

Two ways to be endemic. Alps and Apennines are different functional refugia during climatic cycles

Endemics co-occur because they evolved in situ and persist regionally or because they evolved ex situ and later dispersed to shared habitats, generating evolutionary or ecological endemicity centres, respectively. We investigate whether different endemicity centres can intertwine in the region ranging from Alps to Sicily, by studying their butterfly fauna. We gathered an extensive occurrence data set for butterflies of the study area (27,123 records, 269 species, in cells of 0.5 × 0.5 degrees of latitude-longitude). We applied molecular-based delimitation methods (GMYC model) to 26,557 cytochrome c oxidase subunit 1 (COI) sequences of Western Palearctic butterflies. We identified entities based on molecular delimitations and/or the checklist of European butterflies and objectively attributed occurrences to their most probable entity. We obtained a zoogeographic regionalisation based on the 69 endemics of the area. Using phylogenetic ANOVA we tested if endemics from different centres differ from each other and from nonendemics for key ecological traits and divergence time. Endemicity showed high incidence in the Alps and Southern Italy. The regionalisation separated the Alps from the Italian Peninsula and Sicily. The endemics of different centres showed a high turnover and differed in phylogenetic distances, phenology and distribution traits. Endemics are on average younger than nonendemics and the Peninsula-Sicily endemics also have lower variance in divergence than those from the Alps. The observed variation identifies Alpine endemics as paleoendemics, now occupying an ecological centre, and the Peninsula-Sicily ones as neoendemics, that diverged in the region since the Pleistocene. The results challenge the common view of the Alpine-Apennine area as a single "Italian refugium".

Genetic diversity and structure of the hedgehogs Erinaceus europaeus and Erinaceus roumanicus: evidence for ongoing hybridization in Eastern Europe

Secondary contact zones between related species are key to understanding speciation mechanisms. The Central European sympatry zone of West European (Erinaceus europaeus) and northern white-breasted (Erinaceus roumanicus) hedgehogs is well studied, whereas data on the Eastern European sympatry zone are scarce. We examined the genetic variation in Russian populations using the mitochondrial Cytb gene, TTR intron 1 and 11 microsatellites to assess genetic variability and distribution patterns. In contrast to the Central European sympatry zone, we found evidence of ongoing hybridization between the two species in the sympatry zone of European Russia, where the proportion of individuals with mixed ancestry was c. 20%. Our data indicate bi-directional mtDNA introgression, but with a higher frequency of E. europaeus haplotypes in hybrids. The proportion of pure specimens with introgressed mitotypes is higher in E. roumanicus than in E. europaeus. Nuclear data showed the prevalence of the genetic contribution from E. roumanicus in admixed individuals. Demographic analyses indicated recent population growth in E. europaeus and little change in E. roumanicus, suggesting that E. europaeus colonized East Europe later than E. roumanicus.

3D Geometric Morphometrics Reveals Convergent Character Displacement in the Central European Contact Zone between Two Species of Hedgehogs (Genus Erinaceus)

Simple Summary

Hedgehogs, being insectivores with slow metabolisms, are quite sensitive to temperature and food availability. As a consequence, their ranges have oscillated in relation to past climate changes. Species that have evolved in different regions, but their ranges have shifted and overlapped subsequently, often represent intense competitors as a result of ecological similarities. The present study focuses on this phenomenon in the contact zone in central Europe and adjacent regions, using genetic determination of species and description of size and shape of skull, the morphological structure mirroring many selection pressures related to ecology. While animals living outside of the contact zone show marked differences between the two species, individuals within the contact zone are more alike with a smaller skull size and a convergent jawbone shape. Changes in skull size can be related to inter-species competition and also facilitated by selection pressure, mediated by overpopulated medium-sized predators such as foxes or badgers. Since the function of the lower jaw is mainly connected to feeding, we hypothesize that this pattern is due to the selection to size and shape related to competition for food resources. The present study helps to describe general patterns related to species formation, as well as species responses to anthropogenic environmental changes.

Abstract

Hedgehogs, as medium-sized plantigrade insectivores with low basal metabolic rates and related defensive anti-predator strategies, are quite sensitive to temperature and ecosystem productivity. Their ranges therefore changed dramatically due to Pleistocene climate oscillations, resulting in allopatric speciation and the subsequent formation of secondary contact zones. Such interactions between closely related species are known to generate strong evolutionary forces responsible for niche differentiation. In this connection, here, we detail the results of research on the phenotypic evolution in the two species of hedgehog present in central Europe, as based on genetics and geometric morphometrics in samples along a longitudinal transect that includes the contact zone between the species. While in allopatry, Erinaceus europaeus is found to have a larger skull than E. roumanicus and distinct cranial and mandibular shapes; the members of the two species in sympatry are smaller and more similar to each other, with a convergent shape of the mandible. The relevant data fail to reveal any major role for either hybridisation or clinal variation. We, therefore, hypothesise that competitive pressure exerted on the studied species does not generate divergent selection sufficient for divergent character displacement to evolve, instead giving rise to convergent selection in the face of resource limitation in the direction of smaller skull size. Considering the multi-factorial constraints present in the relevant adaptive landscape, reduction in size could also be facilitated by predator pressure in ecosystems characterised by mesopredator release and other anthropogenic factors. As the function of the animals’ lower jaw is mainly connected with feeding (in contrast to the cranium whose functions are obviously more complex), we interpret the similarity in shape as reflecting local adaptations to overlapping dietary resources in the two species and hence as convergent character displacement.

Evaluating the Impact of Wildlife Shelter Management on the Genetic Diversity of Erinaceus europaeus and E. roumanicus in Their Contact Zone

Simple Summary

Hedgehogs are regularly brought to wildlife shelters. Depending on the area from where animals are accepted, translocation can occur between different regions or populations. In this study, the genetic diversity of wild hedgehog populations was compared with “shelter populations” within central Europe focusing on the western contact zone between both European hedgehog species. Some shelters were hosting both species at the same time, in one this could be shown genetically. Generally, no difference in genetic diversity between shelter individuals and wild populations was found. Two shelters from Innsbruck hosted individuals that probably belong to two subpopulations. This indicates that shelter management-related translocations could facilitate gene flow across a dispersal barrier.

Abstract

Hedgehogs are among the most abundant species to be found within wildlife shelters and after successful rehabilitation they are frequently translocated. The effects and potential impact of these translocations on gene flow within wild populations are largely unknown. In this study, different wild hedgehog populations were compared with artificially created “shelter populations”, with regard to their genetic diversity, in order to establish basic data for future inferences on the genetic impact of hedgehog translocations. Observed populations are located within central Europe, including the species Erinaceus europaeus and E. roumanicus. Shelters were mainly hosting one species; in one case, both species were present syntopically. Apart from one exception, the results did not show a higher genetic diversity within shelter populations, indicating that individuals did not originate from a wider geographical area than individuals grouped into one of the wild populations. Two shelters from Innsbruck hosted individuals that belonged to two potential clusters, as indicated in a distance analysis. When such a structure stems from the effects of landscape elements like large rivers, the shelter management-related translocations might lead to homogenization across the dispersal barrier.

Last Glacial Maximum led to community-wide population expansion in a montane songbird radiation in highland Papua New Guinea

BACKGROUND

Quaternary climate fluctuations are an engine of biotic diversification. Global cooling cycles, such as the Last Glacial Maximum (LGM), are known to have fragmented the ranges of higher-latitude fauna and flora into smaller refugia, dramatically reducing species ranges. However, relatively less is known about the effects of cooling cycles on tropical biota.

RESULTS

We analyzed thousands of genome-wide DNA markers across an assemblage of three closely related understorey-inhabiting scrubwrens (Sericornis and Aethomyias; Aves) from montane forest along an elevational gradient on Mt. Wilhelm, the highest mountain of Papua New Guinea. Despite species-specific differences in elevational preference, we found limited differentiation within each scrubwren species, but detected a strong genomic signature of simultaneous population expansions at 27-29 ka, coinciding with the onset of the LGM.

CONCLUSION

The remarkable synchronous timing of population expansions of all three species demonstrates the importance of global cooling cycles in expanding highland habitat. Global cooling cycles have likely had strongly different impacts on tropical montane areas versus boreal and temperate latitudes, leading to population expansions in the former and serious fragmentation in the latter.

Genetic structure of the European hedgehog (Erinaceus europaeus) in Denmark

OBJECTIVES

Low genetic diversity can lead to reduced average fitness in a population or even extinction. Preserving genetic connectivity across fragmented landscapes is therefore vital to counteract the negative consequences of genetic drift and inbreeding. This study aimed to assess the genetic composition and consequently the conservation status of a nationwide sample of European hedgehogs (Erinaceus europaeus) in Denmark.

METHODS

We applied an adaptation of the genotyping by sequencing (GBS) technique to 178 individuals from six geographically distinct populations. We used a Bayesian clustering method to subdivide individuals into genetically distinct populations. We estimated individual observed (iHO), observed (HO), and unbiased expected (uHE) heterozygosity, inbreeding coefficient (FIS), percentage of polymorphic loci (P%) and tested for deviations from Hardy-Weinberg equilibrium (HWE). We used linear models to test for potential anthropogenic effects on the genetic variability of hedgehogs with iHO, uHE, P% and FIS as response variables, and assessed the demographic history of the population.

RESULTS

The Danish hedgehog population is composed of three genetic clusters. We found a mean P% of 54.44-94.71, a mean uHE of 0.126-0.318 and a mean HO of 0.124-0.293 in the six populations. The FIS was found to be significantly positive for three of the six populations. We detected a large heterogeneity of iHO values within populations, which can be due to inbreeding and/or fragmentation. FIS values decreased with increasing farmland density, but there was no significant association with human population or road density.

CONCLUSIONS

We found a low level of genetic variability and evidence for genetic substructure and low effective population size, which are all consequences of habitat fragmentation. We failed to detect signs of a recent population bottleneck or population increase or decline. However, because the test only identifies recent changes in population size, we cannot reject the possibility of a longer-term decline in the Danish hedgehog population.

Fluctuating fortunes: genomes and habitat reconstructions reveal global climate-mediated changes in bats' genetic diversity

Over the last approximately 2.6 Myr, Earth's climate has been dominated by cyclical ice ages that have profoundly affected species' population sizes, but the impact of impending anthropogenic climate change on species' extinction potential remains a worrying problem. We investigated 11 bat species from different taxonomic, ecological and geographical backgrounds using combined information from palaeoclimatic habitat reconstructions and genomes to analyse biotic impacts of historic climate change. We discover tightly correlated fluctuations between species' historic distribution and effective population size, identify frugivores as particularly susceptible to global warming, pinpoint large insectivores as having overall low effective population size and flag the onset of the Holocene (approx. 10-12 000 years ago) as the period with the generally lowest effective population sizes across the last approximately 1 Myr. Our study shows that combining genomic and palaeoclimatological approaches reveals effects of climatic shifts on genetic diversity and may help predict impacts of future climate change.

Cross-continental phylogeography of two Holarctic Nymphalid butterflies, Boloria eunomia and Boloria selene

Pleistocene glaciations had significant effects on the distribution and evolution of species inhabiting the Holarctic region. Phylogeographic studies concerning the entire region are still rare. Here, we compared global phylogeographic patterns of one boreo-montane and one boreo-temperate butterflies with largely overlapping distribution ranges across the Northern Hemisphere, but with different levels of range fragmentation and food specialization. We reconstructed the global phylogeographic history of the boreo-montane specialist Boloria eunomia (n = 223) and of the boreo-temperate generalist Boloria selene (n = 106) based on mitochondrial and nuclear DNA markers, and with species distribution modelling (SDM). According to the genetic structures obtained, both species show a Siberian origin and considerable split among populations from Nearctic and Palaearctic regions. According to SDMs and molecular data, both butterflies could inhabit vast areas during the moderate glacials. In the case of B. selene, high haplotype diversity and low geographic structure suggest long-lasting interconnected gene flow among populations. A stronger geographic structuring between populations was identified in the specialist B. eunomia, presumably due to the less widespread, heterogeneously distributed food resources, associated with cooler and more humid climatic conditions. Populations of both species show opposite patterns across major parts of North America and in the case of B. eunomia also across Asia. Our data underline the relevance to cover entire distribution ranges to reconstruct the correct phylogeographic history of species.

Application of a SSR-GBS marker system on investigation of European Hedgehog species and their hybrid zone dynamics

By applying second-generation sequencing technologies to microsatellite genotyping, sequence information is produced which can result in high-resolution population genetics analysis populations and increased replicability between runs and laboratories. In the present study, we establish an approach to study the genetic structure patterns of two European hedgehog species Erinaceaus europaeus and E. roumanicus. These species are usually associated with human settlements and are good models to study anthropogenic impacts on the genetic diversity of wild populations. The short sequence repeats genotyping by sequence (SSR-GBS) method presented uses amplicon sequences to determine genotypes for which allelic variants can be defined according to both length and single nucleotide polymorphisms (SNPs). To evaluate whether complete sequence information improved genetic structure definition, we compared this information with datasets based solely on length information. We identified a total of 42 markers which were successfully amplified in both species. Overall, genotyping based on complete sequence information resulted in a higher number of alleles, as well as greater genetic diversity and differentiation between species. Additionally, the structure patterns were slightly clearer with a division between both species and some potential hybrids. There was some degree of genetic structure within species, although only in E. roumanicus was this related to geographical distance. The statistically significant results obtained by SSR-GBS demonstrate that it is superior to electrophoresis-based methods for SSR genotyping. Moreover, the greater reproducibility and throughput with lower effort which can be obtained with SSR-GBS and the possibility to include degraded DNA into the analysis, allow for continued relevance of SSR markers during the genomic era.

Evolutionary and ecological patterns within the South African Bathyergidae: Implications for taxonomy

The family Bathyergidae (comprising six genera) is a group of subterranean rodents endemic to sub-Saharan Africa. Our understanding of the evolution and species richness of the South African bathyergid genera Georychus, Bathyergus and Cryptomys is limited, with the majority of species listed as Least Concern by the IUCN Red List of Threatened Species. Genetic data suggest that several cryptic species may be present in these genera. To explore genetic and ecological distinctiveness, and evaluate taxonomic richness across the ranges of Georychus, Bathyergus and to a lesser degree, Cryptomys, as well as evaluate possible scenarios which have historically influenced evolutionary patterns, we employed four protein coding markers (one mitochondrial and three nuclear) along with distribution wide sampling schemes and large sample sizes. In addition, possible ecological differences among the different intra-generic clades were explored. Genera appear to have originated in the north-eastern interior of South Africa, following novel habitats created through the Post-African I erosion cycle and dramatic changes in climate and phytogeography. In each genus, multiple geographically discrete genetic lineages (clades) are supported by both the mitochondrial and nuclear data. These lineages bear signature of the fragmentation of wider historical distributions through major environmental changes since the middle Miocene (major uplift events, Post-African II erosion cycle, drainage evolution of major river systems, sea-level fluctuations as well as climatic changes and vegetation shifts), thereby leading to long-term isolation. Along with protracted periods of separation, it appears that ecological differences further delimit the lineages in relation to geology, phytogeographic preference, elevation, rainfall and temperature. As such, two lineages in Georychus (Clades 1 and 2) and one lineages in Cryptomys (Clade I) occur at higher elevations above the Great Escarpment (in older deposits harbouring grassland vegetation, with higher rainfall and lower daily temperatures), with the remaining lineages within these genera (Clades 3, 4 and 5 in Georychus and Clades III and IV in Cryptomys) occupying a low-land distribution with contrasting climatic and geological characteristics. Although significant differences in ecological variables were also observed between Bathyergus clades, these were not consistent, given their largely low-land distributions. Our results corroborate and expand previous suggestions that several cryptic species are present within the South African Bathyergidae.

Pleistocene land bridges act as semipermeable agents of avian gene flow in Wallacea

Cyclical periods of global cooling have been important drivers of biotic differentiation throughout the Quaternary. Ice age-induced sea level fluctuations can lead to changing patterns of land connections, both facilitating and disrupting gene flow. In this study, we test if species with differing life histories are differentially affected by Quaternary land connections. We used genome-wide SNPs in combination with mitochondrial gene sequences to analyse levels of divergence and gene flow between two songbird complexes across two Wallacean islands that have been repeatedly connected during glaciations. Although the two bird complexes are similar in ecological attributes, the forest and edge-inhabiting golden whistler Pachycephala pectoralis is comparatively flexible in its diet and niche requirements as compared to the henna-tailed jungle-flycatcher Cyornis colonus, which is largely restricted to the forest interior. Using population-genomic and coalescent approaches, we estimated levels of gene flow, population differentiation and divergence time between the two island populations. We observed higher levels of differentiation, an approximately two to four times deeper divergence time and near-zero levels of gene flow between the two island populations of the more forest-dependent henna-tailed jungle-flycatcher as compared to the more generalist golden whistler. Our results suggest that Quaternary land bridges act as semipermeable agents of gene flow in Wallacea, allowing only certain taxa to connect between islands while others remain isolated. Quaternary land bridges do not accommodate all terrestrial species equally, differing in suitability according to life history and species biology. More generalist species are likely to use Quaternary land connections as a conduit for gene flow between islands whereas island populations of more specialist species may continue to be reproductively isolated even during periods of Quaternary land bridges.

Selection on MHC in a Context of Historical Demographic Change in 2 Closely Distributed Species of Tuco-tucos (Ctenomys australis and C. talarum)

Selection necessarily acts within the same current and historical demographic framework as neutral evolutionary processes, and the outcome of the interplay between these forces may vary according to their relative strength. In this study, we compare the variation at a major histocompatibility complex (MHC) locus (DRB exon 2), typically subject to strong diversifying selection, and mitochondrial diversity (D-loop) across populations encompassing the entire distribution of 2 species of South American subterranean rodents: Ctenomys australis and C. talarum (tuco-tucos). Although these species are parapatric along most of their distribution, historically they have followed distinct demographic trajectories associated with sea level changes during the Quaternary. We surveyed 8 populations of C. australis and 15 of C. talarum, from which we analyzed 70 and 212 D-loop haplotypes and 91 and 346 DRB genotypes, respectively. Both species have gone through a recent demographic expansion; however, the signal of this process only encompasses the entire distribution of one of the species: C. australis. While balancing selection on MHC in C. talarum-enhanced DRB diversity at the local level compared to D-loop, although not promoting divergence among populations, in C. australis local diversifying selection may have driven higher population differentiation at DRB than at D-loop. Our findings reinforce the idea that the relative strength of selection acting on MHC genes varies spatially and temporally within and among species, even between species using the same macrohabitat and exposed to similar immune challenges.

Gene flow during glacial habitat shifts facilitates character displacement in a Neotropical flycatcher radiation

BACKGROUND

Pleistocene climatic fluctuations are known to be an engine of biotic diversification at higher latitudes, but their impact on highly diverse tropical areas such as the Andes remains less well-documented. Specifically, while periods of global cooling may have led to fragmentation and differentiation at colder latitudes, they may - at the same time - have led to connectivity among insular patches of montane tropical habitat with unknown consequences on diversification. In the present study we utilized ~5.5 kb of DNA sequence data from eight nuclear loci and one mitochondrial gene alongside diagnostic morphological and bioacoustic markers to test the effects of Pleistocene climatic fluctuations on diversification in a complex of Andean tyrant-flycatchers of the genus Elaenia.

RESULTS

Population genetic and phylogenetic approaches coupled with coalescent simulations demonstrated disparate levels of gene flow between the taxon chilensis and two parapatric Elaenia taxa predominantly during the last glacial period but not thereafter, possibly on account of downward shifts of montane forest habitat linking the populations of adjacent ridges. Additionally, morphological and bioacoustic analyses revealed a distinct pattern of character displacement in coloration and vocal traits between the two sympatric taxa albiceps and pallatangae, which were characterized by a lack of gene flow.

CONCLUSION

Our study demonstrates that global periods of cooling are likely to have facilitated gene flow among Andean montane Elaenia flycatchers that are more isolated from one another during warm interglacial periods such as the present era. We also identify a hitherto overlooked case of plumage and vocal character displacement, underpinning the complexities of gene flow patterns caused by Pleistocene climate change across the Andes.

Glacial allopatry vs. postglacial parapatry and peripatry: the case of hedgehogs

Although hedgehogs are well-known examples of postglacial recolonisation, the specific processes that shape their population structures have not been examined by detailed sampling and fast-evolving genetic markers in combination with model based clustering methods. This study aims to analyse the impacts of isolation within glacial refugia and of postglacial expansion on the population structure of the Northern White-breasted hedgehog (Erinaceus roumanicus). It also discusses the role of the processes at edges of species distribution in its evolutionary history. The maternally inherited mitochondrial control region and the bi-parentally inherited nuclear microsatellites were used to examine samples within the Central Europe, Balkan Peninsula and adjacent islands. Bayesian coalescent inference and neutrality tests proposed a recent increase in the population size. The most pronounced pattern of population structure involved differentiation of the insular populations in the Mediterranean Sea and the population within the contact zone with E. europaeus in Central Europe. An interspecies hybrid was detected for the first time in Central Europe. A low genetic diversity was observed in Crete, while the highest genetic distances among individuals were found in Romania. The recent population in the post-refugial area related to the Balkan Peninsula shows a complex pattern with pronounced subpopulations located mainly in the Pannonian Basin and at the Adriatic and Pontic coasts. Detailed analyses indicate that parapatry and peripatry may not be the only factors that limit range expansion, but also strong microevolutionary forces that may change the genetic structure of the species. Here we present evidence showing that population differentiation may occur not only during the glacial restriction of the range into the refugia, but also during the interglacial range expansion. Population differentiation at the Balkan Peninsula and adjacent regions could be ascribed to diversification in steppe/forest biomes and complicated geomorphology, including pronounced geographic barriers as Carpathians.

Strong population structure in a species manipulated by humans since the Neolithic: the European fallow deer (Dama dama dama)

Species that have been translocated and otherwise manipulated by humans may show patterns of population structure that reflect those interactions. At the same time, natural processes shape populations, including behavioural characteristics like dispersal potential and breeding system. In Europe, a key factor is the geography and history of climate change through the Pleistocene. During glacial maxima throughout that period, species in Europe with temperate distributions were forced south, becoming distributed among the isolated peninsulas represented by Anatolia, Italy and Iberia. Understanding modern patterns of diversity depends on understanding these historical population dynamics. Traditionally, European fallow deer (Dama dama dama) are thought to have been restricted to refugia in Anatolia and possibly Sicily and the Balkans. However, the distribution of this species was also greatly influenced by human-mediated translocations. We focus on fallow deer to better understand the relative influence of these natural and anthropogenic processes. We compared modern fallow deer putative populations across a broad geographic range using microsatellite and mitochondrial DNA loci. The results revealed highly insular populations, depauperate of genetic variation and significantly differentiated from each other. This is consistent with the expectations of drift acting on populations founded by small numbers of individuals, and reflects known founder populations in the north. However, there was also evidence for differentiation among (but not within) physically isolated regions in the south, including Iberia. In those regions we find evidence for a stronger influence from natural processes than may be expected for a species with such strong, known anthropogenic influence.

Extra-Mediterranean glacial refugia in a Mediterranean faunal element: the phylogeography of the chalk-hill blue Polyommatus coridon (Lepidoptera, Lycaenidae)

Most warm-adapted Central European species are thought to have survived ice ages exclusively in Mediterranean refugia. During recent years, this point of view has been questioned. Therefore, we tested the hypothesis that extra-Mediterranean refugia also played a role in warm-adapted insect species and selected the chalk-hill blue, Polyommatus coridon. We sequenced two mitochondrial loci (COI, CR) in 150 individuals from 30 populations covering nearly the complete range. Minimum spanning networks and other statistical analyses concordantly revealed four genetic lineages with strong phylogeographic signal: a western group in Italy, France and western/central Germany, an eastern lineage in the Balkan Peninsula, the Carpathian Basin and eastern Central Europe, an Alpine group with populations in the Alps and southern Germany and a Pyrenean group. Our results are generally consistent with previous analyses for P. coridon based on allozymes and DNA sequences, but provide additional insights. We propose that these four lineages have evolved during allopatry in different glacial refugia, two in typical Mediterranean refugia (Apennines and Balkan Peninsulas), but two in extra-Mediterranean areas south of the Alps and Pyrenees. This supports survival of warm-adapted organisms in these regions in close geographic proximity to the refugia of high mountain species.

A combination of long term fragmentation and glacial persistence drove the evolutionary history of the Italian wall lizard Podarcis siculus

BACKGROUND

The current distribution of genetic diversity is the result of a vast array of microevolutionary processes, including short-term demographic and ecological mechanisms and long-term allopatric isolation in response to Quaternary climatic fluctuations. We investigated past processes that drove the population differentiation and spatial genetic distribution of the Italian wall lizard Podarcis siculus by means of sequences of mitochondrial cytb (n = 277 from 115 localities) and nuclear mc1r and β-fibint7genes (n = 262 and n = 91, respectively) from all its distribution range. The pattern emerging from the genetic data was compared with current and past (last glacial maximum) species distribution modeling (SDM).

RESULTS

We identified seven deeply divergent parapatric clades which presumably remained isolated in different refugia scattered mainly throughout the Tyrrhenian coast. Conversely, the Adriatic coast showed only two haplogroups with low genetic variability. These results appear to agree with the SDM prediction at the last glacial maximum (LGM) indicating a narrow area of habitat suitability along the Tyrrhenian coast and much lower suitability along the Adriatic one. However, the considerable land exposure of the Adriatic coastline favored a glacial colonization of the Balkan Peninsula.

CONCLUSIONS

Our population-level historical demography showed a common trend consistent with glacial expansions and regional persistence during the last glacial maximum. This complex genetic signature appears to be inconsistent with the expectation of the expansion-contraction model and post-LGM (re)colonizations from southern refugia. Hence it is one of an increasing number of cases in which these assumptions are not met, indicating that long-term fragmentation and pre-LGM events such as glacial persistence were more prominent in shaping genetic variation in this temperate species.

Genetic population structure and relatedness in the narrow-striped mongoose (Mungotictis decemlineata), a social Malagasy carnivore with sexual segregation

Information on the genetic structure of animal populations can allow inferences about mechanisms shaping their social organization, dispersal, and mating system. The mongooses (Herpestidae) include some of the best-studied mammalian systems in this respect, but much less is known about their closest relatives, the Malagasy carnivores (Eupleridae), even though some of them exhibit unusual association patterns. We investigated the genetic structure of the Malagasy narrow-striped mongoose (Mungotictis decemlineata), a small forest-dwelling gregarious carnivore exhibiting sexual segregation. Based on mtDNA and microsatellite analyses, we determined population-wide haplotype structure and sex-specific and within-group relatedness. Furthermore, we analyzed parentage and sibship relationships and the level of reproductive skew. We found a matrilinear population structure, with several neighboring female units sharing identical haplotypes. Within-group female relatedness was significantly higher than expected by chance in the majority of units. Haplotype diversity of males was significantly higher than in females, indicating male-biased dispersal. Relatedness within the majority of male associations did not differ from random, not proving any kin-directed benefits of male sociality in this case. We found indications for a mildly promiscuous mating system without monopolization of females by males, and low levels of reproductive skew in both sexes based on parentages of emergent young. Low relatedness within breeding pairs confirmed immigration by males and suggested similarities with patterns in social mongooses, providing a starting point for further investigations of mate choice and female control of reproduction and the connected behavioral mechanisms. Our study contributes to the understanding of the determinants of male sociality in carnivores as well as the mechanisms of female competition in species with small social units.

Different parasite faunas in sympatric populations of sister hedgehog species in a secondary contact zone

Providing descriptive data on parasite diversity and load in sister species is a first step in addressing the role of host-parasite coevolution in the speciation process. In this study we compare the parasite faunas of the closely related hedgehog species Erinaceus europaeus and E. roumanicus from the Czech Republic where both occur in limited sympatry. We examined 109 hedgehogs from 21 localities within this secondary contact zone. Three species of ectoparasites and nine species of endoparasites were recorded. Significantly higher abundances and prevalences were found for Capillaria spp. and Brachylaemus erinacei in E. europaeus compared to E. roumanicus and higher mean infection rates and prevalences for Hymenolepis erinacei, Physaloptera clausa and Nephridiorhynchus major in E. roumanicus compared to E. europaeus. Divergence in the composition of the parasite fauna, except for Capillaria spp., which seem to be very unspecific, may be related to the complicated demography of their hosts connected with Pleistocene climate oscillations and consequent range dynamics. The fact that all parasite species with different abundances in E. europaeus and E. roumanicus belong to intestinal forms indicates a possible diversification of trophic niches between both sister hedgehog species.

Cryptic speciation and chromosomal repatterning in the South African climbing mice Dendromus (Rodentia, Nesomyidae)

We evaluate the intra- and interspecific diversity in the four South African rodent species of the genus Dendromus. The molecular phylogenetic analysis on twenty-three individuals have been conducted on a combined dataset of nuclear and mitochondrial markers. Moreover, the extent and processes underlying chromosomal variation, have been investigated on three species by mean of G-, C-bands, NORs and Zoo-FISH analysis. The molecular analysis shows the presence of six monophyletic lineages corresponding to D. mesomelas, D. mystacalis and four lineages within D. cfr. melanotis with high divergence values (ranges: 10.6% – 18.3%) that raises the question of the possible presence of cryptic species. The first description of the karyotype for D. mesomelas and D. mystacalis and C- and G- banding for one lineage of D. cfr. melanotis are reported highlighting an extended karyotype reorganization in the genus. Furthermore, the G-banding and Zoo-FISH evidenced an autosome-sex chromosome translocation characterizing all the species and our timing estimates this mutation date back 7.4 mya (Late Miocene). Finally, the molecular clock suggests that cladogenesis took place since the end of Miocene to Plio-Pleistocene, probably due to ecological factors, isolation in refugia followed by differential adaptation to the mesic or dry habitat.

Mitochondrial and nuclear genes-based phylogeography of Arvicanthis niloticus (Murinae) and sub-Saharan open habitats pleistocene history

A phylogeographic study was conducted on the Nile grass rat, Arvicanthis niloticus, a rodent species that is tightly associated with open grasslands from the Sudano-Sahelian regions. Using one mitochondrial (cytochrome b) and one nuclear (intron 7 of Beta Fibrinogen) gene, robust patterns were retrieved that clearly show that (i) the species originated in East Africa concomitantly with expanding grasslands some 2 Ma, and (ii) four parapatric and genetically well-defined lineages differentiated essentially from East to West following Pleistocene bioclimatic cycles. This strongly points towards allopatric genetic divergence within savannah refuges during humid episodes, then dispersal during arid ones; secondary contact zones would have then stabilized around geographic barriers, namely, Niger River and Lake Chad basins. Our results pertinently add to those obtained for several other African rodent as well as non-rodent species that inhabit forests, humid zones, savannahs and deserts, all studies that now allow one to depict a more comprehensive picture of the Pleistocene history of the continent south of the Sahara. In particular, although their precise location remains to be determined, at least three Pleistocene refuges are identified within the West and Central African savannah biome.

A range-wide synthesis and timeline for phylogeographic events in the red fox (Vulpes vulpes)

BACKGROUND

Many boreo-temperate mammals have a Pleistocene fossil record throughout Eurasia and North America, but only few have a contemporary distribution that spans this large area. Examples of Holarctic-distributed carnivores are the brown bear, grey wolf, and red fox, all three ecological generalists with large dispersal capacity and a high adaptive flexibility. While the two former have been examined extensively across their ranges, no phylogeographic study of the red fox has been conducted across its entire Holarctic range. Moreover, no study included samples from central Asia, leaving a large sampling gap in the middle of the Eurasian landmass.

RESULTS

Here we provide the first mitochondrial DNA sequence data of red foxes from central Asia (Siberia), and new sequences from several European populations. In a range-wide synthesis of 729 red fox mitochondrial control region sequences, including 677 previously published and 52 newly obtained sequences, this manuscript describes the pattern and timing of major phylogeographic events in red foxes, using a Bayesian coalescence approach with multiple fossil tip and root calibration points. In a 335 bp alignment we found in total 175 unique haplotypes. All newly sequenced individuals belonged to the previously described Holarctic lineage. Our analyses confirmed the presence of three Nearctic- and two Japan-restricted lineages that were formed since the Mid/Late Pleistocene.

CONCLUSIONS

The phylogeographic history of red foxes is highly similar to that previously described for grey wolves and brown bears, indicating that climatic fluctuations and habitat changes since the Pleistocene had similar effects on these highly mobile generalist species. All three species originally diversified in Eurasia and later colonized North America and Japan. North American lineages persisted through the last glacial maximum south of the ice sheets, meeting more recent colonizers from Beringia during postglacial expansion into the northern Nearctic. Both brown bears and red foxes colonized Japan's northern island Hokkaido at least three times, all lineages being most closely related to different mainland lineages. Red foxes, grey wolves, and brown bears thus represent an interesting case where species that occupy similar ecological niches also exhibit similar phylogeographic histories.

Population biology of establishment in New Zealand hedgehogs inferred from genetic and historical data: conflict or compromise?

The crucial steps in biological invasions, related to the shaping of genetic architecture and the current evolution of adaptations to a novel environment, usually occur in small populations during the phases of introduction and establishment. However, these processes are difficult to track in nature due to invasion lag, large geographic and temporal scales compared with human observation capabilities, the frequent depletion of genetic variance, admixture and other phenomena. In this study, we compared genetic and historical evidence related to the invasion of the West European hedgehog to New Zealand to infer details about the introduction and establishment. Historical information indicates that the species was initially established on the South Island. A molecular assay of populations from Great Britain and New Zealand using mitochondrial sequences and nuclear microsatellite loci was performed based on a set of analyses including approximate Bayesian computation, a powerful approach for disentangling complex population demographies. According to these analyses, the population of the North Island was most similar to that of the native area and showed greatest reduction in genetic variation caused by founder demography and/or drift. This evidence indicated the location of the establishment phase. The hypothesis was corroborated by data on climate and urbanization. We discuss the contrasting results obtained by the molecular and historical approaches in the light of their different explanatory power and the possible biases influencing the description of particular aspects of invasions, and we advocate the integration of the two types of approaches in invasion biology.