Population genetic structure and diversity of the Apennine endemic stream frog,Rana italica- insights on the Pleistocene evolutionary history of the Italian peninsular biota

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.

Distribution and conservation status of threatened endemic amphibians within the Aspromonte mountain region, a hotspot of Mediterranean biodiversity

Amphibian biodiversity loss in recent years has exceeded that of all other groups of vertebrates. In this context, biodiversity hotspots represent priority targets for conservation in amphibian populations. However, little information is available on the distribution and conservation status of amphibian species within most biodiversity hotspots. Here, we characterized the distribution and conservation status of four threatened endemic amphibians (Bombina pachypus, Salamandra salamandra gigliolii, Salamandrina terdigitata, and Rana italica) in the Aspromonte Mountain region, a biodiversity hotspot in southern Italy where the conservation status of amphibians is almost unexplored. We conducted an intensive field survey of 507 potential breeding sites spanning over 2,326 km2. We found that all four species were widespread in the study area. We observed 337 species occurrences: 63 for S. s. gigliolii, 29 for S. terdigitata, 84 for B. pachypus, and 161 for R. italica. Species distribution analysis revealed that S. s. gigliolii and R. italica populations had an extended and homogenous distribution. Conversely, S. terdigitata showed a dispersed pattern, with long distances among breeding sites, and B. pachypus an aggregated pattern, associated with the availability of suitable artificial habitats. On the other hand, we reported a decrease in B. pachypus occurrence in its natural habitats, which reflects a negative trend of its populations. Overall, our results provide an encouraging framework for the conservation of amphibian populations in this area, but highlight the low coverage of threatened amphibian populations in protected areas, highlighting the need for a reassessment of conservation policies and spatial conservation planning for the Aspromonte region.

Phylogeography reveals the origin of the two phenological forms of large blue, Phengaris arion (Lepidoptera: Lycaenidae)

Our main goal was to investigate the phylogeography of the butterfly Phengaris arion to reveal the evolutionary origin of its ‘spring’ and ‘summer’ forms. Molecular analyses based on highly variable microsatellites, together with Wolbachia screening, were carried out on 34 populations in Europe. We found three well-defined genetic lineages of different origins: the Apennine, the central and the eastern. The highly distinct Apennine lineage is limited by the Alps and evaluated as an Evolutionary Significant Unit (ESU). Therefore, the taxon name ligurica, described from the Ligurian coast (Italy), should not be applied to denote the ‘summer form’ of the Pannonian region. The central lineage is limited by the Carpathians and the most eastern ranges of the Alps, and lacks major range fluctuations related to glaciations, although there is evidence for extra-Mediterranean refugia in the Carpathian Basin. The eastern clade could have had refugia in central Asia. Our results exclude the potential allopatric origin of the ‘spring’ and ‘summer’ arion, and support the hypothesis that the existence of the two forms could be a result of local adaptation to the distinctive phenology of host plant flowering which is manifested in the genetic differences between them. Wolbachia infection has been ruled out as a driver of sympatric speciation in P. arion.

European common frog (Rana temporaria) recolonized Switzerland from multiple glacial refugia in northern Italy via trans- and circum-Alpine routes

The high mountain ranges of Western Europe had a profound effect on the biotic recolonization of Europe from glacial refugia. The Alps present a particularly interesting case because they form an absolute barrier to dispersal for most taxa, obstructing recolonization from multiple refugia in northern Italy. Here, we investigate the effect of the European Alps on the phylogeographic history of the European common frog Rana temporaria. Based on partial cytochrome b and COXI sequences from Switzerland, we find two mitochondrial lineages roughly north and south of the Alpine ridge, with contact zones between them in eastern and western Switzerland. The northern haplogroup falls within the previously identified Western European haplogroup, while the southern haplogroup is unique to Switzerland. We find that the lineages diverged ~110 kya, at approximately the onset of the last glacial glaciation; this indicates that they are from different glacial refugia. Phylogenetic analyses suggest that the northern and southern haplogroups colonized Switzerland via trans- and circum-Alpine routes from at least two separate refugia in northern Italy. Our results illustrate how a complex recolonization history of the central European Alps can arise from the semi-permeable barrier created by high mountains.

Stability in the Balkans: phylogeography of the endemic Greek stream frog, Rana graeca

