Genetics of a Difference in Male Cuticular Hydrocarbons between Two Sibling Species, Drosophila Simulans and D. Sechellia

In seven of the eight species of the Drosophila melanogaster group, the predominant cuticular hydrocarbon of males is 7-tricosene, but in the island endemic species D. sechellia it is 6-tricosene. The phylogeny of the group implies that the novel hydrocarbon profile of D. sechellia is a derived character. Genetic analysis of hybrids between D. sechellia and its close relative D. simulans show that each of the five major chromosome arms carries at least one gene affecting the ratio of the two tricosene isomers, with the right arm of the third chromosome having the largest effect. The species difference in this character is therefore polygenic with the effects of the different chromosome arms generally additive, although there is some epistasis among third-chromosome genes. Observations of courtship by males who have been coated with foreign hydrocarbons suggest that a male's hydrocarbon profile may slightly affect the degree of sexual isolation in one of the reciprocal hybridizations between these species, but that this role is small compared to that played by hydrocarbon differences between females.     

Genetics of egg production in Drosophila sechellia

Drosophila sechellia, an island endemic that specializes on a single host plant, has a lower rate of egg production than its generalist sister species D. melanogaster, D. simulans, and D. mauritiana. Earlier work showed that part of this difference in egg production was due to a reduction in the number of ovarioles in D. sechellia relative to its sister species. Here, I extend this earlier work by genetically analyzing the difference in egg production between D. sechellia and D. simulans. In all, 10 genetic markers were used in several interspecific backcrosses to identify chromosome regions that affected the rate of egg production. While previously mapped factors affecting ovariole number appear to impact the rate of egg production, new, non-ovariole factors were also identified. Overall, the difference in egg production between D. sechellia and D. simulans appears to be a polygenic trait. The relationship between these factors and genes involved the adaptation of D. sechellia to its host plant is not yet clear. The data are consistent with the hypothesis that decline in egg production is, in part, a negative pleiotropic effect of genetic changes required for host specialization in D. sechellia, although finer-scale genetic analysis of both traits is needed to truly test this hypothesis.     

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.     

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.     

Quantitative trait loci for cuticular hydrocarbons associated with sexual isolation between Drosophila simulans and D. sechellia

The identification of genes with large effects on sexual isolation and speciation is an important link between classic evolutionary genetics and molecular biology. Few genes that affect sexual isolation and speciation have been identified, perhaps because many traits influencing sexual isolation are complex behaviors. Cuticular hydrocarbons (CHs) of species of the Drosophila melanogaster group play a large role in sexual isolation by functioning as contact pheromones influencing mate recognition. Some of the genes that play key roles in determining species-specific CHs have been identified. We have performed separate quantitative trait locus (QTL) analyses of 7-tricosene (7-T) and 7,11-heptacosadiene (7,11-HD), the two major female CHs differing between D. simulans and D. sechellia. We find that approximately 40% of the phenotypic variance in each CH is associated with two to four chromosomal regions. A region on the right arm of chromosome 3 contains QTL that affect both traits, but other QTL are in distinct chromosomal regions. Epistatic interactions were detected between two pairs of QTL for 7,11-HD such that if either were homozygous for the D. simulans allele, the fly was similar to D. simulans in phenotype, with a low level of 7,11-HD. We discuss the location of these regions with regard to candidate genes for CH production, including those for desaturases.     

Variation in fiber number of a male-specific muscle between Drosophila species: a genetic and developmental analysis

We characterize a newly discovered morphological difference between species of the Drosophila melanogaster subgroup. The muscle of Lawrence (MOL) contains about four to five fibers in D. melanogaster and Drosophila simulans and six to seven fibers in Drosophila mauritiana and Drosophila sechellia. The same number of nuclei per fiber is present in these species but their total number of MOL nuclei differs. This suggests that the number of muscle precursor cells has changed during evolution. Our comparison of MOL development indicates that the species difference appears during metamorphosis. We mapped the quantitative trait loci responsible for the change in muscle fiber number between D. sechellia and D. simulans to two genomic regions on chromosome 2. Our data eliminate the possibility of evolving mutations in the fruitless gene and suggest that a change in the twist might be partly responsible for this evolutionary change.     

Heterochromatic distribution of HeT-A- and TART-like sequences in several Drosophila species

Drosophila melanogaster telomeres contain arrays of two non-LTR retrotransposons called HeT-A and TART. Previous studies have shown that HeT-A- and TART-like sequences are also located at non-telomeric sites in the Y chromosome heterochromatin. By in situ hybridization experiments, we mapped TART sequences in the h16 region of the long arm close to the centromere of the Y chromosome of D. melanogaster. HeT-A sequences were localized in two different regions on the Y chromosome, one very close to the centromere in the short arm (h18-h19) and the other in the long arm (h13-h14). To assess a possible heterochromatic location of TART and HeT-A elements in other Drosophila species, we performed in situ hybridization experiments, using both TART and HeT-A probes, on mitotic and polytene chromosomes of D. simulans, D. sechellia, D. mauritiana, D. yakuba and D. teissieri. We found that TART and HeT-A probes hybridize at specific heterochromatic regions of the Y chromosome in all Drosophila species that we analyzed.     

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.     

The genetics of mating recognition between Drosophila simulans and D. sechellia

