An essential cell cycle regulation gene causes hybrid inviability in Drosophila

See allHide authors and affiliations

Science  18 Dec 2015:
Vol. 350, Issue 6267, pp. 1552-1555
DOI: 10.1126/science.aac7504

Essential genes and species incompatibilities

Crosses between two fruit fly species, Drosophila melanogaster and D. simulans, result in hybrid progeny that are all female. Although some of the genes responsible for this species barrier are known, the full complement of molecular determinants that lead to inviable males has remained mysterious. Phadnis et al. used mutagenesis and a sequencing-based genomic screen to link hybrid inviability to the cell cycle. The inviable males result from an interaction between three genes, one of which is essential, which precluded its identification with standard genetic screens. This strategy to identify speciation genes can be applied to other model and nonmodel systems.

Science, this issue p. 1552


Speciation, the process by which new biological species arise, involves the evolution of reproductive barriers, such as hybrid sterility or inviability between populations. However, identifying hybrid incompatibility genes remains a key obstacle in understanding the molecular basis of reproductive isolation. We devised a genomic screen, which identified a cell cycle–regulation gene as the cause of male inviability in hybrids resulting from a cross between Drosophila melanogaster and D. simulans. Ablation of the D. simulans allele of this gene is sufficient to rescue the adult viability of hybrid males. This dominantly acting cell cycle regulator causes mitotic arrest and, thereby, inviability of male hybrid larvae. Our genomic method provides a facile means to accelerate the identification of hybrid incompatibility genes in other model and nonmodel systems.

Genetic crosses between Drosophila melanogaster females and males from its closest sister species, D. simulans, produce only adult hybrid F1 females (1, 2). These unisexual broods are a result of hybrid F1 male inviability between these species, which manifests during larval stages of development. Despite decades of investigation, the genetic basis of this hybrid F1 male inviability remains incompletely resolved (3, 4). A series of x-ray mutagenesis experiments previously revealed that a complex interaction between the D. melanogaster X-chromosome and dominant alleles from the D. simulans second and third chromosomes is necessary to kill hybrids (5, 6). The isolation of hybrid rescue strains that produce viable hybrid F1 males led to the identification of two causal elements of this hybrid incompatibility: Hybrid male rescue (Hmr) on the D. melanogaster X-chromosome (7, 8) and Lethal hybrid rescue (Lhr) on the D. simulans second chromosome (9, 10). The absence of either Lhrsim or Hmrmel results in viable hybrid males (fig. S1). However, D. melanogaster males that carry transgenic copies of D. simulans Lhr are viable despite carrying both the Hmrmel and Lhrsim incompatible alleles (9). These results suggest that the presence of at least one additional unidentified hybrid incompatibility gene is necessary to cause hybrid male inviability.

Traditional genetic approaches have failed to identify this missing hybrid incompatibility gene for several reasons. First, hybrid sterility and inviability between D. melanogaster and D. simulans hinder recombination-based methods for gene identification. Second, genetic disruptions in D. melanogaster do not assist in identifying this gene because it is a dominantly acting D. simulans factor. Third, the lack of efficient balancer chromosomes in D. simulans prevents the construction and maintenance of mutation-accumulation lines that could help identify this missing incompatibility gene. Finally, all known naturally occurring hybrid rescue alleles are mutations of either Hmr or Lhr; no new rescue alleles have been identified that may correspond to a third gene. Together, these roadblocks have prevented the identification of this missing hybrid incompatibility gene.

Because no null alleles for the missing D. simulans hybrid incompatibility gene have been isolated from natural populations, we speculated that—in contrast to Hmr and Lhr—this gene might be essential for viability. We reasoned that the complex epistatic interaction underlying hybrid F1 male inviability is analogous to a multicomponent toxin; reconstitution of this toxin requires the simultaneous presence of all components. Under this model, hybrid inviability does not occur when even one of the components or hybrid incompatibility genes is missing (e.g., loss of either Lhrsim or Hmrmel rescues hybrid males). Extending this analogy, we sought other genes whose ablation results in viable hybrid males using a simple genomics-based approach (Fig. 1A).

Fig. 1 A genomics screen identifies gfzfsim as a hybrid inviability gene.

(A) We mutagenized D. simulans males (new mutations shown in blue) and crossed them with D. melanogaster females. When a D. simulans sperm carrying a mutation at a hybrid incompatibility gene fertilizes a D. melanogaster egg, a viable “rescue” hybrid F1 male is produced. Sequencing the genomes of multiple “rescue” hybrid males identify the causative hybrid incompatibility gene mutated across these rescue males (shown in red and outlined). (B) Single-fly genome sequencing of all “rescue” hybrid males allows assignment of new mutations (including large deletions in two of the males) to the D. simulans–derived component of hybrid genomes. (C) A single gene, gfzfsim, is mutated across all six “rescue” hybrid F1 males. The x axis represents the number of genes mutated across all six males, any five males, and so on. The y axis represents the number of genes mutated across these males. (D) gfzf encodes two alternative transcripts. The larger transcript encodes FLYWCH zinc finger domains along with a GST domain, whereas the shorter transcript encodes only the GST domain.

