Genetics of sexual isolation in females of the Drosophila simulans species complex.

Genetic analysis of hybrids between Drosophila simulans and D. sechellia shows that sexual isolation in females is caused by at least two genes, one on each major autosome, while the X chromosome has no effect. These results are similar to those of a previous study of hybrids between D. simulans and another sibling species, D. mauritiana. In this latter hybridization, each arm of the second chromosome carries genes causing sexual isolation in females, implying a total divergence of at least three loci. The genetic similarity between the D. simulans/D. mauritiana and D. simulans/D. sechellia hybridizations probably results from independent evolution and not phylogenetic artifacts, because the dominance relationships and behavioural interactions differ between the two hybridizations. The lack of an X-chromosome effect on sexual isolation contrasts with genetic studies of post-zygotic reproductive isolation, which invariably show strong effects of this chromosome.     

The population genetics of the origin and divergence of the Drosophila simulans complex species

The origins and divergence of Drosophila simulans and close relatives D. mauritiana and D. sechellia were examined using the patterns of DNA sequence variation found within and between species at 14 different genes. D. sechellia consistently revealed low levels of polymorphism, and genes from D. sechellia have accumulated mutations at a rate that is approximately 50% higher than the same genes from D. simulans. At synonymous sites, D. sechellia has experienced a significant excess of unpreferred codon substitutions. Together these observations suggest that D. sechellia has had a reduced effective population size for some time, and that it is accumulating slightly deleterious mutations as a result. D. simulans and D. mauritiana are both highly polymorphic and the two species share many polymorphisms, probably since the time of common ancestry. A simple isolation speciation model, with zero gene flow following incipient species separation, was fitted to both the simulans/mauritiana divergence and the simulans/sechellia divergence. In both cases the model fit the data quite well, and the analyses revealed little evidence of gene flow between the species. The exception is one gene copy at one locus in D. sechellia, which closely resembled other D. simulans sequences. The overall picture is of two allopatric speciation events that occurred quite near one another in time.     

Population genetics and phylogenetics of DNA sequence variation at multiple loci within the Drosophila melanogaster species complex

Two regions of the genome, a 1-kbp portion of the zeste locus and a 1.1- kbp portion of the yolk protein 2 locus, were sequenced in six individuals from each of four species: Drosophila melanogaster, D. simulans, D. mauritiana, and D. sechellia. The species and strains were the same as those of a previous study of a 1.9-kbp region of the period locus. No evidence was found for recent balancing or directional selection or for the accumulation of selected differences between species. Yolk protein 2 has a high level of amino acid replacement variation and a low level of synonymous variation, while zeste has the opposite pattern. This contrast is consistent with information on gene function and patterns of codon bias. Polymorphism levels are consistent with a ranking of effective population sizes, from low to high, in the following order: D. sechellia, D. melanogaster, D.mauritiana, and D. simulans. The apparent species relationships are very similar to those suggested by the period locus study. In particular, D. simulans appears to be a large population that is still segregating variation that arose before the separation of D. mauritiana and D. sechellia. It is estimated that the separation of ancestral D. melanogaster from the other species occurred 2.5-3.4 Mya. The separations of D. sechellia and D. mauritiana from ancestral D. simulans appear to have occurred 0.58- 0.86 Mya, with D. mauritiana having diverged from ancestral D. simulans 0.1 Myr more recently than D. sechellia.      

Molecular population genetics of male accessory gland proteins in the Drosophila simulans complex

