Molecular evidence for a species complex in the patagonian lizard Liolaemus bibronii and phylogeography of the closely related Liolaemus gracilis (Squamata: Liolaemini)

The lizard genus Liolaemus is endemic to temperate South America and includes 190 species. Liolaemus bibronii has a large geographic distribution and inhabits a great diversity of habitats, including the Monte, Steppe, and high Andean grassland environments. Liolaemus gracilis has a similar body size and shape to L. bibronii; the two are parapatrically distributed, and L. gracilis is also widely distributed. Here we use the mtDNA cytb sequence data of these two species to investigate lizard phylogeographic patterns in southern South America. L. bibronii is paraphyletic with respect to L. gracilis, Liolaemus ramirezae, Liolaemus robertmertensi and Liolaemus saxatilis; it is composed of many genetically different allopatric haploclades, some of which are reciprocally monophyletic. We also found evidence for introgression between L. bibronii and L. gracilis in the same area that introgression was hypothesized in the Liolaemus darwinii complex. We discuss the distribution of the major haploclades with inferences of their population histories, the concordance of these clades' distributions and histories with other lizard complexes studied with the same markers and methods, and taxonomic implications of these results.     

Phylogeography and introgressive hybridization: chipmunks (genus Tamias) in the northern Rocky Mountains

Abstract. If phylogeographic studies are to be broadly used for assessing population-level processes relevant to speciation and systematics, the ability to identify and incorporate instances of hybridization into the analytical framework is essential. Here, we examine the evolutionary history of two chipmunk species, Tamias ruficaudus and Tamias amoenus , in the northern Rocky Mountains by integrating multivariate morphometrics of bacular (os penis) variation, phylogenetic estimation, and nested clade analysis with regional biogeography. Our results indicate multiple examples of mitochondrial DNA introgression layered within the evolutionary history of these nonsister species. Three of these events are most consistent with recent and/or ongoing asymmetric introgression of mitochondrial DNA across morphologically defined secondary contact zones. In addition, we find preliminary evidence where a fourth instance of nonconcordant characters may represent complete fixation of introgressed mitochondrial DNA via a more ancient hybridization event, although alternative explanations of convergence or incomplete sorting of ancestral polymorphisms cannot be dismissed with these data. The demonstration of hybridization among chipmunks with strongly differentiated bacular morphology contradicts long-standing assumptions that variation within this character is diagnostic of complete reproductive isolation within Tamias . Our results illustrate the utility of phylogeographic analyses for detecting instances of reticulate evolution and for incorporating this and other information in the inference of the evolutionary history of species.     

Mitochondrial introgression and incomplete lineage sorting through space and time: phylogenetics of crotaphytid lizards

We investigate the roles of mitochondrial introgression and incomplete lineage sorting during the phylogenetic history of crotaphytid lizards. Our Bayesian phylogenetic estimate for Crotaphytidae is based on analysis of mitochondrial DNA sequence data for 408 individuals representing the 12 extant species of Crotaphytus and Gambelia. The mitochondrial phylogeny disagrees in several respects with a previously published morphological tree, as well as with conventional species designations, and we conclude that some of this disagreement stems from hybridization-mediated mitochondrial introgression, as well as from incomplete lineage sorting. Unidirectional introgression of Crotaphytus collaris (western collared lizard) mitochondria into C. reticulatus (reticulate collared lizard) populations in the Rio Grande Valley of Texas has resulted in the replacement of ancestral C. reticulatus mitochondria over approximately two-thirds of the total range of the species, a linear distance of approximately 270 km. Introgression of C. collaris mitochondria into C. bicinctores (Great Basin collared lizard) populations in southwestern Arizona requires a more complex scenario because at least three temporally separated and superimposed introgression events appear to have occurred in this region. We propose an "introgression conveyor" model to explain this unique pattern of mitochondrial variation in this region. We show with ecological niche modeling that the predicted geographical ranges of C. collaris, C. bicinctores, and C. reticulatus during glacial maxima could have provided enhanced opportunities for past hybridization. Our analyses suggest that incomplete lineage sorting and/or introgression has further confounded the phylogenetic placements of additional species including C. nebrius, C. vestigium, C. insularis, C. grismeri, and perhaps G. copei. Despite many independent instances of interspecific hybridization among crotaphytid lizards, the species continue to maintain morphological and geographic cohesiveness throughout their ranges.     

Molecular phylogeny of three Mesalina (Reptilia: Lacertidae) species (M. guttulata, M. brevirostris and M. bahaeldini) from North Africa and the Middle East: another case of paraphyly?

