Versatile simulations of admixture and accurate local ancestry inference with mixnmatch and ancestryinfer

Pervasive gene flow despite strong and varied reproductive barriers in swordtails

The evolution of reproductive barriers leads to the formation of new species. However, recent research has demonstrated that hybridization has been pervasive across the tree of life even in the presence of strong barriers. Using swordtail fishes (genus Xiphophorus), an emerging model system, we document overlapping mechanisms that act as barriers to gene flow between Xiphophorus birchmanni and Xiphophorus cortezi by combining genomic sequencing from natural hybrid populations, experimental laboratory crosses, behavioural assays, sperm measures and developmental studies. We show that assortative mating plays a role in maintaining subpopulations with distinct ancestry within natural hybrid populations. Using F2 hybrids we identify several genomic regions that strongly impact hybrid viability. Strikingly, two of these regions underlie genetic incompatibilities in hybrids between X. birchmanni and its sister species Xiphophorus malinche. Our results demonstrate that ancient hybridization has played a role in the origin of this shared genetic incompatibility. Moreover, ancestry mismatch at these incompatible regions has remarkably similar consequences for phenotypes and hybrid survival in X. cortezi × X. birchmanni hybrids as in X. malinche × X. birchmanni hybrids. Our findings identify varied reproductive barriers that shape genetic exchange between naturally hybridizing species and highlight the complex evolutionary outcomes of hybridization.

Opportunities and challenges of local ancestry in genetic association analyses

Recently, admixed populations make up an increasing percentage of the US and global populations, and the admixture is not uniform over space or time or across genomes. Therefore, it becomes indispensable to evaluate local ancestry in addition to global ancestry to improve genetic epidemiological studies. Recent advances in representing human genome diversity, coupled with large-scale whole-genome sequencing initiatives and improved tools for local ancestry inference, have enabled studies to demonstrate that incorporating local ancestry information enhances both genetic association analyses and polygenic risk predictions. Along with the opportunities that local ancestry provides, there exist challenges preventing its full usage in genetic analyses. In this review, we first summarize methods for local ancestry inference and illustrate how local ancestry can be utilized in various analyses, including admixture mapping, association testing, and polygenic risk score construction. In addition, we discuss current challenges in research involving local ancestry, both in terms of the inference itself and its role in genetic association studies. We further pinpoint some future study directions and methodology development opportunities to help more effectively incorporate local ancestry in genetic analyses. It is worth the effort to pursue those future directions and address these analytical challenges because the appropriate use of local ancestry estimates could help mitigate inequality in genomic medicine and improve our understanding of health and disease outcomes.

Repeatable Selection on Large Ancestry Blocks in an Avian Hybrid Zone

Hybrid zones create natural tests of genetic incompatibilities by combining loci from 2 species in the same genetic background in the wild, making them useful for identifying loci involved in both intrinsic and ecological (extrinsic) isolation. Two Swainson's thrush subspecies form a hybrid zone in western North America. These coastal and inland subspecies exhibit dramatic differences in migration routes; their hybrids exhibit poor migratory survival, suggesting that ecological incompatibilities maintain this zone. We used a panel of ancestry informative markers to identify repeated patterns of selection and introgression across 4 hybrid populations that span the entire length of the Swainson's thrush hybrid zone. Two repeatable patterns consistent with selection against incompatibilities-steep genomic clines and few transitions between ancestry states-were found in large genetic blocks on chromosomes 1 and 5. The block on chromosome 1 showed evidence for inland subspecies introgression while the block on chromosome 5 exhibited coastal subspecies introgression. Some regions previously associated with migratory phenotypes, including migratory orientation, or exhibiting misexpression between the subspecies exhibited signatures of selection in the hybrid zone. Both selection and introgression across the genome were shaped by genomic structural features and evolutionary history, with stronger selection and reduced introgression in regions of low recombination, high subspecies differentiation, positive selection within the subspecies, and on macrochromosomes. Cumulatively, these results suggest that linkage among loci interacts with divergent selection and past divergent evolution between species to strengthen barriers to gene flow within hybrid zones.

Fluctuating reproductive isolation and stable ancestry structure in a fine-scaled mosaic of hybridizing Mimulus monkeyflowers

