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.     

Organization and evolution of the alcohol dehydrogenase gene in Drosophila

The alcohol dehydrogenase (Adh) gene was isolated from Drosophila simulans and D. mauritiana, and the DNA sequence of a 4.6-kb region, containing the structural gene and flanking sequence, was determined for each. These sequences were compared with the Adh region of D. melanogaster to characterize changes that occur in the Drosophila genome during evolution and to identify conserved sequences of functional importance. Drosophila simulans and D. mauritiana Adh are organized in a manner similar to that of D. melanogaster Adh, including the presence of two promoters for the single Adh gene. This study identified conserved flanking elements that, in conjunction with other studies, suggest regions that may be involved in the control of Adh expression. Inter- and intraspecies comparisons revealed differences in the kinds of sequence changes that have accumulated. Sequence divergence in and around the Adh gene was used to assess inter- and intraspecies evolutionary relationships. Finally, there appears to be an unrelated structural gene located directly 3' of the Adh transcribed region.      

Nucleotide sequence divergence in the A+T-rich region of mitochondrial DNA in Drosophila simulans and Drosophila mauritiana

We have determined the nucleotide sequences of two regions within the A+T-rich region of mitochondrial DNA (mtDNA) in the siIII type of Drosophila simulans and the maI type of D. mauritiana. The sequences of the two regions in siIII and maI are almost identical. The sequences include elements corresponding to the type I and type II repeats elements and the T-stretches as reported in D. melanogaster; an approximately 340-bp region (A region) adjacent to the tRNA(Ile) gene includes a part of the type II repeat element, and an approximately 440- bp region (B region) includes a central portion of the A+T-rich region between the type I and type II repeat arrays. Each sequence of the two species was compared with those of D. melanogaster and D. yakuba. The sequences of the A region are relatively well conserved among the four species. The alignment of the two sequences of the B region with those of D. melanogaster and D. yakuba requires numerous insertions/deletions. For both regions, nucleotide differences between D. simulans or D. mauritiana and D. melanogaster are similar to those between the two and D. yakuba. The tendency is obvious in a subregion within the type II repeat element in the A region. These findings suggest that the rate of nucleotide substitution in the subregion is accelerated in the lineage leading to D. melanogaster. Loss of functional constraint in the stem-loop-forming sequence is proposed for this acceleration.      

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.     

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.     

Differences in Crossover Frequency and Distribution among Three Sibling Species of Drosophila

Comparisons of the genetic and cytogenetic maps of three sibling species of Drosophila reveal marked differences in the frequency and cumulative distribution of crossovers during meiosis. The maps for two of these species, Drosophila melanogaster and D. simulans, have previously been described, while this report presents new map data for D. mauritiana, obtained using a set of P element markers. A genetic map covering nearly the entire genome was constructed by estimating the recombination fraction for each pair of adjacent inserts. The P-based genetic map of mauritiana is ~1.8 times longer than the standard melanogaster map. It appears that mauritiana has higher recombination along the entire length of each chromosome, but the difference is greatest in centromere-proximal regions of the autosomes. The mauritiana autosomes show little or no centromeric recombinational suppression, a characteristic that is prominent in melanogaster. D. simulans appears to be intermediate both in terms of total map length and intensity of the autosomal centromeric effect. These interspecific differences in recombination have important evolutionary implications for DNA sequence organization and variability. In particular, mauritiana is expected to differ from melanogaster in patterns and amounts of sequence variation and transposon insertions.     

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.      

Sequence Analysis of Active Mariner Elements in Natural Populations of Drosophila Simulans

Active and inactive mariner elements from natural and laboratory populations of Drosophila simulans were isolated and sequenced in order to assess their nucleotide variability and to compare them with previously isolated mariner elements from the sibling species Drosophila mauritiana and Drosophila sechellia. The active elements of D. simulans are very similar among themselves (average 99.7% nucleotide identity), suggesting that the level of mariner expression in different natural populations is largely determined by position effects, dosage effects and perhaps other factors. Furthermore, the D. simulans elements exhibit nucleotide identities of 98% or greater when compared with mariner elements from the sibling species. Parsimony analysis of mariner elements places active elements from the three species into separate groups and suggests that D. simulans is the species from which mariner elements in D. mauritiana and D. sechellia are most likely derived. This result strongly suggests that the ancestral form of mariner among these species was an active element. The two inactive mariner elements sequenced from D. simulans are very similar to the inactive peach element from D. mauritiana. The similarity may result from introgression between D. simulans and D. mauritiana or from selective constraints imposed by regulatory effects of inactive elements.     

