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.     

Multiple nuclear loci reveal patterns of incomplete lineage sorting and complex species history within western mouse lemurs (Microcebus)

Mouse lemurs (genus Microcebus) are nocturnal primates endemic to the island of Madagascar. Until recently, they were classified as two species, one from eastern and one from western Madagascar. Previously published analyses of morphometric and mitochondrial DNA data show strong support for the recognition of more than eight species, however. Here, we test the eight-species hypothesis with DNA data derived from four independently segregating nuclear loci. We find many areas of congruence between the mitochondrial and nuclear data, but incomplete lineage sorting and low mutation rates limit the phylogenetic resolution of the nuclear data. Even so, the nuclear loci unanimously find evidence for three deeply diverged lineages within the mouse lemur radiation: one that is congruent with the mtDNA "southern clade", another that is congruent with the mtDNA "northern clade", and one monospecific branch comprised of the species Microcebus ravelobensis. The latter result in particular emphasizes the need for careful biological study of this species.     

Multilocus phylogeography of a holarctic duck: colonization of north america from eurasia by gadwall (Anas strepera)

More than 100 species of birds have Holarctic distributions extending across Eurasia and North America, and many of them likely achieved these distributions by recently colonizing one continent from the other. Mitochondrial DNA (mtDNA) and five nuclear introns were sequenced to test the direction and timing of colonization for a Holarctic duck, the gadwall (Anas strepera). Three lines of evidence suggest gadwalls colonized North America from Eurasia. First, New World (NW) gadwalls had fewer alleles at every locus and 61% of the allelic richness found in Old World (OW) gadwalls. Second, NW gadwalls had lower mtDNA allelic richness than other NW ducks. Third, coalescent analysis suggested that less than 5% of the ancestral population contributed to NW gadwalls at the time of divergence. Gadwalls likely colonized North America during the Late Pleistocene (approximately 81,000 years ago), but the confidence interval on that estimate was large (8500-450,000 years ago). Intercontinental gene flow and selection also likely contributed to genetic diversity in gadwalls. This study illustrates the use of multiple loci and coalescent analyses for critically testing a priori hypotheses regarding dispersal and colonization and provides an independent datapoint supporting an OW to NW bias in the direction of colonization.     

The causes of mitochondrial DNA gene tree paraphyly in birds

Gene tree paraphyly is a potentially serious problem because many phylogenetic and phylogeographic studies assume species are monophyletic. Funk and Omland (Funk, D.J., Omland, K.E., 2003. Species-level paraphyly and polyphyly: frequency, causes, and consequences, with insights from animal mitochondrial DNA. Annu. Rev. Ecol. Evol. Syst. 34, 397–423) found that a seemingly high proportion of bird species (16.7%) were paraphyletic in their mtDNA gene trees. This could imply that mtDNA is an unreliable or even misleading marker for delimiting species. We expand on Funk and Omland’s survey and identify the causes of species-level paraphyly in birds. We find that in most cases paraphyly is caused by incorrect taxonomy. In such cases, mtDNA serves systematics by exposing and clarifying taxonomic errors. We find the next most common cause of paraphyly to be incomplete lineage sorting due to recent speciation. Here mtDNA gives a consistent picture of evolution, given the timeframe, but it is not useful for delimiting species and other criteria must be employed. There were relatively few clear instances of paraphyly due to hybridization, though there were more cases where incomplete lineage sorting and hybridization could not be distinguished. We ultimately conclude that, far from a hindrance, mtDNA is generally a useful tool that should continue to facilitate delimitation of avian species.     

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.     

The dawn of human matrilineal diversity

