Sexual Isolation Between Drosophila Melanogaster, D. Simulans and D. Mauritiana: Sex and Species Specific Discrimination

The sexual isolation among the related species Drosophila melanogaster, D. simulans and D. mauritiana is asymmetrical. While D. mauritiana males mate well with both D. melanogaster and D. simulans females, females of D. mauritiana discriminate strongly against males of these two species. Similarly, D. simulans males mate with D. melanogaster females but the reciprocal cross is difficult. Interspecific crosses between several populations of the three species were performed to determine if (i) males and females of the same species share a common sexual isolation genetic system, and (ii) males (or females) use the same genetic system to discriminate against females (or males) of the other two species. Results indicate that although differences in male and female isolation depend on the populations tested, the isolation behaviour between a pair of species is highly correlated despite the variations. However, the rank order of the isolation level along the populations was not correlated in both sexes, which suggests that different genes act in male and female sexual isolation. Neither for males nor for females, the isolation behaviour of one species was paralleled in the other two species, which indicates that the genetic systems involved in this trait are species-pair specific. The implications of these results are discussed. This revised version was published online in August 2006 with corrections to the Cover Date.     

Hybrid lethal systems in theDrosophila melanogaster species complex

Lethal phases of the hybrids betweenDrosophila melanogaster and its sibling species,D. simulans are classified into three types: (1) embryonic lethality in hybrids carryingD. simulans cytoplasm andD. melanogaster X chromosome, (2) larval lethality in hybrids not carryingD. simulans X, and (3) temperature-sensitive pupal lethality in hybrids carryingD. simulans X. The same lethal phases are also observed when either of the two other sibling species,D. mauritiana orD. sechellia, is employed for hybridization withD. melanogaster. Here, we describe genetic analyses of each hybrid lethality, and demonstrate that these three types of lethality are independent phenomena. We then propose two models to interpret the mechanisms of each hybrid lethality. The first model is a modification of the conventional X/autosome imbalance hypothesis assuming a lethal gene and a suppressor gene are involved in the larval lethality, while the second model is for embryonic lethality assuming an interaction between a maternal-effect lethal gene and a suppressor gene.     

A New Retroposed Gene in <Emphasis Type="Italic">Drosophila</Emphasis> Heterochromatin Detected by Microarray-Based Comparative Genomic Hybridization

A genomic pattern of new gene origination is often dependent on a genomic method that can efficiently identify a statistically adequate number of recently originated genes. The heterochromatic regions have often been viewed as genomic deserts with low coding potential and thus a low flux of new genes. However, increasing reports revealed unexpected roles of heterochromatic regions in the evolution of genes and genomes. We identified recently retroposed genes that originated in heterochromatic regions in Drosophila, by developing microarray-based comparative genomic hybridization (CGH) with multiple species. This new gene family, named Ifc-2h, originated in the common ancestor of the clade of D. simulans, D. mauritiana, and D. sechellia. The sequence features and phylogenetic distribution indicated that Ifc-2h resulted from the retroposition from its parental gene, Infertile crescent (Ifc), and integrated into heterochromatic region of common ancester of the three sibling species 2 million years ago. Expression analysis revealed that Ifc-2h had developed a new expression pattern by recruiting a putative regulatory element from its target sequence. The distribution of indel variation in Ifc-2h of D. simulans and D. mauritiana revealed a significant sequence constraint, suggesting that the Ifc-2h gene may be functional. These analyses cast fresh insight into the evolution of heterochromatin and the origin of its coding regions. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Dr. Martin Kreitman]     

Contrasting patterns of geographic variation in the cosmopolitan sibling species Drosophila melanogaster and Drosophila simulans

