Isolation and introgression in the Intermountain West: contrasting gene genealogies reveal the complex biogeographic history of the American pika (Ochotona princeps)

Aim We studied the history of colonization, diversification and introgression among major phylogroups in the American pika, Ochotona princeps (Lagomorpha), using comparative and statistical phylogeographic methods. Our goal was to understand how Pleistocene climatic fluctuations have shaped the distribution of diversity at mitochondrial DNA (mtDNA) and nuclear DNA (nDNA) loci in this alpine specialist. Location North America’s Intermountain West. Methods We accumulated mtDNA sequence data (c. 560–1700 bp) from 232 pikas representing 64 localities, and sequenced two nuclear introns (mast cell growth factor, c. 550 bp, n = 148; protein kinase C iota, c. 660 bp, n = 139) from a subset of individuals. To determine the distribution of major mtDNA lineages, we conducted a phylogenetic analysis on the mtDNA sequence data, and we calculated divergence times among the lineages using a Bayesian Markov chain Monte Carlo approach. Relationships among nuclear alleles were explored with minimum spanning networks. Finally, we conducted coalescent simulations of alternative models of population history to test for congruence between nDNA and mtDNA responses to Pleistocene glacial cycles. Results We found that: (1) all individuals could be assigned to one of five allopatric mtDNA lineages; (2) lineages are associated with separate mountain provinces; (3) lineages originated from at least two rounds of differentiation; (4) nDNA and mtDNA markers exhibited overall phylogeographic congruence; and (5) introgression among phylogroups has occurred at nuclear loci since their initial isolation. Main conclusions Pika populations associated with different mountain systems have followed separate but not completely independent evolutionary trajectories through multiple glacial cycles. Range expansion associated with climate cooling (i.e. glaciations) promoted genetic admixture among populations within mountain ranges. It also permitted periodic contact and introgression between phylogroups associated with different mountain systems, the record of which is retained at nDNA but not mtDNA loci. Evidence for different histories at nuclear and mtDNA loci (i.e. periodic introgression versus deep isolation, respectively) emphasizes the importance of multilocus perspectives for reconstructing complete population histories.     

A mirage of cryptic species: Genomics uncover striking mitonuclear discordance in the butterfly Thymelicus sylvestris

Mitochondrial DNA (mtDNA) sequencing has led to an unprecedented rise in the identification of cryptic species. However, it is widely acknowledged that nuclear DNA (nuDNA) sequence data are also necessary to properly define species boundaries. Next generation sequencing techniques provide a wealth of nuclear genomic data, which can be used to ascertain both the evolutionary history and taxonomic status of putative cryptic species. Here, we focus on the intriguing case of the butterfly Thymelicus sylvestris (Lepidoptera: Hesperiidae). We identified six deeply diverged mitochondrial lineages; three distributed all across Europe and found in sympatry, suggesting a potential case of cryptic species. We then sequenced these six lineages using double‐digest restriction‐site associated DNA sequencing (ddRADseq). Nuclear genomic loci contradicted mtDNA patterns and genotypes generally clustered according to geography, i.e., a pattern expected under the assumption of postglacial recolonization from different refugia. Further analyses indicated that this strong mtDNA/nuDNA discrepancy cannot be explained by incomplete lineage sorting, sex‐biased asymmetries, NUMTs, natural selection, introgression or Wolbachia‐mediated genetic sweeps. We suggest that this mitonuclear discordance was caused by long periods of geographic isolation followed by range expansions, homogenizing the nuclear but not the mitochondrial genome. These results highlight T. sylvestris as a potential case of multiple despeciation and/or lineage fusion events. We finally argue, since mtDNA and nuDNA do not necessarily follow the same mechanisms of evolution, their respective evolutionary history reflects complementary aspects of past demographic and biogeographic events.     

Cryptic diversity in the pen shell Atrina pectinata (Bivalvia: Pinnidae): high divergence and hybridization revealed by molecular and morphological data