Mesalina is a widespread lacertid genus occurring throughout the Saharo-Sindian region from North Africa to Pakistan. It has been through a series of taxonomic revisions, but the phylogenetic relationships among the species remain unclear. In this study we estimate the phylogeographic structure of M. guttulata across most of its distributional range and we evaluate the relationships between M. guttulata and the sympatric species M. brevirostris and M. bahaeldini using partial mitochondrial DNA (mtDNA) sequences (cyt b and 16S). M. guttulata and M. brevirostris represent species complexes, whereas M. bahaeldini considered before as M. guttulata is a recently described species with very restricted distribution. Here we present the first evidence that M. guttulata is a paraphyletic taxon with respect to M. bahaeldini, while M. brevirostris proves to be a polytypic species or even a species complex, confirming previous studies. Although mtDNA markers have several properties that make them suitable for phylogeographic studies, they are not free of difficulties. Phylogeographic inferences within and between closely related species can be mislead by introgression and retention of ancestral polymorphism (incomplete lineage sorting). However, the present distribution pattern, the estimated times of divergence and the significant variation in morphology within M. guttulata led us to accept that the paraphyletic pattern observed, is most likely due to inaccurate taxonomy. Our hypothesis is that what has hitherto been considered as intraspecific variation, actually reflects species-level variation. Furthermore, our biogeographic analyses and the estimated time of divergences suggest that the present distribution of M. guttulata was the result of several dispersal and vicariant events, which are associated with historical changes (climatic oscillations and paleogeographic barriers) of late Miocene and Pliocene period.     

Genetic admixture between the Iberian endemic lizards Podarcis bocagei and Podarcis carbonelli: evidence for limited natural hybridization and a bimodal hybrid zone

When recently diverged taxa come into contact, the extent of introgression between them is related to the degree of differentiation that they have achieved. Studying contact zones is therefore essential to understand if differentiated taxa are reproductively isolated and, ultimately, if they are likely to remain distinct. Recent work on Iberian and North African wall lizards ( Podarcis ) has documented the existence of multiple evolutionary units, diagnosable both by genetic markers and morphology, but suggests that gene flow between distinct forms has occurred. Therefore, we were interested in evaluating how species boundaries are maintained in the areas where they meet. In this work, we study the contact zone between Podarcis bocagei and Podarcis carbonelli . We sampled a transect including the only locality where these two lizards are known to occur in syntopy and analysed a battery of 15 unlinked nuclear genetic markers and mitochondrial DNA. We also conducted a preliminary analysis of morphology and fertility. Using model-based clustering approaches, we show that the two species hybridize in the population where they have direct contact, but evidences of introgression are low for nearby populations. Although a significant number of individuals show evidence of admixture, this hybrid zone is clearly bimodal, suggesting strong barriers to gene flow, of which the putative nature are discussed. Interestingly, morphological analyses do not support the existence of intermediate forms among individuals that are admixed genetically. Taken together, these results constitute further evidence validating P. bocagei and P. carbonelli as distinct species.     

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.     

Mitochondrial DNA phylogeography of the Mesoamerican spiny-tailed lizards (Ctenosaura quinquecarinata complex): historical biogeography, species status and conservation

Through the examination of past and present distributions of plants and animals, historical biogeographers have provided many insights on the dynamics of the massive organismal exchange between North and South America. However, relatively few phylogeographic studies have been attempted in the land bridge of Mesoamerica despite its importance to better understand the evolutionary forces influencing this biodiversity 'hotspot'. Here we use mitochondrial DNA sequence data from fresh samples and formalin-fixed museum specimens to investigate the genetic and biogeographic diversity of the threatened Mesoamerican spiny-tailed lizards of the Ctenosaura quinquecarinata complex. Species boundaries and their phylogeographic patterns are examined to better understand their disjunct distribution. Three monophyletic, allopatric lineages are established using mtDNA phylogenetic and nested clade analyses in (i) northern: México, (ii) central: Guatemala, El Salvador and Honduras, and (iii) southern: Nicaragua and Costa Rica. The average sequence divergence observed between lineages varied between 2.0% and 3.7% indicating that they do not represent a very recent split and the patterns of divergence support the recently established nomenclature of C. quinquecarinata, Ctenosaura flavidorsalis and Ctenosaura oaxacana. Considering the geological history of Mesoamerica and the observed phylogeographic patterns of these lizards, major evolutionary episodes of their radiation in Mesoamerica are postulated and are indicative of the regions' geological complexity. The implications of these findings for the historical biogeography, taxonomy and conservation of these lizards are discussed.     

