Incompatibility between X chromosome factor and pericentric heterochromatic region causes lethality in hybrids between Drosophila melanogaster and its sibling species

The Dobzhansky-Muller model posits that postzygotic reproductive isolation results from the evolution of incompatible epistatic interactions between species: alleles that function in the genetic background of one species can cause sterility or lethality in the genetic background of another species. Progress in identifying and characterizing factors involved in postzygotic isolation in Drosophila has remained slow, mainly because Drosophila melanogaster, with all of its genetic tools, forms dead or sterile hybrids when crossed to its sister species, D. simulans, D. sechellia, and D. mauritiana. To circumvent this problem, we used chromosome deletions and duplications from D. melanogaster to map two hybrid incompatibility loci in F(1) hybrids with its sister species. We mapped a recessive factor to the pericentromeric heterochromatin of the X chromosome in D. simulans and D. mauritiana, which we call heterochromatin hybrid lethal (hhl), which causes lethality in F(1) hybrid females with D. melanogaster. As F(1) hybrid males hemizygous for a D. mauritiana (or D. simulans) X chromosome are viable, the lethality of deficiency hybrid females implies that a dominant incompatible partner locus exists on the D. melanogaster X. Using small segments of the D. melanogaster X chromosome duplicated onto the Y chromosome, we mapped a dominant factor that causes hybrid lethality to a small 24-gene region of the D. melanogaster X. We provide evidence suggesting that it interacts with hhl(mau). The location of hhl is consistent with the emerging theme that hybrid incompatibilities in Drosophila involve heterochromatic regions and factors that interact with the heterochromatin.     

Genetic complexity underlying hybrid male sterility in Drosophila

Recent genetic analyses of closely related species of Drosophila have indicated that hybrid male sterility is the consequence of highly complex synergistic effects among multiple genes, both conspecific and heterospecific. On the contrary, much evidence suggests the presence of major genes causing hybrid female sterility and inviability in the less-related species, D. melanogaster and D. simulans. Does this contrast reflect the genetic distance between species? Or, generally, is the genetic basis of hybrid male sterility more complex than that of hybrid female sterility and inviability? To clarify this point, the D. simulans introgression of the cytological region 34D-36A to the D. melanogaster genome, which causes recessive male sterility, was dissected by recombination, deficiency, and complementation mapping. The 450-kb region between two genes, Suppressor of Hairless and snail, exhibited a strong effect on the sterility. Males are (semi-)sterile if this region of the introgression is made homozygous or hemizygous. But no genes in the region singly cause the sterility; this region has at least two genes, which in combination result in male sterility. Further, the males are less fertile when heterozygous with a larger introgression, which suggests that dominant modifiers enhance the effects of recessive genes of male sterility. Such an epistatic view, even in the less-related species, suggests that the genetic complexity is special to hybrid male sterility.     

Nup96-dependent hybrid lethality occurs in a subset of species from the simulans clade of Drosophila

The cross of Drosophila melanogaster females to D. simulans males typically produces lethal F(1) hybrid males. F(1) male lethality is suppressed when the D. simulans Lhr(1) hybrid rescue strain is used. Viability of these F(1) males carrying Lhr(1) is in turn substantially reduced when the hybrids are heterozygous for some mutant alleles of the D. melanogaster Nup96 gene. I show here that similar patterns of Nup96-dependent lethality occur when other hybrid rescue mutations are used to create F(1) males, demonstrating that Nup96 does not reduce hybrid viability by suppressing the Lhr(1) rescue effect. The penetrance of this Nup96-dependent lethality does not correlate with the penetrance of the F(1) hybrid rescue, arguing that these two phenomena reflect genetically independent processes. D. simulans, together with two additional sister species, forms a clade that speciated after the divergence of their common ancestor from D. melanogaster. I report here that Nup96(-) reduces F(1) viability in D. melanogaster hybrids with one of these sister species, D. sechellia, but not with the other, D. mauritiana. These results suggest that Nup96-dependent lethality evolved after the speciation of D. melanogaster from the common ancestor of the simulans clade and is caused by an interaction among Nup96, unknown gene(s) on the D. melanogaster X chromosome, and unknown autosomal gene(s), at least some of which have diverged in D. simulans and D. sechellia but not in D. mauritiana. The genetic properties of Nup96 are also discussed relative to other hybrid lethal genes.     

