Ghost Lineages Highly Influence the Interpretation of Introgression Tests

Testing times: disentangling admixture histories in recent and complex demographies using ancient DNA

Our knowledge of human evolutionary history has been greatly advanced by paleogenomics. Since the 2020s, the study of ancient DNA has increasingly focused on reconstructing the recent past. However, the accuracy of paleogenomic methods in resolving questions of historical and archaeological importance amidst the increased demographic complexity and decreased genetic differentiation remains an open question. We evaluated the performance and behavior of two commonly used methods, qpAdm and the f3-statistic, on admixture inference under a diversity of demographic models and data conditions. We performed two complementary simulation approaches-firstly exploring a wide demographic parameter space under four simple demographic models of varying complexities and configurations using branch-length data from two chromosomes-and secondly, we analyzed a model of Eurasian history composed of 59 populations using whole-genome data modified with ancient DNA conditions such as SNP ascertainment, data missingness, and pseudohaploidization. We observe that population differentiation is the primary factor driving qpAdm performance. Notably, while complex gene flow histories influence which models are classified as plausible, they do not reduce overall performance. Under conditions reflective of the historical period, qpAdm most frequently identifies the true model as plausible among a small candidate set of closely related populations. To increase the utility for resolving fine-scaled hypotheses, we provide a heuristic for further distinguishing between candidate models that incorporates qpAdm model P-values and f3-statistics. Finally, we demonstrate a significant performance increase for qpAdm using whole-genome branch-length f2-statistics, highlighting the potential for improved demographic inference that could be achieved with future advancements in f-statistic estimations.

Introgression and incomplete lineage sorting blurred phylogenetic relationships across the genomes of sclerophyllous oaks from southwest China

Resolving evolutionary relationships among closely related species with interspecific gene flow is challenging. Genome-scale data provide opportunities to clarify complex evolutionary relationships in closely related species and to observe variations in species relationships across the genomes of such species. The Himalayan-Hengduan subalpine oaks have a nearly completely sympatric distribution in southwest China and probably constitute a syngameon. In this study, we mapped resequencing data from different species in this group to the Quercus aquifolioides reference genome to obtain a high-quality filtered single nucleotide polymorphism (SNP) dataset. We also assembled their plastomes. We reconstructed their phylogenetic relationships, explored the level and pattern of introgression among these species and investigated gene tree variation in the genomes of these species using sliding windows. The same or closely related plastomes were found to be shared extensively among different species within a specific geographical area. Phylogenomic analyses of genome-wide SNP data found that most oaks in the Himalayan-Hengduan subalpine clade showed genetic coherence, but several species were found to be connected by introgression. The gene trees obtained using sliding windows showed that the phylogenetic relationships in the genomes of oaks are highly heterogeneous and therefore highly obscured. Our study found that all the oaks of the Himalayan-Hengduan subalpine clade from southwest China form a syngameon. The obscured phylogenetic relationships observed empirically across the genome are best explained by interspecific gene flow in conjunction with incomplete lineage sorting.

Artifactual Orthologs and the Need for Diligent Data Exploration in Complex Phylogenomic Datasets: A Museomic Case Study from the Andean Flora

The Andes mountains of western South America are a globally important biodiversity hotspot, yet there is a paucity of resolved phylogenies for plant clades from this region. Filling an important gap in our understanding of the World's richest flora, we present the first phylogeny of Freziera (Pentaphylacaceae), an Andean-centered, cloud forest radiation. Our dataset was obtained via hybrid-enriched target sequence capture of Angiosperms353 universal loci for 50 of the ca. 75 spp., obtained almost entirely from herbarium specimens. We identify high phylogenomic complexity in Freziera, including the presence of data artifacts. Via by-eye observation of gene trees, detailed examination of warnings from recently improved assembly pipelines, and gene tree filtering, we identified that artifactual orthologs (i.e., the presence of only one copy of a multicopy gene due to differential assembly) were an important source of gene tree heterogeneity that had a negative impact on phylogenetic inference and support. These artifactual orthologs may be common in plant phylogenomic datasets, where multiple instances of genome duplication are common. After accounting for artifactual orthologs as source of gene tree error, we identified a significant, but nonspecific signal of introgression using Patterson's D and f4 statistics. Despite phylogenomic complexity, we were able to resolve Freziera into 9 well-supported subclades whose evolution has been shaped by multiple evolutionary processes, including incomplete lineage sorting, historical gene flow, and gene duplication. Our results highlight the complexities of plant phylogenomics, which are heightened in Andean radiations, and show the impact of filtering data processing artifacts and standard filtering approaches on phylogenetic inference.

Ancient and Recent Hybridization in the Oreochromis Cichlid Fishes

Cichlid fishes of the genus Oreochromis (tilapia) are among the most important fish for inland capture fisheries and global aquaculture. Deliberate introductions of non-native species for fisheries improvement and accidental escapees from farms have resulted in admixture with indigenous species. Such hybridization may be detrimental to native biodiversity, potentially leading to genomic homogenization of populations and the loss of important genetic material associated with local adaptation. By contrast, introgression may fuel diversification when combined with ecological opportunity, by supplying novel genetic combinations. To date, the role of introgression in the evolutionary history of tilapia has not been explored. Here we studied both ancient and recent hybridization in tilapia, using whole genome resequencing of 575 individuals from 23 species. We focused on Tanzania, a natural hotspot of tilapia diversity, and a country where hybridization between exotic and native species in the natural environment has been previously reported. We reconstruct the first genome-scale phylogeny of the genus and reveal prevalent ancient gene flow across the Oreochromis phylogeny. This has likely resulted in the hybrid speciation of one species, O. chungruruensis. We identify multiple cases of recent hybridization between native and introduced species in the wild, linked to the use of non-native species in both capture fisheries improvement and aquaculture. This has potential implications for both conservation of wild populations and the development of the global tilapia aquaculture industry.