Wanderings of Hobo: A Transposon in Drosophila Melanogaster and its Close Relatives

The transposon hobo is present in the genomes of Drosophila melanogaster and Drosophila simulans (and D. mauritiana and probably D. sechellia, based on Southern blots) as full-size elements and internally deleted copies. The full-size melanogaster, simulans and mauritiana hobo elements are 99.9% identical at the DNA sequence level, and internally deleted copies in these species essentially differ only in having deletions. In addition to these, hobo-related sequences are present and detectable with a hobo probe in all these species. Those in D. melanogaster are 86-94% identical to the canonical hobo, but with many indels. We have sequenced one that appears to be inserted in heterochromatin (GenBank Acc. No. AF520587). It is 87.6% identical to the canonical hobo, but quite fragmented by indels, with remnants of other transposons inserted in and near it, and clearly is defunct. Numerous similar elements are found in the sequenced D. melanogaster genome. It has recently been shown that some are fixed in the euchromatic genome, but it is probable that still more reside in heterochromatic regions not included in the D. melanogaster genome database. They are probably all relics of an earlier introduction of hobo into the ancestral species. There appear to have been a minimum of two introductions of hobo into the melanogaster subgroup, and more likely three, two ancient and one quite recent. The recent introduction of hobo was probably followed by transfers between the extant species (whether 'horizontally' or by infrequent interspecific hybridization).     

Intraspecific diversity of nucleotide sequences within the adenine + thymine-rich region of mitochondrial DNA molecules of Drosophila mauritiana, Drosophila melanogaster and Drosophila simulans.

Mitochondrial DNA (mtDNA) molecules from Drosophila mauritiana, D. melanogaster, and D. simulans contain a single adenine + thymine (A+T)-rich region, which is similarly located in all molecules, but varies in size among these species. Using agarose gel electrophoresis and electron microscopy, a difference in occurrence of one EcoRI site, and a difference in size (approximately 0.7 kb) of the A+T-rich regions was found between mtDNA molecules of flies of two female lines of D. mauritiana. In heteroduplexes constructed between these two kinds of mtDNA molecules, two or three regions of strand separation, each comprising single strands of unequal length, were apparent near the center of the A+T-rich region. Using the structural differences between D. mauritiana mtDNA molecules it was demonstrated the mtDNA of this species is maternally inherited. Differences in length of A+T-rich regions were also found between mtDNA molecules of two geographically separated strains of D. melanogaster, and between mtDNA molecules of two geographically separated strains of D. simulans. However, in both cases, in heteroduplexes constructed between mtDNA molecules of different strains of one species, the A+T-rich regions appeared completely paired.     

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.     

One of the copia genes is adjacent to satellite DNA in Drosophila melanogaster.

A method for purifying sequences adjacent to satellite DNA in the heterochromatin of D. melanogaster is described. A cloned DNA segment containing part of a copia gene adjacent to 1.688 g/cm3 satellite DNA has been isolated. The copia genes compose a repeated gene family which codes for abundant cytoplasmic poly(a)-containing RNA (Young and Hogness, 1977; Finnegan et al., 1978). We have identified two major poly (A)-containing RNA species [5.2 and 2.1 kilobases (kb)] produced by the copia gene family. The cloned segment contains copia sequences homologous to the 5' end of RNA within 0.65 kb of the 1.688 satellite DNA sequences. Seven different cloned copia genes from elsewhere in the genome have also been isolated, and a 5.2 kb region present in five of the clones was identified as copia by heteroduplex analysis. In addition, three ususual copies of copia were found: a "partial" copy of the gene (3.7 kb) which has one endpoint in common with the 5.2 kb unit; a copia gene flanked on one side by a 1.6 kb sequence and on the other by the same 1.6 kb sequence in the inverted orientation; and a copia gene flanked only on one side by the same sequence.     

