Sympatric divergence and clinal variation in multiple coloration traits of Ficedula flycatchers

Geographic variation in phenotypes plays a key role in fundamental evolutionary processes such as local adaptation, population differentiation and speciation, but the selective forces behind it are rarely known. We found support for the hypothesis that geographic variation in plumage traits of the pied flycatcher Ficedula hypoleuca is explained by character displacement with the collared flycatcher Ficedula albicollis in the contact zone. The plumage traits of the pied flycatcher differed strongly from the more conspicuous collared flycatcher in a sympatric area but increased in conspicuousness with increasing distance to there. Phenotypic differentiation (P_ST) was higher than that in neutral genetic markers (F_ST), and the effect of geographic distance remained when statistically controlling for neutral genetic differentiation. This suggests that a cline created by character displacement and gene flow explains phenotypic variation across the distribution of this species. The different plumage traits of the pied flycatcher are strongly to moderately correlated, indicating that they evolve non‐independently from each other. The flycatchers provide an example of plumage patterns diverging in two species that differ in several aspects of appearance. The divergence in sympatry and convergence in allopatry in these birds provide a possibility to study the evolutionary mechanisms behind the highly divergent avian plumage patterns.     

Speciation, introgressive hybridization and nonlinear rate of molecular evolution in flycatchers

Evolutionary history of Muscicapidae flycatchers is inferred from nuclear and mitochondrial DNA (mtDNA) sequence comparisons and population genetic analysis of nuclear and mtDNA markers. Phylogenetic reconstruction based on sequences from the two genomes yielded similar trees with respect to the order at which the species split off. However, the genetic distances fitted a nonlinear, polynomial model reflecting diminishing divergence rate of the mtDNA sequences compared to the nuclear DNA sequences. This could be explained by Haldane's rule because genetic isolation might evolve more rapidly on the mitochondrial rather than the nuclear genome in birds. This is because hybrid sterility of the heterogametic sex (females) would predate that of the homogametic sex (males), leading to sex biased introgression of nuclear genes. Analyses of present hybrid zones of pied (Ficedula hypoleuca) and collared flycatchers (F. albicollis) may indicate a slight sexual bias in rate of introgression, but the introgression rates were too low to allow proper statistical analyses. It is suggested, however, that the observed deviation from linearity can be explained by a more rapid mutational saturation of the mtDNA sequences than of the nuclear DNA sequences, as supported by analyses of third codon position transversions at two protein coding mtDNA genes. A phylogeographic scenario for the black and white flycatcher species is suggested based on interpretation of the genetic data obtained. Four species appear to have diverged from a common ancestor relatively simultaneously during the Pleistocene. After the last glaciation period, pied and collared flycatchers expanded their breeding ranges and eventually came into secondary contact in Central and Eastern Europe and on the Baltic Isles.     

Perpendicular axes of differentiation generated by mitochondrial introgression

Differential introgression of mitochondrial vs. nuclear DNA generates discordant patterns of geographic variation and can promote population divergence and speciation. We examined a potential case of mitochondrial introgression leading to two perpendicular axes of differentiation. The Eastern Yellow Robin Eopsaltria australis, a widespread Australian bird, shows a deep mitochondrial split that is perpendicular to north–south nuclear DNA and plumage colour differentiation. We propose a scenario to explain this pattern: (i) first, both nuclear and mitochondrial genomes differentiated in concert during north–south population divergence; (ii) later, their histories disconnected after two mitochondrial introgression events resulting in a deep mitochondrial split perpendicular to the nuclear DNA structure. We explored this scenario by coalescent modelling of ten mitochondrial genes and 400 nuclear DNA loci. Initial mitochondrial and nuclear genome divergences were estimated to have occurred in the early Pleistocene, consistent with the proposed scenario. Subsequent climatic transitions may have driven later mitochondrial introgression. We consider neutral introgression unlikely and instead propose that the evidence is more consistent with adaptive mitochondrial introgression and selection against incompatible mitochondrial‐nuclear combinations. This likely generated an axis of coastal‐inland mitochondrial differentiation in the face of nuclear gene flow, perpendicular to the initial north–south axis of differentiation (reflected in genomewide nuclear DNA and colour variation).     

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.     

