Ancient Hybridization Research Articles

Secondary contact and genomic admixture between rhesus and long-tailed macaques in the Indochina Peninsula.

Understanding the process and consequences of hybridization is one of the major challenges in evolutionary biology. A growing body of literature has reported evidence of ancient hybridization events or natural hybrid zones in primates, including humans; however, we still have relatively limited knowledge about the pattern and history of admixture because there have been little studies that simultaneously achieved genome-scale analysis and a geographically wide sampling of wild populations. Our study applied double-digest restriction site-associated DNA sequencing to samples from the six localities in and around the provisional hybrid zone of rhesus and long-tailed macaques and evaluated population structure, phylogenetic relationships, demographic history, and geographic clines of morphology and allele frequencies. A latitudinal gradient of genetic components was observed, highlighting the transition from rhesus (north) to long-tailed macaque distribution (south) as well as the presence of one northern population of long-tailed macaques exhibiting unique genetic structure. Interspecific gene flow was estimated to have recently occurred after an isolation period, and the migration rate from rhesus to long-tailed macaques was slightly greater than in the opposite direction. Although some rhesus macaque-biased alleles have widely introgressed into long-tailed macaque populations, the inflection points of allele frequencies have been observed as concentrated around the traditionally recognized interspecific boundary where morphology discontinuously changed; this pattern was more pronounced in the X chromosome than in autosomes. Thus, due to geographic separation before secondary contact, reproductive isolation could have evolved, contributing to the maintenance of an interspecific boundary and species-specific morphological characteristics.

Hidden Hybridization and Habitat Differentiation in a Mediterranean Macrophyte, the Euryhaline Genus Ruppia.

In many aquatic plant taxa, classification based on morphology has always been difficult. Molecular markers revealed that the complexity in several of these aquatic taxa could be addressed to recurrent hybridization events and cryptic species diversity. The submerged macrophyte genus Ruppia is one of these aquatic genera with a complex taxonomy due to the absence of clear distinguishable traits and several hybridization events. Two species co-exist throughout Europe, R. maritima and R. spiralis (previously known as R. cirrhosa), but recent molecular studies also found several indications of hybridization, introgression and chloroplast capture between these species. However, the full extent and frequency of hybridization and introgression in this genus has not been studied so far, nor is it clear how these hybrid lineages can co-exist locally with their parental species. In this paper, we wanted to detect whether a single coastal wetland where both species co-exist can act as a Ruppia hybrid zone. As a case study, we chose the Camargue, a Mediterranean coastal wetland that harbors a wide diversity in aquatic habitats, especially in terms of salinity and hydro-regime. We sampled several Ruppia populations within this wetland. To identify each sample and reconstruct the local genetic structure of the two parental species and their hybrids, we used both chloroplast and nuclear microsatellite markers. Afterward, we tested whether different species had different habitat preferences. Our results confirmed that R. maritima and R. spiralis are two strongly divergent species with different reproductive ecologies and different habitat preferences. This prevents frequent hybridization and consequently we could not detect any trace of a recent hybridization event. However, we found several populations of later-generation hybrids, including a population of R. maritima x hybrid backcrosses. The hybrid populations occupy a different habitat and are genetically distinct from their parental species, although they tend to be morphological similar to parental R. maritima. Although local hybridization and introgression in Ruppia is less frequent than we expected, the taxonomy of Ruppia is complicated due to ancient hybridizations and several back-crossings.

Disentangling Sources of Gene Tree Discordance in Phylogenomic Data Sets: Testing Ancient Hybridizations in Amaranthaceae s.l.