An electrophoretic study was carried out to compare the geographic pattern of genetic variation in Drosophila simulans with that of its sibling species, Drosophila melanogaster. An identical set of 32 gene-protein loci was studied in four geographically distant populations of D. simulans and two populations of D. melanogaster, all originating from Europe and Africa. The comparison yielded the following results: (1) tropical populations of D. simulans were, in terms of the number of unique alleles, average heterozygosity per locus, and percentage of loci polymorphic, more variable than conspecific-temperate populations; (2) some loci in both species showed interpopulation differences in allele frequencies that suggest latitudinal clines; and (3) temperate-tropical genetic differentiation between populations was much less in D. simulans than in D. melanogaster. Similar differences between these two species have previously been shown for chromosomal, quantitative, physiological, and middle-repetitive DNA variation. Estimates of N _m (number of migrants per generation) from the spatial distribution of rare alleles suggest that both species have similar levels of interpopulation gene flow. These observations lead us to propose two competing hypotheses: the low level of geographic differentiation in D. simulans is due to its evolutionarily recent worldwide colonization and, alternatively, D. simulans has a narrower niche than D. melanogaster. Geographic variation data on different genetic elements (e.g., mitochondrial DNA, two-dimensional proteins, etc.) are required before these hypotheses can be adequately tested.We thank the Natural Science and Engineering Research Council of Canada for financial support (Grant A0235 to R.S.S.).     

Association Between Levels of Coding Sequence Divergence and Gene Misregulation in <Emphasis Type="Italic">Drosophila</Emphasis> Male Hybrids

Previous studies have shown widespread conservation of gene expression levels between species of the Drosophila melanogaster subgroup as well as a positive correlation between coding sequence divergence and expression level divergence between species. Meanwhile, large-scale misregulation of gene expression level has been described in interspecific sterile hybrids between D. melanogaster, D. simulans, D. mauritiana, and D. sechellia. Using data from gene expression analysis involving D. simulans, D. melanogaster, and their hybrids, we observed a significant positive correlation between protein sequence divergence and gene expression differences between hybrids and their parental species. Furthermore, we demonstrate that underexpressed misregulated genes in hybrids are evolving more rapidly at the protein sequence level than nonmisregulated genes or overexpressed misregulated genes, highlighting the possible effects of sexual and natural selection as male-biased genes and nonessential genes are the principal gene categories affected by interspecific hybrid misregulation. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users. Carlo G. Artieri and Wilfried Haerty contributed equally to this publication.     

The evolutionary genetics of thehobo transposable element in theDrosophila melanogaster complex

Hobo elements are a family of transposable elements found inDrosophila melanogaster and its three sibling species:D. simulans, D. mauritiana andD. sechellia. Studies inD. melanogaster have shown thathobo may be mobilized, and that the genetic effects of such mobilizations included the general features of hybrid dysgenesis: mutations, chromosomal rearrangements and gonadal dysgenis in F1 individuals. At the evolutionary level somehobo-hybridizing sequences have also been found in the other members of themelanogaster subgroup and in many members of the relatedmontium subgroup. Surveys of older collected strains ofD. melanogaster suggest that completehobo elements were absent prior to 50 years ago and that they have recently been introduced into this species by horizontal transfer. In this paper we review our findings and those of others, in order to precisely describe the geographical distribution and the evolutionary history ofhobo in theD. melanogaster complex. Studies of the DNA sequences reveal a different level of divergence between the groupD. melanogaster, D. simulans andD. mauritiana and the fourth speciesD. sechellia. The hypothesis of multiple transfers in the recent past into theD. melanogaster complex from a common outside source is discussed.     

Conservation and change in structural and 5′ flanking sequences of esterase 6 in siblingDrosophila species

Esterase 6 (Est-6/EST6) is the major β-carboxylesterase inD. melanogaster and its siblingsD. simulans andD. mauritiana. It is expressed in several tissues but its major site of expression is the sperm ejaculatory duct of the adult male. Although EST6 activity affects reproductive fitness, there are high levels of electrophoretic and activity polymorphism, at least withinD. melanogaster andD. simulans. Here we present the nucleotide sequences of anEst-6 allele and its flanking regions from each ofD. simulans andD. mauritiana and compare them with the publishedD. melanogaster sequences. As might be expected, replacement sites are significantly less divergent than exon silent sites in all comparisons, suggesting that selection is acting to maintain EST6 structure and function among the three species. Nevertheless, the ratio of the levels of replacement to silent site divergence is still much higher forEst-6 than for seven of ten other genes (including both isozyme-coding loci) for which comparable data have been published for these species. This is consistent with the high levels of EST6 electrophoretic polymorphism withinD. melanogaster andD. simulans and implies that selective constraints against amino acid change are relatively weak for EST6. By contrast, comparisons involving promotor sequences show that the level of divergence in the first 350bp 5′ of the gene is significantly lower than those for four of the six other loci for which comparable data have been published for these species. In particular, there are two perfectly conserved stretches (−1 to −158bp and −219 to −334bp) each over 100bp long included in this 350bp region. Thus the data suggest a relatively low level of selective constraint on the amino acid sequence of EST6 but a relatively high level of constraint on sequences affecting aspects of its expression.     

