Nuclear loci and coalescent methods support ancient hybridization as cause of mitochondrial paraphyly between gadwall and falcated duck (Anas spp.)

Many species have mitochondrial DNA lineages that are phylogenetically intermixed with other species, but studies have rarely tested the cause of such paraphyly. In this study, we tested two hypotheses that could explain mitochondrial paraphyly of Holarctic gadwalls (Anas strepera) with respect to Asian falcated ducks (A. falcata). First, hybridization could have resulted in falcated duck mitochondrial DNA (mtDNA) introgressing into the gadwall gene pool. Second, gadwalls and falcated ducks could have diverged so recently that mtDNA lineages have not sorted to reciprocal monophyly. We used coalescent analyses of three independent loci to distinguish between these two hypotheses. Two lines of evidence support introgression. First, analyses of the three loci combined show that some introgression is necessary to explain current genetic diversity in gadwalls. Second, we generated alternative predictions regarding time since divergence estimated from mtDNA: falcated ducks and gadwalls would have diverged between 65,000 and 700,000 years before present (ybp) under the introgression hypothesis and between 11,000 and 76,000 ybp under the incomplete lineage sorting hypothesis. The two independent nuclear introns indicated that these species diverged between 210,000 and 5,200,000 ybp, which did not overlap the predicted time for incomplete lineage sorting. These analyses also suggested that ancient introgression ( approximately 14,000 ybp) has resulted in the widespread distribution and high frequency of falcated-like mtDNA (5.5% of haplotypes) in North America. This is the first study to use a rigorous quantitative framework to reject incomplete lineage sorting as the cause of mitochondrial paraphyly.     

A tangled tale of two teal: population history of the grey Anas gracilis and chestnut teal A. castanea of Australia

Two Australian species of teal (Anseriformes: Anatidae: Anas ), the grey teal Anas gracilis and the chestnut teal A. castanea , are remarkable for the zero or near-zero divergence recorded between them in earlier surveys of mitochondrial DNA (mtDNA) diversity. We confirmed this result through wider geographical and population sampling as well as nucleotide sampling in the more rapidly evolving mtDNA control region. Any data set where two species share polymorphism as is the case here can be explained by a model of gene flow through hybridization on one hand or by incomplete lineage sorting on the other hand. Ideally, analysis of such shared polymorphism would simultaneously estimate the likelihood of both phenomena. To do this, we used the underlying principle of the IMa package to explore ramifications to understanding population histories of A. gracilis and A. castanea . We cannot reject that hybridization occurs between the two species but an equally or more plausible finding for their nearly zero divergence is incomplete sorting following very recent divergence between the two, probably in the mid-late Pleistocene. Our data add to studies that explore intermediate stages in the evolution of reciprocal monophyly and paraphyletic or polyphyletic relationships in mtDNA diversity among widespread Australian birds.     

Utility of microsatellites and mitochondrial DNA for species delimitation in the spruce budworm (Choristoneura fumiferana) species complex (Lepidoptera: Tortricidae)

Species identifications have been historically difficult in the Choristoneura fumiferana group, an important insect pest complex. We examined the utility of simple sequence repeats (SSRs, also referred to as microsatellites) and mitochondrial DNA (mtDNA) for delimiting and identifying eight currently recognized species sampled across North America. Four of these species formed discrete clusters using SSRs, while only two species were delimited with mtDNA. There was evidence for hybridization or incomplete lineage sorting between several species pairs. An integrative approach, using both phenotypic traits and molecular markers, allowed for the discrimination of more biologically relevant species units than did the use of molecular markers alone. As species are currently identified using putatively adaptive phenotypic traits, the differences observed between recognized species and neutral SSRs or mtDNA suggests that these species (or evolutionary significant units) have diverged via natural selection in spite of some gene flow.     

Detecting hybrid speciation in the presence of incomplete lineage sorting using gene tree incongruence: A model

The application of phylogenetic inference methods, to data for a set of independent genes sampled randomly throughout the genome, often results in substantial incongruence in the single-gene phylogenetic estimates. Among the processes known to produce discord between single-gene phylogenies, two of the best studied in a phylogenetic context are hybridization and incomplete lineage sorting. Much recent attention has focused on the development of methods for estimating species phylogenies in the presence of incomplete lineage sorting, but phylogenetic models that allow for hybridization have been more limited. Here we propose a model that allows incongruence in single-gene phylogenies to be due to both hybridization and incomplete lineage sorting, with the goal of determining the contribution of hybridization to observed gene tree incongruence in the presence of incomplete lineage sorting. Using our model, we propose methods for estimating the extent of the role of hybridization in both a likelihood and a Bayesian framework. The performance of our methods is examined using both simulated and empirical data.     