Camelus knoblochi genome reveals the complex evolutionary history of Old World camels

Extant Old World camels (genus Camelus) contributed to the economic and cultural exchanges between the East and West for thousands of years.1,2 Although many remains have been unearthed,3,4,5 we know neither whether the prevalent hybridization observed between extant Camelus species2,6,7 also occurred between extinct lineages and the ancestors of extant Camelus species nor why some populations became extinct while others survived. To investigate these questions, we generated paleogenomic and stable isotope data from an extinct two-humped camel species, Camelus knoblochi. We find that in the mitochondrial phylogeny, all C. knoblochi form a paraphyletic group that nests within the diversity of modern, wild two-humped camels (Camelus ferus). In contrast, they are clearly distinguished from both wild and domesticated (Camelus bactrianus) two-humped camels on the nuclear level. Moreover, the divergence pattern of the three camel species approximates a trifurcation, because the most common topology is only slightly more frequent than the two other possible topologies. This mito-nuclear phylogenetic discordance likely arose due to interspecific gene flow between all three species, suggesting that interspecific hybridization is not exclusive to modern camels but a recurrent phenomenon throughout the evolutionary history of the genus Camelus. These results suggest that the genomic complexity of Old World camels' evolutionary history is underestimated when considering data from only modern species. Finally, we find that C. knoblochi populations began declining prior to the last glacial maximum and, by integrating palaeoecological evidence and stable isotope data, suggest that this was likely due to failure to adapt to a changing environment.

Introgression Underlies Phylogenetic Uncertainty But Not Parallel Plumage Evolution in a Recent Songbird Radiation

Instances of parallel phenotypic evolution offer great opportunities to understand the evolutionary processes underlying phenotypic changes. However, confirming parallel phenotypic evolution and studying its causes requires a robust phylogenetic framework. One such example is the "black-and-white wagtails," a group of 5 species in the songbird genus Motacilla: 1 species, Motacilla alba, shows wide intra-specific plumage variation, while the 4r others form 2 pairs of very similar-looking species (M. aguimp + M. samveasnae and M. grandis + M. maderaspatensis, respectively). However, the 2 species in each of these pairs were not recovered as sisters in previous phylogenetic inferences. Their relationships varied depending on the markers used, suggesting that gene tree heterogeneity might have hampered accurate phylogenetic inference. Here, we use whole genome resequencing data to explore the phylogenetic relationships within this group, with a special emphasis on characterizing the extent of gene tree heterogeneity and its underlying causes. We first used multispecies coalescent methods to generate a "complete evidence" phylogenetic hypothesis based on genome-wide variants, while accounting for incomplete lineage sorting (ILS) and introgression. We then investigated the variation in phylogenetic signal across the genome to quantify the extent of discordance across genomic regions and test its underlying causes. We found that wagtail genomes are mosaics of regions supporting variable genealogies, because of ILS and inter-specific introgression. The most common topology across the genome, supporting M. alba and M. aguimp as sister species, appears to be influenced by ancient introgression. Additionally, we inferred another ancient introgression event, between M. alba and M. grandis. By combining results from multiple analyses, we propose a phylogenetic network for the black-and-white wagtails that confirms that similar phenotypes evolved in non-sister lineages, supporting parallel plumage evolution. Furthermore, the inferred reticulations do not connect species with similar plumage coloration, suggesting that introgression does not underlie parallel plumage evolution in this group. Our results demonstrate the importance of investing genome-wide patterns of gene tree heterogeneity to help understand the mechanisms underlying phenotypic evolution. [Gene tree heterogeneity; incomplete lineage sorting; introgression; parallel evolution; phylogenomics; plumage evolution; wagtails.].

Detection of Ghost Introgression Requires Exploiting Topological and Branch Length Information

In recent years, the study of hybridization and introgression has made significant progress, with ghost introgression-the transfer of genetic material from extinct or unsampled lineages to extant species-emerging as a key area for research. Accurately identifying ghost introgression, however, presents a challenge. To address this issue, we focused on simple cases involving 3 species with a known phylogenetic tree. Using mathematical analyses and simulations, we evaluated the performance of popular phylogenetic methods, including HyDe and PhyloNet/MPL, and the full-likelihood method, Bayesian Phylogenetics and Phylogeography (BPP), in detecting ghost introgression. Our findings suggest that heuristic approaches relying on site-pattern counts or gene-tree topologies struggle to differentiate ghost introgression from introgression between sampled non-sister species, frequently leading to incorrect identification of donor and recipient species. The full-likelihood method BPP uses multilocus sequence alignments directly-hence taking into account both gene-tree topologies and branch lengths, by contrast, is capable of detecting ghost introgression in phylogenomic datasets. We analyzed a real-world phylogenomic dataset of 14 species of Jaltomata (Solanaceae) to showcase the potential of full-likelihood methods for accurate inference of introgression.

Phylogenomic discordance is driven by wide-spread introgression and incomplete lineage sorting during rapid species diversification within rattlesnakes (Viperidae: Crotalus and Sistrurus)

Phylogenomics allows us to uncover the historical signal of evolutionary processes through time and estimate phylogenetic networks accounting for these signals. Insight from genome-wide data further allows us to pinpoint the contributions to phylogenetic signal from hybridization, introgression, and ancestral polymorphism across the genome. Here we focus on how these processes have contributed to phylogenetic discordance among rattlesnakes (genera Crotalus and Sistrurus), a group for which there are numerous conflicting phylogenetic hypotheses based on a diverse array of molecular datasets and analytical methods. We address the instability of the rattlesnake phylogeny using genomic data generated from transcriptomes sampled from nearly all known species. These genomic data, analyzed with coalescent and network-based approaches, reveal numerous instances of rapid speciation where individual gene trees conflict with the species tree. Moreover, the evolutionary history of rattlesnakes is dominated by incomplete speciation and frequent hybridization, both of which have likely influenced past interpretations of phylogeny. We present a new framework in which the evolutionary relationships of this group can only be understood in light of genome-wide data and network-based analytical methods. Our data suggest that network radiations, like seen within the rattlesnakes, can only be understood in a phylogenomic context, necessitating similar approaches in our attempts to understand evolutionary history in other rapidly radiating species.