Gene tree discordance in large genomic data sets can be caused by evolutionary processes such as incomplete lineage sorting and hybridization, as well as model violation, and errors in data processing, orthology inference, and gene tree estimation. Species tree methods that identify and accommodate all sources of conflict are not available, but a combination of multiple approaches can help tease apart alternative sources of conflict. Here, using a phylotranscriptomic analysis in combination with reference genomes, we test a hypothesis of ancient hybridization events within the plant family Amaranthaceae s.l. that was previously supported by morphological, ecological, and Sanger-based molecular data. The data set included seven genomes and 88 transcriptomes, 17 generated for this study. We examined gene-tree discordance using coalescent-based species trees and network inference, gene tree discordance analyses, site pattern tests of introgression, topology tests, synteny analyses, and simulations. We found that a combination of processes might have generated the high levels of gene tree discordance in the backbone of Amaranthaceae s.l. Furthermore, we found evidence that three consecutive short internal branches produce anomalous trees contributing to the discordance. Overall, our results suggest that Amaranthaceae s.l. might be a product of an ancient and rapid lineage diversification, and remains, and probably will remain, unresolved. This work highlights the potential problems of identifiability associated with the sources of gene tree discordance including, in particular, phylogenetic network methods. Our results also demonstrate the importance of thoroughly testing for multiple sources of conflict in phylogenomic analyses, especially in the context of ancient, rapid radiations. We provide several recommendations for exploring conflicting signals in such situations. [Amaranthaceae; gene tree discordance; hybridization; incomplete lineage sorting; phylogenomics; species network; species tree; transcriptomics.]

Recent hybrids recapitulate ancient hybrid outcomes

Genomic outcomes of hybridization depend on selection and recombination in hybrids. Whether these processes have similar effects on hybrid genome composition in contemporary hybrid zones versus ancient hybrid lineages is unknown. Here we show that patterns of introgression in a contemporary hybrid zone in Lycaeides butterflies predict patterns of ancestry in geographically adjacent, older hybrid populations. We find a particularly striking lack of ancestry from one of the hybridizing taxa, Lycaeides melissa, on the Z chromosome in both the old and contemporary hybrids. The same pattern of reduced L. melissa ancestry on the Z chromosome is seen in two other ancient hybrid lineages. More generally, we find that patterns of ancestry in old or ancient hybrids are remarkably predictable from contemporary hybrids, which suggests selection and recombination affect hybrid genomes in a similar way across disparate time scales and during distinct stages of speciation and species breakdown.

Plastome Evolution in Saxifragaceae and Multiple Plastid Capture Events Involving Heuchera and Tiarella.

Saxifragaceae, a family of over 600 species and approximately 30 genera of herbaceous perennials, is well-known for intergeneric hybridization. Of the main lineages in this family, the Heuchera group represents a valuable model for the analysis of plastid capture and its impact on phylogeny reconstruction. In this study, we investigated plastome evolution across the family, reconstructed the phylogeny of the Heuchera group and examined putative plastid capture between Heuchera and Tiarella. Seven species (11 individuals) representing Tiarella, as well as Mitella and Heuchera, were selected for genome skimming. We assembled the plastomes, and then compared these to six others published for Saxifragaceae; the plastomes were found to be highly similar in overall size, structure, gene order and content. Moreover, ycf15 was lost due to pseudogenization and rpl2 lost its only intron for all the analyzed plastomes. Comparative plastome analysis revealed that size variations of the plastomes are purely ascribed to the length differences of LSC, SSC, and IRs regions. Using nuclear ITS + ETS and the complete plastome, we fully resolved the species relationships of Tiarella, finding that the genus is monophyletic and the Asian species is most closely related to the western North American species. However, the position of the Heuchera species was highly incongruent between nuclear and plastid data. Comparisons of nuclear and plastid phylogenies revealed that multiple plastid capture events have occurred between Heuchera and Tiarella, through putative ancient hybridization. Moreover, we developed numerous molecular markers for Tiarella (e.g., plastid hotspot and polymorphic nuclear SSRs), which will be useful for future studies on the population genetics and phylogeography of this disjunct genus.

Hybrid turbidite-drift channel complexes: An integrated multiscale model