We still have little knowledge concerning the phylogeography of amphibians and reptiles from the Balkan Peninsula compared with the other two Mediterranean peninsulas. This raises concerns for endemic taxa from these peninsulas, because it might interfere with further conservation efforts. Here we focus on the endemic Greek stream frog (Rana graeca) and reconstruct its biogeography and evolutionary history. Using four genetic markers (Cytb, 16S, COI and BDNF) in > 350 sequences covering the whole distribution range, we conducted phylogenetic, demographic and ecological niche analyses, which revealed the phylogeography of this species. Surprisingly, this examination of R. graeca reveals a very shallow level of intraspecific genetic variability through the Balkans, with two main, statistically supported lineages having a partly sympatric distribution. The most variable marker was Cytb, which showed 19 haplotypes in 123 analysed sequences in the whole species distribution area. Here presented genetic data, together with the environmental niche projection and demographic analyses suggest that R. graeca was probably affected only marginally by climatic oscillations, with the Hellenides as the most suitable area for the occurrence of the species in different geological periods. This is consistent with the observed genetic diversity, which is mostly related to these mountains. Although the species shows a certain level of phenotypic variability and ecological preferences, this might be related to species plasticity affected by the micro-climatic conditions in small areas, which merits further research. Comparing phylogeography of other amphibian and reptile species in the Balkans, we showed that the observed pattern represents a new view on the phylogeography of the Balkan herpetofauna.

Reconstructing hotspots of genetic diversity from glacial refugia and subsequent dispersal in Italian common toads (Bufo bufo)

Genetic diversity feeds the evolutionary process and allows populations to adapt to environmental changes. However, we still lack a thorough understanding of why hotspots of genetic diversity are so 'hot'. Here, we analysed the relative contribution of bioclimatic stability and genetic admixture between divergent lineages in shaping spatial patterns of genetic diversity in the common toad Bufo bufo along the Italian peninsula. We combined population genetic, phylogeographic and species distribution modelling (SDM) approaches to map ancestral areas, glacial refugia, and secondary contact zones. We consistently identified three phylogeographic lineages, distributed in northern, central and southern Italy. These lineages expanded from their ancestral areas and established secondary contact zones, before the last interglacial. SDM identified widespread glacial refugia in peninsular Italy, sometimes located under the present-day sea-level. Generalized linear models indicated genetic admixture as the only significant predictor of the levels of population genetic diversity. Our results show that glacial refugia contributed to preserving both levels and patterns of genetic diversity across glacial-interglacial cycles, but not to their formation, and highlight a general principle emerging in Mediterranean species: higher levels of genetic diversity mark populations with substantial contributions from multiple genetic lineages, irrespective of the location of glacial refugia.

Mapping the geographic origin of captive and confiscated Hermann's tortoises: A genetic toolkit for conservation and forensic analyses

The illegal trade has been threatening tortoise populations worldwide for decades. Nowadays, however, DNA typing and forensic genetic approaches allow us to investigate the geographic origin of confiscated animals and to relocate them into the wild, providing that suitable molecular tools and reference data are available. Here we assess the suitability of a small panel of microsatellite markers to investigate patterns of illegal translocations and to assist forensic genetic applications in the endangered Mediterranean land tortoise Testudo hermanni hermanni. Specific allelic ladders were created for each locus and tested on several reference samples. We used the microsatellite panel to (i) increase our understanding of the population genetic structure in wild populations with new data from previously unsampled geographic areas (overall 461 wild individuals from 28 sampling sites); (ii) detect the presence of non-native individuals in wild populations; and (iii) identify the most likely geographic area of origin of 458 confiscated individuals hosted in Italian seizure and recovery centers. Our analysis initially identified six major genetic clusters corresponding to different geographic macro-areas along the Mediterranean range. Long-distance migrants among wild populations, due to translocations, were found and removed from the reference database. Assignment tests allowed us to allocate approximately 70 % of confiscated individuals of unknown origin to one of the six Mediterranean macro-areas. Most of the assigned tortoises belonged to the genetic cluster corresponding to the area where the respective captivity center was located. However, we also found evidence of long-distance origins of confiscated individuals, especially in centers along the Adriatic coast and facing the Balkan regions, a well-known source of illegally traded individuals. Our results clearly show that the microsatellite panel and the reference dataset can play a beneficial role in reintroduction and repatriation projects when confiscated individuals need to be re-assigned to their respective macro-area of origin before release, and can assist future forensic genetic applications in detecting the illegal trade and possession of Testudo hermanni individuals.

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.