Accessory gland proteins are a major component of Drosophila seminal fluid. These proteins have a variety of functions and may be subject to sexual selection and/or antagonistic evolution between the sexes. Most population genetic data from these proteins are from D. melanogaster and D. simulans. Here, we extend the population genetic analysis of Acp genes to the other simulans complex species, D. mauritiana and D. sechellia. We sequenced population samples of seven Acp's from D. mauritiana, D. sechellia, and D. simulans. We investigated the population genetics of these genes on individual simulans complex lineages and compared Acp polymorphism and divergence to polymorphism and divergence from a set of non-Acp loci in the same species. Polymorphism and divergence data from the simulans complex revealed little evidence for adaptive protein evolution at individual loci. However, we observed a dramatically inflated index of dispersion for amino acid substitutions in the simulans complex at Acp genes, but not at non-Acp genes. This pattern of episodic bursts of protein evolution in Acp's provides the strongest evidence to date that the population genetic mechanisms driving Acp divergence are different from the mechanisms driving evolution at most Drosophila genes.     

Courtship song recognition in the Drosophila melanogaster complex: heterospecific songs make females receptive in D. melanogaster, but not in D. sechellia

Abstract. The courtship song emitted by male wing vibration has been regarded as one of the most important signals in sexual isolation in the species of the Drosophila melanogaster complex. Inter- and intraspecific crosses were observed using males whose wings were removed (mute) or females whose aristae were removed (deaf). Females of D. melanogaster, D. simulans , and D. mauritiana mated with heterospecific males in the song-present condition (cross between normal females and winged males) more often than in the no-song condition (cross between normal females and wingless males or between aristaless females and winged males) or they showed no preference between the two conditions. It is possible that in these females heterospecific courtship songs play a role as if they were conspecific. In contrast, the females of D. sechellia mated with D. melanogaster or D. simulans males in the no-song condition more often than in the song-present condition, suggesting that they reject males with heterospecific song. Female mate recognition depending on the courtship song in D. melanogaster, D. simulans , and D. mauritiana is considered to be relatively broader and that in D. sechellia narrower.     

Patterns of Microsatellite Variability in the Drosophila Melanogaster Complex

Forty-seven microsatellite loci were amplified in Drosophila melanogaster, Drosophila simulans, Drosophila mauritiana and Drosophila sechellia. The two cosmopolitan species D. melanogaster and D. simulans were found to be the most variable ones, followed by D. mauritiana and D. sechellia. A model based clustering algorithm was applied to the population samples of D. melanogaster, D. simulans and D. sechellia. No evidence for population substructure was detected within species--most likely due to insufficient power. A Markov chain Monte Carlo method developed for demographic inference based on microsatellites provided unambiguous evidence for population contraction in D. melanogaster, D. simulans and D. sechellia, despite that the D. melanogaster and D. simulans population samples were of non-African origin and represented recently expanded populations.     

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.     

A Reanalysis of Protein Polymorphism in Drosophila Melanogaster, D. Simulans, D. Sechellia and D. Mauritiana: Effects of Population Size and Selection