THE GENETICS OF AN ISOLATING MECHANISM BETWEEN TWO SIBLING SPECIES OF DROSOPHILA

Matings between Drosophila simulans females and males of the sibling species D. mauritiana are of abnormally short duration. These rapid matings interrupt the transfer of sperm, leading to substantial reproductive isolation in interspecific as compared to intraspecific copulations. Genetic analysis of this behavior shows that it is influenced much more by the male than the female genotype, with genes from D. simulans being dominant. In males, shortened copulation is caused by interspecific divergence at a minimum of three loci, with one gene on each of the major chromosomes. This is an underestimate of the true number of loci affecting the trait, which could be much larger. The two autosomes have the largest effect, whereas that of the X chromosome is much smaller. The genetic architecture of copulation duration and the larger effect of male than female genotype suggest that females can detect and discriminate against differences in male genitalia.     

EXPERIMENTAL EVIDENCE FOR MITOCHONDRIAL DNA INTROGRESSION BETWEEN DROSOPHILA SPECIES

Differential introgression of mitochondrial genomes has been used to explain the occurrence in some species of individuals bearing mtDNA from a related species. This situation has been observed for Drosophila mauritiana (endemic to Mauritius) where a high proportion of individuals (88%) carries an mtDNA also found in D. simulans populations from Madagascar and Réunion. Using these two species, experiments were carried out to test for differential mtDNA introgression. A single virgin female from one species (initial frequency 0.03) was introduced into a population of the other. D. simulans mtDNA can, within three generations, almost entirely displace (frequency up to 0.80) D. mauritiana mtDNA. Hybrid male sterility probably curtails to a large degree parallel introgression of nuclear genes. The progress of cytoplasmic introgression is dependent on the degree of inbreeding of the recipient D. mauritiana strains. In reciprocal experiments, introgression was much less likely: few D. mauritiana migrant females are inseminated and their mtDNA frequency always remains very low. The results of these experiments support the hypothesis that a selective advantage of hybrids (probably at the nuclear level) has promoted mtDNA transfer from D. simulans Madagascar or Reunion populations into D. mauritiana through introgressive hybridization.     

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.     

Molecular Phylogeny of the <Emphasis Type="Italic">Drosophila melanogaster</Emphasis> Species Subgroup

Although molecular and phenotypic evolution have been studied extensively in Drosophila melanogaster and its close relatives, phylogenetic relationships within the D. melanogaster species subgroup remain unresolved. In particular, recent molecular studies have not converged on the branching orders of the D. yakuba–D. teissieri and D. erecta–D. orena species pairs relative to the D. melanogaster–D. simulans–D. mauritiana–D. sechellia species complex. Here, we reconstruct the phylogeny of the melanogaster species subgroup using DNA sequence data from four nuclear genes. We have employed vectorette PCR to obtain sequence data for orthologous regions of the Alcohol dehydrogenase (Adh), Alcohol dehydrogenase related (Adhr), Glucose dehydrogenase (Gld), and rosy (ry) genes (totaling 7164 bp) from six melanogaster subgroup species (D. melanogaster, D. simulans, D. teissieri, D. yakuba, D. erecta, and D. orena) and three species from subgroups outside the melanogaster species subgroup [D. eugracilis (eugracilis subgroup), D. mimetica (suzukii subgroup), and D. lutescens (takahashii subgroup)]. Relationships within the D. simulans complex are not addressed. Phylogenetic analyses employing maximum parsimony, neighbor-joining, and maximum likelihood methods strongly support a D. yakuba–D. teissieri and D. erecta–D. orena clade within the melanogaster species subgroup. D. eugracilis is grouped closer to the melanogaster subgroup than a D. mimetica–D. lutescens clade. This tree topology is supported by reconstructions employing simple (single parameter) and more complex (nonreversible) substitution models. Present address (Ryan M. David): University of Texas Health Science Center at San Antonio, 7703 Floyd Curl Drive, San Antonio, TX 78284, USA     

EVOLUTIONARY GENETICS OF TWO SIBLING SPECIES,DROSOPHILA SIMULANS AND D. SECHELLIA