Bidirectional mitochondrial introgression between Korean cobitid fish mediated by hybridogenetic hybrids

Genomic introgression through interspecific hybridization has been observed in some species of the freshwater fish family Cobitidae. Within this family, a Cobitis hankugensis–Iksookimia longicorpa diploid–triploid hybrid species complex on the Korean peninsula is unique in displaying hybridogenesis, a unisexual reproduction mode that allows hybrids to mediate the transfer of mitochondrial DNA (but not nuclear DNA) between the two parent species. However, populations of the parental species in the wild have never been examined for the potential effect of introgression on their genomes. To address the genetic consequences of unisexual hybridization on the parental species, we examined genetic structure of the two parental species, C. hankugensis and I. longicorpa, in three independent natural habitats where they coexist with their hybrid complex using DNA sequence data of one mitochondrial gene and three nuclear genes. We found that mitochondrial introgression between the two species was extensive in all the examined localities, while there was no evidence of nuclear introgression across the species boundary. This result indicates that the hybridogenetic individuals mediate mitochondrial introgression from one species to the other, producing mito‐nuclear mosaic genomes such as C. hankugensis nuclear genomes associated with I. longicorpa mitochondrial DNA and the reverse. The direction and degree of introgression varied among the three localities, but the underlying mechanisms for this observation proved elusive. Introgression might depend on which species serves as the predominant sperm or ovum donor or the environmental conditions of the localities. The present study suggests that introgressive hybridization between pure C. hankugensis and I. longicorpa species is highly likely where the two species co‐occur with hybridogenetic individuals, but the consequence of introgression could be variable due to the history and environmental characteristics of particular populations across the parental species’ ranges.     

Heterospecific female mimicry in Ficedula flycatchers

Mimicry is a widespread phenomenon. Vertebrate visual mimicry often operates in an intraspecific sexual context, with some males resembling conspecific females. Pied flycatcher (Ficedula hypoleuca) dorsal plumage varies from the ancestral black to female‐like brown. Experimental studies have shown that conspecific and heterospecific (collared flycatcher, F. albicollis) individuals of both sexes respond, at least initially, to brown individuals as if they were female. We quantified the perceptual and biochemical differences between brown feathers and found that brown pied flycatcher males are indistinguishable from heterospecific, but not from conspecific, females in both aspects. To our knowledge, this is the first evidence of a visual mimetic signalling system in a sexual context where the model is heterospecific to the mimic. By only mimicking heterospecific females, brown pied flycatcher males can establish territories next to the more dominant collared flycatcher in sympatry, suffer less aggression by darker conspecifics in allopatry and preserve within‐species sexual recognition throughout the breeding range. A closer look at the evolutionary history and ecology of these two species illustrates how such a mimetic system can evolve. Although likely rare, this phenomenon might not be unique to Ficedula flycatchers.     

Deep sympatric mitochondrial divergence without reproductive isolation in the common redstart Phoenicurus phoenicurus

Mitochondrial DNA usually shows low sequence variation within and high sequence divergence among species, which makes it a useful marker for phylogenetic inference and DNA barcoding. A previous study on the common redstart (Phoenicurus phoenicurus) revealed two very different mtDNA haplogroups (5% K2P distance). This divergence is comparable to that among many sister species; however, both haplogroups coexist and interbreed in Europe today. Herein, we describe the phylogeographic pattern of these lineages and test hypotheses for how such high diversity in mtDNA has evolved. We found no evidence for mitochondrial pseudogenes confirming that both haplotypes are of mitochondrial origin. When testing for possible reproductive barriers, we found no evidence for lineage‐specific assortative mating and no difference in sperm morphology, indicating that they are not examples of cryptic species, nor likely to reflect the early stages of speciation. A gene tree based on a short fragment of cytochrome c oxidase subunit 1 from the common redstart and 10 other Phoenicurus species, showed no introgression from any of the extant congenerics. However, introgression from an extinct congeneric cannot be excluded. Sequences from two nuclear introns did not show a similar differentiation into two distinct groups. Mismatch distributions indicated that the lineages have undergone similar demographic changes. Taken together, these results confirm that deeply divergent mitochondrial lineages can coexist in biological species. Sympatric mtDNA divergences are relatively rare in birds, but the fact that they occur argues against the use of threshold mtDNA divergences in species delineation.     

