Simulans Clade Research Articles

A new hybrid incompatibility locus between Drosophila melanogaster and Drosophila sechellia.

Despite the fundamental importance of hybrid incompatibilities to the process of speciation, there are few cases where the evolution and genetic architecture of hybrid incompatibilities are understood. One of the longest studied hybrid incompatibilities causes F1 hybrid male inviability in crosses between Drosophila melanogaster females and males from the Drosophila simulans clade of species-Drosophila simulans, Drosophila mauritiana, and Drosophila sechellia. Here, we discover dramatic differences in the manifestation of this lethal hybrid incompatibility among the D. simulans clade of species. In particular, F1 hybrid males between D. melanogaster and D. sechellia are resistant to hybrid rescue through RNAi knockdown of an essential hybrid incompatibility gene. To understand the genetic basis of this inter-species difference in hybrid rescue, we developed a triple-hybrid mapping method. Our results show that 2 discrete large effect loci and many dispersed small effect changes across the genome underlie D. sechellia aversion to hybrid rescue. The large effect loci encompass a known incompatibility gene Lethal hybrid rescue (Lhr) and previously unknown factor, Sechellia aversion to hybrid rescue (Satyr). These results show that the genetic architecture of F1 hybrid male inviability is overlapping but not identical in the 3 inter-species crosses. Our results raise questions about whether new hybrid incompatibility genes can integrate into an existing hybrid incompatibility thus increasing in complexity over time, or if the continued evolution of genes can gradually strengthen an existing hybrid incompatibility.

Rapid evolutionary diversification of the flamenco locus across simulans clade Drosophila species.

Suppression of transposable elements (TEs) is paramount to maintain genomic integrity and organismal fitness. In D. melanogaster, the flamenco locus is a master suppressor of TEs, preventing the mobilization of certain endogenous retrovirus-like TEs from somatic ovarian support cells to the germline. It is transcribed by Pol II as a long (100s of kb), single-stranded, primary transcript, and metabolized into ~24-32 nt Piwi-interacting RNAs (piRNAs) that target active TEs via antisense complementarity. flamenco is 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 that flamenco has undergone radical transformations in sequence content and even copy number across simulans clade Drosophilid species. Drosophila simulans flamenco has duplicated and diverged, and neither copy exhibits synteny with D. melanogaster beyond the core promoter. Moreover, flamenco organization is highly variable across D. simulans individuals. Next, we find that D. simulans and D. mauritiana flamenco display 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 of flamenco orthologs is testament to the extremely dynamic consequences of TE arms races on genome organization, not only amongst highly related species, but even amongst individuals.

Hybrid Sterility, Genetic Conflict and Complex Speciation: Lessons From the Drosophila simulans Clade Species.

The three fruitfly species of the Drosophila simulans clade— D. simulans, D. mauritiana, and D. sechellia— have served as important models in speciation genetics for over 40 years. These species are reproductively isolated by geography, ecology, sexual signals, postmating-prezygotic interactions, and postzygotic genetic incompatibilities. All pairwise crosses between these species conform to Haldane’s rule, producing fertile F1 hybrid females and sterile F1 hybrid males. The close phylogenetic proximity of the D. simulans clade species to the model organism, D. melanogaster, has empowered genetic analyses of their species differences, including reproductive incompatibilities. But perhaps no phenotype has been subject to more continuous and intensive genetic scrutiny than hybrid male sterility. Here we review the history, progress, and current state of our understanding of hybrid male sterility among the D. simulans clade species. Our aim is to integrate the available information from experimental and population genetics analyses bearing on the causes and consequences of hybrid male sterility. We highlight numerous conclusions that have emerged as well as issues that remain unresolved. We focus on the special role of sex chromosomes, the fine-scale genetic architecture of hybrid male sterility, and the history of gene flow between species. The biggest surprises to emerge from this work are that (i) genetic conflicts may be an important general force in the evolution of hybrid incompatibility, (ii) hybrid male sterility is polygenic with contributions of complex epistasis, and (iii) speciation, even among these geographically allopatric taxa, has involved the interplay of gene flow, negative selection, and positive selection. These three conclusions are marked departures from the classical views of speciation that emerged from the modern evolutionary synthesis.