Drosophila simulans and D. sechellia are sibling species, the former cosmopolitan and the latter restricted to the Seychelles Islands. We used classical genetic analysis to measure the numbers and effects of genes responsible for reproductive isolation and morphological differences in male genitalia between these species. At least five loci are responsible for male sterility in hybrids, with the strongest effects produced by at least two genes on the X chromosome. At least three (and probably four) loci are responsible for the interspecific difference in the size of the posterior process of the male genital arch. These genetic results, as well as the pattern of morphological divergence between the species, show several parallels with the divergence between D. simulans and its other island relative, D. mauritiana. We also present the DNA sequence of a 4.5 kilobase region containing the alcohol dehydrogenase (Adh) locus of D. sechellia, and combine this with previous data to reconstruct the phylogenies of the three species and their more distant relative D. melanogaster. Both D. mauritiana and D. sechellia are very closely related to D. simulans. Although most phylogenies show the two island species to be independent offshoots of the D. simulans lineage (with D. sechellia the more recent), the branch points are too close to make this conclusion unambiguous. The genetic and evolutionary parallels between the simulans/mauritiana and the simulans/sechellia divergences may therefore represent either a striking evolutionary convergence or a close common ancestry of the island species. A comparison of Adh alleles within species shows that the divergence among them may be almost as large as among alleles from different species. We conclude that many of the nucleotide differences among these species actually represent polymorphisms within common ancestors. It may be difficult to build accurate phylogenies using only a single DNA sequence from each species.     

Clonal analysis in hybrids between Drosophila melanogaster and Drosophila simulans

We have analysed the viability of cellular clones induced by mitotic recombination in Drosophila melanogaster/D. simulans hybrid females during larval growth. These clones contain a portion of either melanogaster or simulans genomes in homozygosity. Analysis has been carried out for the X and the second chromosomes, as well as for the 3L chromosome arm. Clones were not found in certain structures, and in others they appeared in a very low frequency. Only in abdominal tergites was a significant number of clones observed, although their frequency was lower than in melanogaster abdomens. The bigger the portion of the genome that is homozygous, the less viable is the recombinant melano-gaster/simulans hybrid clone. The few clones that appeared may represent cases in which mitotic recombination took place in distal chromosome intervals, so that the clones contained a small portion of either melanogaster or simulans chromosomes in homozygosity. Moreover, Lhr, a gene of D. simulans that suppresses the lethality of male and female melanogaster/simulans hybrids, does not suppress the lethality of the recombinant melanogaster/simulans clones. Thus, it appears that there is not just a single gene, but at least one per tested chromosome arm (and maybe more) that cause hybrid lethality. Therefore, the two species, D. melanogaster and D. simulans, have diverged to such a degree that the absence of part of the genome of one species cannot be substituted by the corresponding part of the genome of the other, probably due to the formation of co-adapted gene complexes in both species following their divergent evolution after speciation. The disruption of those coadapted gene complexes would cause the lethality of the recombinant hybrid clones.     

Insertion sites of the transposable elementmariner are fixed in the genome ofDrosophila sechellia

Summary The abundance of the transposable elementmariner differs dramatically in the genomes of the closely related speciesDrosophila simulans, D. mauritiana, D. sechellia, andD. melanogaster. Natural populations ofD. simulans andD. mauritiana have 1–10 and 20–30 copies per diploid genome, respectively, and the insertion sites are polymorphic. The element has not been found inD. melanogaster. In this paper we show thatD. sechellia, a species endemic to the Seychelles Islands, contains only twomariner elements that are at fixed sites in the genome. One element, inserted in chromosome 2R at 51A1–2, contains three deletions and is missing much of the 3 end. The other element, inserted in chromosome 3L at 64A10–11, is the full length of 1286 bp. Although the sequence of the full-length element is polymorphic in populations ofD. sechellia, at least some of the sequences are closely related to elements fromD. simulans andD. mauritiana that are known to be active. However, judging from the progeny of crosses betweenD. sechellia andD. simulans, the biological activity of the full-lengthD. sechellia element appears to be low, either because of the nucleotide sequence of the element or because of its position in the genome, or both.     

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 F₁ 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 F₁ 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 F₁ 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.     

A New Test for Selection Applied to Codon Usage in Drosophila simulans and D. mauritiana