Abstract The interaction of deep-marine bottom currents with episodic, unsteady sediment gravity flows affects global sediment transport, forms climate archives, and controls the evolution of continental slopes. Despite their importance, contradictory hypotheses for reconstructing past flow regimes have arisen from a paucity of studies and the lack of direct monitoring of such hybrid systems. Here, we address this controversy by analyzing deposits, high-resolution seafloor data, and near-bed current measurements from two sites where eastward-flowing gravity flows interact(ed) with northward-flowing bottom currents. Extensive seismic and core data from offshore Tanzania reveal a 1650-m-thick asymmetric hybrid channel levee-drift system, deposited over a period of ∼20 m.y. (Upper Cretaceous to Paleocene). High-resolution modern seafloor data from offshore Mozambique reveal similar asymmetric channel geometries, which are related to northward-flowing near-bed currents with measured velocities of up to 1.4 m/s. Higher sediment accumulation occurs on the downstream flank of channel margins (with respect to bottom currents), with inhibited deposition or scouring on the upstream flank (where velocities are highest). Toes of the drift deposits, consisting of thick laminated muddy siltstone, which progressively step back into the channel axis over time, result in an interfingering relationship with the sandstone-dominated channel fill. Bottom-current flow directions contrast with those of previous models, which lacked direct current measurements or paleoflow indicators. We finally show how large-scale depositional architecture is built through the temporally variable coupling of these two globally important sediment transport processes. Our findings enable more-robust reconstructions of past oceanic circulation and diagnosis of ancient hybrid turbidite-drift systems.

DdRAD analyses reveal a credible phylogenetic relationship of the four main genera of Bambusa-Dendrocalamus-Gigantochloa complex (Poaceae: Bambusoideae)

The Bambusa-Dendrocalamus-Gigantochloa complex (BDG complex) is the most diversified and phylogenetically recalcitrant group of the paleotropical woody bamboos. Species of this complex occur in tropical and subtropical Asia and most of them are of great economic, cultural and ecological value. The lack of resolution achieved through the analyses of previous molecular datasets has long confounded its phylogenetic estimation and generic delimitation. Here, we adopted a ddRAD-seq strategy to investigate phylogenetic relationships of the four main genera (Bambusa, Dendrocalamus, Gigantochloa, and Melocalamus) in the BDG complex. A total of 102 species were sampled, and SNP data were generated. Both MP and ML analyses of the ddRAD-seq data resulted in a well-resolved topology with Gigantochloa and Melocalamus confirmed as monophyletic, and Melocalamus resolved as sister to the rest of the complex. Bambusa and Dendrocalamus were both resolved as paraphyletic. The phylogenetic relationships were mostly supported by morphological evidence including characters of the branch complement, rachilla, lodicules, filaments and stigma. We also generated and assembled complete plastid genomes of 48 representative species. There were conflicts between the plastome and the ddRAD topologies. Our study demonstrated that RAD-seq can be used to reconstruct evolutionary history of lineages such as the bamboos where ancient hybridization and polyploidy play a significant role. The four genera of the BDG complex have a complex evolutionary history which is likely a product of ancient introgression events.

The study of hybridization processes within genus Sparganium L. Subgenus Xanthosparganium holmb. Based on data of next generation sequencing (NGS)

The study represents the results of research of intragenic polymorphism in transcribed spacer ITS1 of the 35S rRNA genes in representatives of subgenus Xanthosparganium genus Sparganium which were obtained by means of locus-specific next generation sequencing on the platform Illumina MiSeq. It was shown that ribotype variations in studied samples generally correspond to the division of this genus into three sections Erecta (subgenus Sparganium), Natantia and Minima (subgenus Xanthosparganium). High level of intragenic polymorphism was revealed in S. hyperboreum, with ribotypes distributed among several groups. Genome of this species includes ribotypes which are typical for other species in subgenus Xanthosparganium. For two investigated S. glomeratum samples, there were no ribotypes similar to such ribotypes in other species of Natantia section. S. glomeratum has got ribotypes identical with S. hyperboreum of Minima section. This feature may be the evidence of ancient intersectional hybridization of these two species. Characteristics of rDNA in S. glomeratum are in favor of putting this species into Minima section. It was suggested that speciation processes within the genus could be based not only on hybridization but also went on in allopatric way. The fist statement is supported by the presence of similar and identical ribotypes in S. emersum, S. longifolium, S. gramineum and S. hyperboreum, the second as it was mentioned by other researchers, is due to close relationship between North American and Eurasian taxa.

Ancient hybridization and mtDNA introgression behind current paternal leakage and heteroplasmy in hybrid zones