Phylogenomics of Psammodynastes and Buhoma (Elapoidea: Serpentes), with the description of a new Asian snake family

Asian mock vipers of the genus Psammodynastes and African forest snakes of the genus Buhoma are two genera belonging to the snake superfamily Elapoidea. The phylogenetic placements of Psammodynastes and Buhoma within Elapoidea has been extremely unstable which has resulted in their uncertain and debated taxonomy. We used ultraconserved elements and traditional nuclear and mitochondrial markers to infer the phylogenetic relationships of these two genera with other elapoids. Psammodynastes, for which a reference genome has been sequenced, were found, with strong branch support, to be a relatively early diverging split within Elapoidea that is sister to a clade consisting of Elapidae, Micrelapidae and Lamprophiidae. Hence, we allocate Psammodynastes to its own family, Psammodynastidae new family. However, the phylogenetic position of Buhoma could not be resolved with a high degree of confidence. Attempts to identify the possible sources of conflict in the rapid radiation of elapoid snakes suggest that both hybridisation/introgression during the rapid diversification, including possible ghost introgression, as well as incomplete lineage sorting likely have had a confounding role. The usual practice of combining mitochondrial loci with nuclear genomic data appears to mislead phylogeny reconstructions in rapid radiation scenarios, especially in the absence of genome scale data.

Anomalous networks under the multispecies coalescent: theory and prevalence

Reticulations in a phylogenetic network represent processes such as gene flow, admixture, recombination and hybrid speciation. Extending definitions from the tree setting, an anomalous network is one in which some unrooted tree topology displayed in the network appears in gene trees with a lower frequency than a tree not displayed in the network. We investigate anomalous networks under the Network Multispecies Coalescent Model with possible correlated inheritance at reticulations. Focusing on subsets of 4 taxa, we describe a new algorithm to calculate quartet concordance factors on networks of any level, faster than previous algorithms because of its focus on 4 taxa. We then study topological properties required for a 4-taxon network to be anomalous, uncovering the key role of [Formula: see text]-cycles: cycles of 3 edges parent to a sister group of 2 taxa. Under the model of common inheritance, that is, when each gene tree coalesces within a species tree displayed in the network, we prove that 4-taxon networks are never anomalous. Under independent and various levels of correlated inheritance, we use simulations under realistic parameters to quantify the prevalence of anomalous 4-taxon networks, finding that truly anomalous networks are rare. At the same time, however, we find a significant fraction of networks close enough to the anomaly zone to appear anomalous, when considering the quartet concordance factors observed from a few hundred genes. These apparent anomalies may challenge network inference methods.

Reconstructing cell lineage trees with genomic barcoding: approaches and applications

In multicellular organisms, developmental history of cell divisions and functional annotation of terminal cells can be organized into a cell lineage tree (CLT). The reconstruction of the CLT has long been a major goal in developmental biology and other related fields. Recent technological advancements, especially those in editable genomic barcodes and single-cell high-throughput sequencing, have sparked a new wave of experimental methods for reconstructing CLTs. Here we review the existing experimental approaches to the reconstruction of CLT, which are broadly categorized as either image-based or DNA barcode-based methods. In addition, we present a summary of the related literature based on the biological insight provided by the obtained CLTs. Moreover, we discuss the challenges that will arise as more and better CLT data become available in the near future. Genomic barcoding-based CLT reconstructions and analyses, due to their wide applicability and high scalability, offer the potential for novel biological discoveries, especially those related to general and systemic properties of the developmental process.

Summary Tests of Introgression Are Highly Sensitive to Rate Variation Across Lineages

The evolutionary implications and frequency of hybridization and introgression are increasingly being recognized across the tree of life. To detect hybridization from multi-locus and genome-wide sequence data, a popular class of methods are based on summary statistics from subsets of 3 or 4 taxa. However, these methods often carry the assumption of a constant substitution rate across lineages and genes, which is commonly violated in many groups. In this work, we quantify the effects of rate variation on the D test (also known as ABBA-BABA test), the D3 test, and HyDe. All 3 tests are used widely across a range of taxonomic groups, in part because they are very fast to compute. We consider rate variation across species lineages, across genes, their lineage-by-gene interaction, and rate variation across gene-tree edges. We simulated species networks according to a birth-death-hybridization process, so as to capture a range of realistic species phylogenies. For all 3 methods tested, we found a marked increase in the false discovery of reticulation (type-1 error rate) when there is rate variation across species lineages. The D3 test was the most sensitive, with around 80% type-1 error, such that D3 appears to more sensitive to a departure from the clock than to the presence of reticulation. For all 3 tests, the power to detect hybridization events decreased as the number of hybridization events increased, indicating that multiple hybridization events can obscure one another if they occur within a small subset of taxa. Our study highlights the need to consider rate variation when using site-based summary statistics, and points to the advantages of methods that do not require assumptions on evolutionary rates across lineages or across genes.

Speciation-by-Extinction

Extinction is a dominant force shaping patterns of biodiversity through time; however its role as a catalyst of speciation through its interaction with intraspecific variation has been overlooked. Here, we synthesize ideas alluded to by Darwin and others into the model of "speciation-by-extinction" in which speciation results from the extinction of intermediate populations within a single geographically variable species. We explore the properties and distinguishing features of speciation-by-extinction with respect to other established speciation models. We demonstrate its plausibility by showing that the experimental extinction of populations within variable species can result in speciation. The prerequisites for speciation-by-extinction, geographically structured intraspecific variation and local extinction, are ubiquitous in nature. We propose that speciation-by-extinction may be a prevalent, but underappreciated, speciation mechanism.