Longitudinal study of Amphibiocystidium sp. infection in a natural population of the Italian stream frog (Rana italica)

Mesomycetozoean-induced infections (order Dermocystida, genus Amphibiocystidium) in European and North American amphibians are causing alarm. To date, the pathogenicity of these parasites in field conditions has been poorly studied, and demographic consequences on amphibian populations have not been explored. In this study, an Amphibiocystidium sp. infection is reported in a natural population of the Italian stream frog (Rana italica) of Central Italy, over a 7-year period from 2008 to 2014. Light and electron microscope examinations, as well as partial 18S rDNA sequence analysis were used to characterize the parasite. Moreover, a capture-mark-recapture study was conducted to assess the frog demographics in response to infection. Negative effects of amphibiocystidiosis on individual survival and population fitness were absent throughout the sampling period, despite the high estimates of disease prevalence. This might have been due to resistance and/or tolerance strategies developed by the frogs in response to the persistence of Amphibiocystidium infection in this system. We hypothesized that in the examined R. italica population, amphibiocystidiosis is an ongoing endemic/epidemic infection. However, ecological and host-specific factors, interacting in a synergistic fashion, might be responsible for variations in the susceptibility to Amphibiocystidium infection of both conspecific populations and heterospecific individuals of R. italica.

Is the Italian stream frog (Rana italica Dubois, 1987) an opportunistic exploiter of cave twilight zone?

Studies on frogs exploiting subterranean environments are extremely scarce, as these Amphibians are usually considered accidental in these environments. However, according to recent studies, some anurans actively select subterranean environments on the basis of specific environmental features, and thus are able to inhabit these environments throughout the year. We present the first study on the abundance and spatial use of the Italian stream frog, Ranaitalica, in subterranean environments. We monthly collected data from 66 cave sectors during a whole year (2013), recording > 120 detections of R.italica. Frogs were more frequently found close to the cave entrance, without significant differences between age classes or sexes. Adults generally were observed being higher up along cave walls compared to juveniles. Frogs abundance was higher in areas showing specific environmental features, such as warm temperature, low incident light and the presence of potential prey. Ranaitalica likely occupies subterranean areas characterized by a combination of microclimatic suitability and prey availability.

Mitochondrial phylogeography of the Iberian endemic frog Rana iberica, with implications for its conservation

Genetic characterization of species using phylogeographic approaches represents a basic reference to understand their evolutionary history as well as to identify conservation priorities to protect areas of particular interest regarding evolutionary potential. Even in well-studied regions such information is lacking for the majority of species, including many endemic species with reduced distribution ranges. We investigate the phylogeographic pattern of the Iberian frog Rana iberica, an endemic amphibian restricted to Central and North-Western Iberian Peninsula. Using mitochondrial sequences, we reconstruct the phylogeographic history of the species to test the effect of Quaternary climate changes on the evolutionary diversification of lineages, that is, the differentiation of mitochondrial lineages and the formation of genetic diversity melting pots, and integrate phylogeographic evidence for future conservation planning. Our results indicate the existence of 3 main mitochondrial lineages differentiated during the Upper Pleistocene. Both historical demographic analyses and climatic niche modeling show a strong effect of glacial climate changes, suggesting recurrent range contractions and expansions. Under such circumstances, differentiation took place most likely by isolation in allopatric interglacial refugia. Secondary lineage admixture in northern Portugal generated a broad mixed zone with highest nucleotide diversity. Given its particular evolutionary potential, its reduced distribution and eventual threats under current climate change scenario, conservation priorities should focus on the isolated lineage from Sierra de Guadalupe.

Quaternary history, population genetic structure and diversity of the cold-adapted Alpine newt Ichthyosaura alpestris in peninsular Italy

Mediterranean peninsulas are major biodiversity hotspots, and cold-adapted species are an important component of this biodiversity. However, cold-adapted species contributed surprisingly little to our knowledge of the intimate links between Quaternary environmental changes, species' responses to these changes, and current patterns of intraspecific biodiversity. Here, we investigated the genetic structure and evolutionary history of a cold-adapted amphibian, the Alpine newt Ichthyosaura alpestris, within the Italian peninsula. Nuclear and mitochondrial markers consistently identified three distinct genetic lineages, whose divergence dates to the Early Pleistocene (1.9 and 0.8 million years ago). Our results show that the Italian peninsula provided multiple Pleistocene refugia to this cold-adapted species, and suggest that allopatric fragmentation followed by secondary admixture have been key events in the formation of its current pattern of genetic diversity. Indeed, estimates of population genetic diversity clearly identified contact populations as those achieving the highest levels of diversity. Such concordance among cold-adapted and temperate species in terms of processes triggering the formation of regional patterns of genetic diversity provides strong support for the hypothesis that gene exchange between divergent lineages, rather than long-term stability of refugial populations, has been the main step toward the formation of hotspots of intraspecific biodiversity.