The quest to explain demographic history during the early part of human evolution has been limited because of the scarce paleoanthropological record from the Middle Stone Age. To shed light on the structure of the mitochondrial DNA (mtDNA) phylogeny at the dawn of Homo sapiens, we constructed a matrilineal tree composed of 624 complete mtDNA genomes from sub-Saharan Hg L lineages. We paid particular attention to the Khoi and San (Khoisan) people of South Africa because they are considered to be a unique relic of hunter-gatherer lifestyle and to carry paternal and maternal lineages belonging to the deepest clades known among modern humans. Both the tree phylogeny and coalescence calculations suggest that Khoisan matrilineal ancestry diverged from the rest of the human mtDNA pool 90,000-150,000 years before present (ybp) and that at least five additional, currently extant maternal lineages existed during this period in parallel. Furthermore, we estimate that a minimum of 40 other evolutionarily successful lineages flourished in sub-Saharan Africa during the period of modern human dispersal out of Africa approximately 60,000-70,000 ybp. Only much later, at the beginning of the Late Stone Age, about 40,000 ybp, did introgression of additional lineages occur into the Khoisan mtDNA pool. This process was further accelerated during the recent Bantu expansions. Our results suggest that the early settlement of humans in Africa was already matrilineally structured and involved small, separately evolving isolated populations.     

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.     

Maternal and paternal lineages in cross-breeding bovine species. Has wisent a hybrid origin?

The tribe Bovini comprises cattle and cattle-like species. Reconstructions of their phylogeny have so far been incomplete and have yielded conflicting conclusions about the relationship of American bison and wisent (European bison). We have compared the sequences of three mitochondrial and two Y-chromosomal DNA segments. Mitochondrial DNA indicates that four distinct maternal lineages diverged after an early split-off of the buffalo species, leading to (1) taurine cattle and zebu, (2) wisent, (3) American bison and yak, and (4) banteng, gaur, and gayal, respectively. At a higher level, lineages (1) and (2) and lineages (3) and (4) are probably associated. In contrast, Y-chromosomal sequences indicate a close association of American and European bison, which is in agreement with their morphological similarity, complete fertility of hybrid offspring, and amplified fragment length polymorphism (AFLP) fingerprints of nuclear DNA. One explanation for the anomalous divergence of the mitochondrial DNA from the two bison species is lineage sorting, which implies that two distinct mitochondrial lineages coexisted in the bison-yak branch until the recent divergence of American bison and wisent. Alternatively, the wisent may have emerged by species hybridization initiated by introgression of bison bulls in another ancestral species. This "transpatric" mode of species formation would be consistent with the recent appearance of the wisent in the fossil record without clearly identifiable ancestors.     

Species tree inference in a recent radiation of orioles (Genus Icterus): multiple markers and methods reveal cytonuclear discordance in the northern oriole group

Recent computational advances provide novel opportunities to infer species trees based on multiple independent loci. Thus, single gene trees no longer need suffice as proxies for species phylogenies. Several methods have been developed to deal with the challenges posed by incomplete and stochastic lineage sorting. In this study, we employed four Bayesian methods to infer the phylogeny of a clade of 11 recently diverged oriole species within the genus Icterus. We obtained well-resolved and mostly congruent phylogenies using a set of seven unlinked nuclear intron loci and sampling multiple individuals per species. Most notably, Bayesian concordance analysis generally agreed well with concatenation; the two methods agreed fully on eight of nine nodes. The coalescent-based method ^∗BEAST further supported six of these eight nodes. The fourth method used, BEST, failed to converge despite exhaustive efforts to optimize the tree search. Overall, the results obtained by new species tree methods and concatenation generally corroborate our findings from previous analyses and data sets. However, we found striking disagreement between mitochondrial and nuclear DNA involving relationships within the northern oriole group. Our results highlight the danger of reliance on mtDNA alone for phylogenetic inference. We demonstrate that in spite of low variability and incomplete lineage sorting, multiple nuclear loci can produce largely congruent phylogenies based on multiple species tree methods, even for very closely-related species.     

Recurrent introgression of mitochondrial DNA among hares (Lepus spp.) revealed by species-tree inference and coalescent simulations