Genomic divergence and introgression between cryptic species of a widespread North American songbird

Analysis of genomic variation among related populations can sometimes reveal distinct species that were previously undescribed due to similar morphological appearances, and close examination of such cases can provide much insight regarding speciation. Genomic data can also reveal the role of reticulate evolution in differentiation and speciation. White-breasted nuthatches (Sitta carolinensis) are widely distributed North American songbirds that are currently classified as a single species but have been suspected to represent a case of cryptic speciation. Previous genetic analyses suggested four divergent groups, but it was unclear whether these represented multiple reproductively isolated species. Using extensive genomic sampling of over 350 white-breasted nuthatches from across North America and a new chromosome-level reference genome, we asked if white-breasted nuthatches are comprised of multiple species and whether introgression has occurred between divergent populations. Genomic variation of over 300,000 loci revealed four highly differentiated populations (Pacific, n = 45; Eastern, n = 23; Rocky Mountains North, n = 138; and Rocky Mountains South, n = 150) with geographic ranges that are adjacent. We observed a moderate degree of admixture between Rocky Mountain populations but only a small number of hybrids between the Rockies and the Eastern population. The rarity of hybrids together with high levels of differentiation between populations is supportive of populations having some level of reproductive isolation. Between populations, we show evidence for introgression from a divergent ghost lineage of white-breasted nuthatches into the Rocky Mountains South population, which is otherwise closely related to Rocky Mountains North. We conclude that white-breasted nuthatches are best considered at least three species and that ghost lineage introgression has contributed to differentiation between the two Rocky Mountain populations. White-breasted nuthatches provide a dramatic case of morphological similarity despite high genomic differentiation, and the varying levels of reproductive isolation among the four groups provide an example of the speciation continuum.

Testing Times: Challenges in Disentangling Admixture Histories in Recent and Complex Demographies

Paleogenomics has expanded our knowledge of human evolutionary history. Since the 2020s, the study of ancient DNA has increased its focus on reconstructing the recent past. However, the accuracy of paleogenomic methods in answering questions of historical and archaeological importance amidst the increased demographic complexity and decreased genetic differentiation within the historical period remains an open question. We used two simulation approaches to evaluate the limitations and behavior of commonly used methods, qpAdm and the f3-statistic, on admixture inference. The first is based on branch-length data simulated from four simple demographic models of varying complexities and configurations. The second, an analysis of Eurasian history composed of 59 populations using whole-genome data modified with ancient DNA conditions such as SNP ascertainment, data missingness, and pseudo-haploidization. We show that under conditions resembling historical populations, qpAdm can identify a small candidate set of true sources and populations closely related to them. However, in typical ancient DNA conditions, qpAdm is unable to further distinguish between them, limiting its utility for resolving fine-scaled hypotheses. Notably, we find that complex gene-flow histories generally lead to improvements in the performance of qpAdm and observe no bias in the estimation of admixture weights. We offer a heuristic for admixture inference that incorporates admixture weight estimate and P-values of qpAdm models, and f3-statistics to enhance the power to distinguish between multiple plausible candidates. Finally, we highlight the future potential of qpAdm through whole-genome branch-length f2-statistics, demonstrating the improved demographic inference that could be achieved with advancements in f-statistic estimations.

False discovery rates of qpAdm-based screens for genetic admixture

Although a broad range of methods exists for reconstructing population history from genome-wide single nucleotide polymorphism data, just a few methods gained popularity in archaeogenetics: principal component analysis (PCA); ADMIXTURE, an algorithm that models individuals as mixtures of multiple ancestral sources represented by actual or inferred populations; formal tests for admixture such as f3-statistics and D/f4-statistics; and qpAdm, a tool for fitting two-component and more complex admixture models to groups or individuals. Despite their popularity in archaeogenetics, which is explained by modest computational requirements and ability to analyze data of various types and qualities, protocols relying on qpAdm that screen numerous alternative models of varying complexity and find "fitting" models (often considering both estimated admixture proportions and p-values as a composite criterion of model fit) remain untested on complex simulated population histories in the form of admixture graphs of random topology. We analyzed genotype data extracted from such simulations and tested various types of high-throughput qpAdm protocols ("rotating" and "non-rotating", with or without temporal stratification of target groups and proxy ancestry sources, and with or without a "model competition" step). We caution that high-throughput qpAdm protocols may be inappropriate for exploratory analyses in poorly studied regions/periods since their false discovery rates varied between 12% and 68% depending on the details of the protocol and on the amount and quality of simulated data (i.e., >12% of fitting two-way admixture models imply gene flows that were not simulated). We demonstrate that for reducing false discovery rates of qpAdm protocols to nearly 0% it is advisable to use large SNP sets with low missing data rates, the rotating qpAdm protocol with a strictly enforced rule that target groups do not pre-date their proxy sources, and an unsupervised ADMIXTURE analysis as a way to verify feasible qpAdm models. Our study has a number of limitations: for instance, these recommendations depend on the assumption that the underlying genetic history is a complex admixture graph and not a stepping-stone model.

Strong postmating reproductive isolation in Mimulus section Eunanus