Three in One--Multiple Faunal Elements within an Endangered European Butterfly Species

Ice ages within Europe forced many species to retreat to refugia, of which three major biogeographic basic types can be distinguished: "Mediterranean", "Continental" and "Alpine / Arctic" species. However, this classification often fails to explain the complex phylogeography of European species with a wide range of latitudinal and altitudinal distribution. Hence, we tested for the possibility that all three mentioned faunal elements are represented within one species. Our data was obtained by scoring 1,307 Euphydryas aurinia individuals (46 European locations) for 17 allozyme loci, and sequencing a subset of 492 individuals (21 sites) for a 626 base pairs COI fragment. Genetic diversity indices, F statistics, hierarchical analyses of molecular variance, individual-based clustering, and networks were used to explore the phylogeographic patterns. The COI fragment represented 18 haplotypes showing a strong geographic structure. All but one allozyme loci analysed were polymorphic with a mean F_ST of 0.20, supporting a pronounced among population structure. Interpretation of both genetic marker systems, using several analytical tools, calls for the recognition of twelve genetic groups. These analyses consistently distinguished different groups in Iberia (2), Italy, Provence, Alps (3), Slovenia, Carpathian Basin, the lowlands of West and Central Europe as well as Estonia, often with considerable additional substructures. The genetic data strongly support the hypothesis that E. aurinia survived the last glaciation in Mediterranean, extra-Mediterranean and perialpine refugia. It is thus a rare example of a model organism that combines attributes of faunal elements from all three of these sources. The observed differences between allozymes and mtDNA most likely result from recent introgression of mtDNA into nuclear allozyme groups. Our results indicate discrepancies with the morphologically-based subspecies models, underlining the need to revise the current taxonomy.

Amphibian molecular ecology and how it has informed conservation

Molecular ecology has become one of the key tools in the modern conservationist's kit. Here we review three areas where molecular ecology has been applied to amphibian conservation: genes on landscapes, within-population processes, and genes that matter. We summarize relevant analytical methods, recent important studies from the amphibian literature, and conservation implications for each section. Finally, we include five in-depth examples of how molecular ecology has been successfully applied to specific amphibian systems.

Characterization of the Skin Microbiota in Italian Stream Frogs (Rana italica) Infected and Uninfected by a Cutaneous Parasitic Disease

In human and wildlife populations, the natural microbiota plays an important role in health maintenance and the prevention of emerging infectious diseases. In amphibians, infectious diseases have been closely associated with population decline and extinction worldwide. Skin symbiont communities have been suggested as one of the factors driving the different susceptibilities of amphibians to diseases. The activity of the skin microbiota of amphibians against fungal pathogens, such as Batrachochytrium dendrobatidis, has been examined extensively, whereas its protective role towards the cutaneous infectious diseases caused by Amphibiocystidium parasites has not yet been elucidated in detail. In the present study, we investigated, for the first time, the cutaneous microbiota of the Italian stream frog (Rana italica) and characterized the microbial assemblages of frogs uninfected and infected by Amphibiocystidium using the Illumina next-generation sequencing of 16S rRNA gene fragments. A total of 629 different OTUs belonging to 16 different phyla were detected. Bacterial populations shared by all individuals represented only one fifth of all OTUs and were dominated by a small number of OTUs. Statistical analyses based on Bray-Curtis distances showed that uninfected and infected specimens had distinct cutaneous bacterial community structures. Phylotypes belonging to the genera Janthinobacterium, Pseudomonas, and Flavobacterium were more abundant, and sometimes almost exclusively present, in uninfected than in infected specimens. These bacterial populations, known to exhibit antifungal activity in amphibians, may also play a role in protection against cutaneous infectious diseases caused by Amphibiocystidium parasites.

Persistence and dispersal in a Southern Hemisphere glaciated landscape: the phylogeography of the spotted snow skink (Niveoscincus ocellatus) in Tasmania

Background

