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.     

Distribution and conservation of sequences homologous to the 1731 retrotransposon in Drosophila

The distribution of 1731 retrotransposon-hybridizing sequences in the family Drosophilidae has been studied using a 1731 probe from Drosophila melanogaster. Squash blot and Southern blot analyses of 42 species reveal that the 1731 sequences are widespread within both the Sophophora and Drosophila subgenera and are also present in the genera Scaptomyza and Zaprionus. Hence the 1731 retrotransposon family appears to have a long evolutionary history in the Drosophilidae genome. Differences of hybridization signal intensity suggested that the 1731 sequence is well conserved only in the three species most closely related to D. melanogaster (D. simulans, D. mauritiana, and D. sechellia). A survey of insertion sites in numerous different populations of the previous four species by in situ hybridization to polytene chromosomes has shown in all cases both chromocentric hybridizations and a low number of sites (0-5) on the chromosomal arms. This number of sites is among the lowest observed in D. melanogaster and D. simulans when 1731 is compared with other retrotransposon families. In addition, we have observed species-specific patterns of the chromocentric hybridization signal, suggesting rapid modifications of the beta-heterochromatin components since the radiation of the melanogaster subgroup.      

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.     

Chromosome localization of the lambda20p1.4 clone of the Drosophila melanogaster nuclear lamina DNA in the melanogaster species subgroup of the genus Drosophila (Sophophora)

Chromosome localization of sequences homologous to lambda 20p1.4 of the Drosophila melanogaster nuclear lamina DNA (nlDNA) was established by in situ hybridization in species of the melanogaster subgroup. DNA of the lambda 20p1.4 clone was shown to be located in the chromocenter in all the species examined. Laboratory strains of D. simulans, D. mauritiana, and D. sechellia exhibited interspecific differences in localization of lambda 20p1.4 nlDNA on chromosome arms. In eight natural populations, intraspecific polymorphism of lambda 20p1.4 nlDNA chromosome localization was shown to be present in D. simulans but absent in D. melanogaster. The possible participation of transposable elements in nlDNA relocation is discussed.     

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.     

Wolbachia Infection and Cytoplasmic Incompatibility in Drosophila Species

Forty-one stocks from 30 Drosophila species were surveyed for Wolbachia infection using PCR technology. D. sechellia and two strains of D. auraria were found to be infected and were tested for the expression of cytoplasmic incompatibility, along with D. ananassae and D. melanogaster strains, which are already known to be infected. D. ananassae and D. melanogaster show levels of incompatibility up to 25%, while D. auraria and D. sechellia exhibit levels of egg mortality ~60%. A dot-blot assay using the dnaA sequence as probe was developed to assess the infection levels in individual males that were used in incompatibility crosses. A positive correlation between bacterial density and cytoplasmic incompatibility was observed. The stocks examined can be clustered into at least two groups, depending on the levels of infection relative to the degree of cytoplasmic incompatibility exhibited. One group, containing D. simulans Hawaii, D. sechellia, and D. auraria, exhibits high levels of cytoplasmic incompatibility relative to levels of infection; all the other species and D. simulans Riverside exhibit significantly lower levels of cytoplasmic incompatibility relative to levels of infection. These data show that, in addition to bacterial density, bacterial and/or host factors also affect the expression of cytoplasmic incompatibility.     

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.     

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.     

Changes in the recombinational environment affect divergence in the yellow gene of Drosophila

The complete coding region of the yellow (y) gene was sequenced in different Drosophila species. In the species of the melanogaster subgroup (D. melanogaster, D. simulans, D. mauritiana, D. yakuba, and D. erecta), this gene is located at the tip of the X chromosome in a region with a strong reduction in recombination rate. In contrast, in D. ananassae (included in the ananassae subgroup of the melanogaster group) and in the obscura group species (D. subobscura, D. madeirensis, D. guanche, and D. pseudoobscura), the y gene is located in regions with normal recombination rates. As predicted by the hitchhiking and background selection models, this change in the recombinational environment affected synonymous divergence in the y-gene-coding region. Estimates of the number of synonymous substitutions per site were much lower between the obscura group species and D. ananassae than between the species of the obscura group and the melanogaster subgroup. In fact, a highly significant increase in the rate of synonymous substitution was detected in all lineages leading to the species of the melanogaster subgroup relative to the D. ananassae lineage. This increase can be explained by a higher fixation rate of mutations from preferred to unpreferred codons (slightly deleterious mutations). The lower codon bias detected in all species of the melanogaster subgroup relative to D. ananassae (or to the obscura group species) would be consistent with this proposal. Therefore, at least in Drosophila, changes in the recombination rate in different lineages might cause deviations of the molecular-clock hypothesis and contribute to the overdispersion of the rate of synonymous substitution. In contrast, the change in the recombinational environment of the y gene has no detectable effect on the rate of amino acid replacement in the Yellow protein.     

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.     

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.     

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.     

Comparative analysis of the localization and mobility of retrotransposons in sibling species Drosophila simulans and Drosophila melanogaster]