Comparison of synonymous and nonsynonymous variation/substitution within and between species at individual genes has become a widely used general approach to detect the effect of selection versus drift. The sibling species group comprised of two cosmopolitan (Drosophila melanogaster and Drosophila simulans) and two island (Drosophila mauritiana and Drosophila sechellia) species has become a model system for such studies. In the present study we reanalyzed the pattern of protein variation in these species, and the results were compared against the patterns of nucleotide variation obtained from the literature, mostly available for melanogaster and simulans. We have mainly focused on the contrasting patterns of variation between the cosmopolitan pair. The results can be summarized as follows: (1) As expected the island species D. mauritiana and D. sechellia showed much less variation than the cosmopolitan species D. melanogaster and D. simulans. (2) The chromosome 2 showed significantly less variation than chromosome 3 and X in all four species which may indicate effects of past selective sweeps. (3) In contrast to its overall low variation, D. mauritiana showed highest variation for X-linked loci which may indicate introgression from its sibling, D. simulans. (4) An average population of D. simulans was as heterozygous as that of D. melanogaster (14.4% v.s. 13.9%) but the difference was large and significant when considering only polymorphic loci (37.2% v.s. 26.1%). (5) The species-wise pooled populations of these two species showed similar results (all loci = 18.3% v.s. 20.0%, polymorphic loci = 47.2% v.s. 37.6%). (6) An average population of D. simulans had more low-frequency alleles than D. melanogaster, and the D. simulans alleles were found widely distributed in all populations whereas the D. melanogaster alleles were limited to local populations. As a results of this, pooled populations of D. melanogaster showed more polymorphic loci than those of D. simulans (48.0% v.s. 32.0%) but the difference was reduced when the comparison was made on the basis of an average population (29.1% v.s. 21.4%). (7) While the allele frequency distributions within populations were nonsignificant in both D. melanogaster and D. simulans, melanogaster had fewer than simulans, but more than expected from the neutral theory, low frequency alleles. (8) Diallelic loci with the second allele with a frequency less than 20% had similar frequencies in all four species but those with the second allele with a frequency higher than 20% were limited to only melanogaster the latter group of loci have clinal (latitudinal) patterns of variation indicative of balancing selection. (9) The comparison of D. simulans/D. melanogaster protein variation gave a ratio of 1.04 for all loci and 1.42 for polymorphic loci, against a ratio of approximately 2-fold difference for silent nucleotide sites. This suggests that the species ratios of protein and silent nucleotide polymorphism are too close to call for selective difference between silent and allozyme variation in D. simulans. In conclusion, the contrasting levels of allozyme polymorphism, distribution of rare alleles, number of diallelic loci and the patterns of geographic differentiation between the two species suggest the role of natural selection in D. melanogaster, and of possibly ancient population structure and recent worldwide migration in D. simulans. Population size differences alone are insufficient as an explanation for the patterns of variation between these two species.     

DNA Sequence Variation at the Period Locus within and among Species of the Drosophila Melanogaster Complex

A 1.9-kilobase region of the period locus was sequenced in six individuals of Drosophila melanogaster and from six individuals of each of three sibling species: Drosophila simulans, Drosophila sechellia and Drosophila mauritiana. Extensive genealogical analysis of 174 polymorphic sites reveals a complex history. It appears that D. simulans, as a large population still segregating very old lineages, gave rise to the island species D. mauritiana and D. sechellia. Rather than considering these speciation events as having produced ``sister' taxa, it seems more appropriate to consider D. simulans a parent species to D. sechellia and D. mauritiana. The order, in time, of these two phylogenetic events remains unclear. D. mauritiana supports a large number of polymorphisms, many of which are shared with D. simulans, and so appears to have begun and persisted as a large population. In contrast, D. sechellia has very little variation and seems to have experienced a severe population bottleneck. Alternatively, the low variation in D. sechellia could be due to recent directional selection and genetic hitchhiking at or near the per locus.     

Horizontal transfer of hobo transposable elements within the Drosophila melanogaster species complex: evidence from DNA sequencing.

The hobo family of transposable elements, one of three transposable-element families that cause hybrid dysgenesis in Drosophila melanogaster, appears to be present in all members of the D. melanogaster species complex: D. melanogaster, D. simulans, D. mauritiana, and D. sechellia. Some hobo-hybridizing sequences are also found in the other members of the melanogaster subgroup and in many members of the related montium subgroup. Surveys of older isofemale lines of D. melanogaster suggest that complete hobo elements were absent prior to 50 years ago and that hobo has recently been introduced into the species by horizontal transfer. To test the horizontal transfer hypothesis, the 2.6-kb XhoI fragments of hobo elements from D. melanogaster, D. simulans, and D. mauritiana were cloned and sequenced. The DNA sequences reveal an extremely low level of divergence and support the conclusion that the active hobo element has been horizontally transferred into or among these species in the recent past.     

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.     

High level of divergence of male-reproductive-tract proteins, between Drosophila melanogaster and its sibling species, D. simulans