Conflicting mitochondrial and nuclear phylogeographic signals and evolution of host-plant shifts in the boreo-montane leaf beetle Chrysomela lapponica

We conducted a phylogeographic study on the cold-adapted leaf beetle Chrysomela lapponica, that feeds on willow or birch, by sampling several populations throughout most of the geographic distribution of the species, and by sequencing for each individual one mitochondrial and two nuclear DNA fragments. Patterns of DNA sequence variation from the mitochondrial and nuclear loci, as displayed in the median-joining networks, appear to display contradicting historical signal: a deep genealogical divergence is observed with the mitochondrial genome between the Alpine population and all other populations found in the Euro-Siberian distribution of the species, that is completely absent with both nuclear loci. We use coalescence simulations of DNA sequence evolution to test the hypothesis that this apparent conflict is compatible with a neutral model of sequence evolution (i.e., to check whether the stochastic nature of the coalescence process can explain these patterns). Because the simulations show that this is highly unlikely, we consider two alternative hypotheses: (1) introgression of the mitochondrial genome of another species and (2) the effect of natural selection. Although introgression is the most plausible explanation, we fail to identify the source species of the introgressed mitochondrial genome among all known species closely related to C. lapponica. We therefore suggest that the putative introgression event is ancient and the source species is either extinct or currently outside the geographic range of C. lapponica explored in this study. The observed DNA sequence variation also suggests that a host-plant shift from willow to birch has occurred recently and independently in each of the three birch-feeding populations. This emphasizes further the relative ease with which these beetles can escape their ancestral host-plant specialization on willow, but shows at the same time that host-plant shifts are highly constrained, as they only occur between willow and birch.     

Molecular evidence for hybridization between two Australian desert skinks, Ctenotus leonhardii and Ctenotus quattuordecimlineatus (Scincidae: Squamata)

Australian scincid lizards in the genus Ctenotus constitute the most diverse vertebrate radiation in Australia. However, the evolutionary processes that have generated this diversity remain elusive, in part because both interspecific phylogenetic relationships and phylogeographic structure within Ctenotus species remain poorly known. Here we use nucleotide sequences from a mitochondrial locus and two nuclear introns to investigate broad-scale phylogeographic patterns within Ctenotus leonhardii and C. quattuordecimlineatus, two geographically widespread species of skinks that were found to have a surprisingly close genetic relationship in a previous molecular phylogenetic study. We demonstrate that the apparent close relationship between these ecologically and phenotypically distinct taxa is attributable to mitochondrial introgression from C. quattuordecimlineatus to C. leonhardii. In the western deserts, Ctenotus leonhardii individuals carry mtDNA lineages that are derived from C. quattuordecimlineatus mtDNA lineages from that geographic region. Coalescent simulations indicate that this pattern is unlikely to have resulted from incomplete lineage sorting, implicating introgressive hybridization as the cause of this regional gene-tree discordance.     

Combining phylogeography with distribution modeling: multiple Pleistocene range expansions in a parthenogenetic gecko from the Australian arid zone

Phylogenetic and geographic evidence suggest that many parthenogenetic organisms have evolved recently and have spread rapidly. These patterns play a critical role in our understanding of the relative merits of sexual versus asexual reproductive modes, yet their interpretation is often hampered by a lack of detail. Here we present a detailed phylogeographic study of a vertebrate parthenogen, the Australian gecko Heteronotia binoei, in combination with statistical and biophysical modeling of its distribution during the last glacial maximum. Parthenogenetic H. binoei occur in the Australian arid zone and have the widest range of any known vertebrate parthenogen. They are broadly sympatric with their sexual counterparts, from which they arose via hybridization. We have applied nested clade phylogeographic, effective migration, and mismatch distribution analyses to mitochondrial DNA (mtDNA) sequences obtained for 319 individuals sampled throughout the known geographic ranges of two parthenogenetic mitochondrial lineages. These analyses provide strong evidence for past range expansion events from west to east across the arid zone, and for continuing eastward range expansion. Parthenogen formation and range expansion events date to the late Pleistocene, with one lineage expanding from the northwest of its present range around 240,000 years ago and the second lineage expanding from the far west around 70,000 years ago. Statistical and biophysical distribution models support these inferences of recent range expansion, with suitable climatic conditions during the last glacial maximum most likely limited to parts of the arid zone north and west of much of the current ranges of these lineages. Combination of phylogeographic analyses and distribution modeling allowed considerably stronger inferences of the history of this complex than either would in isolation, illustrating the power of combining complementary analytical approaches.     