Genetic dissection of hybrid incompatibilities between Drosophila simulans and D. mauritiana. II. Mapping hybrid male sterility loci on the third chromosome

Hybrid male sterility (HMS) is a rapidly evolving mechanism of reproductive isolation in Drosophila. Here we report a genetic analysis of HMS in third-chromosome segments of Drosophila mauritiana that were introgressed into a D. simulans background. Qualitative genetic mapping was used to localize 10 loci on 3R and a quantitative trait locus (QTL) procedure (multiple-interval mapping) was used to identify 19 loci on the entire chromosome. These genetic incompatibilities often show dominance and complex patterns of epistasis. Most of the HMS loci have relatively small effects and generally at least two or three of them are required to produce complete sterility. Only one small region of the third chromosome of D. mauritiana by itself causes a high level of infertility when introgressed into D. simulans. By comparison with previous studies of the X chromosome, we infer that HMS loci are only approximately 40% as dense on this autosome as they are on the X chromosome. These results are consistent with the gradual evolution of hybrid incompatibilities as a by-product of genetic divergence in allopatric populations.     

Intron length evolution in Drosophila

I present data on the evolution of intron lengths among 3 closely related Drosophila species, D. melanogaster, Drosophila simulans, and Drosophila yakuba. Using D. yakuba as an outgroup, I mapped insertion and deletion mutations in 148 introns (spanning approximately 30 kb) to the D. melanogaster and D. simulans lineages. Intron length evolution in the 2 sister species has been different: in D. melanogaster, X-linked introns have increased slightly in size, whereas autosomal ones have decreased slightly in size; in D. simulans, both X-linked and autosomal introns have decreased in size. To understand the possible evolutionary causes of these lineage- and chromosome-specific patterns of intron evolution, I studied insertion-deletion (indel) polymorphism and divergence in D. melanogaster. Small insertion mutations segregate at elevated frequencies and enjoy elevated probabilities of fixation, particularly on the X chromosome. In contrast, there is no detectable X chromosome effect on fixations in D. simulans. These findings suggest X chromosome-specific selection or biased gene conversion-gap repair favoring insertions in D. melanogaster but not in D. simulans. These chromosome- and lineage-specific patterns of indel substitution are not easily explained by existing general population genetic models of intron length evolution. Genomic data from D. melanogaster further suggest that the forces described here affect introns and intergenic regions similarly.     

Mapping and characterization of a 'speciation gene' in Drosophila.

Almost nothing is known about the identity of the genes causing reproductive isolation between species. As a first step towards molecular isolation of a 'speciation gene', I mapped and partly characterized a gene causing hybrid male sterility in Drosophila. This analysis shows that sterility of D. melanogaster males who carry the 'dot' fourth chromosome from D. simulans is due entirely to a very small region of the D. simulans chromosome (including only about 5 salivary gland bands or approximately 250 kb of DNA). Thus the hybrid sterility effect of the D. simulans fourth chromosome is almost surely due to a single gene of very large effect (here named hms, hybrid male sterile). Hms is zygotically acting, and the D. simulans allele of hms is completely recessive. Furthermore, complementation tests suggest that hms is not an allele of any known locus in D. melanogaster.     

A Genetic Basis for the Inviability of Hybrids between Sibling Species of Drosophila

A mutation of Drosophila melanogaster whose only known effect is the rescue of otherwise lethal interspecific hybrids has been characterized. This mutation, Hmr, maps to 1-31.84 (9D1-9E4). Hmr may be the consequence of a P element insertion. It rescues hybrid males from the cross of D. melanogaster females to males of its three sibling species, D. simulans, D. mauritiana and D. sechellia. This rescue is recessive, since hybrid males that carry both Hmr and a duplication expected to be Hmr+ are not rescued. Hmr also rescues the otherwise inviable female hybrids from the cross of compound-X D. melanogaster females to males of its sibling species. This rescue is also recessive, since a compound-X heterozygous for Hmr does not rescue. Another mutation, discovered on the In(1)AB chromosome of D. melanogaster, is also found to rescue normally inviable species hybrids: unlike Hmr, however, In(1)AB rescues hybrid females from the cross of In(1)AB/Y males to sibling females, as well as hybrid males from the cross of In(1)AB females to sibling males. These data are interpreted on the basis of a model for the genetic basis of hybrid inviability of complementary genes.     