During courtship, visual and chemical signals are often exchanged between the sexes. The proper exchange of such signals ensures intraspecific recognition. We have examined the genetic basis of interspecific differences in male mating behaviour and pheromone concentration between Drosophila simulans and D. sechellia by using Drosophila simulans/D. sechellia introgression lines. Our results show a majority of quantitative trait loci (QTLs) explaining variation in both male mating behaviour and pheromone concentration to be located on the third chromosome. One QTL found on the third chromosome explains variation in time needed to start courtship and copulation as well as time spent courting. The position of such QTL (approximately 84A-88B) with effects on courtship and copulation aspects of mating includes the candidate sex determination gene doublesex (84E5-6) and Voila (86E1-2), a gene that affects male courtship in D. melanogaster. One additional third chromosome QTL explained variation in 7-tricosene pheromone concentrations among males. The interval mapping position of this QTL (approximately 68E-76E) did not overlap with the position detected for differences in mating behaviour and the intervals did not include candidate genes previously identified as having an effect on D. melanogaster cuticular hydrocarbon production. We did not detect any directionality of the effect of Drosophila sechellia allele introgressions in male mating recognition.     

The Genetic Basis of Haldane''s Rule and the Nature of Asymmetric Hybrid Male Sterility among Drosophila Simulans, Drosophila Mauritiana and Drosophila Sechellia

Haldane's rule (i.e., the preferential hybrid sterility and inviability of heterogametic sex) has been known for 70 years, but its genetic basis, which is crucial to the understanding of the process of species formation, remains unclear. In the present study, we have investigated the genetic basis of hybrid male sterility using Drosophila simulans, Drosophila mauritiana and Drosophila sechellia. An introgression of D. sechellia Y chromosome into a fairly homogenous background of D. simulans did not show any effect of the introgressed Y on male sterility. The substitution of D. simulans Y chromosome into D. sechellia, and both reciprocal Y chromosome substitutions between D. simulans and D. mauritiana were unsuccessful. Introgressions of cytoplasm between D. simulans and D. mauritiana (or D. sechellia) also did not have any effect on hybrid male sterility. These results rule out the X-Y interaction hypothesis as a general explanation of Haldane's rule in this species group and indicate an involvement of an X-autosome interaction. Models of symmetrical and asymmetrical X-autosome interaction have been developed which explain the Y chromosome substitution results and suggest that evolution of interactions between different genetic elements in the early stages of speciation is more likely to be of an asymmetrical nature. The model of asymmetrical X-autosome interaction also predicts that different sets of interacting genes may be involved in different pairs of related species and can account for the observation that hybrid male sterility in many partially isolated species is often nonreciprocal or unidirectional.     

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.     

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.      

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.     

BAC library construction and BAC end sequencing of five Drosophila species: the comparative map with the D. melanogaster genome

We constructed and characterized arrayed bacterial artificial chromosome (BAC) libraries of five Drosophila species (D. melanogaster, D. simulans, D. sechellia, D. auraria, and D. ananassae), which are genetically well characterized in the studies of meiosis, evolution, population genetics, and developmental biology. The BAC libraries comprise 8,000 to 12,500 clones for each species, estimated to cover the most of the genomes. We sequenced both ends of most of these BAC clones with a success rate of 91%. Of these, 53,701 clones consisting of non-repetitive BAC end sequences (BESs) were mapped with reference of the public D. melanogaster genome sequences. The BES mapping estimated that the BAC libraries of D. auraria and D. ananassae covered 47% and 57% of the D. melanogaster genome, respectively, and those of D. melanogaster, D. sechellia, and D. simulans covered 94-97%. The low coverage by BESs of D. auraria and D. ananassae may be due to the high sequence divergence with D. melanogaster. From the comparative BES mapping, 111 possible breakpoints of chromosomal rearrangements were identified in these four species. The breakpoints of the major chromosome rearrangement between D. simulans and D. melanogaster on the third chromosome were determined within 20 kb in 84E and 30 kb in 93E/F. Corresponding breakpoints were also identified in D. sechellia. The BAC clones described here will be an important addition to the Drosophila genomic resources.     

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.     

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.     

Comparative analysis of transposable elements in the melanogaster subgroup sequenced genomes

Transposable elements (TEs) are indwelling components of genomes, and their dynamics have been a driving force in genome evolution. Although we now have more information concerning their amounts and characteristics in various organisms, we still have little data from overall comparisons of their sequences in very closely-related species. While the Drosophila melanogaster genome has been extensively studied, we have only limited knowledge regarding the precise TE sequences in the genomes of the related species Drosophila simulans, Drosophila sechellia and Drosophila yakuba. In this study we analyzed the number and structure of TE copies in the sequenced genomes of these four species. Our findings show that, unexpectedly, the number of TE insertions in D. simulans is greater than that in D. melanogaster, but that most of the copies in D. simulans are degraded and in small fragments, as in D. sechellia and D. yakuba. This suggests that all three species were invaded by numerous TEs a long time ago, but have since regulated their activity, as the present TE copies are degraded, with very few full-length elements. In contrast, in D. melanogaster, a recent activation of TEs has resulted in a large number of almost-identical TE copies. We have detected variants of some TEs in D. simulans and D. sechellia, that are almost identical to the reference TE sequences in D. melanogaster, suggesting that D. melanogaster has recently been invaded by active TE variants from the other species. Our results indicate that the three species D. simulans, D. sechellia, and D. yakuba seem to be at a different stage of their TE life cycle when compared to D. melanogaster. Moreover, we show that D. melanogaster has been invaded by active TE variants for several TE families likely to come from D. simulans or the ancestor of D. simulans and D. sechellia. The numerous horizontal transfer events implied to explain these results could indicate introgression events between these species.