Hybridization between heterospecific individuals has been documented as playing a direct role in promoting paternal leakage and mitochondrial heteroplasmy in both natural populations and laboratory conditions, by relaxing the egg-sperm recognition mechanisms. Here, we tested the hypothesis that hybridization can lead to mtDNA heteroplasmy also indirectly via mtDNA introgression. By using a phylogenetic approach, we showed in two reproductively isolated beetle species, Ochthebius quadricollis and O. urbanelliae, that past mtDNA introgression occurred between them in sympatric populations. Then, by developing a multiplex allele-specific PCR assay, we showed the presence of heteroplasmic individuals and argue that their origin was through paternal leakage following mating between mtDNA-introgressed and pure conspecific individuals. Our results highlight that mtDNA introgression can contribute to promote paternal leakage, generating genetic novelty in a way that has been overlooked to date. Furthermore, they highlight that the frequency and distribution of mtDNA heteroplasmy can be deeply underestimated in natural populations, as i) the commonly used PCR-Sanger sequencing approach can fail to detect mitochondrial heteroplasmy, and ii) specific studies aimed at searching for it in populations where mtDNA-introgressed and pure individuals co-occur remain scarce, despite the fact that mtDNA introgression has been widely documented in several taxa and populations.

Plastome phylogenomics, biogeography, and clade diversification of Paris (Melanthiaceae)

BackgroundParis (Melanthiaceae) is an economically important but taxonomically difficult genus, which is unique in angiosperms because some species have extremely large nuclear genomes. Phylogenetic relationships within Paris have long been controversial. Based on complete plastomes and nuclear ribosomal DNA (nrDNA) sequences, this study aims to reconstruct a robust phylogenetic tree and explore historical biogeography and clade diversification in the genus.ResultsAll 29 species currently recognized in Paris were sampled. Whole plastomes and nrDNA sequences were generated by the genome skimming approach. Phylogenetic relationships were reconstructed using the maximum likelihood and Bayesian inference methods. Based on the phylogenetic framework and molecular dating, biogeographic scenarios and historical diversification of Paris were explored. Significant conflicts between plastid and nuclear datasets were identified, and the plastome tree is highly congruent with past interpretations of the morphology. Ancestral area reconstruction indicated that Paris may have originated in northeastern Asia and northern China, and has experienced multiple dispersal and vicariance events during its diversification. The rate of clade diversification has sharply accelerated since the Miocene/Pliocene boundary.ConclusionsOur results provide important insights for clarifying some of the long-standing taxonomic debates in Paris. Cytonuclear discordance may have been caused by ancient and recent hybridizations in the genus. The climatic and geological changes since the late Miocene, such as the intensification of Asian monsoon and the rapid uplift of Qinghai-Tibet Plateau, as well as the climatic fluctuations during the Pleistocene, played essential roles in driving range expansion and radiative diversification in Paris. Our findings challenge the theoretical prediction that large genome sizes may limit speciation.

The complex evolutionary history of apricots: Species divergence, gene flow and multiple domestication events.

Domestication is an excellent model to study diversification and this evolutionary process can be different in perennial plants, such as fruit trees, compared to annual crops. Here, we inferred the history of wild apricot species divergence and of apricot domestication history across Eurasia, with a special focus on Central and Eastern Asia, based on microsatellite markers and approximate Bayesian computation. We significantly extended our previous sampling of apricots in Europe and Central Asia towards Eastern Asia, resulting in a total sample of 271 cultivated samples and 306 wild apricots across Eurasia, mainly Prunus armeniaca and Prunus sibirica, with some Prunus mume and Prunus mandshurica. We recovered wild Chinese species as genetically differentiated clusters, with P.sibirica being divided into two clusters, one possibly resulting from hybridization with P.armeniaca. Central Asia also appeared as a diversification centre of wild apricots. We further revealed at least three domestication events, without bottlenecks, that gave rise to European, Southern Central Asian and Chinese cultivated apricots, with ancient gene flow among them. The domestication event in China possibly resulted from ancient hybridization between wild populations from Central and Eastern Asia. We also detected extensive footprints of recent admixture in all groups of cultivated apricots. Our results thus show that apricot is an excellent model for studying speciation and domestication in long-lived perennial fruit trees.

