Negligible nuclear introgression despite complete mitochondrial capture between two species of chipmunks

The idea that species boundaries can be semipermeable to gene flow is now widely accepted but the evolutionary importance of introgressive hybridization remains unclear. Here we examine the genomic contribution of gene flow between two hybridizing chipmunk species, Tamias ruficaudus and T. amoenus. Previous studies have shown that ancient hybridization has resulted in complete fixation of introgressed T. ruficaudus mitochondrial DNA (mtDNA) in some populations of T. amoenus, but the extent of nuclear introgression is not known. We used targeted capture to sequence over 10,500 gene regions from multiple individuals of both species. We found that most of the nuclear genome is sorted between these species and that overall genealogical patterns do not show evidence for introgression. Our analysis rules out all but very minor levels of interspecific gene flow, indicating that introgressive hybridization has had little impact on the overall genetic composition of these species outside of the mitochondrial genome. Given that much of the evidence for introgression in animals has come from mtDNA, our results underscore that unraveling the importance introgressive hybridization during animal speciation will require a genome‐wide perspective that is still absent for many species.     

Phylogeography and lineage‐specific patterns of genetic diversity and molecular evolution in a group of North American skinks

Geography influences the evolutionary trajectory of species by mediating opportunities for hybridization, gene flow, demographic shifts and adaptation. We sought to understand how geography and introgression can generate species‐specific patterns of genetic diversity by examining phylogeographical relationships in the North American skink species Plestiodon multivirgatus and P. tetragrammus (Squamata: Scincidae). Using a multilocus dataset (three mitochondrial genes, four nuclear genes; a total of 3455 bp) we discovered mito‐nuclear discordance, consistent with mtDNA introgression. We further tested for evidence of species‐wide mtDNA introgression by using comparisons of genetic diversity, selection tests and extended Bayesian skyline analyses. Our findings suggest that P. multivirgatus acquired its mitochondrial genome from P. tetragrammus after their initial divergence. This putative species‐wide mitochondrial capture was further evidenced by statistically indistinguishable substitution rates between mtDNA and nDNA in P. multivirgatus. This rate discrepancy was observed in P. multivirgatus but not P. tetragrammus, which has important implications for studies that combine mtDNA and nDNA sequences when inferring time since divergence between taxa. Our findings suggest that by facilitating opportunities for interspecific introgression, geography can alter the course of molecular evolution between recently diverged lineages.     

Spatial distribution of cpDNA and mtDNA haplotypes in a hybrid zone betweenPinus pumila andP. parviflora var.pentaphylla (Pinaceae)

Pinus species exhibit paternal chloroplast inheritance and maternal mitochondrial inheritance. This independent inheritance of two cytoplasmic genomes provides an exceptional environment for discriminating female (seeds) and male (pollen) components of gene flow across hybridizing species. We obtained mitochondrial genetic markers diagnostic toP. parviflora var.pentaphylla andP. pumila by PCR amplification of the intron ofnad1 on mtDNA, and examined the spatial-distribution pattern of the mtDNA haplotypes in a hybrid zone betweenP. parviflora var.pentaphylla andP. pumila in the Tanigawa Mountains of Japan. These data, in conjunction with previous information on cpDNA haplotypes and needle morphology, revealed contrastive patterns of introgression of two cytoplasmic genomes. CpDNA introgression has occurred uni-directionally fromP. parviflora var.pentaphylla toP. pumila. Conversely, mtDNA introgression has occurred in the opposite direction, fromP. pumila toP. parviflora var.pentaphylla. Levels of introgression are roughly equivalent for cpDNA and mtDNA. The contrastive spatial distribution pattern of cpDNA and mtDNA haplotypes could be caused by differential movement of seeds and pollen for interspecific genetic exchange.     

Genetic drift or natural selection? Hybridization and asymmetric mitochondrial introgression in two Caribbean lizards (Anolis pulchellus and Anolis krugi)