Hybridization among taxa impacts a variety of evolutionary processes from adaptation to extinction. We seek to understand both patterns of hybridization across taxa and the evolutionary and ecological forces driving those patterns. To this end, we use whole-genome low-coverage sequencing of 458 wild-grown and 1565 offspring individuals to characterize the structure, stability, and mating dynamics of admixed populations of Mimulus guttatus and Mimulus nasutus across a decade of sampling. In three streams, admixed genomes are common and a M. nasutus organellar haplotype is fixed in M. guttatus, but new hybridization events are rare. Admixture is strongly unidirectional, but each stream has a unique distribution of ancestry proportions. In one stream, three distinct cohorts of admixed ancestry are spatially structured at ~20-50m resolution and stable across years. Mating system provides almost complete isolation of M. nasutus from both M. guttatus and admixed cohorts, and is a partial barrier between admixed and M. guttatus cohorts. Isolation due to phenology is near-complete between M. guttatus and M. nasutus. Phenological isolation is a strong barrier in some years between admixed and M. guttatus cohorts, but a much weaker barrier in other years, providing a potential bridge for gene flow. These fluctuations are associated with differences in water availability across years, supporting a role for climate in mediating the strength of reproductive isolation. Together, mating system and phenology accurately predict fluctuations in assortative mating across years, which we estimate directly using paired maternal and offspring genotypes. Climate-driven fluctuations in reproductive isolation may promote the longer-term stability of a complex mosaic of hybrid ancestry, preventing either complete isolation or complete collapse of species barriers.

Structural genomic variation and behavioral interactions underpin a balanced sexual mimicry polymorphism

How phenotypic diversity originates and persists within populations are classic puzzles in evolutionary biology. While balanced polymorphisms segregate within many species, it remains rare for both the genetic basis and the selective forces to be known, leading to an incomplete understanding of many classes of traits under balancing selection. Here, we uncover the genetic architecture of a balanced sexual mimicry polymorphism and identify behavioral mechanisms that may be involved in its maintenance in the swordtail fish Xiphophorus birchmanni. We find that ∼40% of X. birchmanni males develop a "false gravid spot," a melanic pigmentation pattern that mimics the "pregnancy spot" associated with sexual maturity in female live-bearing fish. Using genome-wide association mapping, we detect a single intergenic region associated with variation in the false gravid spot phenotype, which is upstream of kitlga, a melanophore patterning gene. By performing long-read sequencing within and across populations, we identify complex structural rearrangements between alternate alleles at this locus. The false gravid spot haplotype drives increased allele-specific expression of kitlga, which provides a mechanistic explanation for the increased melanophore abundance that causes the spot. By studying social interactions in the laboratory and in nature, we find that males with the false gravid spot experience less aggression; however, they also receive increased attention from other males and are disdained by females. These behavioral interactions may contribute to the maintenance of this phenotypic polymorphism in natural populations. We speculate that structural variants affecting gene regulation may be an underappreciated driver of balanced polymorphisms across diverse species.

Swordtail fish hybrids reveal that genome evolution is surprisingly predictable after initial hybridization

Over the past 2 decades, biologists have come to appreciate that hybridization, or genetic exchange between distinct lineages, is remarkably common-not just in particular lineages but in taxonomic groups across the tree of life. As a result, the genomes of many modern species harbor regions inherited from related species. This observation has raised fundamental questions about the degree to which the genomic outcomes of hybridization are repeatable and the degree to which natural selection drives such repeatability. However, a lack of appropriate systems to answer these questions has limited empirical progress in this area. Here, we leverage independently formed hybrid populations between the swordtail fish Xiphophorus birchmanni and X. cortezi to address this fundamental question. We find that local ancestry in one hybrid population is remarkably predictive of local ancestry in another, demographically independent hybrid population. Applying newly developed methods, we can attribute much of this repeatability to strong selection in the earliest generations after initial hybridization. We complement these analyses with time-series data that demonstrates that ancestry at regions under selection has remained stable over the past approximately 40 generations of evolution. Finally, we compare our results to the well-studied X. birchmanni × X. malinche hybrid populations and conclude that deeper evolutionary divergence has resulted in stronger selection and higher repeatability in patterns of local ancestry in hybrids between X. birchmanni and X. cortezi.

Pervasive gene flow despite strong and varied reproductive barriers in swordtails