High rate of translocation-based gene birth on the Drosophila Y chromosome

The Y chromosome is a unique genetic environment defined by a lack of recombination and male-limited inheritance. The Drosophila Y chromosome has been gradually acquiring genes from the rest of the genome, with only seven Y-linked genes being gained over the past 63 million years (0.12 gene gains per million years). Using a next-generation sequencing (NGS)-powered genomic scan, we show that gene transfers to the Y chromosome are much more common than previously suspected: at least 25 have arisen across three Drosophila species over the past 5.4 million years (1.67 per million years for each lineage). The gene transfer rate is significantly lower in Drosophila melanogaster than in the Drosophila simulans clade, primarily due to Y-linked retrotranspositions being significantly more common in the latter. Despite all Y-linked gene transfers being evolutionarily recent (<1 million years old), only three showed evidence for purifying selection (ω ≤ 0.14). Thus, although the resulting Y-linked functional gene acquisition rate (0.25 new genes per million years) is double the longer-term estimate, the fate of most new Y-linked genes is defined by rapid degeneration and pseudogenization. Our results show that Y-linked gene traffic, and the molecular mechanisms governing these transfers, can diverge rapidly between species, revealing the Drosophila Y chromosome to be more dynamic than previously appreciated. Our analytical method provides a powerful means to identify Y-linked gene transfers and will help illuminate the evolutionary dynamics of the Y chromosome in Drosophila and other species.

Comparative Analysis of Satellite DNA in the Drosophila melanogaster Species Complex.

Satellite DNAs are highly repetitive sequences that account for the majority of constitutive heterochromatin in many eukaryotic genomes. It is widely recognized that sequences and locations of satellite DNAs are highly divergent even in closely related species, contributing to the hypothesis that satellite DNA differences may underlie speciation. However, due to its repetitive nature, the mapping of satellite DNAs has been mostly left out of recent genomics analyses, hampering the use of molecular genetics techniques to better understand their role in speciation and evolution. Satellite DNAs are most extensively and comprehensively mapped in Drosophila melanogaster, a species that is also an excellent model system with which to study speciation. Yet the lack of comprehensive knowledge regarding satellite DNA identity and location in its sibling species (D. simulans, D. mauritiana, and D. sechellia) has prevented the full utilization of D. melanogaster in studying speciation. To overcome this problem, we initiated the mapping of satellite DNAs on the genomes of the D. melanogaster species complex (D. melanogaster, D. simulans, D. mauritiana, and D. sechellia) using multi-color fluorescent in situ hybridization (FISH) probes. Our study confirms a striking divergence of satellite DNAs in the D. melanogaster species complex, even among the closely related species of the D. simulans clade (D. simulans, D. mauritiana, and D. sechellia), and suggests the presence of unidentified satellite sequences in these species.

Lineage-Specific Evolution of the Complex Nup160 Hybrid Incompatibility Between Drosophila melanogaster and Its Sister Species.

Two genes encoding protein components of the nuclear pore complex Nup160 and Nup96 cause lethality in F2-like hybrid genotypes between Drosophila simulans and Drosophila melanogaster. In particular, D. simulans Nup160 and Nup96 each cause inviability when hemizygous or homozygous in species hybrids that are also hemizygous (or homozygous) for the D. melanogaster X chromosome. The hybrid lethality of Nup160, however, is genetically complex, depending on one or more unknown additional factors in the autosomal background. Here we study the genetics and evolution of Nup160-mediated hybrid lethality in three ways. First, we test for variability in Nup160-mediated hybrid lethality within and among the three species of the D. simulans clade- D. simulans, D. sechellia, and D. mauritiana. We show that the hybrid lethality of Nup160 is fixed in D. simulans and D. sechellia but absent in D. mauritiana. Second, we explore how the hybrid lethality of Nup160 depends on other loci in the autosomal background. We find that D. simulans Nup160-mediated hybrid lethality does not depend on the presence of D. melanogaster Nup96, and we find that D. simulans and D. mauritiana are functionally differentiated at Nup160 as well as at other autosomal factor(s). Finally, we use population genetics data to show that Nup160 has experienced histories of recurrent positive selection both before and after the split of the three D. simulans clade species ∼240,000 years ago. Our genetic results suggest that a hybrid lethal Nup160 allele evolved before the split of the three D. simulans clade species, whereas the other autosomal factor(s) evolved more recently.