We compared male-reproductive-tract polypeptides of Drosophila melanogaster and D. simulans by using two-dimensional gel electrophoresis. Approximately 64% of male-reproductive-tract polypeptides were identical between two randomly chosen isofemale lines from these two species, compared with 83% identity for third-instar imaginal wing-disc polypeptides. Qualitatively similar differences were found between reproductive tracts and imaginal discs when D. sechellia was compared with D. melanogaster and with D. simulans. When genic polymorphism was taken into account, approximately 10% of male- reproductive-tract polypeptides were apparently fixed for different alleles between D. melanogaster and D. simulans; this proportion is the same as that found for soluble enzymes by one-dimensional gel electrophoresis. Strikingly, approximately 20% of male-reproductive- tract polypeptides of either D. melanogaster or D. simulans had no detectable homologue in the other species. We propose that proteins of the Drosophila male reproductive tract may have diverged more extensively between species than have other types of proteins and that much of this divergence may involve large changes in levels of polypeptide expression.      

Rapid evolution of the sex-determining gene,transformer: structural diversity and rate heterogeneity among sibling species of Drosophila

While developmentally regulated genes are generally conserved, transformer (tra), a key locus involved in the regulation of sexual differentiation, is highly diverged between species of Drosophila. With an aim to understand its divergence between sibling species, we investigated tra sequence variation among members of the Drosophila melanogaster species complex, D. melanogaster, D. simulans, D. mauritiana, and D. sechellia. In this species group, tra divergence is rapid yet clocklike and exhibits large differences in protein size. D. melanogaster contains a 13-amino acid tandem duplication, whereas D. sechellia possesses a 72-amino acid tandem duplication representing a 30% increase in total amino acid residues. We also found evidence of a nonrandom distribution of replacement substitutions and heterogeneity in substitution rates using clustering statistics and a codon substitution model. We show that tra's rapid divergence in this species complex is the result of generally lower selective constraints around regions that encode arginine-serine (RS) domains and a significantly higher rate of substitutions around the insertion site of D. sechellia's large duplication. The proximity of rapidly diverged regions to sites of nucleotide insertion suggests that higher local rates of mutation may provide a causal mechanism for TRA's rapid divergence in this subgroup. A comparison of tra orthologs across the genus Drosophila suggest that TRA maintains an assortment of RS domains for proper sex determining function while much of the protein evolves relatively unconstrained.     

Genetics of a Pheromonal Difference Affecting Sexual Isolation between Drosophila Mauritiana and D. Sechellia

Females of the sibling species Drosophila sechellia and D. mauritiana differ in their cuticular hydrocarbons: the predominant compound in D. sechellia is 7,11-heptacosadiene (7,11-HD), while that in D. mauritiana is 7-tricosene (7-T). We investigate the genetic basis of this difference and its involvement in reproductive isolation between the species. Behavioral studies involving hydrocarbon transfer suggest that these compounds play a large role in the sexual isolation between D. mauritiana males and D. sechellia females, while sexual isolation in the reciprocal hybridization results more from differences in female behavior than hydrocarbons. This interspecific difference in hydrocarbon profile is due to evolutionary change at a minimum of six loci, all on the third chromosome. The localization of evolutionary change to the third chromosome has been seen in every other genetic analysis of female hydrocarbon differences in the D. melanogaster group. We suggest that the high 7,11-HD phenotype seen in two species evolved twice independently from ancestors having the high 7-T phenotype, and that the recurrent third-chromosome effects are evolutionary convergences that may be due to a concentration of ``hydrocarbon genes' on that chromosome.     

Wolbachia Infections and the Expression of Cytoplasmic Incompatibility in Drosophila Sechellia and D. Mauritiana