Diversification of the North American Doellingeria-Eucephalus Clade (Astereae: Asteraceae) Inferred from Molecular and Morphological Evidence

Abstract—Doellingeria and Eucephalus form the earliest-diverging clade of the North American Astereae lineage. Phylogenetic analyses of both nuclear and plastid sequence data show that the Doellingeria-Eucephalus clade consists of two main subclades that differ from current circumscriptions of the two genera. Doellingeria is the sister group to E. elegans, and the Doellingeria + E. elegans subclade in turn is sister to the subclade containing all remaining species of Eucephalus. In the plastid phylogeny, the two subclades are deeply divergent, a pattern that is consistent with an ancient hybridization event involving ancestral species of the Doellingeria-Eucephalus clade and an ancestral taxon of a related North American or South American group. Divergence of the two Doellingeria-Eucephalus subclades may have occurred in association with northward migration from South American ancestors. We combine these two genera under the older of the two names, Doellingeria, and propose 12 new combinations (10 species and two varieties) for all species of Eucephalus.

Introgression drives repeated evolution of winter coat color polymorphism in hares

Changing from summer-brown to winter-white pelage or plumage is a crucial adaptation to seasonal snow in more than 20 mammal and bird species. Many of these species maintain nonwhite winter morphs, locally adapted to less snowy conditions, which may have evolved independently. Mountain hares (Lepus timidus) from Fennoscandia were introduced into the Faroe Islands in 1855. While they were initially winter-white, within ∼65 y all Faroese hares became winter-gray, a morph that occurs in the source population at low frequency. The documented population history makes this a valuable model for understanding the genetic basis and evolution of the seasonal trait polymorphism. Through whole-genome scans of differentiation and single-nucleotide polymorphism (SNP) genotyping, we associated winter coat color polymorphism to the genomic region of the pigmentation gene Agouti, previously linked to introgression-driven winter coat color variation in the snowshoe hare (Lepus americanus). Lower Agouti expression in the skin of winter-gray individuals during the autumn molt suggests that regulatory changes may underlie the color polymorphism. Variation in the associated genomic region shows signatures of a selective sweep in the Faroese population, suggesting that positive selection drove the fixation of the variant after the introduction. Whole-genome analyses of several hare species revealed that the winter-gray variant originated through introgression from a noncolor changing species, in keeping with the history of ancient hybridization between the species. Our findings show the recurrent role of introgression in generating winter coat color variation by repeatedly recruiting the regulatory region of Agouti to modulate seasonal coat color change.

Comparative Transcriptomics Provides Insights into Reticulate and Adaptive Evolution of a Butterfly Radiation.

Butterfly eyes are complex organs that are composed of a diversity of proteins and they play a central role in visual signaling and ultimately, speciation, and adaptation. Here, we utilized the whole eye transcriptome to obtain a more holistic view of the evolution of the butterfly eye while accounting for speciation events that co-occur with ancient hybridization. We sequenced and assembled transcriptomes from adult female eyes of eight species representing all major clades of the Heliconius genus and an additional outgroup species, Dryas iulia. We identified 4,042 orthologous genes shared across all transcriptome data sets and constructed a transcriptome-wide phylogeny, which revealed topological discordance with the mitochondrial phylogenetic tree in the Heliconius pupal mating clade. We then estimated introgression among lineages using additional genome data and found evidence for ancient hybridization leading to the common ancestor of Heliconius hortense and Heliconius clysonymus. We estimated the Ka/Ks ratio for each orthologous cluster and performed further tests to demonstrate genes showing evidence of adaptive protein evolution. Furthermore, we characterized patterns of expression for a subset of these positively selected orthologs using qRT-PCR. Taken together, we identified candidate eye genes that show signatures of adaptive molecular evolution and provide evidence of their expression divergence between species, tissues, and sexes. Our results demonstrate: 1) greater evolutionary changes in younger Heliconius lineages, that is, more positively selected genes in the cydno–melpomene–hecale group as opposed to the sara–hortense–erato group, and 2) suggest an ancient hybridization leading to speciation among Heliconius pupal-mating species.