Phylogeny estimation of the radiation of western North American chipmunks (Tamias) in the face of introgression using reproductive protein genes

The causes and consequences of rapid radiations are major unresolved issues in evolutionary biology. This is in part because phylogeny estimation is confounded by processes such as stochastic lineage sorting and hybridization. Because these processes are expected to be heterogeneous across the genome, comparison among marker classes may provide a means of disentangling these elements. Here we use introns from nuclear-encoded reproductive protein genes expected to be resistant to introgression to estimate the phylogeny of the western chipmunks (Tamias: subgenus: Neotamias), a rapid radiation that has experienced introgressive hybridization of mitochondrial DNA (mtDNA). We analyze the nuclear loci using coalescent-based species-tree estimation methods and concatenation to estimate a species tree and we use parametric bootstraps and coalescent simulations to differentiate between phylogenetic error, coalescent stochasticity and introgressive hybridization. Results indicate that the mtDNA gene tree reflects several introgression events that have occurred between taxa of varying levels of divergence and at different time points in the tree. T. panamintinus and T. speciosus appear to be fixed for ancient mitochondrial introgressions from T. minimus. A southern Rocky Mountains clade appears well sorted (i.e., species are largely monophyletic) at multiple nuclear loci, while five of six taxa are nonmonophyletic based on cytochrome b. Our simulations reject phylogenetic error and coalescent stochasticity as causes. The results represent an advance in our understanding of the processes at work during the radiation of Tamias and suggest that sampling reproductive-protein genes may be a viable strategy for phylogeny estimation of rapid radiations in which reproductive isolation is incomplete. However, a genome-scale survey that can statistically compare heterogeneity of genealogical process at many more loci will be necessary to test this conclusion.     

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.     

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.     

Diversification of Hawaiian Cyrtandra (Gesneriaceae) under the influence of incomplete lineage sorting and hybridization

Cyrtandra (Gesneriaceae) is a genus of flowering plants with over 800 species distributed throughout Southeast Asia and the Pacific Islands. On the Hawaiian Islands, 60 named species and over 89 putative hybrids exist, most of which are identified on the basis of morphology. Despite many previous studies on the Hawaiian Cyrtandra lineage, questions regarding the reconciliation of morphology and genetics remain, many of which can be attributed to the relatively young age and evidence of hybridization between species. We utilized targeted enrichment, high‐throughput sequencing, and modern phylogenomics tools to test 31 Hawaiian Cyrtandra samples (22 species, two putative hybrids, four species with two samples each, one species with four samples) and two outgroups for species relationships and hybridization in the presence of incomplete lineage sorting (ILS). Both concatenated and species‐tree methods were used to reconstruct species relationships, and network analyses were conducted to test for hybridization. We expected to see high levels of ILS and putative hybrids intermediate to their parent species. Phylogenies reconstructed from the concatenated and species‐tree methods were highly incongruent, most likely due to high levels of incomplete lineage sorting. Network analyses inferred gene flow within this lineage, but not always between taxa that we expected. Multiple hybridizations were inferred, but many were on deeper branches of the island lineages suggesting a long history of hybridization. We demonstrated the utility of high‐throughput sequencing and a phylogenomic approach using 569 loci to understanding species relationships and gene flow in the presence of ILS.     

Inferring phylogeny despite incomplete lineage sorting

It is now well known that incomplete lineage sorting can cause serious difficulties for phylogenetic inference, but little attention has been paid to methods that attempt to overcome these difficulties by explicitly considering the processes that produce them. Here we explore approaches to phylogenetic inference designed to consider retention and sorting of ancestral polymorphism. We examine how the reconstructability of a species (or population) phylogeny is affected by (a) the number of loci used to estimate the phylogeny and (b) the number of individuals sampled per species. Even in difficult cases with considerable incomplete lineage sorting (times between divergences less than 1 N(e) generations), we found the reconstructed species trees matched the "true" species trees in at least three out of five partitions, as long as a reasonable number of individuals per species were sampled. We also studied the tradeoff between sampling more loci versus more individuals. Although increasing the number of loci gives more accurate trees for a given sampling effort with deeper species trees (e.g., total depth of 10 N(e) generations), sampling more individuals often gives better results than sampling more loci with shallower species trees (e.g., depth = 1 N(e)). Taken together, these results demonstrate that gene sequences retain enough signal to achieve an accurate estimate of phylogeny despite widespread incomplete lineage sorting. Continued improvement in our methods to reconstruct phylogeny near the species level will require a shift to a compound model that considers not only nucleotide or character state substitutions, but also the population genetics processes of lineage sorting. [Coalescence; divergence; population; speciation.].     