The distribution of four retrotransposon families (MDG1, MDG3, MDG4 and copia) on polytene chromosomes of different (from 9 to 15) Drosophila simulans strains is studied. The mean number of MDG1 and copia euchromatic hybridization sites (3 sites for each element) is drastically decreased in D. simulans in comparison with D. melanogaster (24 and 18 sites respectively). The mean number of MDG3 sites of hybridization is 5 in D. simulans against 12 in D. melanogaster. As for MDG4 both species have on the average about 2-3 euchromatic sites. The majority of MDG1 and copia and about a half of MDG3 euchromatic copies are localized in restricted number of sites (hot spots) on D. simulans polytene chromosomes. In D. melanogaster these elements are scattered along the chromosomes though there are some hot spots too. It appears that euchromatic copies of MDG1 and copia are considerably less mobile in D. simulans in contrast to D. melanogaster. Some common hot spots of retrotransposon localization in D. simulans and D. melanogaster were earlier described as intercalary heterochromatin regions in D. melanogaster. The level of interstrain variability of MDG4 hybridization sites is comparable in both species. Comparative blot-analysis of adult and larval salivary gland DNA shows that MDG1 and copia are situated mainly in euchromatic regions of D. melanogaster chromosomes. In D. simulans genome they are located mainly in heterochromatic regions underreplicated in salivary gland polytene chromosomes. There are interspecies differences in the distribution of retrotransposons in beta-heterochromatic chromosome regions.     

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.     

Molecular Evolution of the Metallothionein Gene Mtn in the Melanogaster Species Group: Results from Drosophila Ananassae

Three distinctly different alleles of the metallothionein gene Mtn have been identified in natural Drosophila melanogaster populations: Mtn(.3), Mtn(1), and Dp(Mtn(1)), where the latter designates a tandem duplication of Mtn(1). In Drosophila simulans, only Mtn(.3)-type alleles have been found. It has been suggested that Mtn(.3) is the ancestral allele and demonstrated that a presumed two-step transition from Mtn(.3) to Mtn(1) to Dp(Mtn(1)) is accompanied by an approximate 5-fold increase in RNA levels. We analyzed the evolutionary genetics of the Mtn locus of Drosophila ananassae, a distant relative of D. melanogaster and D. simulans within the melanogaster species group. The Mtn gene of D. ananassae is most similar to Mtn(.3). (i) it is identical with Mtn(.3) at the amino acid level, but differs from Mtn(1) in its terminal codon; (ii) its 3'' UTR contains a characteristic extra DNA segment of about 50 bp which is present in Mtn(.3), but lacking in Mtn(1); (iii) duplications of Mtn were not found in a worldwide sample of 110 wild D. ananassae chromosomes. However, the intron of the Mtn gene in D. ananassae is only 69 bp long, whereas the length of the Mtn(.3) and Mtn(1) introns is 265 bp; and it lacks a polypyrimidine stretch upstream of the 3'' splice site in contrast to the much greater pyrimidine-richness found in the Mtn(.3) and Mtn(1) introns. A short intron (67 bp) was also identified in a D. pseudoobscura Mtn allele, suggesting that the short intron is the ancestral form and that the transition from the short to the long intron occurred within the melanogaster species group. We discuss the significance of this observation with regard to the recently proposed classification of D. melanogaster introns into two groups: short introns (<90 bp) which tend to lack polypyrimidine stretches, and longer ones which have strong 3'' splice signals similar to mammalian introns. A database search revealed that this length dimorphism is an evolutionarily conserved feature of Drosophila introns; transitions from one size class to the other appear to be rare between closely related species (e.g., within the melanogaster subgroup).     

Ancestral polymorphism and recent invasion of transposable elements in Drosophila species

ABSTRACT: BACKGROUND: During the evolutionary history of transposable elements, some processes, such as ancestral polymorphisms and horizontal transfer of sequences between species, can produce incongruences in phylogenies. We investigated the evolutionary history of the transposable elements Bari and 412 in the sequenced genomes of the Drosophila melanogaster group and in the sibling species D. melanogaster and D. simulans using traditional phylogenetic and network approaches. RESULTS: The maximum likelihood (ML) phylogenetic analyses revealed incongruences and unresolved relationships for both the Bari and 412 elements. The DNA transposon Bari within the D. ananassae genome is more closely related to the element of the melanogaster complex than to the sequence in D. erecta, which is inconsistent with the species phylogeny. Divergence analysis and the comparison of the rate of synonymous substitutions per synonymous site of the Bari and host gene sequences explain the incongruence as an ancestral polymorphism inherited stochastically by the derived species. Unresolved relationships were observed in the ML phylogeny of both elements involving D. melanogaster, D. simulans and D. sechellia. A network approach was used to attempt to resolve these relationships. The resulting tree suggests recent transfers of both elements between D. melanogaster and D. simulans. The divergence values of the elements between these species support this conclusion. CONCLUSIONS: We showed that an ancestral polymorphism and recent invasion of genomes due to introgression or horizontal transfer between species occurred during the evolutionary history of the Bari and 412 elements in the melanogaster group. These invasions likely occurred in Africa during the Pleistocene, before the worldwide expansion of D. melanogaster and D. simulans.     

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.