Cryptic species have been increasingly revealed in the marine realm through an analytical approach incorporating multiple lines of evidence (e.g., mtDNA, nuclear genes and morphology). Illustrations of cryptic taxa improve our understanding of species diversity and evolutionary histories within marine animals. The pen shell Atrina pectinata is known to exhibit extensive morphological variations that may harbour cryptic diversity. In this study, we investigated A. pectinata populations along the coast of China and one from Japan to explore possible cryptic diversity and hybridization using a combination of mitochondrial (cytochrome c oxidase subunit I, mtCOI) and nuclear (ribosomal internal transcribed spacer, nrITS) genes as well as morphology. Phylogenetic analyses of mtCOI ‘DNA barcoding gene’ sequences resolved six divergent lineages with intralineage divergences between 0.4% and 0.8%. Interlineage sequence differences ranged from 4.3% to 22.0%, suggesting that six candidate cryptic species are present. The nrITS gene revealed five deep lineages with Kimura 2‐parameter distances of 3.7–30.3%. The five nuclear lineages generally corresponded to mtCOI lineages 1–4 and (5 + 6), suggestive of five distinct evolutionary lineages. Multiple nrITS sequences of significant variance were found within an individual, clearly implying recent hybridization events between/among the evolutionary lineages, which contributed to cytonuclear discordance. Morphologically, five morphotypes matched the five genetic lineages, although the intermediates may well blur the boundaries of different morphotypes. This study demonstrates the importance of combining multiple lines of evidence to explore species cryptic diversity and past evolutionary histories.     

Subterranean termite phylogeography reveals multiple postglacial colonization events in southwestern Europe

A long‐standing goal of evolutionary biology is to understand how paleoclimatic and geological events shape the geographical distribution and genetic structure within and among species. Using a diverse set of markers (cuticular hydrocarbons, mitochondrial and nuclear gene sequences, microsatellite loci), we studied Reticulitermes grassei and R. banyulensis, two closely related termite species in southwestern Europe. We sought to clarify the current genetic structure of populations that formed following postglacial dispersal from refugia in southern Spain and characterize the gene flow between the two lineages over the last several million years. Each marker type separately provided a fragmented picture of the evolutionary history at different timescales. Chemical analyses of cuticular hydrocarbons and phylogenetic analyses of mitochondrial and nuclear genes showed clear separation between the species, suggesting they diverged following vicariance events in the Late Miocene. However, the presence of intermediate chemical profiles and mtDNA introgression in some Spanish colonies suggests ongoing gene flow. The current genetic structure of Iberian populations is consistent with alternating isolation and dispersal events during Quaternary glacial periods. Analyses of population genetic structure revealed postglacial colonization routes from southern Spain to France, where populations underwent strong genetic bottlenecks after traversing the Pyrenees resulting in parapatric speciation.     

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.     

Introgression of mitochondrial DNA among lineages in a hybridogenetic ant

We report a remarkable pattern of incongruence between nuclear and mitochondrial variations in a social insect, the desert ant Cataglyphis hispanica. This species reproduces by social hybridogenesis. In all populations, two distinct genetic lineages coexist; non-reproductive workers develop from hybrid crosses between the lineages, whereas reproductive offspring (males and new queens) are typically produced asexually by parthenogenesis. Genetic analyses based on nuclear markers revealed that the two lineages remain highly differentiated despite constant hybridization for worker production. Here, we show that, in contrast with nuclear DNA, mitochondrial DNA (mtDNA) does not recover the two lineages as monophyletic. Rather, mitochondrial haplotypes cluster according to their geographical origin. We argue that this cytonuclear incongruence stems from introgression of mtDNA among lineages, and review the mechanisms likely to explain this pattern under social hybridogenesis.     

Conflicting patterns of mitochondrial and nuclear DNA diversity in Phylloscopus warblers