Hybridization and gene introgression can occur frequently between closely related taxa, but appear to be rare phenomena among members of the species‐rich West Indian radiation of Anolis lizards. We investigated the pattern and possible mechanism of introgression between two sister species from Puerto Rico, Anolis pulchellus and Anolis krugi, using mitochondrial (ND2) and nuclear (DNAH3, NKTR) DNA sequences. Our findings demonstrated extensive introgression of A. krugi mtDNA (k‐mtDNA) into the genome of A. pulchellus in western Puerto Rico, to the extent that k‐mtDNA has mostly or completely replaced the native mtDNA of A. pulchellus on this part of the island. We proposed two not mutually exclusive scenarios to account for the interspecific matings between A. pulchellus and A. krugi. We inferred that hybridization events occurred independently in several populations, and determined that k‐mtDNA haplotypes harboured in individuals of A. pulchellus can be assigned to four of the five major mtDNA clades of A. krugi. Further, the spatial distribution of k‐mtDNA clades in the two species is largely congruent. Based on this evidence, we concluded that natural selection was the probable driving mechanism for the extensive k‐mtDNA introgression into A. pulchellus. Our two nuclear data sets yielded different results. DNAH3 showed reciprocal monophyly of A. pulchellus and A. krugi, indicating no effect of hybridization on this marker. In contrast, the two species shared nine NKTR alleles, probably due to incomplete lineage sorting. Our study system will provide an excellent opportunity to experimentally assess the behavioural and ecological mechanisms that can lead to hybridization in closely related taxa.     

Nuclear introgression without mitochondrial introgression in two turtle species exhibiting sex‐specific trophic differentiation

Despite the presence of reproductive barriers between species, interspecific gene introgression has been documented in a range of natural systems. Comparing patterns of genetic introgression in biparental versus matrilineal markers can potentially reveal sex‐specific barriers to interspecific gene flow. Hybridization has been documented in the freshwater turtles Graptemys geographica and G. pseudogeographica, whose ranges are largely sympatric. Morphological differentiation between the species is restricted to females, with female G. geographica possessing large heads and jaws compared to the narrow heads of G. pseudogeographica females. If hybrid females are morphologically intermediate, they may be less successful at exploiting parental feeding niches, thereby limiting the introgression of maternally inherited, but not biparental, molecular markers. We paired sequence data with stable isotope analysis and examined sex‐specific genetic introgression and trophic differentiation in sympatric populations of G. geographica and G. pseudogeographica. We observed introgression from G. pseudogeographica into G. geographica at three nuclear loci, but not at the mitochondrial locus. Analysis of ?¹⁵N and ?¹³C was consistent with species differences in trophic positioning in females, but not males. These results suggest that ecological divergence in females may reduce the opportunity for gene flow in this system.     

Evolutionary basis of mitonuclear discordance between sister species of mole salamanders (Ambystoma sp.)

Distinct genetic markers should show similar patterns of differentiation between species reflecting their common evolutionary histories, yet there are increasing examples of differences in the biogeographic distribution of species‐specific nuclear (nuDNA) and mitochondrial DNA (mtDNA) variants within and between species. Identifying the evolutionary processes that underlie these anomalous patterns of genetic differentiation is an important goal. Here, we analyse the putative mitonuclear discordance observed between sister species of mole salamanders (Ambystoma barbouri and A. texanum) in which A. barbouri‐specific mtDNA is found in animals located within the range of A. texanum. We test three hypotheses for this discordance (undetected range expansion, mtDNA introgression, and hybridization) using nuDNA and mtDNA data analysed with methods that varied in the parameters estimated and the timescales measured. Results from a Bayesian clustering technique (structure ), bidirectional estimates of gene flow (migrate ‐n and IMa2) and phylogeny‐based methods (*beast , buck y) all support the conclusion that the discordance is due to geographically restricted mtDNA introgression from A. barbouri into A. texanum. Limited data on species‐specific tooth morphology match this conclusion. Significant differences in environmental conditions exist between sites where A. texanum with and without A. barbouri‐like mtDNA occur, suggesting a possible role for selection in the process of introgression. Overall, our study provides a general example of the value of using complimentary analyses to make inferences of the directionality, timescale, and source of mtDNA introgression in animals.     

When a clonal genome finds its way back to a sexual species: evidence from ongoing but rare introgression in the hybridogenetic water frog complex