MORE ACCURATE PHYLOGENIES INFERRED FROM LOW-RECOMBINATION REGIONS IN THE PRESENCE OF INCOMPLETE LINEAGE SORTING

When speciation events occur in rapid succession, incomplete lineage sorting (ILS) can cause disagreement among individual gene trees. The probability that ILS affects a given locus is directly related to its effective population size (Ne ), which in turn is proportional to the recombination rate if there is strong selection across the genome. Based on these expectations, we hypothesized that low-recombination regions of the genome, as well as sex chromosomes and nonrecombining chromosomes, should exhibit lower levels of ILS. We tested this hypothesis in phylogenomic datasets from primates, the Drosophila melanogaster clade, and the Drosophila simulans clade. In all three cases, regions of the genome with low or no recombination showed significantly stronger support for the putative species tree, although results from the X chromosome differed among clades. Our results suggest that recurrent selection is acting in these low-recombination regions, such that current levels of diversity also reflect past decreases in the effective population size at these same loci. The results also demonstrate how considering the genomic context of a gene tree can assist in more accurate determination of the true species phylogeny, especially in cases where a whole-genome phylogeny appears to be an unresolvable polytomy.

Genome sequencing reveals complex speciation in the Drosophila simulans clade

The three species of the Drosophila simulans clade—the cosmopolitan species, D. simulans, and the two island endemic species, D. mauritiana and D. sechellia—are important models in speciation genetics, but some details of their phylogenetic and speciation history remain unresolved. The order and timing of speciation are disputed, and the existence, magnitude, and timing of gene flow among the three species remain unclear. Here we report on the analysis of a whole-genome four-species sequence alignment that includes all three D. simulans clade species as well as the D. melanogaster reference sequence. The alignment comprises novel, paired short-read sequence data from a single highly inbred line each from D. simulans, D. mauritiana, and D. sechellia. We are unable to reject a species phylogeny with a basal polytomy; the estimated age of the polytomy is 242,000 yr before the present. However, we also find that up to 4.6% of autosomal and 2.2% of X-linked regions have evolutionary histories consistent with recent gene flow between the mainland species (D. simulans) and the two island endemic species (D. mauritiana and D. sechellia). Our findings thus show that gene flow has occurred throughout the genomes of the D. simulans clade species despite considerable geographic, ecological, and intrinsic reproductive isolation. Last, our analysis of lineage-specific changes confirms that the D. sechellia genome has experienced a significant excess of slightly deleterious changes and a dearth of presumed favorable changes. The relatively reduced efficacy of natural selection in D. sechellia is consistent with its derived, persistently reduced historical effective population size.

Rapid Functional Divergence of a Newly Evolved Polyubiquitin Gene in Drosophila and Its Role in the Trade-off between Male Fecundity and Lifespan

The cost of reproduction is a pivotal trade-off with various biological processes during the evolution of organisms. However, the genes and molecular mechanisms underlying the evolution of balancing reproductive capacity and its cost are still largely unknown. Here, we present a comprehensive study on the evolution, expression, and biological functions of a newly evolved pair of X-linked polyubiquitin tandemly duplicated genes, CG32744 and CG11700, of which the duplication event occurred in Drosophila melanogaster lineage after the split from D. simulans clade. We found that CG32744 retains conserved polyubiquitin-coding sequences across Drosophila species and is ubiquitously expressed, whereas CG11700 has accumulated numerous amino acid changes and shows a male-specific expression pattern. Null mutants of CG11700 have a higher male fecundity but shorter lifespan, whereas its overexpression decreases male fecundity. In contrast, the null mutants of the peptide-conserved CG32744 do not exhibit such phenotypes. These results suggest that CG11700 might have experienced neofunctionalization and evolved important functions in the trade-off between male fecundity and lifespan and that CG32744 likely has retained the ancestral function.

Adaptive Evolution of a Novel Drosophila Lectin Induced by Parasitic Wasp Attack

