Genetic and morphological variation in a Mediterranean glacial refugium: evidence from Italian pygmy shrews, Sorex minutus (Mammalia: Soricomorpha)

Unveiling a hotspot of genetic diversity in southern Italy for the endangered Hermann’s tortoise Testudo hermanni

Background

Hotspots of intraspecific genetic diversity represent invaluable resources for species to cope with environmental changes, and their identification is increasingly recognized as a major goal of conservation ecology research. However, even for iconic and endangered species, conservation strategies are often planned without thorough information on the geographic patterns of genetic variation. Here, we investigated the spatial patterns of genetic variation of the endangered Hermann’s tortoise Testudo hermanni in the Italian Peninsula by genotyping 174 individuals at 7 microsatellite loci, with the aim to contribute to planning effective conservation strategies.

Results

Ordination-based and Bayesian clustering analyses consistently identified three main genetic clusters, one spread in the central and northern part of the peninsula, and two restricted to southern Italy and Sicily, respectively. The highest levels of genetic diversity were found in populations of the southern cluster and, in particular, at the northern edges of its distribution (He > 0.6, Ar > 2.8 ), that correspond to areas of putative secondary contact and admixture between distinct lineages. Our results clearly identify a hotspot of genetic diversity for the Hermann’s tortoise in southern Italy.

Conclusion

We inferred the evolutionary history and the spatial patterns of genetic variation of the Hermann’s tortoise in the Italian Peninsula. We identified three main genetic clusters along the peninsula and a hotspot of intraspecific diversity in southern Italy. Our results underline the urgent need for conservation actions to warrant the long-term persistence of viable tortoise populations in this area. Furthrmore, these data add further evidence to the role of southern Italy as a biodiversity hotspot for temperate fauna, claiming for higher consideration of this area in large scale conservation programs.

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.

Are the western water shrew (Sorex navigator) and American water shrew (Sorex palustris) morphologically distinct?

Recent phylogenetic studies applying the mitochondrial cytochrome b (Cytb) gene and nuclear genes revealed divergent cordilleran and boreal lineages in western populations of the Sorex palustris complex that are now recognized as distinct species: Sorex navigator (Baird, 1858) (western water shrew) and Sorex palustris Richardson, 1828 (American water shrew). Their distributions in western Canada are parapatric with a potential contact zone in the foothills of the Rocky Mountains. Using 80 specimens including 69 of known genotype, we assessed morphological differentiation in metric and nonmetric traits among S. navigator, S. palustris, and Sorex bendirii (Merriam, 1894) (marsh shrew), a sister species of S. navigator. Multivariate analyses revealed discrete morpho groups concordant with their genetic lineages and species-level divergence. Three individuals from northwestern British Columbia and one from southwestern Alberta with discordant mitochondrial and nuclear DNA aligned with the S. navigator morpho group showing no evidence for morphological intergradation with S. palustris. Although differentiated in genotype, S. navigator and S. palustris are cryptic species that can only be discriminated from morphology with metric cranio-dental data. We developed diagnostic morphometric criteria to identify specimens lacking genetic data.

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.

Genetic diversity and structure of Dalmatian pyrethrum (Tanacetum cinerariifolium Trevir./Sch./Bip., Asteraceae) within the Balkan refugium

Dalmatian pyrethrum (Tanacetum cinerariifolium Trevir. /Sch./ Bip.) is an outcrossing, perennial insecticidal plant, restricted to the eastern Adriatic coast (Mediterranean). Amplified fragment-length polymorphisms (AFLP) were used to investigate the genetic diversity and structure within and among 20 natural plant populations. The highest level of gene diversity, the number of private alleles and the frequency down-weighted marker values (DW) were found in northern Adriatic populations and gradually decreased towards the southern boundary of the species range. Genetic impoverishment of these southern populations is most likely the result of human-related activities. An analysis of molecular variance (AMOVA) indicated that most of the genetic diversity was attributed to differences among individuals within populations (85.78%), which are expected due to the outcrossing nature of the species. A Bayesian analysis of the population structure identified two dominant genetic clusters. A spatial analysis of the genetic diversity indicated that 5.6% of the genetic differentiation resulted from isolation by distance (IBD), while 12.3% of the genetic differentiation among populations followed the pattern of isolation by environmental distance (IBED). Knowledge of the genetic diversity patterns of the natural populations and the mechanism behind these patterns is required for the exploitation and possible conservation management of this endemic and economically important species.

Microscale vicariance and diversification of Western Balkan caddisflies linked to karstification

