Search

Animal evolution is characterized by frequent turnover of sexually dimorphic traits – new sex- specific characters are gained, and some ancestral sex-specific characters are lost, in many lineages. In insects, sexual differentiation is predominantly cell-autonomous and depends on the expression of the doublesex (dsx) transcription factor. In most cases, cells that transcribe dsx have the potential to undergo sex-specific differentiation, while those that lack dsx expression do not. Consistent with this mode of development, comparative research has shown that the origin of new sex-specific traits can be associated with the origin of new spatial domains of dsx expression. In this report, we examine the opposite situation – a secondary loss of the sex comb, a male-specific grasping structure that develops on the front legs of some drosophilid species. We show that, while the origin of the sex comb is linked to an evolutionary gain of dsx expression in the leg, sex comb loss in a newly identified species of Lordiphosa (Drosophilidae) is associated with a secondary loss of dsx expression. We discuss how the developmental control of sexual dimorphism affects the mechanisms by which sex-specific traits can evolve.

There is increasing interest in the African spiny mouse (Acomys cahirinus) as a model organism because of its ability for regeneration of tissue after injury in skin, muscle, and internal organs such as the kidneys. A high-quality reference genome is needed to better understand these regenerative properties at the molecular level. Here, we present an improved reference genome forA. cahirinusgenerated from long Nanopore sequencing reads. We confirm the quality of our annotations using RNA sequencing data from four different tissues. Our genome is of higher contiguity and quality than previously reported genomes from this species and will facilitate ongoing efforts to better understand the regenerative properties of this organism.

Sequence variation among antigenicvargenes enablesPlasmodium falciparummalaria parasites to evade host immunity. Using long sequence reads from haploid clones from a mutation accumulation experiment, we detectvardiversity inconsistent with simple chromosomal inheritance. We discover putatively circular DNA that is strongly enriched forvargenes, which exist in multiple alleles per locus separated by recombination and indel events. Extrachromosomal DNA likely contributes to rapid antigenic diversification inP. falciparum.

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.

Effective suppression of transposable elements (TEs) is paramount to maintain genomic integrity and organismal fitness. InD. melanogaster,flamencois a master suppressor of TEs, preventing their movement from somatic ovarian support cells to the germline. It is transcribed by Pol II as a long (100s of kb), single-stranded, primary transcript, that is metabolized into Piwi-interacting RNAs (piRNAs) that target active TEs via antisense complementarity.flamencois thought to operate as a trap, owing to its high content of recent horizontally transferred TEs that are enriched in antisense orientation. Using newly-generated long read genome data, which is critical for accurate assembly of repetitive sequences, we find thatflamencohas undergone radical transformations in sequence content and even copy number acrosssimulansclade Drosophilid species.D. simulans flamencohas duplicated and diverged, and neither copy exhibits synteny withD. melanogasterbeyond the core promoter. Moreover,flamencoorganization is highly variable acrossD. simulansindividuals. Next, we find thatD. simulansandD. mauritiana flamencodisplay signatures of a dual-stranded cluster, with ping-pong signals in the testis and/or embryo. This is accompanied by increased copy numbers of germline TEs, consistent with these regions operating as functional dual stranded clusters. Overall, the physical and functional diversity offlamencoorthologs is testament to the extremely dynamic consequences of TE arms races on genome organization, not only amongst highly related species, but even amongst individuals.

We have been assembling the genomes of many Drosophila species. This protocol fork allows for the sequencing of single wild-caught flies using Nanopore MinION sequencers, in contrast to our previous protocol that originally used gDNA extracted from tens or hundreds of flies from an inbred line. While read lengths are not nearly close to what we were seeing before, libraries range from 2-10kb read N50 and typically produce genome assemblies with contig N50 >1 Mbp. The protocol for ultra-long read sequencing from large pools of flies can be found here: dx.doi.org/10.17504/protocols.io.bdfqi3mw For now, please cite the original genomes paper if you use this protocol: https://doi.org/10.7554/eLife.66405 An update that includes many single wild-caught fly genomes is in progress.