Drosophila melanogaster has long been used as a model for the molecular genetics of innate immunity. Such work has uncovered several immune receptors that recognize bacterial and fungal pathogens by binding unique components of their cell walls and membranes. Drosophila also act as hosts to metazoan pathogens such as parasitic wasps, which can infect a majority of individuals in natural populations, but many aspects of their immune responses against these more closely related pathogens are poorly understood. Here, we present data describing the transcriptional induction and molecular evolution of a candidate Drosophila anti-wasp immunity gene, lectin-24A. Lectin-24A has a secretion signal sequence and its lectin domain suggests a function in sugar group binding. Transcript levels of lectin-24A were induced significantly stronger and faster following wasp attack than following wounding or bacterial infection, demonstrating lectin-24A is not a general stress response or defense response gene but is instead part of a specific response against wasps. The major site of lectin-24A transcript production is the fat body, the main humoral immune tissue of flies. Interestingly, lectin-24A is a new gene of the D. melanogaster/Drosophila simulans clade, displaying very little homology to any other Drosophila lectins. Population genetic analyses of lectin-24A DNA sequence data from African and North American populations of D. melanogaster and D. simulans revealed gene length polymorphisms segregating at high frequencies as well as strong evidence of repeated and recent selective sweeps. Thus, lectin-24A is a rapidly evolving new gene that has seemingly developed functional importance for fly resistance against infection by parasitic wasps.

Comparative morphology and phylogeny of Nicotiana section Suaveolentes (Solanaceae) in Australia and the South Pacific

Molecular phylogenies alone have failed to resolve evolutionary relationships of Nicotiana L. section Suaveolentes Goodsp. (Solanaceae), a section largely comprising Australian endemics. Comparative morphology of Suaveolentes is illustrated and characters, together with the chromosome number, coded for phylogenetic analysis. Morphological characters included discrete characters of seeds and trichomes studied by scanning electron and light microscopy, and gap-coded quantitative measurements of flowers and vegetative organs. These data were analysed using both maximum parsimony and splits network, and compared to and combined with molecular analyses based on published nuclear and chloroplast-DNA sequences. Among morphological characters, there was a high level of homoplasy, possibly attributable to convergent evolution, hybridisation and introgression, and the underlying polyploid origin of the group. There was some conflict between the morphological and molecular datasets; however, overall there was a level of concordance that identified a phylogenetic sequence of taxa that reflected a reduction in the chromosome number. With N. africana as the functional outgroup, analyses showed N. fragrans (New Caledonia and Tongatapu) and N. fatuhivensis (Marquesas islands, formerly N. fragrans var. fatuhivensis) to be a basal lineage in the Australian and South Pacific clade of Suaveolentes. N. forsteri (eastern mainland Australia, Lord Howe Island and New Caledonia) is the sister taxon to all other Australian species of Suaveolentes, which form a well supported monophyletic group. Within this Australian endemic group, species with a chromosome count of n = 24 or 23 are early lineages. Two clades with a reduced chromosome number are the ‘N. simulans clade’ (all taxa with n = 20) and the ‘N. suaveolens clade’ (n = 19, 18, 16, or 15). The position of N. cavicola (n = 20 or 23) is equivocal, and it and several other species require further study. Some widespread taxa require greater sampling of populations to test for variation in morphology, DNA sequences and chromosome number for further elucidation of the evolutionary history of Suaveolentes.

Genomic distribution of retrotransposons 297, 1731, copia, mdg1 and roo in the Drosophila melanogaster species subgroup

The intragenomic distribution of five retrotransposon families (297, 1731, copia, mdg1 and roo) in the species of the melanogaster complex was studied by comparing results of the Southern blotting technique in males and females with those of in situ hybridization. The degree of structural polymorphism of each family in the different species was also investigated by restriction enzyme analysis. It was found that genomic distribution is a trait that depends on the family and species. The distribution of roo is mainly euchromatic in the four species and 1731 is heterochromatic, but the distribution of families 297, copia and mdg1 is markedly different in the melanogaster and simulans clades. These families were mainly euchromatic in D. melanogaster but heterochromatic in its sibling species. In the simulans clade most copia and mdg1 elements are located on chromosome Y. Differences in genomic distribution are unrelated with structural conservation. The relation of intragenomic distribution to phylogeny, transpositional activity and the role of the host genome are discussed.

