X-autosome incompatibilities in Drosophila melanogaster: tests of Haldane's rule and geographic patterns within species

Hybrid breakdown in male reproduction between recently diverged Drosophila melanogaster populations has a complex and variable genetic architecture

Secondary contact between formerly isolated populations may result in hybrid breakdown, in which untested allelic combinations in hybrids are maladaptive and limit genetic exchange. Studying early-stage reproductive isolation may yield key insights into the genetic architectures and evolutionary forces underlying the first steps toward speciation. Here, we leverage the recent worldwide expansion of Drosophila melanogaster to test for hybrid breakdown between populations that diverged within the last 13,000 years. We found clear evidence for hybrid breakdown in male reproduction, but not female reproduction or viability, supporting the prediction that hybrid breakdown affects the heterogametic sex first. The frequency of non-reproducing F2 males varied among different crosses involving the same southern African and European populations, as did the qualitative effect of cross direction, implying a genetically variable basis of hybrid breakdown and a role for uniparentally inherited factors. The levels of breakdown observed in F2 males were not recapitulated in backcrossed individuals, consistent with the existence of incompatibilities with at least three partners. Thus, some of the very first steps toward reproductive isolation could involve incompatibilities with complex and variable genetic architectures. Collectively, our findings emphasize this system's potential for future studies on the genetic and organismal basis of early-stage reproductive isolation.

Divergence of X-linked trans regulatory proteins and the misexpression of gene targets in sterile Drosophila pseudoobscura hybrids

BACKGROUND

The genetic basis of hybrid incompatibilities is characterized by pervasive cases of gene interactions. Sex chromosomes play a major role in speciation and X-linked hybrid male sterility (HMS) genes have been identified. Interestingly, some of these genes code for proteins with DNA binding domains, suggesting a capability to act as trans-regulatory elements and disturb the expression of a large number of gene targets. To understand how interactions between trans- and cis-regulatory elements contribute to speciation, we aimed to map putative X-linked trans-regulatory elements and to identify gene targets with disrupted gene expression in sterile hybrids between the subspecies Drosophila pseudoobscura pseudoobscura and D. p. bogotana.

RESULTS

We find six putative trans-regulatory proteins within previously mapped X chromosome HMS loci with sequence changes that differentiate the two subspecies. Among them, the previously characterized HMS gene Overdrive (Ovd) had the largest number of amino acid changes between subspecies, with some substitutions localized within the protein's DNA binding domain. Using an introgression approach, we detected transcriptional responses associated with a sterility/fertility Ovd allele swap. We found a network of 52 targets of Ovd and identified cis-regulatory effects among target genes with disrupted expression in sterile hybrids. However, a combined analysis of polymorphism and divergence in non-coding sequences immediately upstream of target genes found no evidence of changes in candidate regulatory proximal cis-elements. Finally, peptidases were over-represented among target genes.

CONCLUSIONS

We provide evidence of divergence between subspecies within the DNA binding domain of the HMS protein Ovd and identify trans effects on the expression of 52 gene targets. Our results identify a network of trans-cis interactions with possible effects on HMS. This network provides molecular evidence of gene × gene incompatibilities as contributors to hybrid dysfunction.

Estimating the Timing of Multiple Admixture Pulses During Local Ancestry Inference

Admixture, the mixing of genetically distinct populations, is increasingly recognized as a fundamental biological process. One major goal of admixture analyses is to estimate the timing of admixture events. Whereas most methods today can only detect the most recent admixture event, here, we present coalescent theory and associated software that can be used to estimate the timing of multiple admixture events in an admixed population. We extensively validate this approach and evaluate the conditions under which it can successfully distinguish one- from two-pulse admixture models. We apply our approach to real and simulated data of Drosophila melanogaster We find evidence of a single very recent pulse of cosmopolitan ancestry contributing to African populations, as well as evidence for more ancient admixture among genetically differentiated populations in sub-Saharan Africa. These results suggest our method can quantify complex admixture histories involving genetic material introduced by multiple discrete admixture pulses. The new method facilitates the exploration of admixture and its contribution to adaptation, ecological divergence, and speciation.