Besides several exceptions, asexual metazoans are usually viewed as ephemeral sinks for genomes, which become ‘frozen’ in clonal lineages after their emergence from ancestral sexual species. Here, we investigated whether and at what rate the asexuals are able to introgress their genomes back into the parental sexual population, thus more or less importantly affecting the gene pools of sexual species. We focused on hybridogenetic hybrids of western Palaearctic water frogs (Pelophylax esculentus), which originate through hybridization between P. ridibundus and P. lessonae, but transmit only clonal ridibundus genome into their gametes. Although usually mating with P. lessonae, P. esculentus may upon mating with P. ridibundus or another hybrid produce sexually reproducing P. ridibundus offspring with the introgressed ex‐clonal genome. We compared the rate of nuclear amplified fragment length polymorphism (AFLP) and mitochondrial introgression in two types of populations, that is, those where P. ridibundus occurs in isolation and those where it lives with the hybridogens. Although significant differentiation (Φpt) between sexual and clonal ridibundus genomes suggested limited gene flow between sexuals and hybridogens, a non‐negligible (~5%) proportion of P. ridibundus bore introgressed mtDNA and AFLP markers. Whereas transfer of mtDNA was exclusively unidirectional, introgression of nuclear markers was bidirectional. The proportion of introgressed P. ridibundus was highest in syntopic populations with P. esculentus, proving an ongoing and site‐specific interspecific genetic transfer mediated by hybridogenetic hybrids. It turns out that asexual hybrids are not just a sink for genes of sexual species, but may significantly influence the genetic architecture of their sexual counterparts.     

Bidirectional cytoplasmic and nuclear introgression in the New World cottons Gossypium barbadense and G. Hirsutum (Malvaceae)

Patterns of interspecific cytoplasmic (plastid and mitochondrial) and nuclear introgression are typically asymmetrical: cytoplasmic gene flow or “capture” is frequently observed without evidence of nuclear introgression. In contrast, nuclear introgression without concomitant cytoplasmic introgression has rarely been demonstrated. Gossypium barbadense L. and G. hirsutum L. have large indigenous ranges in the New World semiarid tropics, with an extensive area of sympatry in the Caribbean and Central America. Numerous accessions of both species were surveyed for diagnostic cpDNA restriction sites. These data, in conjunction with previous information on nuclear markers, lead to several conclusions: 1) introgression between G. hirsutum and G. barbadense is bidirectional for both nuclear and cytoplasmic genes; 2) patterns of introgression between the two species are not symmetrical—in G. barbadense, introgression of G. hirsutum alleles is largely restricted to modem cultivars and is uncommon in areas of sympatry; in contrast, introgression of G. barbadense alleles into G. hirsutum is relatively common in areas of sympatry and is rare in modem cultivars; 3) nuclear introgression is geographically more widespread and more frequently detected than cytoplasmic introgression. Several mechanisms may underlie the differential patterns of cytoplasmic and nuclear gene flow observed, including differential fitness of infraspecific and interspecific cytonuclear combinations and selection against female function in interspecific backcrosses. Possible explanations for the observed asymmetrical patterns of introgression include differences in population sizes combined with phenological differences that bias interspecific pollen transfer.     

Discordant patterns of introgression across a narrow hybrid zone between two cryptic lineages of an Iberian endemic newt

The study of natural hybrid zones can illuminate aspects of lineage divergence and speciation in morphologically cryptic taxa. We studied a hybrid zone between two highly divergent but morphologically similar lineages (south‐western and south‐eastern) of the Iberian endemic Bosca's newt (Lissotriton boscai) in SW Iberia with a multilocus dataset (microsatellites, nuclear and mitochondrial genes). STRUCTURE and NEWHYBRIDS analyses retrieved few admixed individuals, which classified as backcrosses involving parental individuals of the south‐western lineage. Our results show asymmetric introgression of mtDNA beyond the contact from this lineage into the south‐eastern lineage. Analysis of nongeographic introgression patterns revealed asymmetries in the direction of introgression, but except for mtDNA, we did not find evidence for nonconcordant introgression patterns across nuclear loci. Analysis of a 150‐km transect across the hybrid zone showed broadly coincident cline widths (ca. 3.2–27.9 km), and concordant cline centres across all markers, except for mtDNA that is displaced ca. 60 km northward. Results from ecological niche modelling show that the hybrid zone is in a climatically homogenous area with suitable habitat for the species, suggesting that contact between the two lineages is unlikely to occur further south as their distributions are currently separated by an extensive area of unfavourable habitat. Taken together, our findings suggest the genetic structure of this hybrid zone results from the interplay of historical (biogeographic) and population‐level processes. The narrowness and coincidence of genetic clines can be explained by weak selection against hybrids and reflect a degree of reproductive isolation that is consistent with cryptic speciation.     

Historical male‐mediated introgression in horseshoe bats revealed by multilocus DNA sequence data