Whole‐genome patterns of linkage disequilibrium across flycatcher populations clarify the causes and consequences of fine‐scale recombination rate variation in birds

Recombination rate is heterogeneous across the genome of various species and so are genetic diversity and differentiation as a consequence of linked selection. However, we still lack a clear picture of the underlying mechanisms for regulating recombination. Here we estimated fine‐scale population recombination rate based on the patterns of linkage disequilibrium across the genomes of multiple populations of two closely related flycatcher species (Ficedula albicollis and F. hypoleuca). This revealed an overall conservation of the recombination landscape between these species at the scale of 200 kb, but we also identified differences in the local rate of recombination despite their recent divergence (<1 million years). Genetic diversity and differentiation were associated with recombination rate in a lineage‐specific manner, indicating differences in the extent of linked selection between species. We detected 400–3,085 recombination hotspots per population. Location of hotspots was conserved between species, but the intensity of hotspot activity varied between species. Recombination hotspots were primarily associated with CpG islands (CGIs), regardless of whether CGIs were at promoter regions or away from genes. Recombination hotspots were also associated with specific transposable elements (TEs), but this association appears indirect due to shared preferences of the transposition machinery and the recombination machinery for accessible open chromatin regions. Our results suggest that CGIs are a major determinant of the localization of recombination hotspots, and we propose that both the distribution of TEs and fine‐scale variation in recombination rate may be associated with the evolution of the epigenetic landscape.     

Chlorophyll-binding proteins revisited - a multigenic family of light-harvesting and stress proteins from a brown algal perspective

Background

Mitochondrial introgression may result in the mitochondrial genome of one species being replaced by that of another species without leaving any trace of past hybridization in its nuclear genome. Such introgression can confuse the species genealogy estimates and lead to absurd inferences of species history. We used a phylogenetic approach to explore the potential mitochondrial genome introgression event(s) between two closely related green pond frog species, Pelophylax nigromaculatus and P. plancyi.

Results

DNA sequence data of one mitochondrial and two nuclear genes from an extensive sampling of the two species were collected, and the genealogies of the three genes were constructed and compared. While the two nuclear genes congruently showed mutual reciprocal monophyly of both species, the mitochondrial phylogeny separated a Korean P. nigromaculatus clade, a paraphyletic central China P. plancyi assemblage, and a large well-supported introgression clade. Within the introgression clade, the mitochondrial haplotypes of the two species were mixed together. This reticulated pattern can be most parsimoniously explained by an ancient mitochondrial introgression event from P. plancyi to P. nigromaculatus that occurred at least 1.36 MYA, followed by multiple recent introgression events from P. nigromaculatus back to P. plancyi within the last 0.63 MY. The re-constitution of previously co-adapted genomes in P. plancyi may be responsible for the recent rampant introgression events. The Korean P. nigromaculatus clade likely represents the only surviving "true" mitochondrial lineage of P. nigromaculatus, and the central China P. plancyi assemblage likely represents the "original" P. plancyi mitochondrial lineage. Refugia in the Korean Peninsula and central China may have played a significant role in preserving these ancient lineages.

Conclusions

The majority of individuals in the two species have either introgressed (P. nigromaculatus) or reclaimed (P. plancyi) mitochondrial genomes while no trace of past hybridization in their nuclear genomes was detected. Asymmetrical reproductive ability of hybrids and continuous backcrossing are likely responsible for the observed mitochondrial introgression. This case is unique in that it includes an ancient "forward" introgression and many recent "backward" introgressions, which re-constitutes the original nuclear and mitochondrial genomes of P. plancyi. This hybrid system provides an excellent opportunity to study cyto-nuclear interaction and co-adaptation.     