Deleterious variation shapes the genomic landscape of introgression

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.

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.

Moving Speciation Genetics Forward: Modern Techniques Build on Foundational Studies in Drosophila

The question of how new species evolve has been examined at every level, from macroevolutionary patterns of diversification to molecular population genetic analyses of specific genomic regions between species pairs. Drosophila has been at the center of many of these research efforts. Though our understanding of the speciation process has grown considerably over the past few decades, very few genes have been identified that contribute to barriers to reproduction. The development of advanced molecular genetic and genomic methods provides promising avenues for the rapid discovery of more genes that contribute to speciation, particularly those involving prezygotic isolation. The continued expansion of tools and resources, especially for species other than Drosophila melanogaster, will be most effective when coupled with comparative approaches that reveal the genetic basis of reproductive isolation across a range of divergence times. Future research programs in Drosophila have high potential to answer long-standing questions in speciation. These include identifying the selective forces that contribute to divergence between populations and the genetic basis of traits that cause reproductive isolation. The latter can be expanded upon to understand how the genetic basis of reproductive isolation changes over time and whether certain pathways and genes are more commonly involved.

A Hidden Markov Model Approach for Simultaneously Estimating Local Ancestry and Admixture Time Using Next Generation Sequence Data in Samples of Arbitrary Ploidy

Admixture-the mixing of genomes from divergent populations-is increasingly appreciated as a central process in evolution. To characterize and quantify patterns of admixture across the genome, a number of methods have been developed for local ancestry inference. However, existing approaches have a number of shortcomings. First, all local ancestry inference methods require some prior assumption about the expected ancestry tract lengths. Second, existing methods generally require genotypes, which is not feasible to obtain for many next-generation sequencing projects. Third, many methods assume samples are diploid, however a wide variety of sequencing applications will fail to meet this assumption. To address these issues, we introduce a novel hidden Markov model for estimating local ancestry that models the read pileup data, rather than genotypes, is generalized to arbitrary ploidy, and can estimate the time since admixture during local ancestry inference. We demonstrate that our method can simultaneously estimate the time since admixture and local ancestry with good accuracy, and that it performs well on samples of high ploidy-i.e. 100 or more chromosomes. As this method is very general, we expect it will be useful for local ancestry inference in a wider variety of populations than what previously has been possible. We then applied our method to pooled sequencing data derived from populations of Drosophila melanogaster on an ancestry cline on the east coast of North America. We find that regions of local recombination rates are negatively correlated with the proportion of African ancestry, suggesting that selection against foreign ancestry is the least efficient in low recombination regions. Finally we show that clinal outlier loci are enriched for genes associated with gene regulatory functions, consistent with a role of regulatory evolution in ecological adaptation of admixed D. melanogaster populations. Our results illustrate the potential of local ancestry inference for elucidating fundamental evolutionary processes.

Spatially Heterogeneous Environmental Selection Strengthens Evolution of Reproductively Isolated Populations in a Dobzhansky-Muller System of Hybrid Incompatibility

Within-species hybrid incompatibility can arise when combinations of alleles at more than one locus have low fitness but where possession of one of those alleles has little or no fitness consequence for the carriers. Limited dispersal with small numbers of mate potentials alone can lead to the evolution of clusters of reproductively isolated genotypes despite the absence of any geographical barriers or heterogeneous selection. In this paper, we explore how adding heterogeneous natural selection on the genotypes (e.g., gene environment associations) that are involved in reproductive incompatibility affects the frequency, size and duration of evolution of reproductively isolated clusters. We conducted a simulation experiment that varied landscape heterogeneity, dispersal ability, and strength of selection in a continuously distributed population. In our simulations involving spatially heterogeneous selection, strong patterns of adjacency of mutually incompatible genotypes emerged such that these clusters were truly reproductively isolated from each other, with no reproductively compatible “bridge” individuals in the intervening landscape to allow gene flow between the clusters. This pattern was strong across levels of gene flow and strength of selection, suggesting that even relatively weak selection acting in the context of strong gene flow may produce reproductively isolated clusters that are large and persistent, enabling incipient speciation in a continuous population without geographic isolation.