Is homoplasy or lineage sorting the source of incongruent mtdna and nuclear gene trees in the stiff-tailed ducks (Nomonyx-Oxyura)?

We evaluated the potential effects of homoplasy, ancestral polymorphism, and hybridization as obstacles to resolving phylogenetic relationships within Nomonyx-Oxyura stiff-tailed ducks (Oxyurinae; subtribe Oxyurina). Mitochondrial DNA (mtDNA) control region sequences from 94 individuals supported monophyly of mtDNA haplotypes for each of the six species and provided no evidence of extant incomplete lineage sorting or inter-specific hybridization. The ruddy ducks (O. j. jamaicensis,O. j. andina, O. j. ferruginea) are each others' closest relatives, but the lack of shared haplotypes between O. j. jamaicensis and O. j. ferruginea suggests long-standing historical isolation. In contrast, O. j. andina shares haplotypes with O. j. jamaicensis and O. j. ferruginea, which supports Todd's (1979) and Fjelds?'s (1986) hypothesis that O. j. andina is an intergrade or hybrid subspecies of O. j. jamaicensis and O. j. ferruginea. Control region data and a much larger data set composed of approximately 8800 base pairs of mitochondrial and nuclear sequence for each species indicate that the two New World species, O. vittata and O. jamaicensis, branch basally within Oxyura. A clade of three Old World species (O. australis, O. maccoa, O. leucocephala) is well supported, but different loci and also different characters within the mtDNA data support three different resolutions of the Old World clade, yielding an essentially unresolved trichotomy. Fundamentally different factors limited the resolution of the mtDNA and nuclear gene trees. Gene trees for most nuclear loci were unresolved due to slow rates of mutation and a lack of informative variation, whereas uncertain resolution of the mtDNA gene tree was due to homoplasy. Within the mtDNA, approximately equal numbers of characters supported each of three possible resolutions. Parametric and nonparametric bootstrap analyses suggest that resolution of the mtDNA tree based on ~4300 bp per taxon is uncertain but that complete mtDNA sequences would yield a fully resolved gene tree. A short internode separating O. leucocephala from (O. australis, O. maccoa) in the best mtDNA tree combined with long terminal branches and substantial rate variation among nucleotide sites allowed the small number of changes occurring on the internode to be obscured by homoplasy in a significant portion of simulated data sets. Although most nuclear loci were uninformative, two loci supported a resolution of the Old World clade (O. maccoa, O. leucocephala) that is incongruent with the best mtDNA tree. Thus, incongruence between nuclear and mtDNA trees may be due to random sorting of ancestral lineages during the short internode, homoplasy in the mtDNA data, or both. The Oxyura trichotomy represents a difficult though likely common problem in molecular systematics. Given a short internode, the mtDNA tree has a greater chance of being congruent with the history of speciation because its effective population size (N(e)) is one-quarter that of any nuclear locus, but its resolution is more likely to be obscured by homoplasy. In contrast, gene trees for more slowly evolving nuclear loci will be difficult to resolve due to a lack of substitutions during the internode, and when resolved are more likely to be incongruent with the species history due to the stochastic effects of lineage sorting. We suggest that researchers consider first whether independent gene trees are adequately resolved and then whether those trees are congruent with the species history. In the case of Oxyura, the answer to both questions may be no. Complete mtDNA sequences combined with data from a very large number of nuclear loci may be the only way to resolve such trichotomies.     

The evolutionary history of the coral genus Acropora (Scleractinia, Cnidaria) based on a mitochondrial and a nuclear marker: reticulation, incomplete lineage sorting, or morphological convergence?