Balancers are rearranged chromosomes used in Drosophila melanogaster to maintain deleterious mutations in stable populations, preserve sets of linked genetic elements and construct complex experimental stocks. Here, we assess the phenotypes associated with breakpoint-induced mutations on commonly used third chromosome balancers and show remarkably few deleterious effects. We demonstrate that a breakpoint in p53 causes loss of radiation-induced apoptosis and a breakpoint in Fucosyltransferase A causes loss of fucosylation in nervous and intestinal tissue—the latter study providing new markers for intestinal cell identity and challenging previous conclusions about the regulation of fucosylation. We also describe thousands of potentially harmful mutations shared among X or third chromosome balancers, or unique to specific balancers, including an Ankyrin 2 mutation present on most TM3 balancers, and reiterate the risks of using balancers as experimental controls. We used long-read sequencing to confirm or refine the positions of two inversions with breakpoints lying in repetitive sequences and provide evidence that one of the inversions, In(2L)Cy, arose by ectopic recombination between foldback transposon insertions and the other, In(3R)C, cleanly separates subtelomeric and telomeric sequences and moves the subtelomeric sequences to an internal chromosome position. In addition, our characterization of In(3R)C shows that balancers may be polymorphic for terminal deletions. Finally, we present evidence that extremely distal mutations on balancers can add to the stability of stocks whose purpose is to maintain homologous chromosomes carrying mutations in distal genes. Overall, these studies add to our understanding of the structure, diversity and effectiveness of balancer chromosomes.

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.

Over 100 years of studies in Drosophila melanogaster and related species in the genus Drosophila have facilitated key discoveries in genetics, genomics, and evolution. While high-quality genome assemblies exist for several species in this group, they only encompass a small fraction of the genus. Recent advances in long-read sequencing allow high-quality genome assemblies for tens or even hundreds of species to be efficiently generated. Here, we utilize Oxford Nanopore sequencing to build an open community resource of genome assemblies for 101 lines of 93 drosophilid species encompassing 14 species groups and 35 sub-groups. The genomes are highly contiguous and complete, with an average contig N50 of 10.5 Mb and greater than 97% BUSCO completeness in 97/101 assemblies. We show that Nanopore-based assemblies are highly accurate in coding regions, particularly with respect to coding insertions and deletions. These assemblies, along with a detailed laboratory protocol and assembly pipelines, are released as a public resource and will serve as a starting point for addressing broad questions of genetics, ecology, and evolution at the scale of hundreds of species.

To advance our understanding of adaptation to temporally varying selection pressures, we identified signatures of seasonal adaptation occurring in parallel among Drosophila melanogaster populations. Specifically, we estimated allele frequencies genome-wide from flies sampled early and late in the growing season from 20 widely dispersed populations. We identified parallel seasonal allele frequency shifts across North America and Europe, demonstrating that seasonal adaptation is a general phenomenon of temperate fly populations. Seasonally fluctuating polymorphisms are enriched in large chromosomal inversions, and we find a broad concordance between seasonal and spatial allele frequency change. The direction of allele frequency change at seasonally variable polymorphisms can be predicted by weather conditions in the weeks prior to sampling, linking the environment and the genomic response to selection. Our results suggest that fluctuating selection is an important evolutionary force affecting patterns of genetic variation in Drosophila.

Genome-scale sequence data have invigorated the study of hybridization and introgression, particularly in animals. However, outside of a few notable cases, we lack systematic tests for introgression at a larger phylogenetic scale across entire clades. Here we leverage 155 genome assemblies, from 149 species, to generate a fossil-calibrated phylogeny and conduct multilocus tests for introgression across nine monophyletic radiations within the genus Drosophila. Using complementary phylogenomic approaches, we identify widespread introgression across the evolutionary history of Drosophila. Mapping gene-tree discordance onto the phylogeny revealed that both ancient and recent introgression has occurred across most of the nine clades that we examined. Our results provide the first evidence of introgression occurring across the evolutionary history of Drosophila and highlight the need to continue to study the evolutionary consequences of hybridization and introgression in this genus and across the Tree of Life.

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 fishXiphophorus corteziandX. birchmannifrom collection sites on two rivers. We develop sensitive and accurate local ancestry calling for this system based on low coverage whole genome sequencing. Strikingly, we find that hybrid populations on both rivers consist of two genetically distinct subpopulations: a cluster of nearly pureX. birchmanniindividuals and one of phenotypically intermediate hybrids that derive ~85-90% of their genome fromX. cortezi. Simulations and empirical data suggest that at both sites initial hybridization occurred ~150 generations ago, with little evidence for contemporary gene flow between subpopulations, likely due to strong assortative mating. The patterns of population structure uncovered here mirror those seen in hybridization betweenX. birchmanniand its sister species,X. malinche. Future comparisons will provide a window into the repeatability of the outcomes of hybridization not only across independent hybridization events between the same species but also across distinct species pairs.