The karst areas in the Dinaric region of the Western Balkan Peninsula are a hotspot of freshwater biodiversity. Many investigators have examined diversification of the subterranean freshwater fauna in these karst systems. However, diversification of surface-water fauna remains largely unexplored. We assessed local and regional diversification of surface-water species in karst systems and asked whether patterns of population differentiation could be explained by dispersal-diversification processes or allopatric diversification following karst-related microscale vicariance. We analyzed mitochondrial cytochrome c oxidase subunit I (mtCOI) sequence data of 4 caddisfly species (genus Drusus) in a phylogeographic framework to assess local and regional population genetic structure and Pliocene/Pleistocene history. We used BEAST software to assess the timing of intraspecific diversification of the target species. We compared climate envelopes of the study species and projected climatically suitable areas during the last glacial maximum (LGM) to assess differences in the species climatic niches and infer potential LGM refugia. The haplotype distribution of the 4 species (324 individuals from 32 populations) was characterized by strong genetic differentiation with few haplotypes shared among populations (16%) and deep divergence among populations of the 3 endemic species, even at local scales. Divergence among local populations of endemics often exceeded divergence among regional and continental clades of the widespread D. discolor. Major divergences among regional populations dated to 2.0 to 0.5 Mya. Species distribution model projections and genetic structure suggest that the endemic species persisted in situ and diversified locally throughout multiple Pleistocene climate cycles. The pattern for D. discolor was different and consistent with multiple invasions into the region. Patterns of population genetic structure and diversification were similar for the 3 regional endemic Drusus species and consistent with microscale vicariance after the onset of intensified karstification in the Dinaric region. Karstification may induce microscale vicariance of running surface-water habitats and probably promotes allopatric fragmentation of stream insects at small spatial scales.

The integration of multiple independent data reveals an unusual response to Pleistocene climatic changes in the hard tick Ixodes ricinus

In the last few years, improved analytical tools and the integration of genetic data with multiple sources of information have shown that temperate species exhibited more complex responses to ice ages than previously thought. In this study, we investigated how Pleistocene climatic changes affected the current distribution and genetic diversity of European populations of the tick Ixodes ricinus, an ectoparasite with high ecological plasticity. We first used mitochondrial and nuclear genetic markers to investigate the phylogeographic structure of the species and its Pleistocene history using coalescent-based methods; then we used species distribution modelling to infer the climatic niche of the species at last glacial maximum; finally, we reviewed the literature on the I. ricinus hosts to identify the locations of their glacial refugia. Our results support the scenario that during the last glacial phase, I. ricinus never experienced a prolonged allopatric divergence in separate glacial refugia, but persisted with interconnected populations across Southern and Central Europe. The generalist behaviour in host choice of I. ricinus would have played a major role in maintaining connections between its populations. Although most of the hosts persisted in separate refugia, from the point of view of I. ricinus, they represented a continuity of 'bridges' among populations. Our study highlights the importance of species-specific ecology in affecting responses to Pleistocene glacial-interglacial cycles. Together with other cases in Europe and elsewhere, it contributes to setting new hypotheses on how species with wide ecological plasticity coped with Pleistocene climatic changes.

Multiple lines of evidence for demographic and range expansion of a temperate species (Hyla sarda) during the last glaciation

Many temperate species experienced demographic and range contractions in response to climatic changes during Pleistocene glaciations. In this study, we investigate the evolutionary history of the Tyrrhenian tree frog Hyla sarda, a species inhabiting the Corsica-Sardinia island system (Western Mediterranean basin). We used sequence analysis of two mitochondrial (overall 1229 bp) and three nuclear (overall 1692 bp) gene fragments to assess the phylogeography and demographic history of this species, and species distribution modelling (SDM) to predict its range variation over time. Phylogeographic, historical demographic and SDM analyses consistently indicate that H. sarda does not conform to the scenario generally expected for temperate species but rather underwent demographic and range expansion mostly during the last glacial phase. Palaeogeographic data and SDM analyses suggest that such expansion was driven by the glaciation-induced increase in lowland areas during marine regression. This unusual scenario suggests that at least some temperate species may not have suffered the adverse effects of glacial climate on their population size and range extent, owing to the mitigating effects of other glaciations-induced palaeoenvironmental changes. We discuss previous clues for the occurrence of such a scenario in other species and some possible challenges with its identification. Early phylogeographic literature suggested that responses to the Pleistocene glacial-interglacial cycles were expected to vary among species and regions. Our results point out that such variation may have been greater than previously thought.

Intricate variation patterns in the diploid-polyploid complex of Alyssum montanum-A. repens (Brassicaceae) in the Apennine Peninsula: evidence for long-term persistence and diversification