Secondary contact and local adaptation contribute to genome-wide patterns of clinal variation in Drosophila melanogaster

Populations arrayed along broad latitudinal gradients often show patterns of clinal variation in phenotype and genotype. Such population differentiation can be generated and maintained by both historical demographic events and local adaptation. These evolutionary forces are not mutually exclusive and can in some cases produce nearly identical patterns of genetic differentiation among populations. Here, we investigate the evolutionary forces that generated and maintain clinal variation genome-wide among populations of Drosophila melanogaster sampled in North America and Australia. We contrast patterns of clinal variation in these continents with patterns of differentiation among ancestral European and African populations. Using established and novel methods we derive here, we show that recently derived North America and Australia populations were likely founded by both European and African lineages and that this hybridization event likely contributed to genome-wide patterns of parallel clinal variation between continents. The pervasive effects of admixture mean that differentiation at only several hundred loci can be attributed to the operation of spatially varying selection using an FST outlier approach. Our results provide novel insight into the well-studied system of clinal differentiation in D. melanogaster and provide a context for future studies seeking to identify loci contributing to local adaptation in a wide variety of organisms, including other invasive species as well as temperate endemics.

The Mosaic Ancestry of the Drosophila Genetic Reference Panel and the D. melanogaster Reference Genome Reveals a Network of Epistatic Fitness Interactions

North American populations of Drosophila melanogaster derive from both European and African source populations, but despite their importance for genetic research, patterns of ancestry along their genomes are largely undocumented. Here, I infer geographic ancestry along genomes of the Drosophila Genetic Reference Panel (DGRP) and the D. melanogaster reference genome, which may have implications for reference alignment, association mapping, and population genomic studies in Drosophila. Overall, the proportion of African ancestry was estimated to be 20% for the DGRP and 9% for the reference genome. Combining my estimate of admixture timing with historical records, I provide the first estimate of natural generation time for this species (approximately 15 generations per year). Ancestry levels were found to vary strikingly across the genome, with less African introgression on the X chromosome, in regions of high recombination, and at genes involved in specific processes (e.g., circadian rhythm). An important role for natural selection during the admixture process was further supported by evidence that many unlinked pairs of loci showed a deficiency of Africa–Europe allele combinations between them. Numerous epistatic fitness interactions may therefore exist between African and European genotypes, leading to ongoing selection against incompatible variants. By focusing on hubs in this network of fitness interactions, I identified a set of interacting loci that include genes with roles in sensation and neuropeptide/hormone reception. These findings suggest that admixed D. melanogaster samples could become an important study system for the genetics of early-stage isolation between populations.

Clusters of incompatible genotypes evolve with limited dispersal

Theoretical and empirical studies have shown heterogeneous selection to be the primary driver for the evolution of reproductively isolated genotypes in the absence of geographic barriers. Here, we ask whether limited dispersal alone can lead to the evolution of reproductively isolated genotypes despite the absence of any geographic barriers or heterogeneous selection. We use a spatially-explicit, individual-based, landscape genetics program to explore the influences of dispersal strategies on reproductive isolation. We simulated genetic structure in a continuously distributed population and across various dispersal strategies (ranging from short- to long-range individual movement), as well as potential mate partners in entire population (ranging from 20 to 5000 individuals). We show that short-range dispersal strategies lead to the evolution of clusters of reproductively isolated genotypes despite the absence of any geographic barriers or heterogeneous selection. Clusters of genotypes that are reproductively isolated from other clusters can persist when migration distances are restricted such that the number of mating partners is below about 350 individuals. We discuss how our findings may be applicable to particular speciation scenarios, as well as the need to investigate the influences of heterogeneous gene flow and spatial selection gradients on the emergence of reproductively isolating genotypes.

Species-wide genetic incompatibility analysis identifies immune genes as hot spots of deleterious epistasis