Comparative genomic approaches have been used to identify sites where mutations are under purifying selection and of functional consequence by searching for sequences that are conserved across distantly related species. However, the performance of these approaches has not been rigorously evaluated under population genetic models. Further, short-lived functional elements may not leave a footprint of sequence conservation across many species. We use simulations to study how one measure of conservation, the Genomic Evolutionary Rate Profiling (GERP) score, relates to the strength of selection (Nes). We show that the GERP score is related to the strength of purifying selection. However, changes in selection coefficients or functional elements over time (i.e. functional turnover) can strongly affect the GERP distribution, leading to unexpected relationships between GERP and Nes. Further, we show that for functional elements that have a high turnover rate, adding more species to the analysis does not necessarily increase statistical power. Finally, we use the distribution of GERP scores across the human genome to compare models with and without turnover of sites where mutations are under purifying selection. We show that mutations in 4.51% of the noncoding human genome are under purifying selection and that most of this sequence has likely experienced changes in selection coefficients throughout mammalian evolution. Our work reveals limitations to using comparative genomic approaches to identify deleterious mutations. Commonly used GERP score thresholds miss over half of the noncoding sites in the human genome where mutations are under purifying selection., Author summary One of the most significant and challenging tasks in modern genomics is to assess the functional consequences of a particular nucleotide change in a genome. A common approach to address this challenge prioritizes sequences that share similar nucleotides across distantly related species, with the rationale that mutations at such positions were deleterious and removed from the population by purifying natural selection. Our manuscript shows that one popular measure of sequence conservation, the GERP score, performs well at identifying selected mutations if mutations at a site were under selection across all of mammalian evolution. Changes in selection at a given site dramatically reduces the power of GERP to detect selected mutations in humans. We also combine population genetic models with the distribution of GERP scores at noncoding sites across the human genome to show that the degree of selection at individual sites has changed throughout mammalian evolution. Importantly, we demonstrate that at least 80 Mb of noncoding sequence under purifying selection in humans will not have extreme GERP scores and will likely be missed by modern comparative genomic approaches. Our work argues that new approaches, potentially based on genetic variation within species, will be required to identify deleterious mutations.

We have been assembling the genomes of many Drosophila species. With that in mind, we developed this protocol to keep the cost of sequencing down to 100kbp. Sequencing is halted at about 40-50X depth of coverage (8-10 Gbp for most species). This of course depends on the quality of the sample, quality and quantity of the prepared library, and the frequency at which the flow cell is flushed and reloaded. We typically run 3-4 species per 2 flow cells, usually for ~14-18 Gbp of data per flow cell. This protocol borrows several elements from John Tyson's "Rocky Mountain" protocol and we thank him for several insightful discussions. https://www.protocols.io/view/rocky-mountain-adventures-in-genomic-dna-sample-pr-7euhjew

Admixture with archaic hominins has altered the landscape of genomic variation in modern human populations. Several gene regions have been previously identified as candidates of adaptive introgression (AI) that facilitated human adaptation to specific environments. However, simulation-based studies have suggested that population genetics processes other than adaptive mutations, such as heterosis from recessive deleterious variants private to populations before admixture, can also lead to patterns in genomic data that resemble adaptive introgression. The extent to which the presence of deleterious variants affect the false-positive rate and the power of current methods to detect AI has not been fully assessed. Here, we used extensive simulations to show that recessive deleterious mutations can increase the false positive rates of tests for AI compared to models without deleterious variants. We further examined candidates of AI in modern humans identified from previous studies and show that, although deleterious variants may hinder the performance of AI detection in modern humans, most signals remained robust when deleterious variants are included in the null model. While deleterious variants may have a limited impact on detecting signals of adaptive introgression in humans, we found that at least two AI candidate genes,HYAL2andHLA, are particularly susceptible to high false positive rates due to the recessive deleterious mutations. By quantifying parameters that affect heterosis, we show that the high false positives are largely attributed to the high exon densities together with low recombination rates in the genomic regions, which can further be exaggerated by the population growth in recent human evolution. Although the combination of such parameters is rare in the human genome, caution is still warranted in other species with different genomic composition and demographic histories.