Genetic screens for factors involved in the notum bristle loss of interspecific hybrids between Drosophila melanogaster and D. simulans

Interspecific cross is a powerful means to uncover hidden within- and between-species variation in populations. One example is a bristle loss phenotype of hybrids between Drosophila melanogaster and D. simulans, although both the pure species have exactly the same pattern of bristle formation on the notum. There exists a large amount of genetic variability in the simulans populations with respect to the number of missing bristles in hybrids, and the variation is largely attributable to simulans X chromosomes. Using nine molecular markers, I screened the simulans X chromosome for genetic factors that were responsible for the differences between a pair of simulans lines with high (H) and low (L) missing bristle numbers. Together with duplication-rescue experiments, a single major quantitative locus was mapped to a 13F-14F region. Importantly, this region accounted for most of the differences between H and L lines in three other independent pairs, suggesting segregation of H and L alleles at the single locus in different populations. Moreover, a deficiency screening uncovered several regions with factors that potentially cause the hybrid bristle loss due to epistatic interactions with the other factors.     

Inviability of hybrids between D. melanogaster and D. simulans results from the absence of simulans X not the presence of simulans Y chromosome.

Interspecific crosses between D. melanogaster and D. simulans or its sibling species result in unisexual inviability of the hybrids. Mostly, crosses of D. melanogaster females x D. simulans males produce hybrid females. On the other hand, only hybrid males are viable in the reciprocal crosses. A classical question is the cause of the unisexual hybrid inviability on the chromosomal level. Is it due to the absence of a D. simulans X chromosome or is it due to the presence of a D. simulans Y chromosome? A lack of adequate chromosomal rearrangements available in D. simulans has made it difficult to answer this question. However, it has been assumed that the lethality results from the absence of the D. simulans X rather than the presence of the D. simulans Y. Recently I synthesized the first D. simulans compound-XY chromosome that consists of almost the entire X and Y chromosomes. Males carrying the compound-XY and no free Y chromosome are fertile. By utilizing the compound-XY chromosome, the viability of hybrids with various constitutions of cytoplasm and sex chromosomes has been examined. The results consistently demonstrate that the absence of a D. simulans X chromosome in hybrid genome, and not the presence of the Y chromosome, is a determinant of the hybrid inviability.     

The genetics of inviability and male sterility in hybrids between Anopheles gambiae and An. arabiensis

Male hybrids between Anopheles gambiae and An. arabiensis suffer from hybrid sterility, and inviability effects are sometimes present as well. We examined the genetic basis of these reproductive barriers between the two species, using 21 microsatellite markers. Generally, recessive inviability effects were found on the X chromosome of gambiae that are incompatible with at least one factor on each arabiensis autosome. Inviability is complete when the gambiae and arabiensis inviability factors are hemi- or homozygous. Using a QTL mapping approach, regions that contribute to male hybrid sterility were also identified. The X chromosome has a disproportionately large effect on male hybrid sterility. Additionally, several moderate-to-large autosomal QTL were found in both species. The effect of these autosomal QTL is contingent upon the presence of an X chromosome from the other species. Substantial regions of the autosomes do not contribute markedly to male hybrid sterility. Finally, no evidence for epistatic interactions between conspecific sterility loci was found.     

Genetic dissection of hybrid incompatibilities between Drosophila simulans and D. mauritiana. I. Differential accumulation of hybrid male sterility effects on the X and autosomes

The genetic basis of hybrid incompatibility in crosses between Drosophila mauritiana and D. simulans was investigated to gain insight into the evolutionary mechanisms of speciation. In this study, segments of the D. mauritiana third chromosome were introgressed into a D. simulans genetic background and tested as homozygotes for viability, male fertility, and female fertility. The entire third chromosome was covered with partially overlapping segments. Many segments were male sterile, while none were female sterile or lethal, confirming previous reports of the rapid evolution of hybrid male sterility (HMS). A statistical model was developed to quantify the HMS accumulation. In comparison with previous work on the X chromosome, we estimate that the X has approximately 2.5 times the density of HMS factors as the autosomes. We also estimate that the whole genome contains approximately 15 HMS "equivalents"-i.e., 15 times the minimum number of incompatibility factors necessary to cause complete sterility. Although some caveats for the quantitative estimate of a 2.5-fold density difference are described, this study supports the notion that the X chromosome plays a special role in the evolution of reproductive isolation. Possible mechanisms of a "large X" effect include selective fixation of new mutations that are recessive or partially recessive and the evolution of sex-ratio distortion systems.     