Postmating reproductive isolation can help maintain species boundaries when premating barriers to reproduction are incomplete. The strength and identity of postmating reproductive barriers are highly variable among diverging species, leading to questions about their genetic basis and evolutionary drivers. These questions have been tackled in model systems but are less often addressed with broader phylogenetic resolution. In this study we analyse patterns of genetic divergence alongside direct measures of postmating reproductive barriers in an overlooked group of sympatric species within the model monkeyflower genus, Mimulus. Within this Mimulus brevipes species group, we find substantial divergence among species, including a cryptic genetic lineage. However, rampant gene discordance and ancient signals of introgression suggest a complex history of divergence. In addition, we find multiple strong postmating barriers, including postmating prezygotic isolation, hybrid seed inviability and hybrid male sterility. M. brevipes and M. fremontii have substantial but incomplete postmating isolation. For all other tested species pairs, we find essentially complete postmating isolation. Hybrid seed inviability appears linked to differences in seed size, providing a window into possible developmental mechanisms underlying this reproductive barrier. While geographic proximity and incomplete mating isolation may have allowed gene flow within this group in the distant past, strong postmating reproductive barriers today have likely played a key role in preventing ongoing introgression. By producing foundational information about reproductive isolation and genomic divergence in this understudied group, we add new diversity and phylogenetic resolution to our understanding of the mechanisms of plant speciation.

Defining eukaryotes to dissect eukaryogenesis

The origin of eukaryotes is among the most contentious debates in evolutionary biology, attracting multiple seemingly incompatible theories seeking to explain the sequence in which eukaryotic characteristics were acquired. Much of the controversy arises from differing views on the defining characteristics of eukaryotes. We argue that eukaryotes should be defined phylogenetically, and that doing so clarifies where competing hypotheses of eukaryogenesis agree and how we may test among aspects of disagreement. Some hypotheses make predictions about the phylogenetic origins of eukaryotic genes and are distinguishable on that basis. However, other hypotheses differ only in the order of key evolutionary steps, like mitochondrial endosymbiosis and nuclear assembly, which cannot currently be distinguished phylogenetically. Stages within eukaryogenesis may be made identifiable through the absolute dating of gene duplicates that map to eukaryotic traits, such as in genes of host or mitochondrial origin that duplicated and diverged functionally prior to emergence of the last eukaryotic common ancestor. In this way, it may finally be possible to distinguish heat from light in the debate over eukaryogenesis.

Interploidy Introgression Shaped Adaptation during the Origin and Domestication History of Brassica napus

Polyploidy is recurrent across the tree of life and known as an evolutionary driving force in plant diversification and crop domestication. How polyploid plants adapt to various habitats has been a fundamental question that remained largely unanswered. Brassica napus is a major crop cultivated worldwide, resulting from allopolyploidy between unknown accessions of diploid B. rapa and B. oleracea. Here, we used whole-genome resequencing data of accessions representing the majority of morphotypes and ecotypes from the species B. rapa, B. oleracea, and B. napus to investigate the role of polyploidy during domestication. To do so, we first reconstructed the phylogenetic history of B. napus, which supported the hypothesis that the emergence of B. napus derived from the hybridization of European turnip of B. rapa and wild B. oleracea. These analyses also showed that morphotypes of swede and Siberian kale (used as vegetable and fodder) were domesticated before rapeseed (oil crop). We next observed that frequent interploidy introgressions from sympatric diploids were prominent throughout the domestication history of B. napus. Introgressed genomic regions were shown to increase the overall genetic diversity and tend to be localized in regions of high recombination. We detected numerous candidate adaptive introgressed regions and found evidence that some of the genes in these regions contributed to phenotypic diversification and adaptation of different morphotypes. Overall, our results shed light on the origin and domestication of B. napus and demonstrate interploidy introgression as an important mechanism that fuels rapid diversification in polyploid species.

ANOMALOUS NETWORKS UNDER THE MULTISPECIES COALESCENT: THEORY AND PREVALENCE

Reticulations in a phylogenetic network represent processes such as gene flow, admixture, recombination and hybrid speciation. Extending definitions from the tree setting, an anomalous network is one in which some unrooted tree topology displayed in the network appears in gene trees with a lower frequency than a tree not displayed in the network. We investigate anomalous networks under the Network Multispecies Coalescent Model with possible correlated inheritance at reticulations. Focusing on subsets of 4 taxa, we describe a new algorithm to calculate quartet concordance factors on networks of any level, faster than previous algorithms because of its focus on 4 taxa. We then study topological properties required for a 4-taxon network to be anomalous, uncovering the key role of 32-cycles: cycles of 3 edges parent to a sister group of 2 taxa. Under the model of common inheritance, that is, when each gene tree coalesces within a species tree displayed in the network, we prove that 4-taxon networks are never anomalous. Under independent and various levels of correlated inheritance, we use simulations under realistic parameters to quantify the prevalence of anomalous 4-taxon networks, finding that truly anomalous networks are rare. At the same time, however, we find a significant fraction of networks close enough to the anomaly zone to appear anomalous, when considering the quartet concordance factors observed from a few hundred genes. These apparent anomalies may challenge network inference methods.

Ghost introgression facilitates genomic divergence of a sympatric cryptic lineage in Cycas revoluta

A cryptic lineage is a genetically diverged but morphologically unrecognized variant of a known species. Clarifying cryptic lineage evolution is essential for quantifying species diversity. In sympatric cryptic lineage divergence compared with allopatric divergence, the forces of divergent selection and mating patterns override geographical isolation. Introgression, by supplying preadapted or neutral standing genetic variations, can promote sympatric cryptic lineage divergence via selection. However, most studies concentrated on extant species introgression, ignoring the genetic legacy of introgression from extinct or unsampled lineages ("ghost introgression"). Cycads are an ideal plant for studying the influence of ghost introgression because of their common interspecific gene flow and past high extinction rate. Here, we utilized reference-based ddRADseq to clarify the role of ghost introgression in the evolution of a previously identified sympatric cryptic lineage in Cycas revoluta. After re-evaluating the evolutionary independency of cryptic lineages, the group-wise diverged single-nucleotide polymorphisms among sympatric and allopatric lineages were compared and functionally annotated. Next, we employed an approximate Bayesian computation method for hypothesis testing to clarify the cryptic lineage evolution and ghost introgression effect. SNPs with the genomic signatures of ghost introgression were further annotated. Our results reconfirmed the evolutionary independency of cryptic lineage among C. revoluta and demonstrated that ghost introgression to the noncryptic lineage facilitated their divergence. Gene function related to heat stress and disease resistance implied ecological adaptation of the main extant populations of C. revoluta.