The explosion in population genomic data demands ever more complex modes of analysis, and increasingly these analyses depend on sophisticated simulations. Recent advances in population genetic simulation have made it possible to simulate large and complex models, but specifying such models for a particular simulation engine remains a difficult and error-prone task. Computational genetics researchers currently re-implement simulation models independently, leading to inconsistency and duplication of effort. This situation presents a major barrier to empirical researchers seeking to use simulations for power analyses of upcoming studies or sanity checks on existing genomic data. Population genetics, as a field, also lacks standard benchmarks by which new tools for inference might be measured. Here we describe a new resource,stdpopsim, that attempts to rectify this situation.Stdpopsimis a community-driven open source project, which provides easy access to a growing catalog of published simulation models from a range of organisms and supports multiple simulation engine backends. This resource is available as a well-documented python library with a simple command-line interface. We share some examples demonstrating howstdpopsimcan be used to systematically compare demographic inference methods, and we encourage a broader community of developers to contribute to this growing resource.

The origin, diversification, and secondary loss of sexually dimorphic characters are common in animal evolution. In some cases, structurally and functionally similar traits have evolved independently in multiple lineages. Prominent examples of such traits include the male-specific grasping structures that develop on the front legs of many dipteran insects. In this report, we describe the evolution and development of one of these structures, the male-specific “sex brush”. The sex brush is composed of densely packed, irregularly arranged modified bristles and is found in several distantly related lineages in the family Drosophilidae. Phylogenetic analysis using 250 genes from over 200 species provides modest support for a single origin of the sex brush followed by many secondary losses; however, independent origins of the sex brush cannot be ruled out completely. We show that sex brushes develop in very similar ways in all brush-bearing lineages. The dense packing of brush hairs is explained by the specification of bristle precursor cells at a near-maximum density permitted by the lateral inhibition mechanism, as well as by the reduced size of the surrounding epithelial cells. In contrast to the female and the ancestral male condition, where bristles are arranged in stereotypical, precisely spaced rows, cell migration does not contribute appreciably to the formation of the sex brush. The complex phylogenetic history of the sex brush can make it a valuable model for investigating coevolution of sex-specific morphology and mating behavior.

Summary Genetic studies of rare and endangered species often focus on defining and preserving genetically distinct populations, especially those having unique adaptations [1, 2]. Much less attention is directed at understanding the landscape of deleterious variation, an insidious consequence of geographic isolation and the inefficiency of natural selection to eliminate harmful variants in small populations [3–5]. With population sizes of many vertebrates decreasing and isolation increasing through habitat fragmentation and loss, understanding the extent and nature of deleterious variation in small populations is essential for predicting and enhancing population persistence. The Channel Island fox ( Urocyon littoralis ) is a dwarfed species that inhabits six of California's Channel Islands and is derived from the mainland gray fox ( U. cinereoargenteus ). These isolated island populations have persisted for thousands of years at extremely small population sizes [6, 7] and, consequently, are a model for testing ideas about the accumulation of deleterious variation in small populations under natural conditions. Analysis of complete genome sequence data from island foxes shows a dramatic decrease in genome-wide variation and a sharp increase in the homozygosity of deleterious variants. The San Nicolas Island population has a near absence of variation, demonstrating a unique genetic flatlining that is punctuated by heterozygosity hotspots, enriched for olfactory receptor genes and other genes with high levels of ancestral variation. These findings question the generality of the small-population paradigm that maintains substantial genetic variation is necessary for short- and long-term persistence.