Molecular variation is often used to infer the demographic history of species, but sometimes the complexity of species history can make such inference difficult. The willow warbler, Phylloscopus trochilus, shows substantially less geographical variation than the chiffchaff, Phylloscopus collybita, both in morphology and in mitochondrial DNA (mtDNA) divergence. We therefore predicted that the willow warbler should harbour less nuclear DNA diversity than the chiffchaff. We analysed sequence data obtained from multiple samples of willow warblers and chiffchaffs for the mtDNA cytochrome b gene and four nuclear genes. We confirmed that the mtDNA diversity among willow warblers is low (pi = 0.0021). Sequence data from three nuclear genes (CHD-Z, AFLP-WW1 and MC1R) not linked to the mitochondria demonstrated unexpectedly high nucleotide diversity (pi values of 0.0172, 0.0141 and 0.0038) in the willow warbler, on average higher than the nucleotide diversity for the chiffchaff (pi values of 0.0025, 0.0017 and 0.0139). In willow warblers, Tajima's D analyses showed that the mtDNA diversity, but not the nuclear DNA diversity, has been reduced relative to the neutral expectation of molecular evolution, suggesting the action of a selective sweep affecting the maternally inherited genes. The large nuclear diversity seen within willow warblers is not compatible with processes of neutral evolution occurring in a population with a constant population size, unless the long-term effective population size has been very large (N(e) > 10(6)). We suggest that the contrasting patterns of genetic diversity in the willow warbler may reflect a more complex evolutionary history, possibly including historical demographic fluctuations or historical male-biased introgression of nuclear genes from a differentiated population of Phylloscopus warblers.     

The on‐again,off‐again relationship between mitochondrial genomes and species boundaries

The study of reproductive isolation and species barriers frequently focuses on mitochondrial genomes and has produced two alternative and almost diametrically opposed narratives. On one hand, mtDNA may be at the forefront of speciation events, with co‐evolved mitonuclear interactions responsible for some of the earliest genetic incompatibilities arising among isolated populations. On the other hand, there are numerous cases of introgression of mtDNA across species boundaries even when nuclear gene flow is restricted. We argue that these seemingly contradictory patterns can result from a single underlying cause. Specifically, the accumulation of deleterious mutations in mtDNA creates a problem with two alternative evolutionary solutions. In some cases, compensatory or epistatic changes in the nuclear genome may ameliorate the effects of mitochondrial mutations, thereby establishing coadapted mitonuclear genotypes within populations and forming the basis of reproductive incompatibilities between populations. Alternatively, populations with high mitochondrial mutation loads may be rescued by replacement with a more fit, foreign mitochondrial haplotype. Coupled with many nonadaptive mechanisms of introgression that can preferentially affect cytoplasmic genomes, this form of adaptive introgression may contribute to the widespread discordance between mitochondrial and nuclear genealogies. Here, we review recent advances related to mitochondrial introgression and mitonuclear incompatibilities, including the potential for cointrogression of mtDNA and interacting nuclear genes. We also address an emerging controversy over the classic assumption that selection on mitochondrial genomes is inefficient and discuss the mechanisms that lead lineages down alternative evolutionary paths in response to mitochondrial mutation accumulation.     

GENETIC ARCHITECTURE OF METABOLIC RATE: ENVIRONMENT SPECIFIC EPISTASIS BETWEEN MITOCHONDRIAL AND NUCLEAR GENES IN AN INSECT

The extent to which mitochondrial DNA (mtDNA) variation is involved in adaptive evolutionary change is currently being reevaluated. In particular, emerging evidence suggests that mtDNA genes coevolve with the nuclear genes with which they interact to form the energy producing enzyme complexes in the mitochondria. This suggests that intergenomic epistasis between mitochondrial and nuclear genes may affect whole‐organism metabolic phenotypes. Here, we use crossed combinations of mitochondrial and nuclear lineages of the seed beetle Callosobruchus maculatus and assay metabolic rate under two different temperature regimes. Metabolic rate was affected by an interaction between the mitochondrial and nuclear lineages and the temperature regime. Sequence data suggests that mitochondrial genetic variation has a role in determining the outcome of this interaction. Our genetic dissection of metabolic rate reveals a high level of complexity, encompassing genetic interactions over two genomes, and genotype × genotype × environment interactions. The evolutionary implications of these results are twofold. First, because metabolic rate is at the root of life histories, our results provide insights into the complexity of life‐history evolution in general, and thermal adaptation in particular. Second, our results suggest a mechanism that could contribute to the maintenance of nonneutral mtDNA polymorphism.     

ELUCIDATION OF CRYPTIC DIVERSITY IN A WIDESPREAD NEARCTIC TREEFROG REVEALS EPISODES OF MITOCHONDRIAL GENE CAPTURE AS FROGS DIVERSIFIED ACROSS A DYNAMIC LANDSCAPE