Ephemeral Association Between Gene CG5762 and Hybrid Male Sterility in Drosophila Sibling Species

Interspecies divergence in regulatory pathways may result in hybrid male sterility (HMS) when dominance and epistatic interactions between alleles that are functional within one genome are disrupted in hybrid genomes. The identification of genes contributing to HMS and other hybrid dysfunctions is essential for understanding the origin of new species (speciation). Previously, we identified a panel of male-specific loci misexpressed in sterile male hybrids of Drosophila simulans and D. mauritiana relative to parental species. In the current work, we attempt to dissect the genetic associations between HMS and one of the genes, CG5762, a Drosophila-unique locus characterized by rapid sequence divergence within the genus, presumably driven by positive natural selection. CG5762 is underexpressed in sterile backcross males compared with their fertile brothers. In CG5762 heterozygotes, the D. mauritiana allele is consistently overexpressed on both the D. simulans and D. mauritiana backcross genomic background, suggesting a cis-acting regulation factor. There is a significant association between heterozygosity and HMS in hybrid males from early but not later backcross generations. Microsatellite markers spanning CG5762 fail to associate with HMS. These observations lead to a conclusion that CG5762 is not a causative factor of HMS. Although genetic linkage between CG5762 and a neighboring causative introgression cannot be ruled out, it seems that the pattern is most consistent with CG5762 participating in epistatic interactions that are disrupted in flies with HMS.     

Finding Hybrid Incompatibilities Using Genome Sequences from Hybrid Populations

Natural hybrid zones offer a powerful framework for understanding the genetic basis of speciation in progress because ongoing hybridization continually creates unfavorable gene combinations. Evidence indicates that postzygotic reproductive isolation is often caused by epistatic interactions between mutations in different genes that evolved independently of one another (hybrid incompatibilities). We examined the potential to detect epistatic selection against incompatibilities from genome sequence data using the site frequency spectrum (SFS) of polymorphisms by conducting individual-based simulations in SLiM. We found that the genome-wide SFS in hybrid populations assumes a diagnostic shape, with the continual input of fixed differences between source populations via migration inducing a mass at intermediate allele frequency. Epistatic selection locally distorts the SFS as non-incompatibility alleles rise in frequency in a manner analogous to a selective sweep. Building on these results, we present a statistical method to identify genomic regions containing incompatibility loci that locates departures in the local SFS compared with the genome-wide SFS. Cross-validation studies demonstrate that our method detects recessive and codominant incompatibilities across a range of scenarios varying in the strength of epistatic selection, migration rate, and hybrid zone age. Our approach takes advantage of whole genome sequence data, does not require knowledge of demographic history, and can be applied to any pair of nascent species that forms a hybrid zone.     

Hybrid Dysgenesis in Drosophila Simulans Lines Transformed with Autonomous P Elements

The molecular and phenotypic analysis of several previously described P element-transformed lines of Drosophila simulans was extended in order to determine whether they had the potential to produce a syndrome of P-M hybrid dysgenesis analogous to the one in Drosophila melanogaster. The transformed line with the highest number of P elements at the beginning of the analysis, DsP pi-5C, developed strong P activity potential and P element regulation, properties characteristic of D. melanogaster P strains. The subsequent analysis of sublines derived from 34 single pair matings of DsP pi-5C revealed that they were heterogeneous with respect to both their P element complements and P activity potentials, but similar with respect to their regulatory capabilities. The subline with the highest P activity, DsP pi-5C-27, was subsequently used as a reference P strain in the genetic analysis of the D. simulans transformants. In these experiments, the reciprocal cross effect was observed with respect to both gonadal sterility and male recombination. As in D. melanogaster, the induction of gonadal sterility in D. simulans was shown to be temperature-dependent. Molecular analysis of DsP pi-5C-27 revealed that it has approximately 30 P elements per genome, at least some of which are defective. The number of potentially complete P elements in its genome is similar to the number in the D. melanogaster P strain, Harwich-77. Overall our analysis indicates that P-transformed lines of D. simulans are capable of expressing the major features of P-M hybrid dysgenesis previously demonstrated in D. melanogaster and that P elements appear to behave in a similar way in the two sibling species.     