Hybridisation has shaped a recent radiation of grass-feeding aphids

BACKGROUND

Aphids are common crop pests. These insects reproduce by facultative parthenogenesis involving several rounds of clonal reproduction interspersed with an occasional sexual cycle. Furthermore, clonal aphids give birth to live young that are already pregnant. These qualities enable rapid population growth and have facilitated the colonisation of crops globally. In several cases, so-called "super clones" have come to dominate agricultural systems. However, the extent to which the sexual stage of the aphid life cycle has shaped global pest populations has remained unclear, as have the origins of successful lineages. Here, we used chromosome-scale genome assemblies to disentangle the evolution of two global pests of cereals-the English (Sitobion avenae) and Indian (Sitobion miscanthi) grain aphids.

RESULTS

Genome-wide divergence between S. avenae and S. miscanthi is low. Moreover, comparison of haplotype-resolved assemblies revealed that the S. miscanthi isolate used for genome sequencing is likely a hybrid, with one of its diploid genome copies closely related to S. avenae (~ 0.5% divergence) and the other substantially more divergent (> 1%). Population genomics analyses of UK and China grain aphids showed that S. avenae and S. miscanthi are part of a cryptic species complex with many highly differentiated lineages that predate the origins of agriculture. The complex consists of hybrid lineages that display a tangled history of hybridisation and genetic introgression.

CONCLUSIONS

Our analyses reveal that hybridisation has substantially contributed to grain aphid diversity, and hence, to the evolutionary potential of this important pest species. Furthermore, we propose that aphids are particularly well placed to exploit hybridisation events via the rapid propagation of live-born "frozen hybrids" via asexual reproduction, increasing the likelihood of hybrid lineage formation.

Echoes of ancient introgression punctuate stable genomic lineages in the evolution of figs

Studies investigating the evolution of flowering plants have long focused on isolating mechanisms such as pollinator specificity. Some recent studies have proposed a role for introgressive hybridization between species, recognizing that isolating processes such as pollinator specialization may not be complete barriers to hybridization. Occasional hybridization may therefore lead to distinct yet reproductively connected lineages. We investigate the balance between introgression and reproductive isolation in a diverse clade using a densely sampled phylogenomic study of fig trees (Ficus, Moraceae). Codiversification with specialized pollinating wasps (Agaonidae) is recognized as a major engine of fig diversity, leading to about 850 species. Nevertheless, some studies have focused on the importance of hybridization in Ficus, highlighting the consequences of pollinator sharing. Here, we employ dense taxon sampling (520 species) throughout Moraceae and 1,751 loci to investigate phylogenetic relationships and the prevalence of introgression among species throughout the history of Ficus. We present a well-resolved phylogenomic backbone for Ficus, providing a solid foundation for an updated classification. Our results paint a picture of phylogenetically stable evolution within lineages punctuated by occasional local introgression events likely mediated by local pollinator sharing, illustrated by clear cases of cytoplasmic introgression that have been nearly drowned out of the nuclear genome through subsequent lineage fidelity. The phylogenetic history of figs thus highlights that while hybridization is an important process in plant evolution, the mere ability of species to hybridize locally does not necessarily translate into ongoing introgression between distant lineages, particularly in the presence of obligate plant-pollinator relationships.

Phylogenomics reveals widespread hybridization and polyploidization in Henckelia (Gesneriaceae)

BACKGROUND AND AIMS

Hybridization has long been recognized as an important process for plant evolution and is often accompanied by polyploidization, another prominent force in generating biodiversity. Despite its pivotal importance in evolution, the actual prevalence and distribution of hybridization across the tree of life remain unclear.

METHODS

We used whole-genome shotgun (WGS) sequencing and cytological data to investigate the evolutionary history of Henckelia, a large genus in the family Gesneriaceae with a high frequency of suspected hybridization and polyploidization events. We generated WGS sequencing data at about 10× coverage for 26 Chinese Henckelia species plus one Sri Lankan species. To untangle the hybridization history, we separately extracted whole plastomes and thousands of single-copy nuclear genes from the sequencing data, and reconstructed phylogenies based on both nuclear and plastid data. We also explored sources of both genealogical and cytonuclear conflicts and identified signals of hybridization and introgression within our phylogenomic dataset using several statistical methods. Additionally, to test the polyploidization history, we evaluated chromosome counts for 45 populations of the 27 Henckelia species studied.

KEY RESULTS

We obtained well-supported phylogenetic relationships using both concatenation- and coalescent-based methods. However, the nuclear phylogenies were highly inconsistent with the plastid phylogeny, and we observed intensive discordance among nuclear gene trees. Further analyses suggested that both incomplete lineage sorting and gene flow contributed to the observed cytonuclear and genealogical discordance. Our analyses of introgression and phylogenetic networks revealed a complex history of hybridization within the genus Henckelia. In addition, based on chromosome counts for 27 Henckelia species, we found independent polyploidization events occurred within Henckelia after different hybridization events.

CONCLUSIONS

Our findings demonstrated that hybridization and polyploidization are common in Henckelia. Furthermore, our results revealed that H. oblongifolia is not a member of the redefined Henckelia and they suggested several other taxonomic treatments in this genus.

A community genomics approach to natural hybridization