Summary The recovery and persistence of rare and endangered species are often threatened by genetic factors, such as the accumulation of deleterious mutations, loss of adaptive potential, and inbreeding depression [ 1 ]. Island foxes (Urocyon littoralis), the dwarfed descendants of mainland gray foxes (Urocyon cinereoargenteus), have inhabited California’s Channel Islands for >9,000 years [ 2 , 3 , 4 ]. Previous genomic analyses revealed that island foxes have exceptionally low levels of diversity and elevated levels of putatively deleterious variation [ 5 ]. Nonetheless, all six populations have persisted for thousands of generations, and several populations rebounded rapidly after recent severe bottlenecks [ 6 , 7 ]. Here, we combine morphological and genomic data with population-genetic simulations to determine the mechanism underlying the enigmatic persistence of these foxes. First, through analysis of genomes from 1929 to 2009, we show that island foxes have remained at small population sizes with low diversity for many generations. Second, we present morphological data indicating an absence of inbreeding depression in island foxes, confirming that they are not afflicted with congenital defects common to other small and inbred populations. Lastly, our population-genetic simulations suggest that long-term small population size results in a reduced burden of strongly deleterious recessive alleles, providing a mechanism for the absence of inbreeding depression in island foxes. Importantly, the island fox illustrates a scenario in which genetic restoration through human-assisted gene flow could be a counterproductive or even harmful conservation strategy. Our study sheds light on the puzzle of island fox persistence, a unique success story that provides a model for the preservation of small populations.

While it is appreciated that population size changes can impact patterns of deleterious variation in natural populations, less attention has been paid to how population admixture affects the dynamics of deleterious variation. Here we use population genetic simulations to examine how admixture impacts deleterious variation under a variety of demographic scenarios, dominance coefficients, and recombination rates. Our results show that gene flow between populations can temporarily reduce the genetic load of smaller populations, especially if deleterious mutations are recessive. Additionally, when fitness effects of new mutations are recessive, between-population differences in the sites at which deleterious variants exist creates heterosis in hybrid individuals. This can lead to an increase in introgressed ancestry, particularly when recombination rates are low. Under certain scenarios, introgressed ancestry can increase from an initial frequency of 5% to 30-75% and fix at many loci, even in the absence of beneficial mutations. Further, deleterious variation and admixture can generate correlations between the frequency of introgressed ancestry and recombination rate or exon density, even in the absence of other types of selection. The direction of these correlations is determined by the specific demography and whether mutations are additive or recessive. Therefore, it is essential that null models include both demography and deleterious variation before invoking reproductive incompatibilities or adaptive introgression to explain unusual patterns of genetic variation.

A long-term debate in evolutionary biology is the extent to which reproductive isolation is a necessary element of speciation. Hybridizing plants in general are cited as evidence against this notion, and oaks specifically have been used as the classic example of species maintenance without reproductive isolation. Here, we use thousands of SNPs generated by RAD sequencing to describe the phylogeny of a set of sympatric white oak species in California and then test whether these species exhibit pervasive interspecific gene exchange. Using RAD sequencing, we first constructed a phylogeny of ten oak species found in California. Our phylogeny revealed that seven scrub oak taxa occur within one clade that diverged from a common ancestor with Q. lobata, that they comprise two subclades, and they are not monophyletic but include the widespread tree oak Q. douglasii. Next, we searched for genomic patterns of allele sharing consistent with gene flow between long-divergent tree oaks with scrub oaks. Specifically, we utilized the D-statistic as well as model-based inference to compare the signature of shared alleles between two focal tree species (Q. lobata and Q. engelmannii) with multiple scrub species within the two subclades. We found that introgression is not equally pervasive between sympatric tree and scrub oak species. Instead, gene flow commonly occurs from scrub oaks to recently sympatric Q. engelmannii, but less so from scrub oaks to long-sympatric Q. lobata. This case study illustrates the influence of ancient introgression and impact of reproductive isolating mechanisms in preventing indiscriminate interspecific gene exchange.

The distribution of fitness effects (DFE) has considerable importance in population genetics. To date, estimates of the DFE come from studies using a small number of individuals. Thus, estimates of the proportion of moderately to strongly deleterious new mutations may be unreliable because such variants are unlikely to be segregating in the data. Additionally, the true functional form of the DFE is unknown, and estimates of the DFE differ significantly between studies. Here we present a flexible and computationally tractable method, called Fit∂a∂i, to estimate the DFE of new mutations using the site frequency spectrum from a large number of individuals. We apply our approach to the frequency spectrum of 1300 Europeans from the Exome Sequencing Project ESP6400 data set, 1298 Danes from the LuCamp data set, and 432 Europeans from the 1000 Genomes Project to estimate the DFE of deleterious nonsynonymous mutations. We infer significantly fewer (0.38–0.84 fold) strongly deleterious mutations with selection coefficient |s| > 0.01 and more (1.24–1.43 fold) weakly deleterious mutations with selection coefficient |s| < 0.001 compared to previous estimates. Furthermore, a DFE that is a mixture distribution of a point mass at neutrality plus a gamma distribution fits better than a gamma distribution in two of the three data sets. Our results suggest that nearly neutral forces play a larger role in human evolution than previously thought.

