Phylogeography and historical demography of the Italian treefrog, Hyla intermedia, reveals multiple refugia, population expansions and secondary contacts within peninsular Italy

The introduction of three cryptic tree frog species in the Dutch coastal dunes challenges conservation paradigms

Invasive species are considered one of the main drivers of the sixth mass extinction. Conservation solutions depend on whether a species is also indigenous to the country it invades (i.e., beyond its native range). In the case of invasive cryptic species, genetic tools are required to establish their identity. We illustrate these issues with the human-mediated colonization of the Dutch coastal dunes by Hyla tree frogs. Although previously assumed to concern the indigenous common tree frog H. arborea, European tree frogs comprise a complex of allopatric cryptic species, meaning the taxonomic identity of introduced Dutch populations warrants investigation. We sequence mtDNA for 173 individuals from native and introduced populations across the Netherlands and compare our dataset with hundreds of Hyla haplotypes previously barcoded in the Western Palearctic. Two of the dune populations carry an mtDNA haplotype of the native species H. arborea that occurs naturally elsewhere in the Netherlands. In contrast, mtDNA assigned to the eastern tree frog H. orientalis was detected in all three other dune populations. In one of these populations mtDNA of the Italian tree frog H. intermedia was also found. Not one, but three species of tree frogs have thus been introduced to the Dutch coastal dunes, only one of them being native to the Netherlands. This situation causes a conservation conundrum as some introduced populations are lawfully protected but could pose a threat to local biodiversity. Regarding the ‘true’ exotic tree frog species, all conservation options should be considered.

Biogeographic and demographic history of the Mediterranean snakes Malpolon monspessulanus and Hemorrhois hippocrepis across the Strait of Gibraltar

BACKGROUND

The contribution of North Africa to the assembly of biodiversity within the Western Palaearctic is still poorly documented. Since the Miocene, multiple biotic exchanges occurred across the Strait of Gibraltar, underlying the high biogeographic affinity between the western European and African sides of the Mediterranean basin. We investigated the biogeographic and demographic dynamics of two large Mediterranean-adapted snakes across the Strait and assess their relevance to the origin and diversity patterns of current European and North African populations.

RESULTS

We inferred phylogeographic patterns and demographic history of M. monspessulanus and H. hippocrepis, based on range-wide multilocus data, combined with fossil data and species distribution modelling, under present and past bioclimatic envelopes. For both species we identified endemic lineages in the High Atlas Mountains (Morocco) and in eastern Iberia, suggesting their persistence in Europe during the Pleistocene. One lineage is shared between North Africa and southern Iberia and likely spread from the former to the latter during the sea-level low stand of the last glacial stage. During this period M. monspessulanus shows a sudden demographic expansion, associated with increased habitat suitability in North Africa. Lower habitat suitability is predicted for both species during interglacial stages, with suitable areas restricted to coastal and mountain ranges of Iberia and Morocco. Compiled fossil data for M. monspessulanus show a continuous fossil record in Iberia at least since the Pliocene and throughout the Pleistocene.

CONCLUSIONS

The previously proposed hypothesis of Pleistocene glacial extinction of both species in Europe is not supported based on genetic data, bioclimatic envelopes models, and the available fossil record. A model of range retraction to mountain refugia during arid periods and of glacial expansion (demographic and spatial) associated to an increase of Mediterranean habitats during glacial epochs emerges as a general pattern for mesic vertebrates in North Africa. Moreover, the phylogeographic pattern of H. hippocrepis conforms to a well-established biogeographic partition between western and eastern Maghreb.

Phylogeography and ecological niche modeling reveal evolutionary history of Leiolepis ocellata (Squamata, Leiolepidae)