The aim of this research was to identify the effects of Pleistocene climate change on the distribution of fauna in Tasmania, and contrast this with biotic responses in other temperate regions in the Northern and Southern Hemisphere that experienced glacial activity during this epoch. This was achieved by examining the phylogeographic patterns in a widely distributed Tasmanian endemic reptile, Niveoscincus ocellatus. 204 individuals from 29 populations across the distributional range of N. ocellatus were surveyed for variation at two mitochondrial genes (ND2, ND4), and two nuclear genes (β-globin, RPS8). Phylogenetic relationships were reconstructed using a range of methods (maximum parsimony, Bayesian inference and haplotype networks), and the demographic histories of populations were assessed (AMOVA, Tajima’s D, Fu’s Fs, mismatch distributions, extended Bayesian skyline plots, and relaxed random walk analyses).

Results

There was a high degree of mitochondrial haplotype diversity (96 unique haplotypes) and phylogeographic structure, where spatially distinct groups were associated with Tasmania’s Northeast and a large area covering Southeast and Central Tasmania. Phylogeographic structure was also present within each major group, but the degree varied regionally, being highest in the Northeast. Only the Southeastern group had a signature of demographic expansion, occurring during the Pleistocene but post-dating the Last Glacial Maximum. In contrast, nuclear DNA had low levels of variation and a lack of phylogeographic structure, and further loci should be surveyed to corroborate the mitochondrial inferences.

Conclusions

The phylogeographic patterns of N. ocellatus indicate Pleistocene range and demographic expansion in N. ocellatus, particularly in the Southeast and Central areas of Tasmania. Expansion in Central and Southeastern areas appears to have been more recent in both demographic and spatial contexts, than in Northeast Tasmania, which is consistent with inferences for other taxa of greater stability and persistence in Northeast Tasmania during the Last Glacial Maximum. These phylogeographic patterns indicate contrasting demographic histories of populations in close proximity to areas directly affected by glaciers in the Southern Hemisphere during the LGM.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0397-y) contains supplementary material, which is available to authorized users.

Influence of ecological and geological features on rangewide patterns of genetic structure in a widespread passerine

Geological and ecological features restrict dispersal and gene flow, leading to isolated populations. Dispersal barriers can be obvious physical structures in the landscape; however microgeographic differences can also lead to genetic isolation. Our study examined dispersal barriers at both macro- and micro-geographical scales in the black-capped chickadee, a resident North American songbird. Although birds have high dispersal potential, evidence suggests dispersal is restricted by barriers. The chickadee's range encompasses a number of physiological features which may impede movement and lead to divergence. Analyses of 913 individuals from 34 sampling sites across the entire range using 11 microsatellite loci revealed as many as 13 genetic clusters. Populations in the east were largely panmictic whereas populations in the western portion of the range showed significant genetic structure, which often coincided with large mountain ranges, such as the Cascade and Rocky Mountains, as well as areas of unsuitable habitat. Unlike populations in the central and southern Rockies, populations on either side of the northern Rockies were not genetically distinct. Furthermore, Northeast Oregon represents a forested island within the Great Basin; genetically isolated from all other populations. Substructuring at the microgeographical scale was also evident within the Fraser Plateau of central British Columbia, and in the southeast Rockies where no obvious physical barriers are present, suggesting additional factors may be impeding dispersal and gene flow. Dispersal barriers are therefore not restricted to large physical structures, although mountain ranges and large water bodies do play a large role in structuring populations in this study.

What triggers the rising of an intraspecific biodiversity hotspot? Hints from the agile frog

Hotspots of genetic diversity are regions of utmost importance for species survival and conservation, and their intimate link with the geographic location of glacial refugia has been well established. Nonetheless, the microevolutionary processes underlying the generation of hotspots in such regions have only recently become a fervent field of research. We investigated the phylogeographic and population genetic structure of the agile frog, Rana dalmatina, within its putative refugium in peninsular Italy. We found this region to harbour far more diversity, phylogeographic structure, and lineages of ancient origin than that by the rest of the species' range in Europe. This pattern appeared to be well explained by climate-driven microevolutionary processes that occurred during both glacial and interglacial epochs. Therefore, the inferred evolutionary history of R. dalmatina in Italy supports a view of glacial refugia as 'factories' rather than as repositories of genetic diversity, with significant implications for conservation strategies for hotspots.

Paleoclimatic modeling and phylogeography of least killifish, Heterandria formosa: insights into Pleistocene expansion-contraction dynamics and evolutionary history of North American Coastal Plain freshwater biota