The distribution of fitness effects (DFE) of new mutations plays a fundamental role in evolutionary genetics. However, the extent to which the DFE differs across species has yet to be systematically investigated. Furthermore, the biological mechanisms determining the DFE in natural populations remain unclear. Here, we show that theoretical models emphasizing different biological factors at determining the DFE, such as protein stability, back-mutations, species complexity, and mutational robustness make distinct predictions about how the DFE will differ between species. Analyzing amino acid-changing variants from natural populations in a comparative population genomic framework, we find that humans have a higher proportion of strongly deleterious mutations than Drosophila melanogaster. Furthermore, when comparing the DFE across yeast, Drosophila, mice, and humans, the average selection coefficient becomes more deleterious with increasing species complexity. Last, pleiotropic genes have a DFE that is less variable than that of nonpleiotropic genes. Comparing four categories of theoretical models, only Fisher's geometrical model (FGM) is consistent with our findings. FGM assumes that multiple phenotypes are under stabilizing selection, with the number of phenotypes defining the complexity of the organism. Our results suggest that long-term population size and cost of complexity drive the evolution of the DFE, with many implications for evolutionary and medical genomics.

While it is appreciated that population size changes can impact patterns of deleterious variation in natural populations, less attention has been paid to how gene flow affects and is affected by the dynamics of deleterious variation. Here we use population genetic simulations to examine how gene flow impacts deleterious variation under a variety of demographic scenarios, mating systems, dominance coefficients, and recombination rates. Our results show that admixture between populations can temporarily reduce the genetic load of smaller populations and cause increases in the frequency of introgressed ancestry, especially if deleterious mutations are recessive. Additionally, when fitness effects of new mutations are recessive, between-population differences in the sites at which deleterious variants exist creates heterosis in hybrid individuals. Together, these factors lead to an increase in introgressed ancestry, particularly when recombination rates are low. Under certain scenarios, introgressed ancestry can increase from an initial frequency of 5% to 30–75% and fix at many loci, even in the absence of beneficial mutations. Further, deleterious variation and admixture can generate correlations between the frequency of introgressed ancestry and recombination rate or exon density, even in the absence of other types of selection. The direction of these correlations is determined by the specific demography and whether mutations are additive or recessive. Therefore, it is essential that null models of admixture include both demography and deleterious variation before invoking other mechanisms to explain unusual patterns of genetic variation., Author summary Individuals from distinct populations sometimes will produce fertile offspring and will exchange genetic material in a process called hybridization. Genomes of hybrid individuals often show non-random patterns of hybrid ancestry across the genome, where some regions have a high frequency of ancestry from the second population and other regions have less. Typically, this pattern has been attributed to adaptive introgression, where beneficial genetic variants are passed from one population to the other, or to genomic incompatibilities between these distinct species. However, other mechanisms could lead to these heterogeneous patterns of ancestry in hybrids. Here we use simulations to investigate whether deleterious mutations affect the patterns of introgressed ancestry across genomes. We show that when ancestry from a larger population is added to a smaller population, the ancestry from the larger population dramatically increases in frequency because it carries fewer deleterious mutations. This occurs even in the absence of beneficial mutations in either population. Additionally, we show that differences in sex chromosome evolution relative to autosomes, or differences in mating system, can affect patterns of introgression in similar ways. Our study argues that deleterious mutations should be included in population genetic models used to identify unusual regions of the genome that appear to be under selection in hybrids.

It has been hypothesized that the greater proportion of Neandertal ancestry in East Asians than in Europeans is due to the fact that purifying selection is less effective at removing weakly deleterious Neandertal alleles from East Asian populations. Using simulations of a broad range of models of selection and demography, we have shown that this hypothesis cannot account for the higher proportion of Neandertal ancestry in East Asians than in Europeans. Instead, more complex demographic scenarios, most likely involving multiple pulses of Neandertal admixture, are required to explain the data.