Contrasting patterns of X-linked and autosomal nucleotide variation in Drosophila melanogaster and Drosophila simulans

Surveys of molecular variation in Drosophila melanogaster and Drosophila simulans have suggested that diversity outside of Africa is a subset of that within Africa. It has been argued that reduced levels of diversity in non-African populations reflect a population bottleneck, adaptation to temperate climates, or both. Here, I summarize the available single-nucleotide polymorphism data for both species. A simple "out of Africa" bottleneck scenario is consistent with geographic patterns for loci on the X chromosome but not with loci on the autosomes. Interestingly, there is a trend toward lower nucleotide diversity on the X chromosome relative to autosomes in non-African populations of D. melanogaster, but the opposite trend is seen in African populations. In African populations, autosomal inversion polymorphisms in D. melanogaster may contribute to reduced autosome diversity relative to the X chromosome. To elucidate the role that selection might play in shaping patterns of variability, I present a summary of within- and between-species patterns of synonymous and replacement variation in both species. Overall, D. melanogaster autosomes harbor an excess of amino acid replacement polymorphisms relative to D. simulans. Interestingly, range expansion from Africa appears to have had little effect on synonymous-to-replacement polymorphism ratios.     

Cis-by-Trans regulatory divergence causes the asymmetric lethal effects of an ancestral hybrid incompatibility gene

The Dobzhansky and Muller (D-M) model explains the evolution of hybrid incompatibility (HI) through the interaction between lineage-specific derived alleles at two or more loci. In agreement with the expectation that HI results from functional divergence, many protein-coding genes that contribute to incompatibilities between species show signatures of adaptive evolution, including Lhr, which encodes a heterochromatin protein whose amino acid sequence has diverged extensively between Drosophila melanogaster and D. simulans by natural selection. The lethality of D. melanogaster/D. simulans F1 hybrid sons is rescued by removing D. simulans Lhr, but not D. melanogaster Lhr, suggesting that the lethal effect results from adaptive evolution in the D. simulans lineage. It has been proposed that adaptive protein divergence in Lhr reflects antagonistic coevolution with species-specific heterochromatin sequences and that defects in LHR protein localization cause hybrid lethality. Here we present surprising results that are inconsistent with this coding-sequence-based model. Using Lhr transgenes expressed under native conditions, we find no evidence that LHR localization differs between D. melanogaster and D. simulans, nor do we find evidence that it mislocalizes in their interspecific hybrids. Rather, we demonstrate that Lhr orthologs are differentially expressed in the hybrid background, with the levels of D. simulans Lhr double that of D. melanogaster Lhr. We further show that this asymmetric expression is caused by cis-by-trans regulatory divergence of Lhr. Therefore, the non-equivalent hybrid lethal effects of Lhr orthologs can be explained by asymmetric expression of a molecular function that is shared by both orthologs and thus was presumably inherited from the ancestral allele of Lhr. We present a model whereby hybrid lethality occurs by the interaction between evolutionarily ancestral and derived alleles.     

GENETIC ARCHITECTURE AND POSTZYGOTIC REPRODUCTIVE ISOLATION: EVOLUTION OF BATESON–DOBZHANSKY–MULLER INCOMPATIBILITIES IN A POLYGENIC MODEL