Distribution of hobo transposable elements in the genus Drosophila

This study describes the distribution of hobo-hybridizing sequences in the genus Drosophila. Southern blot analysis of 134 species revealed that hobo sequences are limited to the melanogaster and montium subgroups of the melanogaster-species group. Of the hobo-bearing species, only D. melanogaster and two of its sibling species, D. simulans and D. mauritiana, were found to contain potentially complete hobo elements. The distribution of hobo sequences is one of the narrowest distributions thus far described for any Drosophila transposable element.     

hobo transposable elements in Drosophila melanogaster and D. simulans.

Genomic patterns of occurrence of the transposable element hobo are polymorphic in the sibling species Drosophila melanogaster and D. simulans. Most tested strains of both species have apparently complete (3.0 kb) and smaller hobo elements (H lines), but in both species some strains completely lack such canonical hobo elements (E lines). The occurrence of H and E lines in D. simulans as well as in D. melanogaster implies that an hypothesis of recent introduction in the latter species is inadequate to explain the phylogenetic occurrence of hobo. Particular internally deleted elements, the approximately 1.5 kb Th1 and Th2 elements, are abundant in many lines of D. melanogaster, and an analogous 1.1 kb internally deleted element, h del sim, is abundant in most lines of D. simulans. Besides the canonical hobo sequences, both species (and their sibling species D. sechellia and D. mauritiana) have many hobo-hybridizing sequences per genome that do not appear to be closely related to the canonical hobo sequence.     

Evolution of the A+T-rich region of mitochondrial DNA in the melanogaster species subgroup of Drosophila

We determined the nucleotide sequences of two regions in the A+T-rich region of mitochondrial DNA (mtDNA) in the siI and siII types of D. simulans, the maII type of D. mauritiana, and D. sechellia. The sequences were aligned with those of the corresponding regions of siIII of D. simulans and maI of D. mauritiana, D. melanogaster, and D. yakuba. The type I and type II elements and the T-stretches were detected in all eight of the mtDNA types compared, indicating that the three elements are essential in the A+T-rich region of this species subgroup. The alignment revealed several short repetitive sequences and relatively large deletions in the central portions of the region. In the highly conserved sequence elements in the type II elements, the substitution rates were not uniform among lineages and acceleration in the substitution rate might have been due to loss of functional constraint in the stem–loop-forming sequences predicted in the type II elements. Patterns of nucleotide substitutions observed in the A+T-rich region were further compared with those in the coding regions and in the intergenic regions of mtDNA. Substitutions between A and T were particularly repressed in the highly conserved sequence elements and in the intergenic regions compared with those in the A+T-rich region excluding the highly conserved sequence elements and in the fourfold degenerate sites in the coding regions. The functional and structural characteristics of the A+T-rich region that might be involved in this substitutional bias are discussed.     

Extensive diversity among Drosophila species with respect to nucleotide sequences within the adenine + thymine-rich region of mitochondrial DNA molecules.

Mitochondrial DNA (mtDNA) molecules from species of the genus Drosophila contain a region exceptionally rich in adenine + thymine (A+T). Using agarose gel electrophoresis and electron microscopy, we determined that in the mtDNA molecules of D. melanogaster, D. simulans, D. mauritiana, D. yakuba, D. takahashii, and D. virilis, the A+T-rich regions, which are 5.1, 4.8, 4.6, 1.1, 2.2, and 1.0 kilobase pairs in size, respectively, are at homologous locations relative to various common EcoRI and HindIII cleavage sites. Under conditions highly permissive for base pairing (35% formamide), heteroduplexes were constructed between EcoRI fragments and whole circular molecules of mtDNAs of the above mentioned six species in a variety of combinations. Complete pairing of molecules outside the A+T-rich region was found in all heteroduplexes examined. However, in contrast, A+T-rich regions of the different species failed to pair in all but those combinations of mtDNAs involving the three most closely related species. In heteroduplexes between D. melanogaster and D. simulans, and between D. melanogaster and D. mauritiana mtDNAs, up to 35% of the A+T-rich regions appeared double-stranded. These data indicate that much more extensive divergence of sequences has occurred in A+T-rich regions than in other regions of Drosophila mtDNA molecules.     

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.     

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.