Instances of hybridization between mammalian taxa in the wild are rarely documented. To test for introgression between sibling species of horseshoe bat (Rhinolophus yunanensis and R. pearsoni) and two subspecies of the latter (R. p. pearsoni and R. p. chinensis), we sequenced two mtDNA and two ncDNA markers in individuals sampled from multiple localities within their overlapping ranges. The interspecific mtDNA gene tree corresponded to the expected taxonomic divisions, and coalescent‐based analyses suggested divergence occurred around 4 MYA. However, these relationships strongly conflicted with those recovered from two independent nuclear gene trees, in which R. yunanensis clustered with R. p. pearsoni to the exclusion of R. p. chinensis. This geographically widespread discordance is best explained by large‐scale historical introgression of ncDNA from R. yunanensis to R. pearsoni by male‐mediated exchange in mixed species colonies during Pleistocene glacial periods, when ranges may have contracted and overlapped more than at present. Further species tree–gene tree conflicts were detected between R. p. pearsoni and R. p. chinensis, also indicating past and/or current introgression in their overlapping regions. However, here the patterns point to asymmetric mtDNA introgression without ncDNA introgression. Analyses of coalescence times indicate this exchange has occurred subsequent to the divergence of these subspecies from their common ancestor. Our work highlights the importance of using multiple data sets for reconstructing phylogeographic histories and resolving taxonomic relationships.     

Limited gene exchange between two sister species of leaf beetles within a hybrid zone in the Alps

Analysing genomic variation within and between sister species is a first step towards understanding species boundaries. We focused on two sister species of cold‐resistant leaf beetles, Gonioctena quinquepunctata and G. intermedia, whose ranges overlap in the Alps. A previous study of DNA sequence variation had revealed multiple instances of mitochondrial genome introgression in this region, suggesting recent hybridization between the two species. To evaluate the extent of gene exchange resulting from these hybridization events, we sampled individuals of both species inside and outside the hybrid zone and analysed genomic variation among them using RAD‐seq markers. Individual levels of introgression in the nuclear genome were estimated first by defining species‐specific SNPs (displaying a fixed difference between species) a priori and second by using model‐based methods. Both types of analyses indicated little gene exchange, if any, between species at the level of the nuclear genome. Whereas the first method suggested slightly more gene flow, we argue that it has likely overestimated introgression in the phylogeographic context of this study. We conclude that strong intrinsic barriers prevent genetic exchange at the level of the nuclear genome between the two species. The apparent discrepancy observed between introgression occurring in the nuclear and mitochondrial genomes could be explained by selection acting in favour of the latter. Also, these results have consequences for the phylogeographic study of each species, since we can assume that genetic diversity in the overlapping portion of their ranges is not the product of introgression.     

On shortcomings of using mtDNA sequence divergence for the systematics of hares (genus Lepus): An example from cape hares

Despite the well-known fact that evolutionary patterns of single genes or sequences are not necessarily paralleled by organismic evolution, using mitochondrial sequence divergence for inferring phylogenetic relationships among hare species is not uncommon. Earlier studies reported interspecific introgression in the genus Lepus that may slur systematic conclusions drawn exclusively from mtDNA data. We examined pairwise divergence in partial mtHV-1 sequences and nuclear gene pools based on microsatellite and allozyme loci separately in a South and North African population of cape hares, Lepus capensis, and did not find significant correspondence on the individual level within either population. Also, the former population had a significantly higher level of mtHV-1 sequence divergence compared to the latter, but levels of individual nuclear gene pool divergence did not differ significantly between the two populations. Hence, the absence of correspondence between mtHV-1 sequence divergence and nuclear gene pool divergence among individuals within a population might be independent of the population-specific level of differentiation among individuals. Our results complement earlier findings in hares, and strongly recommend to include nuclear gene pool evidence for systematic inferences within this genus, even under the absence of mitochondrial introgression.     

The role of mitochondrial introgression in illuminating the evolutionary history of Nearctic treefrogs