This study examines molecular relationships across a wide range of species in the mass spawning scleractinian coral genus Acropora. Molecular phylogenies were obtained for 28 species using DNA sequence analyses of two independent markers, a nuclear intron and the mtDNA putative control region. Although the compositions of the major clades in the phylogenies based on these two markers were similar, there were several important differences. This, in combination with the fact that many species were not monophyletic, suggests either that introgressive hybridization is occurring or that lineage sorting is incomplete. The molecular tree topologies bear little similarity to the results of a recent cladistic analysis based on skeletal morphology and are at odds with the fossil record. We hypothesize that these conflicting results may be due to the same morphology having evolved independently more than once in Acropora and/or the occurrence of extensive interspecific hybridization and introgression in combination with morphology being determined by a small number of genes. Our results indicate that many Acropora species belong to a species complex or syngameon and that morphology has little predictive value with regard to syngameon composition. Morphological species in the genus often do not correspond to genetically distinct evolutionary units. Instead, species that differ in timing of gamete release tend to constitute genetically distinct clades.     

Eight independent nuclear genes support monophyly of the plovers: The role of mutational variance in gene trees

Molecular phylogenies of Charadriiformes based on mtDNA genes and one to three nuclear loci do not support the traditional placement of Pluvialis in the plovers (Charadriidae), assigning it instead to oystercatchers, stilts, and avocets (Haematopodidae and Recurvirostridae). To investigate this hypothesis of plover paraphyly, the relationships among Pluvialis and closely related families were revisited by sequencing two individuals of all taxa except Peltohyas for eight independent single copy nuclear protein-coding loci selected for their informativeness at this phylogenetic depth. The species tree estimated jointly with the gene trees in the coalescent programme (*)BEAST strongly supported plover monophyly, as did Bayesian analysis of the concatenated matrix. The data sets that supported plover paraphyly in Baker et al. (2007) and Fain and Houde (2007) reflect two to four independent gene histories, and thus discordance with the plover monophyly species tree might have arisen by chance through stochastic mutational variance. For the plovers we conclude there is no conclusive evidence of coalescent variance from ancient incomplete lineage sorting across the interior branch leading to Pluvialis in the species tree. Rather, earlier studies seem have been misled by faster evolving mtDNA genes with high mutational variance, and a few nuclear genes that had low resolving power at the Pluvialis sister group level. These findings are of general relevance in avian phylogenetics, as they show that careful attention needs to be paid to the number and the phylogenetic informativeness of genes required to obtain accurate estimates of the species tree, especially where there is mutational heterogeneity in gene trees.     

Coalescent simulations reveal hybridization and incomplete lineage sorting in Mediterranean Linaria

We examined the phylogenetic history of Linaria with special emphasis on the Mediterranean sect. Supinae (44 species). We revealed extensive highly supported incongruence among two nuclear (ITS, AGT1) and two plastid regions (rpl32-trnL(UAG), trnS-trnG). Coalescent simulations, a hybrid detection test and species tree inference in *BEAST revealed that incomplete lineage sorting and hybridization may both be responsible for the incongruent pattern observed. Additionally, we present a multilabelled *BEAST species tree as an alternative approach that allows the possibility of observing multiple placements in the species tree for the same taxa. That permitted the incorporation of processes such as hybridization within the tree while not violating the assumptions of the *BEAST model. This methodology is presented as a functional tool to disclose the evolutionary history of species complexes that have experienced both hybridization and incomplete lineage sorting. The drastic climatic events that have occurred in the Mediterranean since the late Miocene, including the Quaternary-type climatic oscillations, may have made both processes highly recurrent in the Mediterranean flora.     

Genetic introgression, incomplete lineage sorting and faulty taxonomy create multiple cases of polyphyly in a montane clade of tyrant-flycatchers (Elaenia, Tyrannidae)

Species-level paraphyly and polyphyly are pervasive phenomena in modern phylogenetic research and can be due to a number of factors. We explore a complicated pattern of nuclear and mitochondrial polyphyly in montane Neotropical Elaenia flycatchers. Using a combination of phylogenetic and population genetic methods, we demonstrate that no single factor is sufficient to account for this pattern of polyphyly, and that it is likely based on an interplay of three different factors: (i) faulty taxonomy which has led to the recognition of two polyphyletic species that are better classified as four biological species; (ii) a late Pleistocene hybridization event that resulted in two morphologically and ecologically distinct species sharing extremely similar mitochondrial DNA but distinct nuclear DNA profiles; and (iii) incomplete lineage sorting in a nuclear marker that results in a polyphyletic placement of species that are otherwise well-differentiated in mitochondrial DNA, morphology and ecology. Additionally, we demonstrate that the two clades of montane Elaenia exhibit a reverse pattern of mitochondrial and nuclear diversity, with high mitochondrial and low nuclear genetic diversity in one clade and vice versa in the other clade. A possible cause for this pattern is differences in population histories, with large panmictic population structures being conducive to the retention of ancient nuclear polymorphisms in Elaenia albiceps chilensis .     

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.     