We investigate the evolutionary history of the wide‐ranging Nearctic treefrog Hyla arenicolor through the integration of extensive range‐wide sampling, phylogenetic analyses of multilocus genetic data, and divergence dating. Previous phylogeographic studies of this frog documented a potential signature of introgressive hybridization from an ecologically and morphologically divergent sister species. Based on our Bayesian phylogenetic analyses of mitochondrial DNA, we inferred strong phylogeographic structure in H. arenicolor as indicated by seven well‐supported clades, five of which correspond to well‐defined biogeographic regions. Clades from the Balsas Basin and southwestern Central Mexican Plateau in Mexico, and the Grand Canyon of Arizona, group with the morphologically, behaviorally, and ecologically divergent mountain treefrogs in the H. eximia group, rendering H. arenicolor as paraphyletic. The phylogenetic position of at least two of these three H. arenicolor clades within the H. eximia group, however, is most likely the result of several episodes of introgressive hybridization and subsequent mitochondrial gene capture separated in time and space, as supported by evidence from the nuclear genes. Hyla arenicolor from the Balsas Basin appear to be deeply divergent from other H. arenicolor and represent a distinctly different species. Results suggests that introgressive hybridization events, both ancient and contemporary, coupled with late Neogene vicariance and Pleistocene climate‐driven range shifts, have all played a role in the historical diversification of H. arenicolor.     

Phylogeny and temporal diversification of darters (Percidae: Etheostomatinae)

Discussions aimed at resolution of the Tree of Life are most often focused on the interrelationships of major organismal lineages. In this study, we focus on the resolution of some of the most apical branches in the Tree of Life through exploration of the phylogenetic relationships of darters, a species-rich clade of North American freshwater fishes. With a near-complete taxon sampling of close to 250 species, we aim to investigate strategies for efficient multilocus data sampling and the estimation of divergence times using relaxed-clock methods when a clade lacks a fossil record. Our phylogenetic data set comprises a single mitochondrial DNA (mtDNA) gene and two nuclear genes sampled from 245 of the 248 darter species. This dense sampling allows us to determine if a modest amount of nuclear DNA sequence data can resolve relationships among closely related animal species. Darters lack a fossil record to provide age calibration priors in relaxed-clock analyses. Therefore, we use a near-complete species-sampled phylogeny of the perciform clade Centrarchidae, which has a rich fossil record, to assess two distinct strategies of external calibration in relaxed-clock divergence time estimates of darters: using ages inferred from the fossil record and molecular evolutionary rate estimates. Comparison of Bayesian phylogenies inferred from mtDNA and nuclear genes reveals that heterospecific mtDNA is present in approximately 12.5% of all darter species. We identify three patterns of mtDNA introgression in darters: proximal mtDNA transfer, which involves the transfer of mtDNA among extant and sympatric darter species, indeterminate introgression, which involves the transfer of mtDNA from a lineage that cannot be confidently identified because the introgressed haplotypes are not clearly referable to mtDNA haplotypes in any recognized species, and deep introgression, which is characterized by species diversification within a recipient clade subsequent to the transfer of heterospecific mtDNA. The results of our analyses indicate that DNA sequences sampled from single-copy nuclear genes can provide appreciable phylogenetic resolution for closely related animal species. A well-resolved near-complete species-sampled phylogeny of darters was estimated with Bayesian methods using a concatenated mtDNA and nuclear gene data set with all identified heterospecific mtDNA haplotypes treated as missing data. The relaxed-clock analyses resulted in very similar posterior age estimates across the three sampled genes and methods of calibration and therefore offer a viable strategy for estimating divergence times for clades that lack a fossil record. In addition, an informative rank-free clade-based classification of darters that preserves the rich history of nomenclature in the group and provides formal taxonomic communication of darter clades was constructed using the mtDNA and nuclear gene phylogeny. On the whole, the appeal of mtDNA for phylogeny inference among closely related animal species is diminished by the observations of extensive mtDNA introgression and by finding appreciable phylogenetic signal in a modest sampling of nuclear genes in our phylogenetic analyses of darters.     