Inferring the evolutionary history of lineages often becomes difficult when gene histories are in conflict with each other. Introgression, for example, can cause DNA sequences from one species to be more similar to sequences of a different species and lead to incongruence amongst gene trees. However, incorporating congruent and incongruent locus‐specific phylogenetic estimates with the geographical distribution of lineages may provide valuable insight into evolutionary processes important to speciation. In this study, we investigated mitochondrial introgression within the Hyla eximia group to better understand its role in illuminating the evolutionary history and phylogeography of these treefrogs. We reconstructed and compared the matrilineal history of the Hyla eximia group with estimates of evolutionary history inferred from nuclear genes. We tested for introgression within the mitochondrial and nuclear genes using a posterior predictive checking approach. Reconstructions of the species tree based on the mitochondrial DNA (mtDNA) and nuclear DNA data were strongly discordant. Introgression between lineages was widespread in the mtDNA data set (145 occurrences amongst 11 of the 16 lineages), but uncommon in the nuclear genes (12 occurrences amongst four of the 16 lineages). Nonetheless, the geographical structuring of mtDNA within species provides valuable information on biogeographical areas, ancient areas of hybridization, and unique histories of lineages within the H. eximia group. These results suggest that the combination of nuclear, mitochondrial, and spatial information can provide a more complete picture of ‘how evolutionary history played out’, particularly in cases where mitochondrial introgression is known to occur. © 2014 The Linnean Society of London     

Mitonuclear discordance is caused by rampant mitochondrial introgression in Neodiprion (Hymenoptera: Diprionidae) sawflies

We investigate the pervasiveness of hybridization and mitochondrial introgression in Neodiprion Rohwer (Hymenoptera; Diprionidae), a Holarctic genus of conifer-feeding sawflies. A phylogenetic analysis of the lecontei species group revealed extensive discordance between a contiguous mitochondrial region spanning three genes (COI, tRNA-leucine, and COII) and three nuclear loci (EF1alpha, CAD, and an anonymous nuclear locus). Bayesian tests of monophyly and Shimodaira-Hasegawa (SH) tests of topological congruence were consistent with mitochondrial introgression; however, these patterns could also be explained by lineage sorting (i.e., deep coalescence). Therefore, to explicitly test the mitochondrial introgression hypothesis, we used a novel application of coalescent-based isolation with migration (IM) models to measure interspecific gene flow at each locus. In support of our hypothesis, mitochondrial gene flow was consistently higher than nuclear gene flow across 120 pairwise species comparisons (P < 1 x 10(-12)). We combine phylogenetic and coalescent evidence to identify likely cases of recent and ancient introgression in Neodiprion, and based on these observations, we hypothesize that shared hosts and/or pheromones facilitate hybridization, whereas disparate abundances between hybridizing species promote mitochondrial introgression. Our results carry implications for phylogenetic analysis, and we advocate the separation of high and low gene flow regions to inform analyses of hybridization and speciational history, respectively.     

Species‐specific habitat preferences do not shape the structure of a crested newt hybrid zone (Triturus cristatus x T. carnifex)

Reproductive isolation barriers maintain the integrity of species by preventing interspecific gene flow. They involve temporal, habitat or behavioral isolation acting before fertilization, and postzygotic isolation manifested as hybrid mortality or sterility. One of the approaches of how to study reproductive isolation barriers is through the analysis of hybrid zones. In this paper, we describe the structure of a hybrid zone between two crested newt species (Triturus cristatus and T. carnifex) in the southern part of the Czech Republic using morphological, microsatellite, and mitochondrial (mtDNA) markers. Specifically, we tested the hypothesis that the structure of the hybrid zone is maintained by species‐specific habitat preferences. Comparing the genetic structure of populations with geographical and ecological parameters, we found that the hybrid zone was structured primarily geographically, with T. cristatus‐like populations occurring in the northeast and T. carnifex‐like populations in the southwest. Despite T. cristatus tending to occur in deeper ponds and T. carnifex on localities with more shading, the effect of both ecological parameters on the structure of the zone was minimal. Next, we corroborated that T. carnifex individuals and some hybrids possess mtDNA of T. dobrogicus, whose nuclear background was not detected in the studied hybrid zone. Hybridization between T. carnifex and T. dobrogicus (resulting in unidirectional mtDNA introgression) had to predate subsequent formation of the hybrid zone between T. cristatus and T. carnifex. Populations of crested newts in the southern part of the Czech Republic thus represent a genetic mosaic of nuclear and mitochondrial genomes of three species.     

Testing hypotheses for genealogical discordance in a rainforest lizard