Mitonuclear discordance in wolf spiders: Genomic evidence for species integrity and introgression

Systematists and taxonomists have benefited greatly from the emergence of molecular methods. Species identification has become straightforward through DNA barcoding and the rapid build‐up of massive DNA barcode reference libraries. In animals, mitonuclear discordance can significantly complicate the process of species identification and delimitation. The causes of mitonuclear discordance are either biological (e.g., introgression, incomplete lineage sorting, horizontal gene transfer androgenesis) or induced by operational factors (e.g., human error with specimen misidentification or incorrect species delimitation). Moreover, endosymbionts may play an important role in promoting fixation of mitochondrial genomes. Here, we study the mitonuclear discordance of wolf spiders species (Lycosidae) (independent cases from Alopecosa aculeata and Pardosa pullata groups) that share identical COI DNA barcodes. We approached the case utilizing double‐digest restriction site‐associated DNA sequencing (ddRADseq) to obtain and analyse genomic‐scale data. Our results suggest that the observed cases of mitonuclear discordance are not due to operational reasons but result from biological processes. Further analysis indicated introgression and that incomplete lineage sorting is unlikely to have been responsible for the observed discrepancy. Additional survey of endosymbionts provided ideas on further research and their role in shaping mitochondrial DNA distribution patterns. Thus, ddRADseq grants an efficient way to study the taxonomy of problematic groups with insight into underlying evolutionary processes.     

Coalescent histories on phylogenetic networks and detection of hybridization despite incomplete lineage sorting

Analyses of the increasingly available genomic data continue to reveal the extent of hybridization and its role in the evolutionary diversification of various groups of species. We show, through extensive coalescent-based simulations of multilocus data sets on phylogenetic networks, how divergence times before and after hybridization events can result in incomplete lineage sorting with gene tree incongruence signatures identical to those exhibited by hybridization. Evolutionary analysis of such data under the assumption of a species tree model can miss all hybridization events, whereas analysis under the assumption of a species network model would grossly overestimate hybridization events. These issues necessitate a paradigm shift in evolutionary analysis under these scenarios, from a model that assumes a priori a single source of gene tree incongruence to one that integrates multiple sources in a unifying framework. We propose a framework of coalescence within the branches of a phylogenetic network and show how this framework can be used to detect hybridization despite incomplete lineage sorting. We apply the model to simulated data and show that the signature of hybridization can be revealed as long as the interval between the divergence times of the species involved in hybridization is not too small. We reanalyze a data set of 106 loci from 7 in-group Saccharomyces species for which a species tree with no hybridization has been reported in the literature. Our analysis supports the hypothesis that hybridization occurred during the evolution of this group, explaining a large amount of the incongruence in the data. Our findings show that an integrative approach to gene tree incongruence and its reconciliation is needed. Our framework will help in systematically analyzing genomic data for the occurrence of hybridization and elucidating its evolutionary role.     

Ancient hybridization and phenotypic novelty within Lake Malawi's cichlid fish radiation

Does hybridization play a broad innovative role in evolution? Many studies have shown hybrid origins of individual species, particularly in major adaptive radiations, but this may be a consequence, rather than a cause, of the existence of many closely related species. Cases of hybridization in the early stages of major adaptive radiations are comparatively rare. Here, we report phylogenetic evidence for ancient introgression between distinct lineages of the species-rich Lake Malawi haplochromine cichlid fishes. Mitochondrial DNA (mtDNA) sequences indicated surprisingly close relationships between the shallow-water rocky habitat "Mbuna" species and a group of dark-adapted "Deep-Benthic" species specialized for feeding in low-light conditions (dawn/dusk, under overhangs, and deep water). By contrast, analyses of nuclear amplified fragment length polymorphism data demonstrated that these Deep-Benthic cichlids were more closely related to shallow water "Shallow-Benthic" soft-sediment feeders, a group that shares similar head and body morphology. A coalescent-based computer simulation indicated that the mtDNA similarity of rocky habitat Mbuna species and dark-adapted Deep-Benthic species was due to hybridization rather than incomplete lineage sorting. Comparisons of morphology indicated that some Deep-Benthic species possessed novel morphology not present in other Lake Malawi species groups. Thus, these analyses support the hypothesis that ancient hybridization occurred within the Lake Malawi cichlid radiation, that the event occurred before the radiation of a species group adapted to low-light benthic habitats, and that this group went on to dominate the deep-water regions of Lake Malawi. The results of this study contribute to a growing literature consistent with a creative role of hybridization in the evolution of species diversity and adaptive radiations.     