Intraspecific genetic incompatibilities prevent the assembly of specific alleles into single genotypes and influence genome- and species-wide patterns of sequence variation. A common incompatibility in plants is hybrid necrosis, characterized by autoimmune responses due to epistatic interactions between natural genetic variants. By systematically testing thousands of F1 hybrids of Arabidopsis thaliana strains, we identified a small number of incompatibility hot spots in the genome, often in regions densely populated by nucleotide-binding domain and leucine-rich repeat (NLR) immune receptor genes. In several cases, these immune receptor loci interact with each other, suggestive of conflict within the immune system. A particularly dangerous locus is a highly variable cluster of NLR genes, DM2, which causes multiple independent incompatibilities with genes that encode a range of biochemical functions, including NLRs. Our findings suggest that deleterious interactions of immune receptors limit the combinations of favorable disease resistance alleles accessible to plant genomes.

Genetic background and GxE interactions modulate the penetrance of a naturally occurring wing mutation in Drosophila melanogaster

Many genes involved in producing complex traits are incompletely penetrant. One such example is vesiculated, an X-linked gene in Drosophila melanogaster that results in wing defects. To examine the genetic architecture of a complex trait (wings containing vesicles), we placed a naturally occurring variant into multiple autosomal backgrounds and quantified penetrance and expressivity at a range of developmental temperatures. We found significant epistasis, genotype-by-environment interactions, and maternal effects. Sex and temperature effects were modulated by genetic background. The severity of wing phenotypes also varied across different genetic backgrounds, and expressivity was positively correlated with penetrance. We also found evidence of naturally segregating suppressors of vesiculated. These suppressors were present on both the second and third chromosomes, and complex interactions were observed. Taken together, these findings indicate that multiple genetic and environmental factors modulate the phenotypic effects of a naturally occurring vesiculated allele.

The genetics of sex chromosomes: evolution and implications for hybrid incompatibility

Heteromorphic sex chromosomes, where one sex has two different types of sex chromosomes, face very different evolutionary consequences than do autosomes. Two important features of sex chromosomes arise from being present in only one copy in one of the sexes: dosage compensation and the meiotic silencing of sex chromosomes. Other differences arise because sex chromosomes spend unequal amounts of time in each sex. Thus, the impact of evolutionary processes (mutation, selection, genetic drift, and meiotic drive) differs substantially between each sex chromosome, and between the sex chromosomes and the autosomes. Sex chromosomes also play a disproportionate role in Haldane's rule and other important patterns related to hybrid incompatibility, and thus speciation. We review the consequences of sex chromosomes on hybrid incompatibility. A theme running through this review is that epigenetic processes, notably those related to chromatin, may be more important to the evolution of sex chromosomes and the evolution of hybrid incompatibility than previously recognized.

The population genetics of X-autosome synthetic lethals and steriles

Epistatic interactions are widespread, and many of these interactions involve combinations of alleles at different loci that are deleterious when present in the same individual. The average genetic environment of sex-linked genes differs from that of autosomal genes, suggesting that the population genetics of interacting X-linked and autosomal alleles may be complex. Using both analytical theory and computer simulations, we analyzed the evolutionary trajectories and mutation-selection balance conditions for X-autosome synthetic lethals and steriles. Allele frequencies follow a set of fundamental trajectories, and incompatible alleles are able to segregate at much higher frequencies than single-locus expectations. Equilibria exist, and they can involve fixation of either autosomal or X-linked alleles. The exact equilibrium depends on whether synthetic alleles are dominant or recessive and whether fitness effects are seen in males, females, or both sexes. When single-locus fitness effects and synthetic incompatibilities are both present, population dynamics depend on the dominance of alleles and historical contingency (i.e., whether X-linked or autosomal mutations occur first). Recessive synthetic lethality can result in high-frequency X-linked alleles, and dominant synthetic lethality can result in high-frequency autosomal alleles. Many X-autosome incompatibilities in natural populations may be cryptic, appearing to be single-locus effects because one locus is fixed. We also discuss the implications of these findings with respect to standing genetic variation and the origins of Haldane's rule.