BACKGROUND

Climatic and sea-level fluctuations throughout the last Pleistocene glacial cycle (~130-0 ka) profoundly influenced present-day distributions and genetic diversity of Northern Hemisphere biotas by forcing range contractions in many species during the glacial advance and allowing expansion following glacial retreat ('expansion-contraction' model). Evidence for such range dynamics and refugia in the unglaciated Gulf-Atlantic Coastal Plain stems largely from terrestrial species, and aquatic species Pleistocene responses remain relatively uninvestigated. Heterandria formosa, a wide-ranging regional endemic, presents an ideal system to test the expansion-contraction model within this biota. By integrating ecological niche modeling and phylogeography, we infer the Pleistocene history of this livebearing fish (Poeciliidae) and test for several predicted distributional and genetic effects of the last glaciation.

RESULTS

Paleoclimatic models predicted range contraction to a single southwest Florida peninsula refugium during the Last Glacial Maximum, followed by northward expansion. We inferred spatial-population subdivision into four groups that reflect genetic barriers outside this refuge. Several other features of the genetic data were consistent with predictions derived from an expansion-contraction model: limited intraspecific divergence (e.g. mean mtDNA p-distance = 0.66%); a pattern of mtDNA diversity (mean Hd = 0.934; mean π = 0.007) consistent with rapid, recent population expansion; a lack of mtDNA isolation-by-distance; and clinal variation in allozyme diversity with higher diversity at lower latitudes near the predicted refugium. Statistical tests of mismatch distributions and coalescent simulations of the gene tree lent greater support to a scenario of post-glacial expansion and diversification from a single refugium than to any other model examined (e.g. multiple-refugia scenarios).

CONCLUSIONS

Congruent results from diverse data indicate H. formosa fits the classic Pleistocene expansion-contraction model, even as the genetic data suggest additional ecological influences on population structure. While evidence for Plio-Pleistocene Gulf Coast vicariance is well described for many freshwater species presently codistributed with H. formosa, this species demography and diversification departs notably from this pattern. Species-specific expansion-contraction dynamics may therefore have figured more prominently in shaping Coastal Plain evolutionary history than previously thought. Our findings bolster growing appreciation for the complexity of phylogeographical structuring within North America's southern refugia, including responses of Coastal Plain freshwater biota to Pleistocene climatic fluctuations.

Widespread occurrence of Batrachochytrium dendrobatidis in contemporary and historical samples of the endangered Bombina pachypus along the Italian peninsula

Batrachochytrium dendrobatidis is considered a main driver of the worldwide declines and extinctions of amphibian populations. Nonetheless, fundamental questions about its epidemiology, including whether it acts mainly as a “lone killer” or in conjunction with other factors, remain largely open. In this paper we analysed contemporary and historical samples of the endangered Apennine yellow-bellied toad (Bombina pachypus) along the Italian peninsula, in order to assess the presence of the pathogen and its spreading dynamics. Once common throughout its range, B. pachypus started to decline after the mid-1990s in the northern and central regions, whereas no declines have been observed so far in the southern region. We show that Batrachochytrium dendrobatidis is currently widespread along the entire peninsula, and that this was already so at least as early as the late 1970s, that is, well before the beginning of the observed declines. This temporal mismatch between pathogen occurrence and host decline, as well as the spatial pattern of the declines, suggests that the pathogen has not acted as a “lone killer”, but in conjunction with other factors. Among the potentially interacting factors, we identified two as the most probable, genetic diversity of host populations and recent climate changes. We discuss the plausibility of this scenario and its implications on the conservation of B. pachypus populations.

Radically different phylogeographies and patterns of genetic variation in two European brown frogs, genus Rana