The Bateson–Dobzhansky–Muller model predicts that postzygotic isolation evolves due to the accumulation of incompatible epistatic interactions, but few studies have quantified the relationship between genetic architecture and patterns of reproductive divergence. We examined how the direction and magnitude of epistatic interactions in a polygenic trait under stabilizing selection influenced the evolution of hybrid incompatibilities. We found that populations evolving independently under stabilizing selection experienced suites of compensatory allelic changes that resulted in genetic divergence between populations despite the maintenance of a stable, high‐fitness phenotype. A small number of loci were then incompatible with multiple alleles in the genetic background of the hybrid and the identity of these incompatibility loci changed over the evolution of the populations. For F₁ hybrids, reduced fitness evolved in a window of intermediate strengths of epistatic interactions, but F₂ and backcross hybrids evolved reduced fitness across weak and moderate strengths of epistasis due to segregation variance. Strong epistatic interactions constrained the allelic divergence of parental populations and prevented the development of reproductive isolation. Because many traits with varying genetic architectures must be under stabilizing selection, our results indicate that polygenetic drift is a plausible hypothesis for the evolution of postzygotic reproductive isolation.     

Loss of notum macrochaetae as an interspecific hybrid anomaly between Drosophila melanogaster and D. simulans.

With the aim of revealing genetic variation accumulated among closely related species during the course of evolution, this study focuses on loss of macrochaetae on the notum as one of the developmental anomalies seen in interspecific hybrids between Drosophila melanogaster and its closely related species. Interspecific hybrids between a line of D. melanogaster and D. simulans isofemale lines exhibited a wide range in the number of missing bristles. By contrast, D. mauritiana and D. sechellia lines showed almost no reduction in bristle number in hybrids with D. melanogaster. Genetic analysis showed that the D. simulans X chromosome confers a large effect on hybrid bristle loss, although X-autosome interaction may be involved. This suggests that at least one genetic factor contributing to hybrid anomalies arose recently on a D. simulans X chromosome. Moreover, the results indicate sex dependency: the male hybrids were more susceptible to bristle loss than the female hybrids were. Use of cell type markers suggests that the defect does not lie in cell fate decisions during bristle development, but in the maintenance of neural fate and/or differentiation of the descendants of sensory mother cells.     

The impact of shared ancestral variation on hybrid male lethality--a 16 codon indel in the Drosophila simulans Lhr gene

The Lethal hybrid rescue (Lhr) gene causes hybrid male lethality in crosses between Drosophila simulans and D. melanogaster. Lhr(2) is a D. simulans allele, which rescues hybrid males. It has been recently proposed that a 16 codon insertion, which distinguishes the D. melanogaster and the canonical D. simulans allele, and is lacking in Lhr(2), may be responsible for the functional divergence of D. melanogaster and D. simulans Lhr alleles. Here, we show that the Lhr(2) allele lacking the insertion represents an ancestral polymorphism segregating at a moderate frequency in D. simulans. Crosses of D. melanogaster females to males from two D. simulans strains carrying this deletion showed a severe deficiency of viable hybrid males. Our results suggest that the absence of this insertion alone is not sufficient to explain functional differences between D. melanogaster and D. simulans Lhr alleles.     

Genetics of Hybrid Male Sterility between Drosophila Sibling Species: A Complex Web of Epistasis Is Revealed in Interspecific Studies

To study the genetic differences responsible for the sterility of their male hybrids, we introgressed small segments of an X chromosome from Drosophila simulans into a pure Drosophila mauritiana genetic background, then assessed the fertility of males carrying heterospecific introgressions of varying size. Although this analysis examined less than 20% of the X chromosome (roughly 5% of the euchromatic portion of the D. simulans genome), and the segments were introgressed in only one direction, a minimum of four factors that contribute to hybrid male sterility were revealed. At least two of the factors exhibited strong epistasis: males carrying either factor alone were consistently fertile, whereas males carrying both factors together were always sterile. Distinct spermatogenic phenotypes were observed for sterile introgressions of different lengths, and it appeared that an interaction between introgressed segments also influenced the stage of spermatogenic defect. Males with one category of introgression often produced large quantities of motile sperm and were observed copulating, but never inseminated females. Evidently these two species have diverged at a large number of loci which have varied effects on hybrid male fertility. By extrapolation, we estimate that there are at least 40 such loci on the X chromosome alone. Because these species exhibit little DNA-sequence divergence at arbitrarily chosen loci, it seems unlikely that the extensive functional divergence observed could be due mainly to random genetic drift. Significant epistasis between conspecific genes appears to be a common component of hybrid sterility between recently diverged species of Drosophila. The linkage relationships of interacting factors could shed light on the role played by epistatic selection in the dynamics of the allele substitutions responsible for reproductive barriers between species.