Incongruence among mitochondrial,chloroplast and nuclear gene trees in <Emphasis Type="Italic">Pinus</Emphasis> subgenus <Emphasis Type="Italic">Strobus</Emphasis> (Pinaceae)

Introgression has been considered to be one of main factors leading to phylogenetic incongruence among different datasets at lower taxonomic levels. In the plants of Pinaceae, the mtDNA, cpDNA, and nuclear DNA (nrDNA) may have different evolutionary histories through introgression because they are inherited maternally, paternally and biparentally, respectively. We compared mtDNA, cpDNA, and two low-copy nrDNA phylogenetic trees in the genus Pinus subgenus Strobus, in order to detect unknown past introgression events in this group. nrDNA trees were mostly congruent with the cpDNA tree, and supported the recent sectional and subsectional classification system. In contrast, mtDNA trees split the members of sect. Quinquefoliae into two groups that were not observed in the other gene trees. The factors constituting incongruence may be divided into the following two categories: the different splits within subsect. Strobus, and the non-monophyly of subsect. Gerardianae. The former was hypothesized to have been caused by the past introgression of cpDNA, mtDNA or both between Eurasian and North American species through Beringia. The latter was likely caused by the chimeric structure of the mtDNA sequence of P. bungeana, which might have originated through past hybridization, or through a horizontal transfer event and subsequent recombination. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Morphology and nuclear markers reveal extensive mitochondrial introgressions in the Iberian Wall Lizard species complex

Mitochondrial markers are still often used alone to identify evolutionary units, despite widespread evidence for processes such as incomplete lineage sorting or introgressive hybridization that may blur past population history. The combination of mitochondrial DNA data with other sources of information (morphology, nuclear genes) is a powerful tool to reveal when and why mitochondrial markers are potentially misleading. In this study, we evaluate the performance of mtDNA markers to unravel the evolutionary history of Spanish lizards from the Podarcis hispanicus species complex. We first uncover several cases of discordance between morphological and mitochondrial data in delimitation of taxa. To assess the origin of these discordances, we analysed the same populations using several independent nuclear loci. Both morphological and nuclear markers identified the same three evolutionary units in the region, while mitochondrial data revealed four deeply divergent lineages. We suggest here that the most likely scenario to explain this discordance is ancient mitochondrial introgression originating from a fourth evolutionary unit presently absent from the study area. Notably, this resulted in a complete replacement of the original lineage in a large part of the distribution of one of the taxa investigated. We discuss the potential evolutionary scenarios leading to this complete mitochondrial replacement and suggest why the previous studies have failed to recover the correct history of this species complex.     

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.     

Differential introgression across newt hybrid zones: Evidence from replicated transects

Genomic heterogeneity of divergence between hybridizing species may reflect heterogeneity of introgression, but also processes unrelated to hybridization. Heterogeneous introgression and its repeatability can be directly tested in natural hybrid zones by examining multiple transects. Here, we studied hybrid zones between the European newts Lissotriton montandoni and two lineages of Lissotriton vulgaris, with replicate transects within each zone. Over 1,000 nuclear genes located on a linkage map and mitochondrial DNA were investigated using geographical and genomic clines. Overall, the five transects were all similar, showing hallmarks of strong reproductive isolation: bimodal distribution of genotypes in central populations and narrow allele frequency clines. However, the extent of introgression differed between the zones, possibly as a consequence of their different ages, as suggested by the analysis of heterozygosity runs in diagnostic markers. In three transects genomic signatures of small‐scale (~2 km) zone movements were detected. We found limited overlap of cline outliers between transects, and only weak evidence of stronger differentiation of introgression between zones than between transects within zones. Introgression was heterogeneous across linkage groups, with patterns of heterogeneity similar between transects and zones. Predefined candidates for increased or reduced introgression exhibited only a subtle tendency in the expected direction, suggesting that interspecific differentiation is not a reliable indicator for the strength of introgression. These hierarchically sampled hybrid zones of apparently different ages show how introgression unfolds with time and offer an excellent opportunity to dissect the dynamics of hybridization and architecture of reproductive isolation at advanced stages of speciation.     

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.     