One of the mechanisms that can lead to the formation of new species occurs through the evolution of reproductive barriers. However, recent research has demonstrated that hybridization has been pervasive across the tree of life even in the presence of strong barriers. Swordtail fishes (genus Xiphophorus) are an emerging model system for studying the interface between these barriers and hybridization. We document overlapping mechanisms that act as barriers between closely related species, X. birchmanni and X. cortezi, by combining genomic sequencing from natural hybrid populations, artificial crosses, behavioral assays, sperm performance, and developmental studies. We show that strong assortative mating plays a key role in maintaining subpopulations with distinct ancestry in natural hybrid populations. Lab experiments demonstrate that artificial F1 crosses experience dysfunction: crosses with X. birchmanni females were largely inviable and crosses with X. cortezi females had a heavily skewed sex ratio. Using F2 hybrids we identify several genomic regions that strongly impact hybrid viability. Strikingly, two of these regions underlie genetic incompatibilities in hybrids between X. birchmanni and its sister species X. malinche. Our results demonstrate that ancient hybridization has played a role in the origin of this shared genetic incompatibility. Moreover, ancestry mismatch at these incompatible regions has remarkably similar consequences for phenotypes and hybrid survival in X. cortezi × X. birchmanni hybrids as in X. malinche × X. birchmanni hybrids. Our findings identify varied reproductive barriers that shape genetic exchange between naturally hybridizing species and highlight the complex evolutionary outcomes of hybridization.

A lethal mitonuclear incompatibility in complex I of natural hybrids

The evolution of reproductive barriers is the first step in the formation of new species and can help us understand the diversification of life on Earth. These reproductive barriers often take the form of hybrid incompatibilities, in which alleles derived from two different species no longer interact properly in hybrids1-3. Theory predicts that hybrid incompatibilities may be more likely to arise at rapidly evolving genes4-6 and that incompatibilities involving multiple genes should be common7,8, but there has been sparse empirical data to evaluate these predictions. Here we describe a mitonuclear incompatibility involving three genes whose protein products are in physical contact within respiratory complex I of naturally hybridizing swordtail fish species. Individuals homozygous for mismatched protein combinations do not complete embryonic development or die as juveniles, whereas those heterozygous for the incompatibility have reduced complex I function and unbalanced representation of parental alleles in the mitochondrial proteome. We find that the effects of different genetic interactions on survival are non-additive, highlighting subtle complexity in the genetic architecture of hybrid incompatibilities. Finally, we document the evolutionary history of the genes involved, showing signals of accelerated evolution and evidence that an incompatibility has been transferred between species via hybridization.

Experimental evolution of hybrid populations to identify Dobzhansky-Muller incompatibility loci

Epistatic interactions between loci that reduce fitness in interspecies hybrids, Dobzhansky-Muller incompatibilities (DMIs), contribute genetically to the inviability and infertility within hybrid populations. It remains a challenge, however, to identify the loci that contribute to DMIs as causes of reproductive isolation between species. Here, we assess through forward simulation the power of evolve-and-resequence (E&R) experimental evolution of hybrid populations to map DMI loci. We document conditions under which such a mapping strategy may be most feasible and demonstrate how mapping power is sensitive to biologically relevant parameters such as one-way versus two-way incompatibility type, selection strength, recombination rate, and dominance interactions. We also assess the influence of parameters under direct control of an experimenter, including duration of experimental evolution and number of replicate populations. We conclude that an E&R strategy for mapping DMI loci, and other cases of epistasis, can be a viable option under some circumstances for study systems with short generation times like Caenorhabditis nematodes.

Genome evolution is surprisingly predictable after initial hybridization

Over the past two decades, evolutionary biologists have come to appreciate that hybridization, or genetic exchange between distinct lineages, is remarkably common - not just in particular lineages but in taxonomic groups across the tree of life. As a result, the genomes of many modern species harbor regions inherited from related species. This observation has raised fundamental questions about the degree to which the genomic outcomes of hybridization are repeatable and the degree to which natural selection drives such repeatability. However, a lack of appropriate systems to answer these questions has limited empirical progress in this area. Here, we leverage independently formed hybrid populations between the swordtail fish Xiphophorus birchmanni and X. cortezi to address this fundamental question. We find that local ancestry in one hybrid population is remarkably predictive of local ancestry in another, demographically independent hybrid population. Applying newly developed methods, we can attribute much of this repeatability to strong selection in the earliest generations after initial hybridization. We complement these analyses with time-series data that demonstrates that ancestry at regions under selection has remained stable over the past ~40 generations of evolution. Finally, we compare our results to the well-studied X. birchmanni×X. malinche hybrid populations and conclude that deeper evolutionary divergence has resulted in stronger selection and higher repeatability in patterns of local ancestry in hybrids between X. birchmanni and X. cortezi.