Various stocks of Drosophila mauritiana and D. sechellia were found to be infected with Wolbachia, a Rickettsia-like bacterium that is known to cause cytoplasmic incompatibility and other reproductive abnormalities in arthropods. Testing for the expression of cytoplasmic incompatibility in these two species showed partial incompatibility in D. sechellia but no expression of incompatibility in D. mauritiana. To determine whether absence of cytoplasmic incompatibility in D. mauritiana was due to either the bacterial or host genome, we transferred bacteria from D. mauritiana into an uninfected strain of D. simulans, a host species known to express high levels of incompatibility with endogenous Wolbachia. We also performed the reciprocal transfer of the natural D. simulans Riverside infection into a tetracycline-treated stock of D. mauritiana. In each case, the ability to express incompatibility was unaltered by the different host genetic background. These experiments indicate that in D. simulans and D. mauritiana expression of the cytoplasmic incompatibility phenotype is determined by the bacterial strain and that D. mauritiana harbors a neutral strain of Wolbachia.     

Evidence for Complex Genic Interactions between Conspecific Chromosomes Underlying Hybrid Female Sterility in the Drosophila Simulans Clade

F(1) 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 F(2) 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.     

Genetics of Differences in Pheromonal Hydrocarbons between Drosophila Melanogaster and D. Simulans

Females of Drosophila melanogaster and its sibling species D. simulans have very different cuticular hydrocarbons, with the former bearing predominantly 7,11-heptacosadiene and the latter 7-tricosene. This difference contributes to reproductive isolation between the species. Genetic analysis shows that this difference maps to only the third chromosome, with the other three chromosomes having no apparent effect. The D. simulans alleles on the left arm of chromosome 3 are largely recessive, allowing us to search for the relevant regions using D. melanogaster deficiencies. At least four nonoverlapping regions of this arm have large effects on the hydrocarbon profile, implying that several genes on this arm are responsible for the species difference. Because the right arm of chromosome 3 also affects the hydrocarbon profile, a minimum of five genes appear to be involved. The large effect of the thrid chromosome on hydrocarbons has also been reported in the hybridization between D. simulans and its closer relative D. sechellia, implying either an evolutionary convergence or the retention in D. sechellia of an ancestral sexual dimorphism.     

Hybrid Lethal Systems in the Drosophila Melanogaster Species Complex. II. the Zygotic Hybrid Rescue (Zhr) Gene of D. Melanogaster

Hybrid females from Drosophila simulans females X Drosophila melanogaster males die as embryos while hybrid males from the reciprocal cross die as larvae. We have recovered a mutation in melanogaster that rescues the former hybrid females. It was located on the X chromosome at a position close to the centromere, and it was a zygotically acting gene, in contrast with mhr (maternal hybrid rescue) in simulans that rescues the same hybrids maternally. We named it Zhr (Zygotic hybrid rescue). The gene also rescues hybrid females from embryonic lethals in crosses of Drosophila mauritiana females X D. melanogaster males and of Drosophila sechellia females X D. melanogaster males. Independence of the hybrid embryonic lethality and the hybrid larval lethality suggested in a companion study was confirmed by employing two rescue genes, Zhr and Hmr (Hybrid male rescue), in doubly lethal hybrids. A model is proposed to explain the genetic mechanisms of hybrid lethalities as well as the evolutionary pathways.     

The Effects of Interspecific Y Chromosome Replacements on Hybrid Sterility within the Drosophila Simulans Clade

We attempted to introgress Y chromosomes between three sibling species of Drosophila: D. simulans, D. sechellia and D. mauritiana. Four D. sechellia Y chromosomes were introgressed into D. simulans without loss of fertility whereas the four reciprocal introgressions (D. simulans Y introgressed into D. sechellia) all result in sterility. Both reciprocal Y introgressions of D. simulans and D. mauritiana (four of each) also result in sterility. Compared with D. simulans males, the males with the D. sechellia Y chromosome in D. simulans background had lower productivity but only after multiple matings with virgin females. These males also were inferior compared with pure species males in sperm displacement and/or remating ability. The two different Y genotype males, however, were comparable in viability, longevity and mating success in female choice tests. We also use our results to estimate the effective number of autosomal loci interacting with X-linked genes to produce hybrid male sterility.