Resolving evolutionary relationships among six closely related taxa of the horseshoe bats (Rhinolophus) with targeted resequencing data

Recently diverged taxa are often characterised by high rates of introgressive hybridization and incomplete lineage sorting, both of which can complicate phylogenetic reconstructions of species histories. Here we use a sequence capture approach to obtain genome-wide data to resolve the evolutionary relationships, and infer the extent and timescale of hybridization and introgression events, among six recently diverged taxa of the horseshoe bat species complexes Rhinolophus sinicus and R. thomasi. We show that two different methods of species tree reconstruction applied to a set of ~1500 nuclear loci all recover species trees with similar topologies, differing from the previous phylogeny based on two nuclear loci. By comparing the tree topology obtained from the nuclear loci with that inferred from the mitochondrial genome, we observed at least three cases of conflict, each of which likely results from past introgression. Of these, the occurrence of a highly similar mitogenome sequence shared by individuals of two taxa in a sympatric region points to very recent mtDNA introgression. The other cases are characterised by greater divergence and strong phylogeographic structure among putative introgressed individuals and their source populations, and thus likely reflect more ancient hybridization events. These results also suggest that two of the subspecies (R. s. septentrionalis and the undescribed taxon R. s. ssp) are likely to represent full species, warranting full taxonomic descriptions. This work adds a growing number of studies showing the potential problems of relying solely on mitochondrial sequences, or a limited number of loci, to infer phylogenetic relationships among recently diverged taxa.

Ancient Hybridization and Adaptive Introgression of an Invadolysin Gene in Schistosome Parasites.

Introgression among parasite species has the potential to transfer traits of biomedical importance across species boundaries. The parasitic blood fluke Schistosoma haematobium causes urogenital schistosomiasis in humans across sub-Saharan Africa. Hybridization with other schistosome species is assumed to occur commonly, because genetic crosses between S. haematobium and livestock schistosomes, including S. bovis, can be staged in the laboratory, and sequencing of mtDNA and rDNA amplified from microscopic miracidia larvae frequently reveals markers from different species. However, the frequency, direction, age, and genomic consequences of hybridization are unknown. We hatched miracidia from eggs and sequenced the exomes from 96 individual S. haematobium miracidia from infected patients from Niger and the Zanzibar archipelago. These data revealed no evidence for contemporary hybridization between S. bovis and S. haematobium in our samples. However, all Nigerien S. haematobium genomes sampled show hybrid ancestry, with 3.3–8.2% of their nuclear genomes derived from S. bovis, providing evidence of an ancient introgression event that occurred at least 108–613 generations ago. Some S. bovis-derived alleles have spread to high frequency or reached fixation and show strong signatures of directional selection; the strongest signal spans a single gene in the invadolysin gene family (Chr. 4). Our results suggest that S. bovis/S. haematobium hybridization occurs rarely but demonstrate profound consequences of ancient introgression from a livestock parasite into the genome of S. haematobium, the most prevalent schistosome species infecting humans.

The poplar pangenome provides insights into the evolutionary history of the genus

The genus Populus comprises a complex amalgam of ancient and modern species that has become a prime model for evolutionary and taxonomic studies. Here we sequenced the genomes of 10 species from five sections of the genus Populus, identified 71 million genomic variations, and observed new correlations between the single-nucleotide polymorphism–structural variation (SNP–SV) density and indel–SV density to complement the SNP–indel density correlation reported in mammals. Disease resistance genes (R genes) with heterozygous loss-of-function (LOF) were significantly enriched in the 10 species, which increased the diversity of poplar R genes during evolution. Heterozygous LOF mutations in the self-incompatibility genes were closely related to the self-fertilization of poplar, suggestive of genomic control of self-fertilization in dioecious plants. The phylogenetic genome-wide SNPs tree also showed possible ancient hybridization among species in sections Tacamahaca, Aigeiros, and Leucoides. The pangenome resource also provided information for poplar genetics and breeding.