Leiolepis ocellata is a lizard species distributing in topographically diverse habitats in northern Thailand. To explore its evolutionary history, 113 samples of L. ocellata were collected from 11 localities covering its distributional range in northern Thailand, and sequenced for mtDNA fragments (Cyt b and ND2). Pairwise comparisons across sampling localities yielded significant genetic differentiation (F ST and Jost's D) but no clear pattern of isolation by distance could be demonstrated based on the Mantel test. Phylogenetic and network analyses highlighted six haplogroups. Their divergence times were estimated to occur during the Pleistocene, much more recent than major orogenic events affecting northern Thailand. Instead, the results suggested that lineage divergences, of particularly eastern and western haplogroups of the region, coincided with the major rivers in the region (Yom river and Ping river, respectively), indicating vicariance in response to riverine barriers. Furthermore, ecological niche modeling suggested an expansion of suitable habitats of L. ocellata, when LGM-liked conditions. This expansion potentially facilitated their dispersal among adjacent localities leading to lineage diversification and genetic admixture, after the riverine divergence.

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.

The effect of phylogeographic history on species boundaries: a comparative framework in Hyla tree frogs

Because it is indicative of reproductive isolation, the amount of genetic introgression across secondary contact zones is increasingly considered in species delimitation. However, patterns of admixture at range margins can be skewed by the regional dynamics of hybrid zones. In this context, we posit an important role for phylogeographic history: hybrid zones located within glacial refugia (putatively formed during the Late-Pleistocene) should be better defined than those located in post-glacial or introduced ranges (putatively formed during the Holocene and the Anthropocene). We test this hypothesis in a speciation continuum of tree frogs from the Western Palearctic (Hyla), featuring ten identified contacts between species spanning Plio-Pleistocene to Miocene divergences. We review the rich phylogeographic literature of this group and examine the overlooked transition between H. arborea and H. molleri in Western France using a multilocus dataset. Our comparative analysis supports a trend that contacts zones resulting from post-glacial expansions and human translocations feature more extensive introgression than those established within refugial areas. Integrating the biogeographic history of incipient species, i.e. their age since first contact together with their genetic divergence, thus appears timely to draw sound evolutionary and taxonomic inferences from patterns of introgression across hybrid zones.

Identification of novel haplotypes and interpretation of gene flow of mitochondrial DNA control region of Eurasian otter (Lutra lutra) for the effective conservation

The number and distribution of Eurasian otters have declined during twentieth century due to human activity and water pollution. The global conservation status of Eurasian otter is presently ‘Near Threatened (NT)’ and strictly protected by being listed on the international legislation and conventions. A number of studies using the mitochondrial DNA (mtDNA) control region (CR) have been conducted in order to effectively apply conservation and reintroduction programs, especially in Europe. However, aside from Europe, there have been few studies concerning genetic diversity and phylogeny of Eurasian otters. Therefore, in this study, we sequenced partial mtDNA CR sequences (232 bp) from five South Korean Eurasian otters and analyzed 27 otters originating from parts of northeast Asia (South Korea, China, Japan and Russia (Sakhalin)), and Europe. Out of 232 bp partial mtDNA CR sequences, 13 polymorphic sites (5.6%) were identified and 4 novel mtDNA CR haplotypes (Lut16–19) were discovered from 12 Eurasian otters originating from northeast Asian region. In this study, a comprehensive analysis of genetic diversity and population structure of Eurasian otter between Europe and northeast Asia continents were conducted. Of these, different past demographic histories in Pleistocene period might have largely impacted the genetic structure of each population differently. In addition, low degree of gene flow, isolation by distance (IBD) pattern from geographically wide distanced dataset and analysis of molecular variance (AMOVA) also represented distinct genetic characteristics of Eurasian otter between Europe and northeast Asia.

Chasing the phantom: biogeography and conservation of Vipera latastei-monticola in the Maghreb (North Africa)

The Maghreb region (North Africa) constitutes a major component of the Mediterranean Basin biodiversity hotspot. During the last centuries, a consistent human population growth has led to an unprecedented rate of habitat transformation and loss in the region and thus, threatening its biodiversity. The Western Mediterranean viper Vipera latastei-monticola inhabits humid and subhumid areas in the main mountain ranges of the Maghreb, facing such threatening factors; however, its elusive character and rarity hindered data collection for distinct biological purposes. Here, we study the biogeographical patterns and conservation status of the Maghrebian V. latastei-monticola resulting from recent sampling campaigns in Morocco and Tunisia. We update species distribution, and integrate phylogeographic and ecological niche modelling analyses at both species and lineage level to identify suitable areas, and to evaluate the impact of anthropogenic transformation and level of protection of their suitable space. We identified four highly divergent mitochondrial lineages, including a new lineage endemic to the Western High Atlas, with allopatric distributions and restricted to mountain ranges, supporting the role of mountains as past climatic refugia. Despite the remoteness of suitable areas, we report widespread habitat degradation and identify the low effectiveness of the current protected areas system in preserving the species and lineages range. Our study shows the urgent need to apply management actions for the long-term conservation of this vulnerable species and suggests a revaluation of the specific status of V. monticola, as these populations likely represent an ecotype of V. latastei.

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.

Evolutionary melting pots: a biodiversity hotspot shaped by ring diversifications around the Black Sea in the Eastern tree frog (Hyla orientalis)

Hotspots of intraspecific genetic diversity, which are of primary importance for the conservation of species, have been associated with glacial refugia, that is areas where species survived the Quaternary climatic oscillations. However, the proximate mechanisms generating these hotspots remain an open issue. Hotspots may reflect the long-term persistence of large refugial populations; alternatively, they may result from allopatric differentiation between small and isolated populations, that later admixed. Here, we test these two scenarios in a widely distributed species of tree frog, Hyla orientalis, which inhabits Asia Minor and southeastern Europe. We apply a fine-scale phylogeographic survey, combining fast-evolving mitochondrial and nuclear markers, with a dense sampling throughout the range, as well as ecological niche modelling, to understand what shaped the genetic variation of this species. We documented an important diversity centre around the Black Sea, composed of multiple allopatric and/or parapatric diversifications, likely driven by a combination of Pleistocene climatic fluctuations and complex regional topography. Remarkably, this diversification forms a ring around the Black Sea, from the Caucasus through Anatolia and eastern Europe, with terminal forms coming into contact and partially admixing in Crimea. Our results support the view that glacial refugia generate rather than host genetic diversity and can also function as evolutionary melting pots of biodiversity. Moreover, we report a new case of ring diversification, triggered by a large, yet cohesive dispersal barrier, a very rare situation in nature. Finally, we emphasize the Black Sea region as an important centre of intraspecific diversity in the Palearctic with implications for conservation.

Contrasting evolutionary histories of the legless lizards slow worms (Anguis) shaped by the topography of the Balkan Peninsula

Background

Genetic architecture of a species is a result of historical changes in population size and extent of distribution related to climatic and environmental factors and contemporary processes of dispersal and gene flow. Population-size and range contractions, expansions and shifts have a substantial effect on genetic diversity and intraspecific divergence, which is further shaped by gene-flow limiting barriers. The Balkans, as one of the most important sources of European biodiversity, is a region where many temperate species persisted during the Pleistocene glaciations and where high topographic heterogeneity offers suitable conditions for local adaptations of populations. In this study, we investigated the phylogeographical patterns and demographic histories of four species of semifossorial slow-worm lizards (genus Anguis) present in the Balkan Peninsula, and tested the relationship between genetic diversity and topographic heterogeneity of the inhabited ranges.

Results

We inferred phylogenetic relationships, compared genetic structure and historical demography of slow worms using nucleotide sequence variation of mitochondrial DNA. Four Anguis species with mostly parapatric distributions occur in the Balkan Peninsula. They show different levels of genetic diversity. A signature of population growth was detected in all four species but with various courses in particular populations. We found a strong correlation between genetic diversity of slow-worm populations and topographic ruggedness of the ranges (mountain systems) they inhabit. Areas with more rugged terrain harbour higher genetic diversity.

Conclusions

Phylogeographical pattern of the genus Anguis in the Balkans is concordant with the refugia-within-refugia model previously proposed for both several other taxa in the region and other main European Peninsulas. While slow-worm populations from the southern refugia mostly have restricted distributions and have not dispersed much from their refugial areas, populations from the extra-Mediterranean refugia in northern parts of the Balkans have colonized vast areas of eastern, central, and western Europe. Besides climatic historical events, the heterogeneous topography of the Balkans has also played an important role in shaping genetic diversity of slow worms.

Electronic supplementary material

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

Fine Dissection of Human Mitochondrial DNA Haplogroup HV Lineages Reveals Paleolithic Signatures from European Glacial Refugia

Genetic signatures from the Paleolithic inhabitants of Eurasia can be traced from the early divergent mitochondrial DNA lineages still present in contemporary human populations. Previous studies already suggested a pre-Neolithic diffusion of mitochondrial haplogroup HV*(xH,V) lineages, a relatively rare class of mtDNA types that includes parallel branches mainly distributed across Europe and West Asia with a certain degree of structure. Up till now, variation within haplogroup HV was addressed mainly by analyzing sequence data from the mtDNA control region, except for specific sub-branches, such as HV4 or the widely distributed haplogroups H and V. In this study, we present a revised HV topology based on full mtDNA genome data, and we include a comprehensive dataset consisting of 316 complete mtDNA sequences including 60 new samples from the Italian peninsula, a previously underrepresented geographic area. We highlight points of instability in the particular topology of this haplogroup, reconstructed with BEAST-generated trees and networks. We also confirm a major lineage expansion that probably followed the Late Glacial Maximum and preceded Neolithic population movements. We finally observe that Italy harbors a reservoir of mtDNA diversity, with deep-rooting HV lineages often related to sequences present in the Caucasus and the Middle East. The resulting hypothesis of a glacial refugium in Southern Italy has implications for the understanding of late Paleolithic population movements and is discussed within the archaeological cultural shifts occurred over the entire continent.

Timeframe of speciation inferred from secondary contact zones in the European tree frog radiation (Hyla arborea group)

Background

Hybridization between incipient species is expected to become progressively limited as their genetic divergence increases and reproductive isolation proceeds. Amphibian radiations and their secondary contact zones are useful models to infer the timeframes of speciation, but empirical data from natural systems remains extremely scarce. Here we follow this approach in the European radiation of tree frogs (Hyla arborea group). We investigated a natural hybrid zone between two lineages (Hyla arborea and Hyla orientalis) of Mio-Pliocene divergence (~5 My) for comparison with other hybrid systems from this group.

Results

We found concordant geographic distributions of nuclear and mitochondrial gene pools, and replicated narrow transitions (~30 km) across two independent transects, indicating an advanced state of reproductive isolation and potential local barriers to dispersal. This result parallels the situation between H. arborea and H. intermedia, which share the same amount of divergence with H. orientalis. In contrast, younger lineages show much stronger admixture at secondary contacts.

Conclusions

Our findings corroborate the negative relationship between hybridizability and divergence time in European tree frogs, where 5 My are necessary to achieve almost complete reproductive isolation. Speciation seems to progress homogeneously in this radiation, and might thus be driven by gradual genome-wide changes rather than single speciation genes. However, the timescale differs greatly from that of other well-studied amphibians. General assumptions on the time necessary for speciation based on evidence from unrelated taxa may thus be unreliable. In contrast, comparative hybrid zone analyses within single radiations such as our case study are useful to appreciate the advance of speciation in space and time.

Electronic supplementary material

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

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.

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.

Cryptic diversity among Western Palearctic tree frogs: postglacial range expansion, range limits, and secondary contacts of three European tree frog lineages (Hyla arborea group)

We characterize divergence times, intraspecific diversity and distributions for recently recognized lineages within the Hyla arborea species group, based on mitochondrial and nuclear sequences from 160 localities spanning its whole distribution. Lineages of H. arborea, H. orientalis, H. molleri have at least Pliocene age, supporting species level divergence. The genetically uniform Iberian H. molleri, although largely isolated by the Pyrenees, is parapatric to H. arborea, with evidence for successful hybridization in a small Aquitanian corridor (southwestern France), where the distribution also overlaps with H. meridionalis. The genetically uniform H. arborea, spread from Crete to Brittany, exhibits molecular signatures of a postglacial range expansion. It meets different mtDNA clades of H. orientalis in NE-Greece, along the Carpathians, and in Poland along the Vistula River (there including hybridization). The East-European H. orientalis is strongly structured genetically. Five geographic mitochondrial clades are recognized, with a molecular signature of postglacial range expansions for the clade that reached the most northern latitudes. Hybridization with H. savignyi is suggested in southwestern Turkey. Thus, cryptic diversity in these Pliocene Hyla lineages covers three extremes: a genetically poor, quasi-Iberian endemic (H. molleri), a more uniform species distributed from the Balkans to Western Europe (H. arborea), and a well-structured Asia Minor-Eastern European species (H. orientalis).

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.

Strong reproductive barriers in a narrow hybrid zone of West-Mediterranean green toads (Bufo viridis subgroup) with Plio-Pleistocene divergence

BACKGROUND

One key question in evolutionary biology deals with the mode and rate at which reproductive isolation accumulates during allopatric speciation. Little is known about secondary contacts of recently diverged anuran species. Here we conduct a multi-locus field study to investigate a contact zone between two lineages of green toads with an estimated divergence time of 2.7 My, and report results from preliminary experimental crosses.

RESULTS

The Sicilian endemic Bufo siculus and the Italian mainland-origin B. balearicus form a narrow hybrid zone east of Mt. Etna. Despite bidirectional mtDNA introgression over a ca. 40 km North-South cline, no F1 hybrids could be found, and nuclear genomes display almost no admixture. Populations from each side of the contact zone showed depressed genetic diversity and very strong differentiation (FST = 0.52). Preliminary experimental crosses point to a slightly reduced fitness in F1 hybrids, a strong hybrid breakdown in backcrossed offspring (F1 x parental, with very few reaching metamorphosis) and a complete and early mortality in F2 (F1 x F1).

CONCLUSION

Genetic patterns at the contact zone are molded by drift and selection. Local effective sizes are reduced by the geography and history of the contact zone, B. balearicus populations being at the front wave of a recent expansion (late Pleistocene). Selection against hybrids likely results from intrinsic genomic causes (disruption of coadapted sets of genes in backcrosses and F2-hybrids), possibly reinforced by local adaptation (the ranges of the two taxa roughly coincide with the borders of semiarid and arid climates). The absence of F1 in the field might be due to premating isolation mechanisms. Our results, show that these lineages have evolved almost complete reproductive isolation after some 2.7 My of divergence, contrasting sharply with evidence from laboratory experiments that some anuran species may still produce viable F1 offspring after > 20 My of divergence.