Gene Expression Disruptions of Organism versus Organ in Drosophila Species Hybrids

Hybrid dysfunctions, such as sterility, may result in part from disruptions in the regulation of gene expression. Studies of hybrids within the Drosophila simulans clade have reported genes expressed above or below the expression observed in their parent species, and such misexpression is associated with male sterility in multigenerational backcross hybrids. However, these studies often examined whole bodies rather than testes or had limited replication using less-sensitive but global techniques. Here, we use a new RNA isolation technique to re-examine hybrid gene expression disruptions in both testes and whole bodies from single Drosophila males by real-time quantitative RT-PCR. We find two early-spermatogenesis transcripts are underexpressed in hybrid whole-bodies but not in assays of testes alone, while two late-spermatogenesis transcripts seem to be underexpressed in both whole-bodies and testes alone. Although the number of transcripts surveyed is limited, these results provide some support for a previous hypothesis that the spermatogenesis pathway in these sterile hybrids may be disrupted sometime after the expression of the early meiotic arrest genes.

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.

Strong Purifying Selection on the Odysseus Gene in Two Clades of Sibling Species of the Drosophila montium Species Subgroup

The Odysseus (OdsH) gene was duplicated from its ancestral neuron-expressed gene, unc-4, and then evolved very rapidly under strong positive Darwinian selection as a speciation gene causing hybrid-male sterility between closely related species of the Drosophila simulans clade. Has OdsH also experienced similar positive selection between Drosophila sibling species other than those of the simulans clade? We cloned and sequenced OdsH and unc-4 from two clades of the Drosophila montium species subgroup, the Drosophila lini and the Drosophila kikkawai clades. The ratios of Ka/Ks for OdsH were remarkably low between sibling species of these two clades, suggesting that OdsH has been subjected to strong purifying selection in these two clades.

Drosophila melanogaster and D. simulans rescue strains produce fit offspring, despite divergent centromere-specific histone alleles.

The interaction between rapidly evolving centromere sequences and conserved kinetochore machinery appears to be mediated by centromere-binding proteins. A recent theory proposes that the independent evolution of centromere-binding proteins in isolated populations may be a universal cause of speciation among eukaryotes. In Drosophila the centromere-specific histone, Cid (centromere identifier), shows extensive sequence divergence between D. melanogaster and the D. simulans clade, indicating that centromere machinery incompatibilities may indeed be involved in reproductive isolation and speciation. However, it is presently unclear whether the adaptive evolution of Cid was a cause of the divergence between these species, or merely a product of postspeciation adaptation in the separate lineages. Furthermore, the extent to which divergent centromere identifier proteins provide a barrier to reproduction remains unknown. Interestingly, a small number of rescue lines from both D. melanogaster and D. simulans can restore hybrid fitness. Through comparisons of cid sequence between nonrescue and rescue strains, we show that cid is not involved in restoring hybrid viability or female fertility. Further, we demonstrate that divergent cid alleles are not sufficient to cause inviability or female sterility in hybrid crosses. Our data do not dispute the rapid divergence of cid or the coevolution of centromeric components in Drosophila; however, they do suggest that cid underwent adaptive evolution after D. melanogaster and D. simulans diverged and, consequently, is not a speciation gene.

CYTOLOGICAL CHARACTERIZATION OF PREMEIOTIC VERSUS POSTMEIOTIC DEFECTS PRODUCING HYBRID MALE STERILITY AMONG SIBLING SPECIES OF THE DROSOPHILA MELANOGASTER COMPLEX.