Introgression from Lepus europaeus to L. timidus in Russia revealed by mitochondrial single nucleotide polymorphisms and nuclear microsatellites

Hybridisation among wild mammal populations may lead to introgression of genes and genomes over the species barrier. In Sweden, in northern Europe, and on the Iberian Peninsula in southern Europe, mitochondrial DNA from L. timidus occurs among L. europaeus specimens, presumably as a result of interspecific hybridisation. In Russia, the species are believed to hybridise as well, but no investigations have confirmed introgression. Here we develop species diagnostic single nucleotide polymorphisms in the mitochondrial genomes and combine them with analysis of nuclear microsatellite markers to investigate hybridisation and introgression in 71 Lepus specimens from Russia. A total of 58 specimens are typical representatives of either species. An additional nine specimens have slightly intermediate genotypes, potentially as a result of introgression of nuclear genes. Finally, we find three specimens with L. europaeus mitochondrial genome and apparent L. timidus nuclear genome. This indicates that the reciprocal transfer of mtDNA occur among Russian populations of these species. Our observation differs from previous observations of mtDNA introgression in Sweden and Iberia, and provides further support for a reticulated mode of introgression within the genus Lepus.     

The puzzling phylogeography of the haplochromine cichlid fish Astatotilapia burtoni

Astatotilapia burtoni is a member of the “modern haplochromines,” the most species‐rich lineage within the family of cichlid fishes. Although the species has been in use as research model in various fields of research since almost seven decades, including developmental biology, neurobiology, genetics and genomics, and behavioral biology, little is known about its spatial distribution and phylogeography. Here, we examine the population structure and phylogeographic history of A. burtoni throughout its entire distribution range in the Lake Tanganyika basin. In addition, we include several A. burtoni laboratory strains to trace back their origin from wild populations. To this end, we reconstruct phylogenetic relationships based on sequences of the mitochondrial DNA (mtDNA) control region (d‐loop) as well as thousands of genomewide single nucleotide polymorphisms (SNPs) derived from restriction‐associated DNA sequencing. Our analyses reveal high population structure and deep divergence among several lineages, however, with discordant nuclear and mtDNA phylogenetic inferences. Whereas the SNP‐based phylogenetic hypothesis uncovers an unexpectedly deep split in A. burtoni, separating the populations in the southern part of the Lake Tanganyika basin from those in the northern part, analyses of the mtDNA control region suggest deep divergence between populations from the southwestern shoreline and populations from the northern and southeastern shorelines of Lake Tanganyika. This phylogeographic pattern and mitochondrial haplotype sharing between populations from the very North and the very South of Lake Tanganyika can only partly be explained by introgression linked to lake‐level fluctuations leading to past contact zones between otherwise isolated populations and large‐scale migration events.     

Invasion from the cold past: extensive introgression of mountain hare (Lepus timidus) mitochondrial DNA into three other hare species in northern Iberia

Mitochondrial DNA introgression from Lepus timidus into Lepus granatensis and Lepus europaeus was recently reported in Iberia, although L. timidus presumably retreated from this region at the end of the last ice age. Here we assess the extent of this ancient mtDNA introgression by RFLP analysis of 695 specimens representing the three hare species present in Iberia. The introgressed L. timidus lineage was found in 23 of the 37 populations sampled. It is almost fixed in L. europaeus across its Iberian range in the Pyrenean foothills, and in L. granatensis, which occupies the rest of the peninsula, it is predominant in the north and gradually disappears further south. We also found it in Lepus castroviejoi, a species endemic to Cantabria. Multiple hybridizations and, potentially, a selective advantage for the L. timidus lineage can explain the remarkable taxonomic and geographical range of this mitochondrial introgression.     

Phylogeography and demographic history of the neotropical otter (Lontra longicaudis)