We mutagenized 55,000 D. simulans males by feeding adults with ethyl methane sulfonate (EMS) and crossed these males to D. melanogaster females. All resulting progeny inherit one mutagenized complement of the D. simulans genome and one intact complement from D. melanogaster. When D. simulans sperm carrying null mutations at any F1 hybrid incompatibility gene fertilize D. melanogaster eggs, the resulting hybrid male progeny are predicted to be viable. This strategy allows us to survey mutations in all D. simulans genes that may be involved in the F1 hybrid incompatibility, even those in essential genes; however, haploinsufficient genes (i.e., genes that require two copies for viability) would not be sampled.

We recovered 32 viable hybrid F1 males from these crosses (compared with >300,000 hybrid F1 females). Of these, 26 males were the result of a nondisjunction event that led to their inheriting a D. simulans X-chromosome (11); these males were viable, as shown previously (2). Further tests confirmed that we had recovered six independently produced rescue hybrid F1 males, each of which is viable because of mutations at a locus different from Lhrsim (11).

Because rescue hybrid F1 males isolated from these crosses are sterile, they cannot be used in genetic crosses to map the causal gene. Instead, we performed high-throughput sequencing to obtain whole-genome sequences of each of the six independently derived rescue hybrid males and of both parental strains. We then compared the D. simulans–derived component of the genomes of rescue hybrid F1 males to the unmutagenized D. simulans parental strain. This allowed us to identify all new mutations in each of the rescue males (11) (table S1 and fig. S2). As expected, most of the EMS-induced mutations were point substitutions (Fig. 1B). However, we also identified two large partially overlapping deletions, which mapped to the D. simulans–derived chromosomal arm 3R (Fig. 1B and fig. S3). Each of the six rescue males carried between 600 and 1200 new mutations, as expected, on the basis of the random mutagenesis strategy. Only one D. simulans gene, however, was disrupted across all six rescue hybrid males (Fig. 1C). This gene was Suppressor of Killer-of-prune [Su(Kpn)]–glutathione-S-transferase–containing FLYWCH zinc finger protein (gfzf) (we refer to this as gfzf) (12, 13).

The gfzf gene encodes two alternative transcripts. The longer transcript encodes a polypeptide with four FLYWCH zinc finger domains and one glutathione-S-transferase (GST) domain, whereas the shorter transcript encodes a polypeptide with only the GST domain. The D. simulans allele of gfzf (i.e., gfzfsim) incurred unique mutations (two nonsense, one frameshift, two deletions, and one missense mutation in a highly conserved residue) in each of the six rescue hybrid F1 males (Fig. 1D and table S2). Four of these mutations only disrupt the longer of the two alternate transcripts encoded by gfzf (Fig. 1D and table S2). These results suggest that the longer gfzfsim transcript is involved in hybrid incompatibility. None of the rescue hybrid F1 males we collected had mutations in the Lhr gene, which suggested that our genetic screen did not achieve saturation. We attribute this to the fact that the coding sequence of Lhrsim [1188 base pairs (bp)] is smaller and may present a smaller mutagenesis target than gfzfsim (3117 bp).

Consistent with our predictions (6), gfzfsim resides on the D. simulans third chromosome and is essential for viability (13). To circumvent the difficulty of testing the contribution of an essential gene in hybrid inviability, we depleted (knocked down) the expression of the gfzfsim longer transcript in F1 hybrids using RNA interference (RNAi) knockdown constructs (pValium20- gfzfsim) that target only gfzfsim but not gfzfmel (11) (figs. S4 and S5). We produced transgenic D. melanogaster strains that carry these constructs under the control of the inducible promoter upstream activating sequence (UAS), inserted on the D. melanogaster X-chromosome (fig. S6).

We crossed these transgenic flies to a heterozygous D. melanogaster strain carrying a CyO balancer and a ubiquitously expressing GAL4 driver P{Act5C-GAL4} on the second chromosome. This cross produces two types of daughters. The first set inherits the CyO balancer chromosome but not the Actin5C-GAL4 driver and, therefore, does not express the knockdown construct. When these D. melanogaster females are crossed to D. simulans males, the resulting hybrid F1 males are inviable, as expected, because they do not express the knockdown construct (“RNAi off,” Fig. 2A). The second set of daughters inherits one copy each of the Act5C-GAL4 driver and the RNAi construct. In crosses between these D. melanogaster females and D. simulans males, one out of four hybrid F1 sons inherit both the RNAi construct and the Act5C-GAL4 driver and therefore robustly express the RNAi construct. We found that these hybrid F1 males are viable (“RNAi on,” Fig. 2A and figs. S7 to S9). Thus, knocking down the expression of only the long transcript of gfzfsim in hybrid F1 males is sufficient to reverse hybrid inviability. These results confirm that gfzfsim is the missing hybrid incompatibility gene.