Understanding recent speciation history requires merging phylogenetic and population genetics approaches, taking into account the persistence of ancestral polymorphism and possible introgression. The emergence of a clear phylogeny of hares (genus Lepus) has been hampered by poor genomic sampling and possible occurrence of mitochondrial DNA (mtDNA) introgression from the arctic/boreal Lepus timidus into several European temperate and possibly American boreal species. However, no formal test of introgression, taking also incomplete lineage sorting into account, has been done. Here, to clarify the yet poorly resolved species phylogeny of hares and test hypotheses of mtDNA introgression, we sequenced 14 nuclear DNA and 2 mtDNA fragments (8205 and 1113 bp, respectively) in 50 specimens from 11 hare species from Eurasia, North America, and Africa. By applying an isolation-with-migration model to the nuclear data on subsets of species, we find evidence for very limited gene flow from L. timidus into most temperate European species, and not into the American boreal ones. Using a multilocus coalescent-based method, we infer the species phylogeny, which we find highly incongruent with mtDNA phylogeny using parametric bootstrap. Simulations of mtDNA evolution under the speciation history inferred from nuclear genes did not support the hypothesis of mtDNA introgression from L. timidus into the American L. townsendii but did suggest introgression from L. timidus into 4 temperate European species. One such event likely resulted in the complete replacement of the aboriginal mtDNA of L. castroviejoi and of its sister species L. corsicanus. It is remarkable that mtDNA introgression in hares is frequent, extensive, and always from the same donor arctic species. We discuss possible explanations for the phenomenon in relation to the dynamics of range expansions and species replacements during the climatic oscillations of the Pleistocene.     

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.     

Genetic admixture supports an ancient hybrid origin of the endangered Hawaiian duck

Speciation is regarded primarily as a bifurcation from an ancestral species into two distinct taxonomic units, but gene flow can create complex signals of phylogenetic relationships, especially among different loci. We evaluated several hypotheses that could account for phylogenetic discord between mitochondrial DNA (mtDNA) and nuclear DNA (nuDNA) within Hawaiian duck (Anas wyvilliana), including stochastic lineage sorting, mtDNA capture and widespread genomic introgression. Our results best support the hypothesis that the contemporary Hawaiian duck is descended from an ancient hybridization event between the mallard (Anas platyrhynchos) and Laysan duck (Anas laysanensis). Whereas mtDNA clearly shows a sister relationship between Hawaiian duck and mallard, nuDNA is consistent with a genetic mosaic with nearly equal contributions from Laysan duck and mallard. In addition, coalescent analyses suggest that gene flow from either mallard or Laysan duck, depending on the predefined tree topology, is necessary to explain contemporary genetic diversity in Hawaiian ducks, and these estimates are more consistent with ancient, rather than contemporary, hybridization. Time since divergence estimates suggest that the genetic admixture event occurred around the Pleistocene–Holocene boundary, which is further supported by circumstantial evidence from the Hawaiian subfossil record. Although the extent of reproductive isolation from either putative parental taxon is not currently known, these species are phenotypically, genetically and ecologically different, and they meet primary criteria used in avian taxonomy for species designation. Thus, the available data are consistent with an admixed origin and support the hypothesis that the Hawaiian duck may represent a young hybrid species.     

Contrasting patterns of genetic structure in two species of the coral trout Plectropomus (Serranidae) from east and west Australia: introgressive hybridisation or ancestral polymorphisms

Inter-specific genetic relationships among regional populations of two species of grouper (Plectropomus maculatus and Plectropomus leopardus) were examined using mitochondrial and nuclear markers. mtDNA revealed contrasting regional inter-specific patterns whilst nuclear markers revealed contrasting patterns among markers, irrespective of region. In eastern Australia (EA) the species form a single mtDNA lineage, but the two species are reciprocally monophyletic in Western Australia (WA). This supports previous evidence for hybridisation between these species on the east coast. WA P. leopardus forms a sister relationship with the EA P. leopardus-maculatus clade while WA P. maculatus is more basal and sister to the P. leopardus lineages, indicating mtDNA does not suffer from incomplete lineage sorting for these species. In contrast, one of three nuclear markers (locus 7-90TG) differentiated the species into two reciprocally monophyletic clades, with no evidence of hybridisation or ancestral polymorphism. The remaining two nuclear markers (2-22 and ETS-2) did not separate these two species, while distinguishing other plectropomid species, suggesting incomplete lineage sorting at these nuclear loci. These results together with coalescence analyses suggest that P. leopardus females have hybridised historically with P. maculatus males and that P. maculatus mitochondria were displaced through introgressive hybridisation and fixation in the P. maculatus founder population on the Great Barrier Reef. The contrasting regional patterns of mtDNA structure may be attributed to Quaternary sea-level changes and shelf width differences driving different reef configurations on each coast. These reef configurations have provided opportunities for local scale interaction and reproduction among species on the narrower EA continental shelves, but not on the broader WA continental shelves.     