Hybridization is a complicated, oft-misunderstood process. Once deemed unnatural and uncommon, hybridization is now recognized as ubiquitous among species. But hybridization rates within and among communities are poorly understood despite the relevance to ecology, evolution and conservation. To clarify, we examined hybridization across 75 freshwater fish communities within the Ozarks of the North American Interior Highlands (USA) by single nucleotide polymorphism (SNP) genotyping 33 species (N = 2865 individuals; double-digest restriction site-associated DNA sequencing (ddRAD)). We found evidence of hybridization (70 putative hybrids; 2.4% of individuals) among 18 species-pairs involving 73% (24/33) of study species, with the majority being concentrated within one family (Leuciscidae/minnows; 15 species; 66 hybrids). Interspecific genetic exchange-or introgression-was evident from 24 backcrossed individuals (10/18 species-pairs). Hybrids occurred within 42 of 75 communities (56%). Four selected environmental variables (species richness, protected area extent, precipitation (May and annually)) exhibited 73-78% accuracy in predicting hybrid occurrence via random forest classification. Our community-level assessment identified hybridization as spatially widespread and environmentally dependent (albeit predominantly within one diverse, omnipresent family). Our approach provides a more holistic survey of natural hybridization by testing a wide range of species-pairs, thus contrasting with more conventional evaluations.

Exploring Conflicts in Whole Genome Phylogenetics: A Case Study Within Manakins (Aves: Pipridae)

Some phylogenetic problems remain unresolved even when large amounts of sequence data are analyzed and methods that accommodate processes such as incomplete lineage sorting are employed. In addition to investigating biological sources of phylogenetic incongruence, it is also important to reduce noise in the phylogenomic dataset by using appropriate filtering approach that addresses gene tree estimation errors. We present the results of a case study in manakins, focusing on the very difficult clade comprising the genera Antilophia and Chiroxiphia. Previous studies suggest that Antilophia is nested within Chiroxiphia, though relationships among Antilophia+Chiroxiphia species have been highly unstable. We extracted more than 11,000 loci (ultra-conserved elements and introns) from whole genomes and conducted analyses using concatenation and multispecies coalescent methods. Topologies resulting from analyses using all loci differed depending on the data type and analytical method, with 2 clades (Antilophia+Chiroxiphia and Manacus+Pipra+Machaeopterus) in the manakin tree showing incongruent results. We hypothesized that gene trees that conflicted with a long coalescent branch (e.g., the branch uniting Antilophia+Chiroxiphia) might be enriched for cases of gene tree estimation error, so we conducted analyses that either constrained those gene trees to include monophyly of Antilophia+Chiroxiphia or excluded these loci. While constraining trees reduced some incongruence, excluding the trees led to completely congruent species trees, regardless of the data type or model of sequence evolution used. We found that a suite of gene metrics (most importantly the number of informative sites and likelihood of intralocus recombination) collectively explained the loci that resulted in non-monophyly of Antilophia+Chiroxiphia. We also found evidence for introgression that may have contributed to the discordant topologies we observe in Antilophia+Chiroxiphia and led to deviations from expectations given the multispecies coalescent model. Our study highlights the importance of identifying factors that can obscure phylogenetic signal when dealing with recalcitrant phylogenetic problems, such as gene tree estimation error, incomplete lineage sorting, and reticulation events. [Birds; c-gene; data type; gene estimation error; model fit; multispecies coalescent; phylogenomics; reticulation].

Deep reticulation: the long legacy of hybridization in vascular plant evolution

Hybridization has long been recognized as a fundamental evolutionary process in plants but, until recently, our understanding of its phylogenetic distribution and biological significance across deep evolutionary scales has been largely obscure. Over the past decade, genomic and phylogenomic datasets have revealed, perhaps not surprisingly, that hybridization, often associated with polyploidy, has been common throughout the evolutionary history of plants, particularly in various lineages of flowering plants. However, phylogenomic studies have also highlighted the challenges of disentangling signals of ancient hybridization from other sources of genomic conflict (in particular, incomplete lineage sorting). Here, we provide a critical review of ancient hybridization in vascular plants, outlining well-documented cases of ancient hybridization across plant phylogeny, as well as the challenges unique to documenting ancient versus recent hybridization. We provide a definition for ancient hybridization, which, to our knowledge, has not been explicitly attempted before. Further documenting the extent of deep reticulation in plants should remain an important research focus, especially because published examples likely represent the tip of the iceberg in terms of the total extent of ancient hybridization. However, future research should increasingly explore the macroevolutionary significance of this process, in terms of its impact on evolutionary trajectories (e.g. how does hybridization influence trait evolution or the generation of biodiversity over long time scales?), as well as how life history and ecological factors shape, or have shaped, the frequency of hybridization across geologic time and plant phylogeny. Finally, we consider the implications of ubiquitous ancient hybridization for how we conceptualize, analyze, and classify plant phylogeny. Networks, as opposed to bifurcating trees, represent more accurate representations of evolutionary history in many cases, although our ability to infer, visualize, and use networks for comparative analyses is highly limited. Developing improved methods for the generation, visualization, and use of networks represents a critical future direction for plant biology. Current classification systems also do not generally allow for the recognition of reticulate lineages, and our classifications themselves are largely based on evidence from the chloroplast genome. Updating plant classification to better reflect nuclear phylogenies, as well as considering whether and how to recognize hybridization in classification systems, will represent an important challenge for the plant systematics community.

Inference of Gene Flow between Species under Misspecified Models

Genomic sequence data provide a rich source of information about the history of species divergence and interspecific hybridization or introgression. Despite recent advances in genomics and statistical methods, it remains challenging to infer gene flow, and as a result, one may have to estimate introgression rates and times under misspecified models. Here we use mathematical analysis and computer simulation to examine estimation bias and issues of interpretation when the model of gene flow is misspecified in analysis of genomic datasets, for example, if introgression is assigned to the wrong lineages. In the case of two species, we establish a correspondence between the migration rate in the continuous migration model and the introgression probability in the introgression model. When gene flow occurs continuously through time but in the analysis is assumed to occur at a fixed time point, common evolutionary parameters such as species divergence times are surprisingly well estimated. However, the time of introgression tends to be estimated towards the recent end of the period of continuous gene flow. When introgression events are assigned incorrectly to the parental or daughter lineages, introgression times tend to collapse onto species divergence times, with introgression probabilities underestimated. Overall, our analyses suggest that the simple introgression model is useful for extracting information concerning between-specific gene flow and divergence even when the model may be misspecified. However, for reliable inference of gene flow it is important to include multiple samples per species, in particular, from hybridizing species.