Fig. 2 Knockdown of gfzfsim rescues cell proliferation defects and restores hybrid male viability.

(A) No hybrid males are recovered in crosses where the pValium20- gfzfsim RNAi construct is not expressed (no GAL4 driver, “RNAi OFF”). In crosses between D. simulans males and D. melanogaster females carrying one copy each of pValium20-gfzfsim and a ubiquitously expressing Actin5C-GAL4 driver, one out of four possible hybrid male progeny inherit both pValium20-gfzfsim and the Actin5C-GAL4 driver (“RNAi ON”) and produce viable F1 hybrid male progeny. P values were calculated using Fisher’s exact test. (B) RNAi knockdown of gfzfsim by a T80-GAL4 driver, more specific to larval neuroblasts and imaginal discs, successfully restores the viability of F1 male hybrids. (C) EdU staining shows the diminutive larval brains and cell proliferation defects in “inviable” hybrid males compared with viable F1 hybrid female larvae. These cell proliferation defects are also partially rescued in hybrid males upon gfzfsim knockdown.

In contrast to our results with gfzfsim knockdown, we found that disrupting gfzfmel does not rescue the viability of F1 hybrid males (table S3). Thus, allelic differences between gfzfmel and gfzfsim are important for hybrid inviability as are both Hmr and Lhr, with the limitation of comparing results between gfzfsim knockdown and gfzfmel disruption. Because positive selection likely resulted in the divergent functional properties of Hmr and Lhr orthologs (9, 14), we tested whether gfzf had also been subject to positive selection. We obtained gfzf sequences from 9 D. melanogaster and 13 D. simulans strains (table S4). Using a McDonald-Kreitman test and an outgroup species D. yakuba, we found that an excess of fixed nonsynonymous changes had occurred leading up to the hybrid inviability–associated gfzfsim, especially in the FLYWCH zinc finger domains (fig. S10). In contrast, we found no evidence for positive selection along the D. melanogaster lineage.

Although our results demonstrate the role of gfzfsim in hybrid male inviability between D. melanogaster and D. simulans, previous studies have found that gfzfmel also affects dominant genetic incompatibility between strains of D. melanogaster (13). Indeed, gfzfmel plays an essential role in potentiating inviability seen in crosses between D. melanogaster females homozygous for the eye color mutation prune (pn), and D. melanogaster males carrying Killer-of-prune (Kpn) (13) [hence, the name Su(Kpn), or Suppressor of Killer of prune]. The essential, dominant role of gfzf in lethal incompatibilities within and between species suggests that there may be limited genetic paths to the evolution of dominant lethal incompatibilities.

The ability of gfzf in mediating dominant lethal incompatibilities may stem from its role in the DNA damage–induced G2-M cell cycle–checkpoint mechanism, where it can potentiate the dE2F2/RBF pathway to block cell proliferation (15, 16). In contrast, gfzf has also shown to be required for cell proliferation by transcriptionally regulating the RAS/MAPK pathway (17). Despite its essential role in both cell cycle arrest and regulation of cell proliferation, the precise molecular function of gfzf is still uncharacterized. Moreover, the biological consequence of gfzfsim activity on hybrid male viability is unknown. Nevertheless, the developmental timing and consequences of either gfzf deficiency or gfzf-mediated dominant lethality are suggestive of a common mechanism that manifests in the larval-pupal transition. Larval tissues in Drosophila mostly consist of polyploid cells, whereas the larval nervous system and imaginal discs are composed of diploid cells. During pupation, the polyploid tissues are degraded, and the diploid imaginal discs proliferate to produce the adult body form. Individuals that lack proper imaginal discs can survive and continue to grow as larvae but die during the larval-pupal transition. Notably, homozygous gfzfmel null mutants, gfzfmel-pn-Kpn males, and gfzfsim-expressing hybrid males, all lack imaginal discs and die as larvae (13, 18, 19). This phenotype of lethality during the larval-pupal transition along with an absence of imaginal discs is diagnostic of dysfunction in cell cycle–regulation mechanisms (20).