PREMISE OF THE STUDY

The Apennine Peninsula, one of the three peninsulas of the European Mediterranean, is an important hotspot of genetic and species diversity, but studies devoted to plant evolution are still very scarce in this region. We studied the diploid-polyploid complex of Alyssum montanum-A. repens, focusing on Apennine and adjacent southwestern Alpine populations from southeastern France, with the aim of examining their taxonomic position and evolutionary patterns.

METHODS

We employed AFLP markers and cpDNA sequences, along with cytotype determination using flow cytometry, and a multivariate morphometric approach.

KEY RESULTS

The Italian and French populations formed two well-delimited groupings within the studied complex that were, in contrast to previous taxonomic treatments, clearly distinct from A. montanum. Populations from southeastern France represent A. orophilum, a previously described but abandoned species. Those from central and southern Italy correspond to A. diffusum, exhibiting high, geographically structured variation (central Apennines, Gargano, and southern Apennines/Calabria). This pattern coincides with hotspot refugial regions, in congruence with the "refugia-within-refugia" hypothesis, and is reflected here in the recognition of three subspecies within A. diffusum.

CONCLUSIONS

We provide evidence for the presence of Mediterranean refugia for the studied Alyssum montanum-A. repens complex located in central and southern Italy, which, however, did not contribute to the postglacial colonization of Central Europe. Past extinctions, genetic bottlenecks, and recent expansion were inferred in Central Europe, while long-term accumulation of diversity as well as polyploidization occurred in the Apennines.

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.

Palaeoclimatic events, dispersal and migratory losses along the Afro-European axis as drivers of biogeographic distribution in Sylvia warblers

BACKGROUND

The Old World warbler genus Sylvia has been used extensively as a model system in a variety of ecological, genetic, and morphological studies. The genus is comprised of about 25 species, and 70% of these species have distributions at or near the Mediterranean Sea. This distribution pattern suggests a possible role for the Messinian Salinity Crisis (from 5.96-5.33 Ma) as a driving force in lineage diversification. Other species distributions suggest that Late Miocene to Pliocene Afro-tropical forest dynamics have also been important in the evolution of Sylvia lineages. Using a molecular phylogenetic hypothesis and other methods, we seek to develop a biogeographic hypothesis for Sylvia and to explicitly assess the roles of these climate-driven events.

RESULTS

We present the first strongly supported molecular phylogeny for Sylvia. With one exception, species fall into one of three strongly supported clades: one small clade of species distributed mainly in Africa and Europe, one large clade of species distributed mainly in Africa and Asia, and another large clade with primarily a circum-Mediterranean distribution. Asia is reconstructed as the ancestral area for Sylvia. Long-distance migration is reconstructed as the ancestral character state for the genus, and sedentary behavior subsequently evolved seven times.

CONCLUSION

Molecular clock calibration suggests that Sylvia arose in the early Miocene and diverged into three main clades by 12.6 Ma. Divergence estimates indicate that the Messinian Salinity Crisis had a minor impact on Sylvia. Instead, over-water dispersals, repeated loss of long-distance migration, and palaeo-climatic events in Africa played primary roles in Sylvia divergence and distribution.

Birth of a hotspot of intraspecific genetic diversity: notes from the underground

Hotspots of intraspecific diversity have been observed in most species, often within areas of putative Pleistocene refugia. They have thus mostly been viewed as the outcome of prolonged stability of large populations within the refugia. However, recent evidence has suggested that several other microevolutionary processes could also be involved in their formation. Here, we investigate the contribution of these processes to current range-wide patterns of genetic diversity in the Italian endemic mole Talpa romana, using both nuclear (30 allozyme loci) and mitochondrial markers (cytochrome b sequences). Southern populations of this species showed an allozyme variation that is amongst the highest observed in small mammals (most populations had an expected heterozygosity of 0.10 or above), which was particularly unexpected for a subterranean species. Population genetic, phylogeographic and historical demographic analyses indicated that T. romana populations repeatedly underwent allopatric differentiations followed by secondary admixture within the refugial range in southern Italy. A prolonged demographic stability was reliably inferred from the mitochondrial DNA data only for a population group located north and east of the Calabrian peninsula, showing comparatively lower levels of allozyme variability, and lacking evidence of secondary admixture with other groups. Thus, our results point to the admixture between differentiated lineages as the main cause of the higher levels of diversity of refugial populations. When compared with the Pleistocene evolutionary history recently inferred for species from both the same and other geographic regions, these results suggest the need for a reappraisal of the role of gene exchange in the formation of intraspecific hotspots of genetic diversity.