In many organisms, synonymous codon usage is biased by a history of natural selection. However, codon bias, itself, does not indicate that selection is ongoing; it may be a vestige of past selection. Simple statistical tests have been devised to infer ongoing selection on codon usage by comparing the derived state frequency spectra at polymorphic sites segregating either derived preferred codons or derived unpreferred codons; if selection is effective, the frequency of derived states should be higher in the former. We propose a new test that uses the inferred degree of preference, essentially calculating the correlation of derived state frequency and the difference in preference between the derived and the ancestral states; the correlation should be positive if selection is effective. When implementing the test, derived and ancestral states can be assigned by parsimony or on the basis of relative probability. In either case, statistical significance is estimated by a simple permutation test. We explored the statistical power of the test by sampling polymorphism data from 14 loci in 16 strains of D. simulans, finding that the test retains 80% power even when quite a few of the data are discarded. The power of the test likely reflects better use of multiple features of the data, combining population frequencies of polymorphic variants and quantitative estimates of codon preferences. We also applied this novel test to 14 newly sequenced loci in five strains of D. mauritiana, showing for the first time ongoing selection on codon usage in this species.     

Interspecific transfer of P elements by crosses between Drosophila simulans and Drosophila mauritiana

Interspecific crosses were carried out between P element-transformed strains of D. simulans and a strain of D. mauritiana, a species devoid of this transposable element family. Four lines were established from hybrid females backcrossed with D. mauritiana males for four generations, and then maintained by intra-line mass mating. In situ hybridization of polytene chromosomes and southern blots showed that full-length and deleted P elements were present in all of the lines after 15 generations. We conclude that at least some of the P elements observed in two lines result from their transposition into D. mauritiana genome. Gonadal sterility, induced at 29°C in D. melanogaster by P elements also occurred with these two latter lines.     

RACE FORMATION,SPECIATION, AND INTROGRESSION WITHIN DROSOPHILA SIMULANS,D. MAURITIANA,AND D. SECHELLIA INFERRED FROM MITOCHONDRIAL DNA ANALYSIS

Mitochondrial DNA cleavage maps from three chromosomally homosequential species Drosophila simulans, D. mauritiana, and D. sechellia, were established for 12 restriction enzymes. One isofemale strain was studied in D. sechellia (se), 13 in D. simulans, and 17 in D. mauritiana: in the last two species, respectively, three (siI, II, and III) and two (maI and II) cleavage morphs were found. The evolutionary relationships based on mtDNA cleavage map comparisons show that the maI and se mtDNAs are internal branches of the phylogenetic tree of the D. simulans mtDNA. D. mauritiana and D. sechellia species appear to be derived from a population of D. simulans which carried an ancestral form of the current siI mtDNA type. In addition, two cleavage morphs (siIII [only present in D. simulans from Madagascar] and maI) appeared to be identical, although found in different species. We present a speculative interpretation of data on biogeography and hybridization which is consistent with the hypothesis of a recent introgression of mitochondrial DNA of D. simulans from Madagascar into D. mauritiana.     

Identification of One Intron Loss and Phylogenetic Evolution of Dfak Gene in the Drosophila melanogaster Species Group

Intron loss and its evolutionary significance have been noted in Drosophila. The current study provides another example of intron loss within a single-copy Dfak gene in Drosophila. By using polymerase chain reaction (PCR), we amplified about 1.3 kb fragment spanning intron 5–10, located in the position of Tyr kinase (TyK) domain of Dfak gene from Drosophila melanogaster species group, and observed size difference among the amplified DNA fragments from different species. Further sequencing analysis revealed that D. melanogaster and D. simulans deleted an about 60 bp of DNA fragment relative to other 7 Drosophila species, such as D. elegans, D. ficusphila, D. biarmipes, D. takahashii, D. jambulina, D. prostipennis and D. pseudoobscura, and the deleted fragment located precisely in the position of one intron. The data suggested that intron loss might have occurred in the Dfak gene evolutionary process of D. melanogaster and D. simulans of Drosophila melanogaster species group. In addition, the constructed phylogenetic tree based on the Dfak TyK domains clearly revealed the evolutionary relationships between subgroups of Drosophila melanogaster species group, and the intron loss identified from D. melanogaster and D. simulans provides a unique diagnostic tool for taxonomic classification of the melanogaster subgroup from other group of genus Drosophila.