Ghost lineages can invalidate or even reverse findings regarding gene flow

Introgression, endosymbiosis, and gene transfer, i.e., horizontal gene flow (HGF), are primordial sources of innovation in all domains of life. Our knowledge on HGF relies on detection methods that exploit some of its signatures left on extant genomes. One of them is the effect of HGF on branch lengths of constructed phylogenies. This signature has been formalized in statistical tests for HGF detection and used for example to detect massive adaptive gene flows in malaria vectors or to order evolutionary events involved in eukaryogenesis. However, these studies rely on the assumption that ghost lineages (all unsampled extant and extinct taxa) have little influence. We demonstrate here with simulations and data reanalysis that when considering the more realistic condition that unsampled taxa are legion compared to sampled ones, the conclusion of these studies become unfounded or even reversed. This illustrates the necessity to recognize the existence of ghosts in evolutionary studies.

Impact of Ghost Introgression on Coalescent-based Species Tree Inference and Estimation of Divergence Time

The species studied in any evolutionary investigation generally constitute a small proportion of all the species currently existing or that have gone extinct. It is therefore likely that introgression, which is widespread across the tree of life, involves "ghosts," i.e., unsampled, unknown, or extinct lineages. However, the impact of ghost introgression on estimations of species trees has rarely been studied and is poorly understood. Here, we use mathematical analysis and simulations to examine the robustness of species tree methods based on the multispecies coalescent model to introgression from a ghost or extant lineage. We found that many results originally obtained for introgression between extant species can easily be extended to ghost introgression, such as the strongly interactive effects of incomplete lineage sorting (ILS) and introgression on the occurrence of anomalous gene trees (AGTs). The relative performance of the summary species tree method (ASTRAL) and the full-likelihood method (*BEAST) varies under different introgression scenarios, with the former being more robust to gene flow between non-sister species whereas the latter performing better under certain conditions of ghost introgression. When an outgroup ghost (defined as a lineage that diverged before the most basal species under investigation) acts as the donor of the introgressed genes, the time of root divergence among the investigated species generally was overestimated, whereas ingroup introgression, as commonly perceived, can only lead to underestimation. In many cases of ingroup introgression that may or may not involve ghost lineages, the stronger the ILS, the higher the accuracy achieved in estimating the time of root divergence, although the topology of the species tree is more prone to be biased by the effect of introgression.

Population genetics reveals divergent lineages and ongoing hybridization in a declining migratory fish species complex

Deciphering the effects of historical and recent demographic processes responsible for the spatial patterns of genetic diversity and structure is a key objective in evolutionary and conservation biology. Using population genetic analyses, we investigated the demographic history, the contemporary genetic diversity and structure, and the occurrence of hybridization and introgression of two species of anadromous fish with contrasting life history strategies and which have undergone recent demographic declines, the allis shad (Alosa alosa) and the twaite shad (Alosa fallax). We genotyped 706 individuals from 20 rivers and 5 sites at sea in Southern Europe at thirteen microsatellite markers. Genetic structure between populations was lower for the nearly semelparous species A. alosa, which disperses greater distances compared to the iteroparous species, A. fallax. Individuals caught at sea were assigned at the river level for A. fallax and at the region level for A. alosa. Using an approximate Bayesian computation framework, we inferred that the most likely long term historical divergence scenario between both species and lineages involved historical separation followed by secondary contact accompanied by strong population size decline. Accordingly, we found evidence for contemporary hybridization and bidirectional introgression due to gene flow between both species and lineages. Moreover, our results support the existence of at least one distinct species in the Mediterrannean sea: A. agone in Golfe du Lion area, and another divergent lineage in Corsica. Overall, our results shed light on the interplay between historical and recent demographic processes and life history strategies in shaping population genetic diversity and structure of closely related species. The recent demographic decline of these species' populations and their hybridization should be carefully considered while implementing conservation programs.

Phylogenomic analyses in Phrymaceae reveal extensive gene tree discordance in relationships among major clades

PREMISE

Phylogenomic datasets using genomes and transcriptomes provide rich opportunities beyond resolving bifurcating phylogenetic relationships. Monkeyflower (Phrymaceae) is a model system for evolutionary ecology. However, it lacks a well-supported phylogeny as a basis for a stable taxonomy and for macroevolutionary comparisons.

METHODS

We sampled 24 genomes and transcriptomes in Phrymaceae and closely related families, including eight newly sequenced transcriptomes. We reconstructed the phylogeny using IQ-TREE and ASTRAL, evaluated gene tree discordance using PhyParts, Quartet Sampling, and a cloudogram, and carried out reticulation analyses using PhyloNet and HyDe. We searched for whole genome duplication (WGD) events using chromosome numbers, synonymous distances, and gene duplication events as evidence.

RESULTS

Most gene trees support the monophyly of Phrymaceae and each of its tribes. Most gene trees also support tribe Mimuleae being sister to Phrymeae + Diplaceae + Leucocarpeae, with extensive gene tree discordance among the latter three. Despite the discordance, the monophyly of Mimulus s.l. is rejected, and no individual reticulation event among the Phrymaceae tribes is well-supported. Reticulation likely occurred among Erythranthe bicolor and closely related species. No ancient WGD was detected in Phrymaceae. Instead, small-scale duplications are among potential drivers of macroevolutionary diversification of Phrymaceae.

CONCLUSIONS

We show that analysis of reticulate evolution is sensitive to taxon sampling and methods used. We also demonstrate that phylogenomic datasets using genomes and transcriptomes present rich opportunities to investigate gene family evolution and genome duplication events involved in lineage diversification and adaptation.