link-ancestors: fast simulation of local ancestry with tree sequence software

Summary

It is challenging to simulate realistic tracts of genetic ancestry on a scale suitable for simulation-based inference. We present an algorithm that enables this information to be extracted efficiently from tree sequences produced by simulations run with msprime and SLiM.

Availability and implementation

A C-based implementation of the link-ancestors algorithm is in tskit (https://tskit.dev/tskit/docs/stable/). We also provide a user-friendly wrapper for link-ancestors in tspop, a Python-based utility package.

Localizing Post-Admixture Adaptive Variants with Object Detection on Ancestry-Painted Chromosomes

Gene flow between previously differentiated populations during the founding of an admixed or hybrid population has the potential to introduce adaptive alleles into the new population. If the adaptive allele is common in one source population, but not the other, then as the adaptive allele rises in frequency in the admixed population, genetic ancestry from the source containing the adaptive allele will increase nearby as well. Patterns of genetic ancestry have therefore been used to identify post-admixture positive selection in humans and other animals, including examples in immunity, metabolism, and animal coloration. A common method identifies regions of the genome that have local ancestry "outliers" compared with the distribution across the rest of the genome, considering each locus independently. However, we lack theoretical models for expected distributions of ancestry under various demographic scenarios, resulting in potential false positives and false negatives. Further, ancestry patterns between distant sites are often not independent. As a result, current methods tend to infer wide genomic regions containing many genes as under selection, limiting biological interpretation. Instead, we develop a deep learning object detection method applied to images generated from local ancestry-painted genomes. This approach preserves information from the surrounding genomic context and avoids potential pitfalls of user-defined summary statistics. We find the method is robust to a variety of demographic misspecifications using simulated data. Applied to human genotype data from Cabo Verde, we localize a known adaptive locus to a single narrow region compared with multiple or long windows obtained using two other ancestry-based methods.

Imbalanced segregation of recombinant haplotypes in hybrid populations reveals inter- and intrachromosomal Dobzhansky-Muller incompatibilities

Dobzhansky-Muller incompatibilities (DMIs) are a major component of reproductive isolation between species. DMIs imply negative epistasis and are exposed when two diverged populations hybridize. Mapping the locations of DMIs has largely relied on classical genetic mapping. Approaches to date are hampered by low power and the challenge of identifying DMI loci on the same chromosome, because strong initial linkage of parental haplotypes weakens statistical tests. Here, we propose new statistics to infer negative epistasis from haplotype frequencies in hybrid populations. When two divergent populations hybridize, the variance in heterozygosity at two loci decreases faster with time at DMI loci than at random pairs of loci. When two populations hybridize at near-even admixture proportions, the deviation of the observed variance from its expectation becomes negative for the DMI pair. This negative deviation enables us to detect intermediate to strong negative epistasis both within and between chromosomes. In practice, the detection window in hybrid populations depends on the demographic scenario, the recombination rate, and the strength of epistasis. When the initial proportion of the two parental populations is uneven, only strong DMIs can be detected with our method unless migration prevents parental haplotypes from being lost. We use the new statistics to infer candidate DMIs from three hybrid populations of swordtail fish. We identify numerous new DMI candidates, some of which are inferred to interact with several loci within and between chromosomes. Moreover, we discuss our results in the context of an expected enrichment in intrachromosomal over interchromosomal DMIs.

Predictability and parallelism in the contemporary evolution of hybrid genomes

Hybridization between species is widespread across the tree of life. As a result, many species, including our own, harbor regions of their genome derived from hybridization. Despite the recognition that this process is widespread, we understand little about how the genome stabilizes following hybridization, and whether the mechanisms driving this stabilization tend to be shared across species. Here, we dissect the drivers of variation in local ancestry across the genome in replicated hybridization events between two species pairs of swordtail fish: Xiphophorus birchmanni × X. cortezi and X. birchmanni × X. malinche. We find unexpectedly high levels of repeatability in local ancestry across the two types of hybrid populations. This repeatability is attributable in part to the fact that the recombination landscape and locations of functionally important elements play a major role in driving variation in local ancestry in both types of hybrid populations. Beyond these broad scale patterns, we identify dozens of regions of the genome where minor parent ancestry is unusually low or high across species pairs. Analysis of these regions points to shared sites under selection across species pairs, and in some cases, shared mechanisms of selection. We show that one such region is a previously unknown hybrid incompatibility that is shared across X. birchmanni × X. cortezi and X. birchmanni × X. malinche hybrid populations.

A dense linkage map for a large repetitive genome: discovery of the sex-determining region in hybridizing fire-bellied toads (Bombina bombina and Bombina variegata)

Genomic analysis of hybrid zones offers unique insights into emerging reproductive isolation and the dynamics of introgression. Because hybrid genomes consist of blocks inherited from one or the other parental taxon, linkage information is essential. In most cases, the spectrum of local ancestry tracts can be efficiently uncovered from dense linkage maps. Here, we report the development of such a map for the hybridizing toads, Bombina bombina and Bombina variegata (Anura: Bombinatoridae). Faced with the challenge of a large (7–10 Gb), repetitive genome, we set out to identify a large number of Mendelian markers in the nonrepetitive portion of the genome that report B. bombina vs B. variegata ancestry with appropriately quantified statistical support. Bait sequences for targeted enrichment were selected from a draft genome assembly, after filtering highly repetitive sequences. We developed a novel approach to infer the most likely diplotype per sample and locus from the raw read mapping data, which is robust to over-merging and obviates arbitrary filtering thresholds. Validation of the resulting map with 4755 markers underscored the large-scale synteny between Bombina and Xenopus tropicalis. By assessing the sex of late-stage F2 tadpoles from histological sections, we identified the sex-determining region in the Bombina genome to 7 cM on LG5, which is homologous to X. tropicalis chromosome 5, and inferred male heterogamety. Interestingly, chromosome 5 has been repeatedly recruited as a sex chromosome in anurans with XY sex determination.

Two new hybrid populations expand the swordtail hybridization model system

Natural hybridization events provide unique windows into the barriers that keep species apart as well as the consequences of their breakdown. Here, we characterize hybrid populations formed between the northern swordtail fish Xiphophorus cortezi and Xiphophorus birchmanni from collection sites on two rivers. We use simulations and new genetic reference panels to develop sensitive and accurate local ancestry calling in this novel system. Strikingly, we find that hybrid populations on both rivers consist of two genetically distinct subpopulations: a cluster of pure X. birchmanni individuals and one of phenotypically intermediate hybrids that derive ∼85-90% of their genome from X. cortezi. Simulations suggest that initial hybridization occurred ∼150 generations ago at both sites, with little evidence for contemporary gene flow between subpopulations. This population structure is consistent with strong assortative mating between individuals of similar ancestry. The patterns of population structure uncovered here mirror those seen in hybridization between X. birchmanni and its sister species, Xiphophorus malinche, indicating an important role for assortative mating in the evolution of hybrid populations. Future comparisons will provide a window into the shared mechanisms driving the outcomes of hybridization not only among independent hybridization events between the same species but also across distinct species pairs.

The genomic consequences of hybridization

In the past decade, advances in genome sequencing have allowed researchers to uncover the history of hybridization in diverse groups of species, including our own. Although the field has made impressive progress in documenting the extent of natural hybridization, both historical and recent, there are still many unanswered questions about its genetic and evolutionary consequences. Recent work has suggested that the outcomes of hybridization in the genome may be in part predictable, but many open questions about the nature of selection on hybrids and the biological variables that shape such selection have hampered progress in this area. We synthesize what is known about the mechanisms that drive changes in ancestry in the genome after hybridization, highlight major unresolved questions, and discuss their implications for the predictability of genome evolution after hybridization.

A beginner's guide to low-coverage whole genome sequencing for population genomics

Low-coverage whole genome sequencing (lcWGS) has emerged as a powerful and cost-effective approach for population genomic studies in both model and non-model species. However, with read depths too low to confidently call individual genotypes, lcWGS requires specialized analysis tools that explicitly account for genotype uncertainty. A growing number of such tools have become available, but it can be difficult to get an overview of what types of analyses can be performed reliably with lcWGS data, and how the distribution of sequencing effort between the number of samples analyzed and per-sample sequencing depths affects inference accuracy. In this introductory guide to lcWGS, we first illustrate how the per-sample cost for lcWGS is now comparable to RAD-seq and Pool-seq in many systems. We then provide an overview of software packages that explicitly account for genotype uncertainty in different types of population genomic inference. Next, we use both simulated and empirical data to assess the accuracy of allele frequency and genetic diversity estimation, detection of population structure, and selection scans under different sequencing strategies. Our results show that spreading a given amount of sequencing effort across more samples with lower depth per sample consistently improves the accuracy of most types of inference, with a few notable exceptions. Finally, we assess the potential for using imputation to bolster inference from lcWGS data in non-model species, and discuss current limitations and future perspectives for lcWGS-based population genomics research. With this overview, we hope to make lcWGS more approachable and stimulate its broader adoption.

Systematic Review on Local Ancestor Inference From a Mathematical and Algorithmic Perspective

Genotypic data provide deep insights into the population history and medical genetics. The local ancestry inference (LAI) (also termed local ancestry deconvolution) method uses the hidden Markov model (HMM) to solve the mathematical problem of ancestry reconstruction based on genomic data. HMM is combined with other statistical models and machine learning techniques for particular genetic tasks in a series of computer tools. In this article, we surveyed the mathematical structure, application characteristics, historical development, and benchmark analysis of the LAI method in detail, which will help researchers better understand and further develop LAI methods. Firstly, we extensively explore the mathematical structure of each model and its characteristic applications. Next, we use bibliometrics to show detailed model application fields and list articles to elaborate on the historical development. LAI publications had experienced a peak period during 2006-2016 and had kept on moving in the following years. The efficiency, accuracy, and stability of the existing models were evaluated by the benchmark. We find that phased data had higher accuracy in comparison with unphased data. We summarize these models with their distinct advantages and disadvantages. The Loter model uses dynamic programming to obtain a globally optimal solution with its parameter-free advantage. Aligned bases can be used directly in the Seqmix model if the genotype is hard to call. This research may help model developers to realize current challenges, develop more advanced models, and enable scholars to select appropriate models according to given populations and datasets.

Inferring Adaptive Introgression Using Hidden Markov Models

Adaptive introgression-the flow of adaptive genetic variation between species or populations-has attracted significant interest in recent years and it has been implicated in a number of cases of adaptation, from pesticide resistance and immunity, to local adaptation. Despite this, methods for identification of adaptive introgression from population genomic data are lacking. Here, we present Ancestry_HMM-S, a hidden Markov model-based method for identifying genes undergoing adaptive introgression and quantifying the strength of selection acting on them. Through extensive validation, we show that this method performs well on moderately sized data sets for realistic population and selection parameters. We apply Ancestry_HMM-S to a data set of an admixed Drosophila melanogaster population from South Africa and we identify 17 loci which show signatures of adaptive introgression, four of which have previously been shown to confer resistance to insecticides. Ancestry_HMM-S provides a powerful method for inferring adaptive introgression in data sets that are typically collected when studying admixed populations. This method will enable powerful insights into the genetic consequences of admixture across diverse populations. Ancestry_HMM-S can be downloaded from https://github.com/jesvedberg/Ancestry_HMM-S/.

The Genetic Architecture of Variation in the Sexually Selected Sword Ornament and Its Evolution in Hybrid Populations

Biologists since Darwin have been fascinated by the evolution of sexually selected ornaments, particularly those that reduce viability. Uncovering the genetic architecture of these traits is key to understanding how they evolve and are maintained. Here, we investigate the genetic architecture and evolutionary loss of a sexually selected ornament, the "sword" fin extension that characterizes many species of swordtail fish (Xiphophorus). Using sworded and swordless sister species of Xiphophorus, we generated a mapping population and show that the sword ornament is polygenic-with ancestry across the genome explaining substantial variation in the trait. After accounting for the impacts of genome-wide ancestry, we identify one major-effect quantitative trait locus (QTL) that explains ~5% of the overall variation in the trait. Using a series of approaches, we narrow this large QTL interval to several likely candidate genes, including genes involved in fin regeneration and growth. Furthermore, we find evidence of selection on ancestry at one of these candidates in four natural hybrid populations, consistent with selection against the sword in these populations.