TESTING MITOCHONDRIAL CAPTURE AND DEEP COALESCENCE IN AMAZONIAN CICHLID FISHES (CICHLIDAE: CICHLA)

Hybridization and introgression have important consequences in evolution, such as increasing the genetic diversity and adaptive potential of a species. One of their most conspicuous footprints is discordance among gene trees or between genes and phenotypes. However, most studies that report introgression fail to disprove the null hypothesis that genetic incongruence may result from stochastic sorting of ancestral allelic polymorphisms. In the case of ancient introgression, these two processes may be especially difficult to distinguish topologically, but they make different predictions about the patterns of coalescence among loci. Here we apply three methods, molecular dating, multispecies coalescent models, and gene tree simulation under coalescence, to compare these two hypotheses that explain the polyphyletic mtDNA of the butterfly peacock bass, Cichla orinocensis. In comparison with a species tree based on 20 unlinked nuclear loci, we determined that mtDNA divergences were too recent to be explained by ancestral polymorphism. Similarly, coalescent species tree branches were significantly shorter when putative introgressed mtDNA was incorporated, and simulations showed the mtDNA topology to be unlikely under lineage sorting only. We conclude that introgression approximately 1.5 million years ago resulted in capture by C. orinocensis of an mtDNA lineage ancestral to the modern subspecies C. oc. monoculus.     

ASYMMETRIC INTROGRESSION OF MITOCHONDRIAL DNA AMONG EUROPEAN POPULATIONS OF BLUE MUSSELS (MYTILUS SPP.)

Abstract.—Mytilus edulis and M. galloprovincialis are two blue mussel species that coexist in western Europe. Previously, we reported that M. galloprovincialis populations contain female and male haplotypes that are fixed in M. edulis populations as well as unique haplotypes. This study assesses whether paraphyly for these species is due to introgression or incomplete lineage extinction. The lineage extinction hypothesis predicts that the shared mtDNA haplotypes in M. galloprovincialis will be significantly diverged from those in M. edulis and form distinct sequence clades. In contrast, the introgression hypothesis proposes that M. edulis haplotypes have only recently been introduced into M. galloprovincialis through hybridization with relatively little divergence accumulating between the shared RFLP haplotypes. We examined 80 mtl6S gene sequences for both the maternal and paternal mtDNA lineages from mussels sampled from various European populations and found strong support for the introgression hypothesis. In addition, we found that M. edulis mtDNA haplotypes appear to be introgressing into mussel populations in the Baltic Sea, which have predominantly M. trossulus nuclear genotypes, indicating that introgressive hybridization is prevalent among European mussel populations.     

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.     

DNA-based species delimitation in algae

Given the problems of species delimitation in algae using morphology or sexual compatibility, molecular data are becoming the standard for delimiting species and testing their traditional boundaries. The idea that species are separately evolving metapopulation lineages, along with theoretical progress in phylogenetic and population genetic analyses, has led to the development of new methods of species delimitation. We review these recent developments in DNA-based species delimitation methods, and discuss how they have changed and continue to change our understanding of algal species boundaries. Although single-locus approaches have proven effective for a first rapid and large-scale assessment of species diversity, species delimitation based on single gene trees falls short due to gene tree–species tree incongruence, caused by confounding processes like incomplete lineage sorting, trans-species polymorphism, hybridization and introgression. Data from unlinked loci and multi-species coalescent methods, which combine principles from phylogenetics and population genetics, may now be able to account for these complicating factors. Several of these methods also provide statistical support regarding species boundaries, which is important because speciation is a process and therefore uncertainty about precise species boundaries is inevitable in recently diverged lineages.     

Isolation‐driven divergence: speciation in a widespread North American songbird (Aves: Certhiidae)