Mitochondrial capture and incomplete lineage sorting in the diversification of balitorine loaches (Cypriniformes,Balitoridae) revealed by mitochondrial and nuclear genes

Tang, Q.‐Y., Liu, S.‐Q., Yu, D., Liu, H.‐Z. & Danley, P.D. (2012) Mitochondrial capture and incomplete lineage sorting in the diversification of balitorine loaches (Cypriniformes, Balitoridae) revealed by mitochondrial and nuclear genes. —Zoologica Scripta, 41, 233–247. Understanding the diversification of species is a central goal of evolutionary biological studies. One powerful tool to investigate the speciation process is molecular systematics. Here, we use molecular methods to investigate the evolution of balitorine loaches belonging to two genera, Lepturichthys and Jinshaia. Both genera contain only two species (Lepturichthys fimbriata, Lepturichthys dolichopterus and Jinshaia sinensis and Jinshaia abbreviata), all of which are endemic to China. These species share many morphological and ecological characters and exhibit overlapping distributions in the Upper Yangtze River. In this study, we used two mitochondrial genes (Cytb and COI) and one nuclear gene (RAG1) to investigate the phylogenetic relationships within and between these two genera. Phylogenetic analyses and network construction based on mitochondrial and nuclear genes consistently supported the monophyly of Jinshaia. In contrast, the mitochondrial and nuclear genes yielded conflicting results in Lepturichthys. The phylogenetic analyses of mitochondrial sequences identify two distinct Lepturichthys lineages, Lepturichthys A and Lepturichthys B. Lepturichthys A includes most of L. fimbriata individuals from the Upper Yangtze River and is the sister group to all Jinshaia species. Lepturichthys B consists of the remaining L. fimbriata individuals from the Upper and Middle Yangtze River, and all L. dolichopterus individuals from the Minjiang River in Southeastern China. However, the analysis of the nuclear sequence indicates that the genus Lepturichthys is monophyletic and is only distantly related to Jinshaia. This incongruence suggests that introgressive hybridization might have occurred between L. fimbriata (Lepturichthys A) and Jinshaia species. As a result of this hybridization event, L. fimbriata captured the mitochondrial genome of the sympatric Jinshaia species. This capture event appears to have occurred at least 1.74 million years ago. Additionally, L. fimbriata appears to be paraphyletic; the nuclear data indicated that L. dolichopterus forms a monophyletic clade nested within L. fimbriata. Because L. dolichopterus and L. fimbriata are allopatric and hybridization may not be possible, we suggest that the observed paraphyly of L. fimbriata is a product of incomplete lineage sorting. In addition, the reciprocal monophyly of J. sinensis and J. abbreviata could not be resolved. This may be the result of interspecific hybridization as these species occur sympatrically. However, incomplete lineage sorting may have caused the observed topology of the Jinshaia species. The data presented here illustrate the complex evolutionary history of the balitorine loach species: intergeneric hybridization and interspecific hybridization have likely occurred in this lineage. In addition, possible incomplete lineage sorting may further obscure the evolutionary history of this group. The complex relationships of the balitorine loaches provide a rich evolutionary system to study the creation of sympatric and sister species.     

Bayesian Inference of Reticulate Phylogenies under the Multispecies Network Coalescent

The multispecies coalescent (MSC) is a statistical framework that models how gene genealogies grow within the branches of a species tree. The field of computational phylogenetics has witnessed an explosion in the development of methods for species tree inference under MSC, owing mainly to the accumulating evidence of incomplete lineage sorting in phylogenomic analyses. However, the evolutionary history of a set of genomes, or species, could be reticulate due to the occurrence of evolutionary processes such as hybridization or horizontal gene transfer. We report on a novel method for Bayesian inference of genome and species phylogenies under the multispecies network coalescent (MSNC). This framework models gene evolution within the branches of a phylogenetic network, thus incorporating reticulate evolutionary processes, such as hybridization, in addition to incomplete lineage sorting. As phylogenetic networks with different numbers of reticulation events correspond to points of different dimensions in the space of models, we devise a reversible-jump Markov chain Monte Carlo (RJMCMC) technique for sampling the posterior distribution of phylogenetic networks under MSNC. We implemented the methods in the publicly available, open-source software package PhyloNet and studied their performance on simulated and biological data. The work extends the reach of Bayesian inference to phylogenetic networks and enables new evolutionary analyses that account for reticulation.