Cryptic genetic diversity in Rattus of the San Francisco Bay region, California

Invasive species can have complex invasion histories, harbor cryptic levels of diversity, and pose taxonomic problems for pest management authorities. Roof rats, Rattus rattus sensu lato, are common invasive pests of the San Francisco Bay Area in California, USA. They are a significant health risk and pest management efforts impose a large financial investment from public institutions and private individuals. Recent molecular genetic and taxonomic studies of black rats in their native range in Asia have shown that the species is a complex of two karyotypic forms and four mitochondrial genetic lineages that may represent four distinct species. We used mtDNA sequences and nuclear microsatellite variation to identify which mitochondrial lineages of the R. rattus group are present in the San Francisco Bay Area and to test for gene flow among them. We recovered specimens with mtDNA sequences representing two of the major mtDNA lineages of the R. rattus group. Microsatellite variation, however, was not structured in concordance with mtDNA lineages, suggesting a more complex history involving hybridization and introgression between these lineages. Although Aplin et al. (2011) and Lack et al. (2012) reported R. rattus Lineage II in North America, this is the first detailed examination of possible gene flow amongst lineages in this region.     

PERCHED AT THE MITO‐NUCLEAR CROSSROADS: DIVERGENT MITOCHONDRIAL LINEAGES CORRELATE WITH ENVIRONMENT IN THE FACE OF ONGOING NUCLEAR GENE FLOW IN AN AUSTRALIAN BIRD

Relationships among multilocus genetic variation, geography, and environment can reveal how evolutionary processes affect genomes. We examined the evolution of an Australian bird, the eastern yellow robin Eopsaltria australis, using mitochondrial (mtDNA) and nuclear (nDNA) genetic markers, and bioclimatic variables. In southeastern Australia, two divergent mtDNA lineages occur east and west of the Great Dividing Range, perpendicular to latitudinal nDNA structure. We evaluated alternative scenarios to explain this striking discordance in landscape genetic patterning. Stochastic mtDNA lineage sorting can be rejected because the mtDNA lineages are essentially distinct geographically for > 1500 km. Vicariance is unlikely: the Great Dividing Range is neither a current barrier nor was it at the Last Glacial Maximum according to species distribution modeling; nuclear gene flow inferred from coalescent analysis affirms this. Female philopatry contradicts known female‐biased dispersal. Contrasting mtDNA and nDNA demographies indicate their evolutionary histories are decoupled. Distance‐based redundancy analysis, in which environmental temperatures explain mtDNA variance above that explained by geographic position and isolation‐by‐distance, favors a nonneutral explanation for mitochondrial phylogeographic patterning. Thus, observed mito‐nuclear discordance accords with environmental selection on a female‐linked trait, such as mtDNA, mtDNA–nDNA interactions or genes on W‐chromosome, driving mitochondrial divergence in the presence of nuclear gene flow.     

Introgression and selection shaped the evolutionary history of sympatric sister‐species of coral reef fishes (genus: Haemulon)

Closely related marine species with large overlapping ranges provide opportunities to study mechanisms of speciation, particularly when there is evidence of gene flow between such lineages. Here, we focus on a case of hybridization between the sympatric sister‐species Haemulon maculicauda and H. flaviguttatum, using Sanger sequencing of mitochondrial and nuclear loci, as well as 2422 single nucleotide polymorphisms (SNPs) obtained via restriction site‐associated DNA sequencing (RADSeq). Mitochondrial markers revealed a shared haplotype for COI and low divergence for CytB and CR between the sister‐species. On the other hand, complete lineage sorting was observed at the nuclear loci and most of the SNPs. Under neutral expectations, the smaller effective population size of mtDNA should lead to fixation of mutations faster than nDNA. Thus, these results suggest that hybridization in the recent past (0.174–0.263 Ma) led to introgression of the mtDNA, with little effect on the nuclear genome. Analyses of the SNP data revealed 28 loci potentially under divergent selection between the two species. The combination of mtDNA introgression and limited nuclear DNA introgression provides a mechanism for the evolution of independent lineages despite recurrent hybridization events. This study adds to the growing body of research that exemplifies how genetic divergence can be maintained in the presence of gene flow between closely related species.     

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.