The Neotropical otter (Lontra longicaudis) is a medium-sized semiaquatic carnivore with a broad distribution in the Neotropical region. Despite being apparently common in many areas, it is one of the least known otters, and genetic studies on this species are scarce. Here, we have investigated its genetic diversity, population structure, and demographic history across a large portion of its geographic range by analyzing 1471 base pairs (bp) of mitochondrial DNA from 52 individuals. Our results indicate that L. longicaudis presents high levels of genetic diversity and a consistent phylogeographic pattern, suggesting the existence of at least 4 distinct evolutionary lineages in South America. The observed phylogeographic partitions are partially congruent with the subspecies classification previously proposed for this species. Coalescence-based analyses indicate that Neotropical otter mitochondrial DNA lineages have shared a rather recent common ancestor, approximately 0.5 Ma, and have subsequently diversified into the observed phylogroups. A consistent scenario of recent population expansion was identified in Eastern South America based on several complementary analyses of historical demography. The results obtained here provide novel insights on the evolutionary history of this largely unknown Neotropical mustelid and should be useful to design conservation and management policies on behalf of this species and its habitats.     

Species boundaries and phylogeography of the "Euscorpius carpathicus complex" (Scorpiones: Euscorpiidae) in Italy

Euscorpius tergestinus (C.L. Koch, 1837), Euscorpius concinnus (C.L. Koch, 1837) and Euscorpius sicanus (C.L. Koch, 1837), three presumed closely related species belonging to the "carpathicus group", occur in the Italian peninsula with a largely parapatric distribution and some zones of range overlap. These areas of sympatry represent interesting opportunities to investigate species boundaries in natural populations. Here we report on a study exploring genetic variation in sympatric populations of the three species from central Tuscany. Additional collecting sites, from different localities across Italy, were also included in the analysis in order to explore the phylogeographic structure of the group. Species boundaries and evolutionary relationships were examined by sequence comparison of mitochondrial 16S rRNA and nuclear ITS-1 rRNA gene fragments. DNA sequence data show no evidence of genetic introgression between different evolutionary lineages from the area of range overlap, suggesting the absence of either past or ongoing inter-specific gene flow. It is therefore probable that reproductive barriers exist, preventing gene pools from amalgamating. Furthermore, our results support the recent morphological distinction of E. tergestinus, as traditionally classified, into two different species: E. tergestinus and E. concinnus. Both mitochondrial and nuclear sequence data clearly indicate that the two taxa represent well-supported and deeply divergent lineages. Euscorpius sicanus seems to represent a monophyletic taxon, but the high genetic variability observed within this taxon calls for future investigation. The present distribution patterns across the Italian peninsula were mainly interpreted as the consequence of climatic oscillations.     

Species trees for the tree swallows (Genus Tachycineta): An alternative phylogenetic hypothesis to the mitochondrial gene tree

The New World swallow genus Tachycineta comprises nine species that collectively have a wide geographic distribution and remarkable variation both within- and among-species in ecologically important traits. Existing phylogenetic hypotheses for Tachycineta are based on mitochondrial DNA sequences, thus they provide estimates of a single gene tree. In this study we sequenced multiple individuals from each species at 16 nuclear intron loci. We used gene concatenated approaches (Bayesian and maximum likelihood) as well as coalescent-based species tree inference to reconstruct phylogenetic relationships of the genus. We examined the concordance and conflict between the nuclear and mitochondrial trees and between concatenated and coalescent-based inferences. Our results provide an alternative phylogenetic hypothesis to the existing mitochondrial DNA estimate of phylogeny. This new hypothesis provides a more accurate framework in which to explore trait evolution and examine the evolution of the mitochondrial genome in this group.     

Quaternary range and demographic expansion of Liolaemus darwinii (Squamata: Liolaemidae) in the Monte Desert of Central Argentina using Bayesian phylogeography and ecological niche modelling

Until recently, most phylogeographic approaches have been unable to distinguish between demographic and range expansion processes, making it difficult to test for the possibility of range expansion without population growth and vice versa. In this study, we applied a Bayesian phylogeographic approach to reconstruct both demographic and range expansion in the lizard Liolaemus darwinii of the Monte Desert in Central Argentina, during the Late Quaternary. Based on analysis of 14 anonymous nuclear loci and the cytochrome b mitochondrial DNA gene, we detected signals of demographic expansion starting at ~55 ka based on Bayesian Skyline and Skyride Plots. In contrast, Bayesian relaxed models of spatial diffusion suggested that range expansion occurred only between ~95 and 55 ka, and more recently, diffusion rates were very low during demographic expansion. The possibility of population growth without substantial range expansion could account for the shared patterns of demographic expansion during the Last Glacial Maxima (OIS 2 and 4) in fish, small mammals and other lizards of the Monte Desert. We found substantial variation in diffusion rates over time, and very high rates during the range expansion phase, consistent with a rapidly advancing expansion front towards the southeast shown by palaeo‐distribution models. Furthermore, the estimated diffusion rates are congruent with observed dispersal rates of lizards in field conditions and therefore provide additional confidence to the temporal scale of inferred phylogeographic patterns. Our study highlights how the integration of phylogeography with palaeo‐distribution models can shed light on both demographic and range expansion processes and their potential causes.     