GENOME DIVERGENCE AND THE GENETIC ARCHITECTURE OF BARRIERS TO GENE FLOW BETWEEN LYCAEIDES IDAS AND L. MELISSA

Genome divergence during speciation is a dynamic process that is affected by various factors, including the genetic architecture of barriers to gene flow. Herein we quantitatively describe aspects of the genetic architecture of two sets of traits, male genitalic morphology and oviposition preference, that putatively function as barriers to gene flow between the butterfly species Lycaeides idas and L. melissa. Our analyses are based on unmapped DNA sequence data and a recently developed Bayesian regression approach that includes variable selection and explicit parameters for the genetic architecture of traits. A modest number of nucleotide polymorphisms explained a small to large proportion of the variation in each trait, and average genetic variant effects were nonnegligible. Several genetic regions were associated with variation in multiple traits or with trait variation within‐ and among‐populations. In some instances, genetic regions associated with trait variation also exhibited exceptional genetic differentiation between species or exceptional introgression in hybrids. These results are consistent with the hypothesis that divergent selection on male genitalia has contributed to heterogeneous genetic differentiation, and that both sets of traits affect fitness in hybrids. Although these results are encouraging, we highlight several difficulties related to understanding the genetics of speciation.     

SIGNIFICANT ROLE FOR HISTORICAL EFFECTS IN THE EVOLUTION OF REPRODUCTIVE ISOLATION: EVIDENCE FROM PATTERNS OF INTROGRESSION BETWEEN THE CYPRINID FISHES,LUXILUS CORNUTUS AND LUXILUS CHRYSOCEPHALUS

Samples of Luxilus cornutus, Luxilus chrysocephalus, and their hybrids were collected along hypothesized routes of dispersal from Pleistocene refugia to examine the significance of geographic variation in patterns of introgression between these species. Patterns of allozyme and mitochondrial DNA (mtDNA) variation were generally consistent with those from previous studies. Tests of Hardy-Weinberg equilibrium revealed significant deficiencies of heterozygotes in all samples, indicating some form of reproductive isolation. Mitochondrial DNAs of each species were not equally represented in F₁ hybrids; however, this bias was eliminated when the two largest samples were excluded from the analysis. Backcross hybrids exhibited biased mtDNA introgression, as samples from Lake Erie (eastern) and Lake Michigan (western) drainages showed significant excesses of mtDNAs from L. chrysocephalus and L. cornutus, respectively, relative to frequencies of diagnostic allozyme markers. The extent and direction of allozyme and mtDNA introgression was quantified by calculating isolation index values from morphologically “pure” individuals of each species from each locality. Analysis of variance of these measures identified limited introgression of allozyme variants with no geographic pattern, but significant differences in direction of mtDNA introgression between drainages (i.e., postglacial dispersal route). Association between patterns of mtDNA introgression and dispersal route across the latitudinal width of the contact zone is best explained by genetic divergence during past isolation of ancestral populations from these drainages. These results identify a significant role for historical effects in the evolution of reproductive isolation and the process of speciation.     

Speciation in the Rana chensinensis species complex and its relationship to the uplift of the Qinghai-Tibetan Plateau

Speciation remains a fundamental issue in biology. Herein, we report an investigation into speciation in the Rana chensinensis species complex using DNA sequence data from one mitochondrial and five nuclear genes. A phylogenetic analysis of the data revealed four major clades in the complex, and each of them was found to likely represent a species, including one cryptic species. Ecological niche models were generated from 19 climatic variables for three of the four major clades, which were represented by widespread sampling, including R. chensinensis, Rana kukunoris and the potential cryptic species. Each clade is associated with a unique ecological unit, and this indicates that ecological divergence probably drove speciation. Ecological divergence is likely related to the late Cenozoic orogenesis of the Qinghai–Tibetan Plateau. In addition, gene flow between species was detected but only in peripheral portions of the ranges of the four major clades, thus likely had little influence on the speciation processes. Discordances between mitochondrial and nuclear genes were also found; the nominal species, R. chensinensis, contains multiple maternal clades, suggesting potential mitochondrial introgression between R. chensinensis and R. kukunoris.     