In accordance with Haldane's rule, hybridizations between species of the Drosophila simulans clade produce fertile females but sterile males. In this study, a comprehensive characterization was undertaken on the six types of F1 males that were the result of the crosses between D. simulans, D. sechellia, and D. mauritiana. With the use of light and electron microscopy, it was shown that while each particular hybrid genotype exhibited a specific sterility phenotype, these phenotypes fell into two distinct classes. The two hybrid genotypes that possessed D. mauritiana X-chromosomes contained spermatogenic defects that caused arrests in premeiotic spermatogenic stages. The other four F1 hybrids possessed postmeiotic spermatogenic defects. Nonsynchronous cell divisions, underdeveloped mitochondrial derivative-axonemal associations, and microtubule abnormalities were common to all of these hybrids. Each particular postmeiotically defective hybrid genotype demonstrated characteristically distinct profiles in sperm bundle number in addition to characteristic spermiogenic arrests in the furthest developed spermatids. These results in species hybrids contrast with the absence of significant differences in spermatogenic characters between species of this clade. In addition, by utilizing an attached-X cross, we investigated the influence of maternal effects and cytoplasmic factors on the sterility of D. simulans F1 hybrids and found none. However, we discovered a strain of D. simulans (2119) that caused a large shift in sterility from postmeiotic to premeiotic when crossed to D. sechellia. This suggests that D. simulans is polymorphic for genes involving premeiotic and postmeiotic sterility and that the two types of sterilities between species may have a simple genetic basis.

Discord between the phylogenies inferred from molecular versus functional data: uneven rates of functional evolution or low levels of gene flow?

According to measures of molecular divergence, the three species of the Drosophila simulans clade are closely related to and essentially equidistant from each other. We introgressed 10% of the D. sechellia X chromosome into a pure D. simulans genetic background and found that males carrying this introgressed region were consistently fertile; in contrast, males carrying the same segment from D. mauritiana are sterile and suffer from incompatibilities at a minimum of four loci. Together with other recent results, these data suggest that D. simulans and D. sechellia are much more closely related to each other than either is to D. mauritiana. How can we reconcile the phylogeny inferred from the density of hybrid sterility genes with that inferred from molecular divergence? If the molecular phylogeny is correct, the discrepancy might be explained by uneven rates of functional evolution, resulting in the uneven accumulation of substitutions with corresponding negative effects in hybrids. If the functional phylogeny is correct, then low levels of gene flow across nascent species boundaries, particularly for loci not tightly linked to a hybrid sterility gene, may have erased the original pattern of lineage splitting. We propose tests that will allow us to discriminate between these hypotheses.

The genetics of reproductive isolation in the Drosophila simulans clade: X vs. autosomal effects and male vs. female effects.

A strong effect of homozygous autosomal regions on reproductive isolation was found for crosses between the species in the Drosophila simulans clade. Second chromosome regions were introgressed from D. mauritiana and D. sechellia into D. simulans and tested for their homozygous effects on hybrid male and hybrid female sterility and inviability. Most introgressions are fertile as heterozygotes, yet produce sterile male offspring when made homozygous. The density of homozygous autosomal factors contributing to hybrid male sterility is comparable to the density of X chromosome factors for this level of resolution. Female sterility was also revealed, yet the disparity between male and female levels of sterility was great, with male sterility being up to 23 times greater than female sterility. Complete hybrid inviability was also associated with some regions of the second chromosome, yet there were no strong sex differences. In conclusion, we find no evidence to support a strong X chromosome bias in the evolution of hybrid sterility or inviability but do find a very strong sex bias in the evolution of hybrid sterility. In light of these findings, we reevaluate the current models proposed to explain the genetic pattern of reproductive isolation.

Evidence for complex genic interactions between conspecific chromosomes underlying hybrid female sterility in the Drosophila simulans clade.

F1 hybrid females between the sibling species Drosophila simulans, Drosophila mauritiana and Drosophila sechellia are completely fertile. However, we have found that female sterility can be observed in F2 backcross females who are homozygous for D. simulans X chromosomes and homozygous for autosomal regions from either D. mauritiana or D. sechellia. Our results indicate that neither D. mauritiana autosome (2 or 3) can cause complete female sterility in a D. simulans background. The simultaneous presence of homozygous regions from both the second and third chromosomes of D. mauritiana, however, causes nearly complete female sterility which cannot be accounted for by their individual effects. The two autosomes of D. sechellia may show a similar pattern. From the same crosses, we also obtained evidence against a role for cytoplasmic or maternal effects in causing hybrid male sterility between these species. Taken with the results presented elsewhere, these observations suggest that epistatic interactions between conspecific genes in a hybrid background may be the prevalent mode of hybrid sterility between recently diverged species.