Finding complexity in complexes: Assessing the causes of mitonuclear discordance in a problematic species complex of Mesoamerican toads

Mitonuclear discordance is a frequently encountered pattern in phylogeographic studies and occurs when mitochondrial and nuclear DNA display conflicting signals. Discordance among these genetic markers can be caused by several factors including confounded taxonomies, gene flow, and incomplete lineage sorting. In this study, we present a strong case of mitonuclear discordance in a species complex of toads (Bufonidae: Incilius coccifer complex) found in the Chortís Block of Central America. To determine the cause of mitonuclear discordance in this complex, we used spatially explicit genetic data to test species limits and relationships, characterize demographic history, and quantify gene flow. We found extensive mitonuclear discordance among the three recognized species within this group, especially in populations within the Chortís Highlands of Honduras. Our data reveal nuclear introgression within the Chortís Highlands populations that was most probably driven by cyclical range expansions due to climatic fluctuations. Though we determined introgression occurred within the nuclear genome, our data suggest that it is not the key factor in driving mitonuclear discordance in the entire species complex. Rather, due to a lack of discernible geographic pattern between mitochondrial and nuclear DNA, as well as a relatively recent divergence time of this complex, we concluded that mitonuclear discordance has been caused by incomplete lineage sorting. Our study provides a framework to test sources of mitonuclear discordance and highlights the importance of using multiple marker types to test species boundaries in cryptic species.     

The extent and meaning of hybridization and introgression between Siberian spruce (Picea obovata) and Norway spruce (Picea abies): cryptic refugia as stepping stones to the west?

Boreal species were repeatedly exposed to ice ages and went through cycles of contraction and expansion while sister species alternated periods of contact and isolation. The resulting genetic structure is consequently complex, and demographic inferences are intrinsically challenging. The range of Norway spruce (Picea abies) and Siberian spruce (Picea obovata) covers most of northern Eurasia; yet their geographical limits and histories remain poorly understood. To delineate the hybrid zone between the two species and reconstruct their joint demographic history, we analysed variation at nuclear SSR and mitochondrial DNA in 102 and 88 populations, respectively. The dynamics of the hybrid zone was analysed with approximate Bayesian computation (ABC) followed by posterior predictive structure plot reconstruction and the presence of barriers across the range tested with estimated effective migration surfaces. To estimate the divergence time between the two species, nuclear sequences from two well‐separated populations of each species were analysed with ABC. Two main barriers divide the range of the two species: one corresponds to the hybrid zone between them, and the other separates the southern and northern domains of Norway spruce. The hybrid zone is centred on the Urals, but the genetic impact of Siberian spruce extends further west. The joint distribution of mitochondrial and nuclear variation indicates an introgression of mitochondrial DNA from Norway spruce into Siberian spruce. Overall, our data reveal a demographic history where the two species interacted frequently and where migrants originating from the Urals and the West Siberian Plain recolonized northern Russia and Scandinavia using scattered refugial populations of Norway spruce as stepping stones towards the west.     

Testing phylogeographic predictions on an active volcanic island: Brachyderes rugatus (Coleoptera: Curculionidae) on La Palma (Canary Islands)

Volcanic islands with well-characterized geological histories can provide ideal templates for generating and testing phylogeographic predictions. Many studies have sought to utilize these to investigate patterns of colonization and speciation within groups of closely related species across a number of islands. Here we focus attention within a single volcanic island with a well-characterized geological history to develop and test phylogeographic predictions. We develop phylogeographic predictions within the island of La Palma of the Canary Islands and test these using 69 haplotypes from 570 base pairs of mitochondrial DNA cytochrome oxidase II sequence data for 138 individuals of Brachyderes rugatus rugatus, a local endemic subspecies of curculionid beetle occurring throughout the island in the forests of Pinus canariensis. Although geological data do provide some explanatory power for the phylogeographic patterns found, our network-based analyses reveal a more complicated phylogeographic history than initial predictions generated from data on the geological history of the island. Reciprocal illumination of geological and phylogeographic history is also demonstrated with previous geological speculation gaining phylogeographic corroboration from our analyses.