Evidence for strong assortative mating,limited gene flow,and strong differentiation across the blue‐footed/Peruvian booby hybrid zone in northern Peru

Hybrid zones represent natural laboratories in which the processes of divergence and genetic isolation can be examined. The generation and maintenance of a hybrid zone requires mispairing and successful reproduction between organisms that differ in one or more heritable traits. Understanding the dynamics of hybridization between two species requires an understanding of the extent to which they have diverged genetically, the frequency of misparing and hybrid production, and the extent of introgression. Three hundred and twenty one blue‐footed Sula nebouxii and Peruvian S. variegata boobies from the eastern tropical Pacific Ocean were analyzed using 19 putatively neutral genetic markers to evaluate interspecific differentiation, to classify morphological hybrids using Bayesian assignments, and to characterize hybridization using cline theory and Bayesian assignments. The species were well differentiated at mitochondrial and nuclear microsatellites, the hybrid zone was bimodal (contained a high frequency of each parental species but a low frequency of hybrids), and morphologically intermediate individuals were most likely F1 hybrids resulting from mating between female Peruvian boobies and male blue‐footed boobies. Clines in allele frequency could be constrained to share a common geographic centre but could not be constrained to share a common width. Peruvian and blue‐footed boobies hybridize infrequently, potentially due to strong premating reproductive isolation; however, backcrossing appears to facilitate introgression from blue‐footed to Peruvian boobies in this hybrid system.     

Insight into the roles of selection in speciation from genomic patterns of divergence and introgression in secondary contact in venomous rattlesnakes

Investigating secondary contact of historically isolated lineages can provide insight into how selection and drift influence genomic divergence and admixture. Here, we studied the genomic landscape of divergence and introgression following secondary contact between lineages of the Western Diamondback Rattlesnake (Crotalus atrox) to determine whether genomic regions under selection in allopatry also contribute to reproductive isolation during introgression. We used thousands of nuclear loci to study genomic differentiation between two lineages that have experienced recent secondary contact following isolation, and incorporated sampling from a zone of secondary contact to identify loci that are resistant to gene flow in hybrids. Comparisons of patterns of divergence and introgression revealed a positive relationship between allelic differentiation and resistance to introgression across the genome, and greater‐than‐expected overlap between genes linked to lineage‐specific divergence and loci that resist introgression. Genes linked to putatively selected markers were related to prominent aspects of rattlesnake biology that differ between populations of Western Diamondback rattlesnakes (i.e., venom and reproductive phenotypes). We also found evidence for selection against introgression of genes that may contribute to cytonuclear incompatibility, consistent with previously observed biased patterns of nuclear and mitochondrial alleles suggestive of partial reproductive isolation due to cytonuclear incompatibilities. Our results provide a genome‐scale perspective on the relationships between divergence and introgression in secondary contact that is relevant for understanding the roles of selection in maintaining partial isolation of lineages, causing admixing lineages to not completely homogenize.     

The role of western Mediterranean islands in the evolutionary diversification of the spotted flycatcher Muscicapa striata,a long‐distance migratory passerine species