We reconstruct range-wide phylogeographies of two widespread and largely co-occurring Western Palearctic frogs, Rana temporaria and R. dalmatina. Based on tissue or saliva samples of over 1000 individuals, we compare a variety of genetic marker systems, including mitochondrial DNA, single-copy protein-coding nuclear genes, microsatellite loci, and single nucleotide polymorphisms (SNPs) of transcriptomes of both species. The two focal species differ radically in their phylogeographic structure, with R. temporaria being strongly variable among and within populations, and R. dalmatina homogeneous across Europe with a single strongly differentiated population in southern Italy. These differences were observed across the various markers studied, including microsatellites and SNP density, but especially in protein-coding nuclear genes where R. dalmatina had extremely low heterozygosity values across its range, including potential refugial areas. On the contrary, R. temporaria had comparably high range-wide values, including many areas of probable postglacial colonization. A phylogeny of R. temporaria based on various concatenated mtDNA genes revealed that two haplotype clades endemic to Iberia form a paraphyletic group at the base of the cladogram, and all other haplotypes form a monophyletic group, in agreement with an Iberian origin of the species. Demographic analysis suggests that R. temporaria and R. dalmatina have genealogies of roughly the same time to coalescence (TMRCA ~3.5 mya for both species), but R. temporaria might have been characterized by larger ancestral and current effective population sizes than R. dalmatina. The high genetic variation in R. temporaria can therefore be explained by its early range expansion out of Iberia, with subsequent cycles of differentiation in cryptic glacial refugial areas followed by admixture, while the range expansion of R. dalmatina into central Europe is a probably more recent event.

Refugia within refugia as a key to disentangle the genetic pattern of a highly variable species: the case of Rana temporaria Linnaeus, 1758 (Anura, Ranidae)

Two distinct lineages of Rana temporaria are known in the Palaearctic region, but it is uncertain whether this species persisted in one or more Pleistocene refugia. We resolved the phylogeographic history and genetic variability of R. temporaria in the Italian peninsula, a 'traditional' Pleistocene refugium, and related our findings to patterns described for other European populations. We sequenced the mitochondrial markers Cox I and cytochrome b. Phylogenetic reconstruction only indicated the presence of haplotypes belonging to the Western lineage in the Italian peninsula. Overall, the genetic variability of Italian populations was higher than other European populations, which shared haplotypes with the Alpine populations. We demonstrated subdivision into five main Italian sublineages, which was associated with a geographical structure of populations in two divergent groups. In particular, one Apennine group might have resulted from bottlenecks during the last interglacials ages. In contrast, Alpine populations were recently diverged and showed incomplete lineage sorting. Our data indicate that the Italian peninsula served as refugium for the Western lineage of R. temporaria. Dispersion towards Central Europe probably started only from the western slope of the Alps via a rapid leading edge expansion. The identified structure is partially congruent with traditional peripheral refugia identified for plants. This evolutionary scenario does not support any taxonomic distinction at the subspecific level for R. temporaria.

One species, three Pleistocene evolutionary histories: phylogeography of the Italian crested newt, Triturus carnifex

Phylogeographic patterns of temperate species from the Mediterranean peninsulas have been investigated intensively. Nevertheless, as more phylogeographies become available, either unique patterns or new lines of concordance continue to emerge, providing new insights on the evolution of regional biotas. Here, we investigated the phylogeography and evolutionary history of the Italian crested newt, Triturus carnifex, through phylogenetic, molecular dating and population structure analyses of two mitochondrial gene fragments (ND2 and ND4; overall 1273 bp). We found three main mtDNA lineages having parapatric distribution and estimated divergence times between Late Pliocene and Early Pleistocene. One lineage (S) was widespread south of the northern Apennine chain and was further geographically structured into five sublineages, likely of Middle Pleistocene origin. The second lineage (C) was widespread throughout the Padano-Venetian plain and did not show a clear phylogeographic structure. The third lineage (N) was observed in only two populations located on western Croatia/Slovenia. Results of analysis of molecular variance suggested that partitioning populations according to the geographic distribution of these lineages and sublineages explains 76% of the observed genetic variation. The phylogeographic structure observed within T. carnifex and divergence time estimates among its lineages, suggest that responses to Pleistocene environmental changes in this single species have been as diverse as those found previously among several codistributed temperate species combined. Consistent with the landscape heterogeneity, physiographic features, and palaeogeographical evolution of its distribution range, these responses encompass multiple refugia along the Apennine chain, lowland refugia in large peri-coastal plains, and a 'cryptic' northern refugium.

Phylogeography and species biogeography of montane Great Basin stoneflies