Lineage, or true ‘species’, trees may differ from gene trees because of stochastic processes in molecular evolution leading to gene‐tree heterogeneity. Problems with inferring species trees because of excessive incomplete lineage sorting may be exacerbated in lineages with rapid diversification or recent divergences necessitating the use of multiple loci and individuals. Many recent multilocus studies that investigate divergence times identify lineage splitting to be more recent than single‐locus studies, forcing the revision of biogeographic scenarios driving divergence. Here, we use 21 nuclear loci from regional populations to re‐evaluate hypotheses identified in an mtDNA phylogeographic study of the Brown Creeper (Certhia americana), as well as identify processes driving divergence. Nuclear phylogeographic analyses identified hierarchical genetic structure, supporting a basal split at approximately 32°N latitude, splitting northern and southern populations, with mixed patterns of genealogical concordance and discordance between data sets within the major lineages. Coalescent‐based analyses identify isolation, with little to no gene flow, as the primary driver of divergence between lineages. Recent isolation appears to have caused genetic bottlenecks in populations in the Sierra Madre Oriental and coastal mountain ranges of California, which may be targets for conservation concerns.     

A multilocus phylogeny of the Sulidae (Aves: Pelecaniformes)

Gene trees will often differ from the true species history, the species tree, as a result of processes such as incomplete lineage sorting. New methods such as Bayesian Estimation of the Species Tree (BEST) use the multispecies coalescent to model lineage sorting, and directly infer the species tree from multilocus DNA sequence data. The Sulidae (Aves: Pelecaniformes) is a family of ten booby and gannet species with a global distribution. We sequenced five nuclear intron loci and one mitochondrial locus to estimate a species tree for the Sulidae using both BEST and by concatenating nuclear loci. We also used fossil calibrated strict and relaxed molecular clocks in BEAST to estimate divergence times for major nodes in the sulid phylogeny. Individual gene trees showed little phylogenetic conflict but varied in resolution. With the exception of the mitochondrial gene tree, no gene tree was completely resolved. On the other hand, both the BEST and concatenated species trees were highly resolved, strongly supported, and topologically consistent with each other. The three sulid genera (Morus, Sula, Papasula) were monophyletic and the relationships within genera were mostly consistent with both a previously estimated mtDNA gene tree and the mtDNA gene tree estimated here. However, our species trees conflicted with the mtDNA gene trees in the relationships among the three genera. Most notably, we find that the endemic and endangered Abbott's booby (Papasula abbotti) is likely basal to all other members of the Sulidae and diverged from them approximately 22 million years ago.     

Mito‐nuclear discord in six congeneric lineages of Holarctic ducks (genus Anas)

Many species have Holarctic distributions that extend across Europe, Asia and North America. Most genetics research on these species has examined only mitochondrial (mt) DNA, which has revealed wide variance in divergence between Old World (OW) and New World (NW) populations, ranging from shallow, unstructured genealogies to deeply divergent lineages. In this study, we sequenced 20 nuclear introns to test for concordant patterns of OW–NW differentiation between mtDNA and nuclear (nu) DNA for six lineages of Holarctic ducks (genus Anas). Genetic differentiation for both marker types varied widely among these lineages (idiosyncratic population histories), but mtDNA and nuDNA divergence within lineages was not significantly correlated. Moreover, compared with the association between mtDNA and nuDNA divergence observed among different species, OW–NW nuDNA differentiation was generally lower than mtDNA divergence, at least for lineages with deeply divergent mtDNA. Furthermore, coalescent estimates indicated significantly higher rates of gene flow for nuDNA than mtDNA for four of the six lineages. Thus, Holarctic ducks show prominent mito‐nuclear discord between OW and NW populations, and we reject differences in sorting rates as the sole cause of the within‐species discord. Male‐mediated intercontinental gene flow is likely a leading contributor to this discord, although selection could also cause increased mtDNA divergence relative to weak nuDNA differentiation. The population genetics of these ducks contribute to growing evidence that mtDNA can be an unreliable indicator of stage of speciation and that more holistic approaches are needed for species delimitation.