We investigated the evolutionary history of the spotted flycatcher Muscicapa striata, a long distance migratory passerine having a widespread range, using mitochondrial markers and nuclear introns. Our mitochondrial results reveal the existence of one insular lineage restricted to the western Mediterranean islands (Balearics, Corsica, Sardinia) and possibly to the Tyrrhenian coast of Italy that diverged from the mainland lineages around 1 Mya. Mitochondrial genetic distance between insular and mainland lineages is around 3.5%. Limited levels of shared nuclear alleles among insular and mainland populations further support the genetic distinctiveness of insular spotted flycatchers with respect to their mainland counterparts. Moreover, lack of mitochondrial haplotypes sharing between Balearic birds (M. s. balearica) and Corso‐Sardinian birds (M. s. tyrrhenica) suggest the absence of recent matrilineal gene flow between these two insular subspecies. Accordingly, we suggest that insular spotted flycatchers could be treated as one polytypic species (Muscicapa tyrrhenica) that differs from M. striata in morphology, migration, mitochondrial and nuclear DNA and comprises two subspecies (the nominate and M. t. balearica) that diverged recently phenotypically and in mitochondrial DNA and but still share the same nuclear alleles. This study provides an interesting case‐study illustrating the crucial role of western Mediterranean islands in the evolution of a passerine showing high dispersal capabilities. Our genetic results highlight the role of glacial refugia of these islands that allowed initial allopatric divergence of insular populations. We hypothesize that differences in migratory and breeding phenology may prevent any current gene flow between insular and mainland populations of the spotted flycatcher that temporarily share the same insular habitats during the spring migration.     

Genetic differentiation and hybridization between greater and lesser spotted eagles (Accipitriformes:Aquila clanga, A. pomarina)

Greater and lesser spotted eagles (Aquila clanga, A. pomarina) are two closely related forest eagles overlapping in breeding range in east-central Europe. In recent years a number of mixed pairs have been observed, some of which fledged hybrid young. Here we use mitochondrial (control region) DNA sequences and AFLP markers to estimate genetic differentiation and possible gene flow between these species. In a sample of 83 individuals (61 pomarina, 20 clanga, 2 F1-hybrids) we found 30 mitochondrial haplotypes which, in a phylogenetic network, formed two distinct clusters differing on average by 3.0% sequence divergence. The two species were significantly differentiated both at the mitochondrial and nuclear (AFLP) genetic level. However, five individuals with pomarina phenotype possessed clanga-type mtDNA, suggesting occasional gene flow. Surprisingly, AFLP markers indicated that these mismatched birds (originating from Germany, E Poland and Latvia) were genetically intermediate between the samples of individuals in which mtDNA haplotype and phenotype agreed. This indicates that mismatched birds were either F1 or recent back-cross hybrids. Mitochondrial introgression was asymmetrical (no pomarina haplotype found in clanga so far), which may be due to assortative mating by size. Gene flow of nuclear markers was estimated to be about ten times stronger than for mtDNA, indicating a sex-bias in hybrid fertility in accordance with Haldanes rule. Hybridization between the two species may be more frequent and may occur much further west than hitherto assumed. This is supported by the recent discovery of a mixed pair producing at least one fledgling in NE Germany.     

Testing hypotheses of mitochondrial gene‐tree paraphyly: unravelling mitochondrial capture of the Streak‐breasted Scimitar Babbler (Pomatorhinus ruficollis) by the Taiwan Scimitar Babbler (Pomatorhinus musicus)

Species‐level paraphyly inferred from mitochondrial gene trees is a prevalent phenomenon in taxonomy and systematics, but there are several potential causes that are not easily explained by currently used methods. This study aimed to test the underlying causes behind the observed paraphyly of Streak‐breasted Scimitar Babbler (Pomatorhinus ruficollis) via statistical analyses of four mitochondrial (mtDNA) and nine nuclear (nuDNA) genes. Mitochondrial gene trees show paraphyly of P. ruficollis with respect to the Taiwan Scimitar Babbler (Pomatorhinus musicus), but nuclear genealogies support a sister‐group relationship. Predictive coalescent simulations imply several hypothetical explanations, the most likely being mitochondrial capture of P. ruficollis by P. musicus for the observed cyto‐nuclear incongruence. Further approximate Bayesian computation suggests a unidirectional introgression model with substantial level of gene flow from P. ruficollis to P. musicus during their initial divergence during the Late Pleistocene. This specific observation frames several potential causes for incongruent outcomes of mitochondrial and nuclear introgression in general, and on the whole, our results underscore the strength of multiple independent loci for species delimitation and importance of testing hypotheses that explain disparate causes of mitochondrial gene‐tree paraphyly.     

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.