Sky islands are ideal systems for determining the effects of climatic oscillations on species distributions and genetic structure. Our study focused on montane stonefly populations in the Great Basin of western North America. We used niche-based distribution modelling, phylogeography and traditional species-based biogeography to test several hypotheses as follows: (i) genetic differentiation among Doroneuria baumanni populations will be independent of hydrologic connectivity (headwater model); (ii) Sky islands were colonized when habitat was more continuous and populations likely experienced multiple expansions and contractions; (iii) Colonization events were coincident with the late Pleistocene and Holocene; and (iv) Shared topography and climate history will result in concordant patterns of genetic differentiation in D. baumanni and occurrences of 32 stonefly species across the region. Overall, Φ(ST) 's and coalescent-based estimates of migration were consistent with the headwater model. Maximum likelihood and Bayesian gene trees identified three major nonoverlapping east-west clades. Distribution modelling indicated more suitable habitat in the Great Basin during the Last Glacial Maximum than at present, but none during the last interglacial period. Demographic analyses showed evidence of population expansion in one of the three major east-west clades. Intra-clade divergence times (60,000-183,000ybp) were well within the late Pleistocene while among-clade divergence times (499.000-719,000ybp) were deeper. Genetic differentiation in D. baumanni and distributions of stonefly species were significantly concordant. These results imply that climatic oscillations have played major roles in shaping the genetic structure and distributions of Great Basin stoneflies, but that divergence among clades occurred much earlier than our late Pleistocence/early Holocene predictions.

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.

Gone with the plate: the opening of the Western Mediterranean basin drove the diversification of ground-dweller spiders

BACKGROUND

The major islands of the Western Mediterranean--Corsica, Sardinia, and the Balearic Islands--are continental terrenes that drifted towards their present day location following a retreat from their original position on the eastern Iberian Peninsula about 30 million years ago. Several studies have taken advantage of this well-dated geological scenario to calibrate molecular rates in species for which distributions seemed to match this tectonic event. Nevertheless, the use of external calibration points has revealed that most of the present-day fauna on these islands post-dated the opening of the western Mediterranean basin. In this study, we use sequence information of the cox1, nad1, 16S, L1, and 12S mitochondrial genes and the 18S, 28S, and h3 nuclear genes, along with relaxed clock models and a combination of biogeographic and fossil external calibration points, to test alternative historical scenarios of the evolutionary history of the ground-dweller spider genus Parachtes (Dysderidae), which is endemic to the region.

RESULTS

We analyse 49 specimens representing populations of most Parachtes species and close relatives. Our results reveal that both the sequence of species formation in Parachtes and the estimated divergence times match the geochronological sequence of separation of the main islands, suggesting that the diversification of the group was driven by Tertiary plate tectonics. In addition, the confirmation that Parachtes diversification matches well-dated geological events provides a model framework to infer substitution rates of molecular markers. Divergence rates estimates ranged from 3.5% My(-1) (nad1) to 0.12% My(-1) (28S), and the average divergence rate for the mitochondrial genes was 2.25% My(-1), very close to the "standard" arthropod mitochondrial rate (2.3% My(-1)).

CONCLUSIONS

Our study provides the first unequivocal evidence of terrestrial endemic fauna of the major western Mediterranean islands, whose origin can be traced back to the Oligocene separation of these islands from the continent. Moreover, our study provides useful information on the divergence rate estimates of the most commonly used genes for phylogenetic inference in non-model arthropods.

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.

Population dynamic of the extinct European aurochs: genetic evidence of a north-south differentiation pattern and no evidence of post-glacial expansion

Background

The aurochs (Bos primigenius) was a large bovine that ranged over almost the entirety of the Eurasian continent and North Africa. It is the wild ancestor of the modern cattle (Bos taurus), and went extinct in 1627 probably as a consequence of human hunting and the progressive reduction of its habitat. To investigate in detail the genetic history of this species and to compare the population dynamics in different European areas, we analysed Bos primigenius remains from various sites across Italy.

Results

Fourteen samples provided ancient DNA fragments from the mitochondrial hypervariable region. Our data, jointly analysed with previously published sequences, support the view that Italian aurochsen were genetically similar to modern bovine breeds, but very different from northern/central European aurochsen. Bayesian analyses and coalescent simulations indicate that the genetic variation pattern in both Italian and northern/central European aurochsen is compatible with demographic stability after the last glaciation. We provide evidence that signatures of population expansion can erroneously arise in stable aurochsen populations when the different ages of the samples are not taken into account.

Conclusions

Distinct groups of aurochsen probably inhabited Italy and northern/central Europe after the last glaciation, respectively. On the contrary, Italian and Fertile Crescent aurochsen likely shared several mtDNA sequences, now common in modern breeds. We argue that a certain level of genetic homogeneity characterized aurochs populations in Southern Europe and the Middle East, and also that post-glacial recolonization of northern and central Europe advanced, without major demographic expansions, from eastern, and not southern, refugia.