We hypothesized that the gfzfsim-associated hybrid male lethality was due to cell proliferation defects in hybrid larvae. We therefore drove the expression of our gfzfsim knockdown construct in hybrid F1 males with a T80-GAL4 driver, which is expressed ubiquitously in early embryonic stages but specifically in the nervous system and imaginal discs in late larval stages (21). We found that T80-GAL4 mediated gfzfsim knockdown robustly rescued the viability of hybrid F1 males to adulthood (Fig. 2B and fig. S11). This result suggests that the primary defect in hybrid F1 males produced in D. melanogaster × D. simulans crosses may be in diploid tissue proliferation. Indeed, previous studies on larval brains have shown both cell cycle arrest and profound mitotic defects in hybrid F1 male larvae. These larvae also display diminutive imaginal discs and reduced larval brain sizes due to cell cycle arrest (22, 23). Using 5-ethynyl-2′-deoxyuridine (EdU) to track DNA synthesis in proliferating cells, we found that cell proliferation is restored in larval brains from hybrid F1 males upon gfzfsim knockdown, which indicates a relief from cell cycle arrest (Fig. 2C and fig. S12). Thus, gfzfsim knockdown relieves both cell cycle arrest and hybrid F1 male inviability. Together, these results support the model that gfzf is a cell cycle regulator of diploid tissues in larvae. Furthermore, they implicate the arrest of cell proliferation as the cause of hybrid F1 male inviability at this late-larval stage of Drosophila development.

Although Hmr and Lhr physically interact with each other (9), there is no evidence of a direct physical interaction between gfzf and either Hmr or Lhr. Both Hmr and Lhr proteins localize to centromeres and pericentric heterochromatin, where they play a role in mitotic chromosome segregation (24) and the suppression of transposable elements (25). These findings have led to a model in which incompatibility between Lhrsim and Hmrmel and their expression levels may cause dysfunction at centromeres or pericentric heterochromatin (24). Although the molecular nature of this dysfunction is still unclear, we speculate that the direct engagement of gfzfsim arrests the proliferation of dysfunctional diploid imaginal discs, which leads to hybrid inviability. Under this scenario, ablation of Lhrsim or Hmrmel removes the primary dysfunction, whereas ablation of gfzfsim removes the cell cycle arrest. Alternatively, gfzfsim may act indirectly by contributing to the sensitization of the hybrid genetic background, making it susceptible to the defects caused by the incompatibility between Lhrsim and Hmrmel that leads to hybrid inviability. In both scenarios, removal of any one of these three genes would restore hybrid viability. Thus, the same checkpoints that normally ensure the correction of mitotic errors may be also responsible for the inviability of hybrid males in the D. melanogaster × D. simulans interspecies cross.

The discovery of hybrid rescue genes, with mutations that reverse hybrid sterility or inviability, has meaningfully advanced our understanding of the genetic mechanisms that underlie the evolution of reproductive isolation during or following speciation. The identification of gfzf, in particular, emphasizes the role of cell cycle–regulation mechanisms in the evolution of hybrid incompatibilities (22, 23) and the complex epistatic interactions that underlie dominant hybrid incompatibilities in F1 hybrids. Our genomics-based approach may also allow mapping of genes that underlie hybrid incompatibilities and other phenotypes even when they lie within chromosomal inversions, which often impedes their precise genetic identification. Although this method requires that there be a single incompatibility separating two species, recently diverged species are likely to meet this criterion (26). Our approach may help accelerate the discovery of genes and genetic mechanisms underlying hybrid dysfunction in multiple taxa, shedding light on how reproductive isolation evolves.

Supplementary Materials

Materials and Methods

Tables S1 to S4

Figs. S1 to S12

References (2732)


  1. Materials and methods are available as supplementary materials on Science Online.
  2. Acknowledgments: We thank T. Levin, M. Levine, M. Patel, Naina Phadnis, B. Ross, and S. Zanders for their comments on the manuscript and members of the Malik, Peichel, and Phadnis labs for discussions. This work was supported by a Howard Hughes Medical Institute–Life Sciences Research Foundation (HHMI-LSRF) postdoctoral fellowship (N.P.); a Mario Cappecchi–endowed assistant professorship (N.P.); an NIH Developmental Biology Training Grant 5T32 HD0741 (J.C.C.); graduate research fellowship DGE-071824 from the NSF (J.O.K.); grants from the NIH: R01 GM115914 (N.P.), HG006283 (J.S.), and R01 GM074108 (H.S.M.); and especially a grant from the Mathers Foundation (H.S.M.). J.S. and H.S.M. are investigators of HHMI. The data reported in this paper are tabulated in the supplementary materials and archived at the following databases: primary read data at the Sequence Read Archive under accession number PRJNA290665 and gfzf sequences at GenBank under accession numbers KU064780 to KU064795.
View Abstract

Stay Connected to Science

Navigate This Article