Recombination-Aware Phylogenomics Reveals the Structured Genomic Landscape of Hybridizing Cat Species.

Current phylogenomic approaches implicitly assume that the predominant phylogenetic signal within a genome reflects the true evolutionary history of organisms, without assessing the confounding effects of postspeciation gene flow that can produce a mosaic of phylogenetic signals that interact with recombinational variation. Here, we tested the validity of this assumption with a phylogenomic analysis of 27 species of the cat family, assessing local effects of recombination rate on species tree inference and divergence time estimation across their genomes. We found that the prevailing phylogenetic signal within the autosomes is not always representative of the most probable speciation history, due to ancient hybridization throughout felid evolution. Instead, phylogenetic signal was concentrated within regions of low recombination, and notably enriched within large X chromosome recombination cold spots that exhibited recurrent patterns of strong genetic differentiation and selective sweeps across mammalian orders. By contrast, regions of high recombination were enriched for signatures of ancient gene flow, and these sequences inflated crown-lineage divergence times by ∼40%. We conclude that existing phylogenomic approaches to infer the Tree of Life may be highly misleading without considering the genomic architecture of phylogenetic signal relative to recombination rate and its interplay with historical hybridization.

An ancient hybridization event reconciles mito-nuclear discordance among spiral-horned antelopes

AbstractThe spiral-horned antelopes (genus Tragelaphus) are among the most phenotypically diverse of all large mammals, and evolved in Africa during an adaptive radiation that began in the late Miocene, around 6 million years ago. Tragelaphus was able to exploit the habitat heterogeneity created by Plio-Pleistocene paleoclimatic fluctuations and tectonic processes to eventually occupy almost every habitat type in present day sub-Saharan Africa. The smallest of the spiral-horned antelopes, the bushbuck (T. scriptus), is also widely distributed across Africa, but is genetically divided into polyphyletic Scriptus and Sylvaticus mitochondrial (mt)DNA superlineages that inhabit opposite halves of the continent, suggesting the convergent evolution of independent bushbuck species. In this study, we provide a species tree reconstruction for the genus Tragelaphus and show that Scriptus and Sylvaticus are reciprocally monophyletic at nuclear DNA loci, comprising a single species across its African range. Given that mtDNA will sort into species-specific lineages more quickly than nuclear DNA, only an ancient interspecific hybridization event between a female from a now-extinct Tragelaphus species and a proto-Scriptus bushbuck male can reconcile the mito-nuclear incongruence. This extinct species diverged from the nyala (T. angasii) in the Pliocene about 4.1 million years ago. This study adds to an increasing body of evidence that suggests interspecific hybridization may be more common than previously thought.

An invariants-based method for efficient identification of hybrid species from large-scale genomic data

BackgroundCoalescent-based species tree inference has become widely used in the analysis of genome-scale multilocus and SNP datasets when the goal is inference of a species-level phylogeny. However, numerous evolutionary processes are known to violate the assumptions of a coalescence-only model and complicate inference of the species tree. One such process is hybrid speciation, in which a species shares its ancestry with two distinct species. Although many methods have been proposed to detect hybrid speciation, only a few have considered both hybridization and coalescence in a unified framework, and these are generally limited to the setting in which putative hybrid species must be identified in advance.ResultsHere we propose a method that can examine genome-scale data for a large number of taxa and detect those taxa that may have arisen via hybridization, as well as their potential “parental” taxa. The method is based on a model that considers both coalescence and hybridization together, and uses phylogenetic invariants to construct a test that scales well in terms of computational time for both the number of taxa and the amount of sequence data. We test the method using simulated data for up 20 taxa and 100,000bp, and find that the method accurately identifies both recent and ancient hybrid species in less than 30 s. We apply the method to two empirical datasets, one composed of Sistrurus rattlesnakes for which hybrid speciation is not supported by previous work, and one consisting of several species of Heliconius butterflies for which some evidence of hybrid speciation has been previously found.ConclusionsThe proposed method is powerful for detecting hybridization for both recent and ancient hybridization events. The computations required can be carried out rapidly for a large number of sequences using genome-scale data, and the method is appropriate for both SNP and multilocus data.