Genealogical discordance, or when different genes tell distinct stories although they evolved under a shared history, often emerges from either coalescent stochasticity or introgression. In this study, we present a strong case of mito‐nuclear genealogical discordance in the Australian rainforest lizard species complex of Saproscincus basiliscus and S. lewisi. One of the lineages that comprises this complex, the Southern S. basiliscus lineage, is deeply divergent at the mitochondrial genome but shows markedly less divergence at the nuclear genome. By placing our results in a comparative context and reconstructing the lineages' demography via multilocus and coalescent‐based approximate Bayesian computation methods, we test hypotheses for how coalescent variance and introgression contribute to this pattern. These analyses suggest that the observed genealogical discordance likely results from introgression. Further, to generate such strong discordance, introgression probably acted in concert with other factors promoting asymmetric gene flow between the mitochondrial and nuclear genomes, such as selection or sex‐biased dispersal. This study offers a framework for testing sources of genealogical discordance and suggests that historical introgression can be an important force shaping the genetic diversity of species and their populations.     

Current and historical hybridization with differential introgression among three species of cyprinid fishes (genus <Emphasis Type="Italic">Cyprinella</Emphasis>)

Hybridization is common among freshwater fishes, particular among the Cyprinidae. We used two mitochondrial genes and one nuclear gene to characterize hybridization among two species pairs of Cyprinella in southwestern North America. Genalogical patterns revealed that C. lutrensis and C. venusta are currently hybridizing in several localities producing apparent F₁, F₂ and backcross generations, yet there was no evidence for introgression outside of local hybrid zones. Alternatively, mitochondrial haplotypes from C. lutrensis appear to have introgressed into a C. lepida population in the Nueces River completely replacing the native C. lepida haplotype. There was no evidence of introgression of nuclear DNA and there does not appear to be ongoing hybridization. The population of C. lepida from the nearby Frio River exhibits no evidence of hybridization with C. lutrensis. Thus, contact between C. lutrensis and C. venusta results in the formation of localized hybrid swarms, while contact between C. lutrensis and C. lepida has resulted in complete mitochondrial introgression in the Nueces River or no apparent hybridization in the Frio River. The three different outcomes of contact between these species illustrate the variable nature of interspecific reproductive interactions and provide an excellent system in which to better understand the factors influencing hybridization among freshwater fishes.     

Wolbachia have made it twice: Hybrid introgression between two sister species of Eurema butterflies

Wolbachia, cytoplasmically inherited endosymbionts of arthropods, are known to hijack their host reproduction in various ways to increase their own vertical transmission. This may lead to the selective sweep of associated mitochondria, which can have a large impact on the evolution of mitochondrial lineages. In Japan, two different Wolbacahia strains (wCI and wFem) are found in two sister species of pierid butterflies, Eurema mandarina and Eurema hecabe. In both species, females infected with wCI (C females) produce offspring with a nearly 1:1 sex ratio, while females infected with both wCI and wFem (CF females) produce all‐female offspring. Previous studies have suggested the historical occurrence of hybrid introgression in C individuals between the two species. Furthermore, hybrid introgression in CF individuals is suggested by the distinct mitochondrial lineages between C females and CF females of E. mandarina. In this study, we performed phylogenetic analyses based on nuclear DNA and mitochondrial DNA markers of E. hecabe with previously published data on E. mandarina. We found that the nuclear DNA of this species significantly diverged from that of E. mandarina. By contrast, mitochondrial DNA haplotypes comprised two clades, mostly reflecting Wolbachia infection status rather than the individual species. Collectively, our results support the previously suggested occurrence of two independent historical events wherein the cytoplasms of CF females and C females moved between E. hecabe and E. mandarina through hybrid introgression.     

The crested newt Triturus cristatus recolonized temperate Eurasia from an extra‐Mediterranean glacial refugium

We assess the role of the Carpathians as an extra‐Mediterranean glacial refugium for the crested newt Triturus cristatus. We combine a multilocus phylogeography (one mitochondrial protein‐coding gene, three nuclear introns, and one major histocompatibility complex gene) with species distribution modelling (projected on current and Last Glacial Maximum climate layers). All genetic markers consistently show extensive genetic variation within and genetic depletion outside the Carpathians. The species distribution model suggests that most of the current range was unsuitable at the Last Glacial Maximum, but a small suitable area remained in the Carpathians. Triturus cristatus dramatically expanded its postglacial range, colonizing much of temperate Eurasia from a glacial refugium in the Carpathians. Within the Carpathians, T. cristatus persisted in multiple geographically discrete regions, providing further support for a Carpathian ‘refugia within refugia’ scenario. © 2015 The Linnean Society of London, Biological Journal of the Linnean Society, 2015, 114 , 574–587.