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.     

Origin and direction of replication in mitochondrial DNA molecules from the genus Drosophila.

Mitochondrial DNA (mtDNA) obtained from ovaries of Drosophila simulans, D. mauritiana, D. takahashii, D. yakuba and D. virilis was examined by electron microscopy. From a consideration of the structural properties of replicative intermediates, it was concluded that in mtDNA molecules of each species, synthesis on one strand can be up to 97% complete before synthesis on the complementary strand is initiated. MtDNA molecules of each species contain a single A+T-rich region which shows species-specific size variation from 1.0 kb (D. virilis) to 4.8 kb (D. simulans), and maps at the same position in all molecules relative to three common EcoRI sites. The structural properties of complex forms, interpreted as having originated from replicative intermediates, and produced by either partial denaturation or EcoRI digestion, are consistent with the hypothesis that replication is initiated within the A+T-rich region and proceeds unidirectionally around the molecule towards the nearest common EcoRI site. The replication origin is located near the center of the A+T-rich region in D. simulans and D. mauritiana, but lies closer to that end of the A+T-rich region which is distal to the nearest common EcoRI site in D. takahashii, D. yakuba and D. virilis.     

Circular DNA Molecules in the Genus Drosophila

The satellite DNA's from the embryos of five species of Drosophila (D. melanogaster, D. simulans, D. nasuta, D. virilis and D. hydei) have been analyzed for the presence of closed circular duplex DNA molecules, as determined by CsCl-EBr gradients. Circular DNA molecules were found in every species but D. melanogaster. Analyses of cell fractions from adult Drosophila and organ fractions from Drosophila larvae show that fractions containing mitochondria are highly enriched in these molecules.     

Determination of some biochemical and structural features of alcohol dehydrogenases from Drosophila simulans and Drosophila virilis. Comparison of their properties with the Drosophila melanogaster Adhs enzyme.

The biochemical properties of the enzyme alcohol dehydrogenase of two different Drosophila species, Drosophila simulans and Drosophila virilis, were studied and compared with those of Drosophila melanogaster Adhs enzyme. All of them consist of two identical subunits of molecular weight 27800 and share significant similarities in function. The substrate specificities of these enzymes were characterized and Km(app.) and Vmax.(app.) values were calculated. All these alcohol dehydrogenases show greater affinity for secondary rather than for primary alcohols. The amino acid compositions of the three enzymes were determined, and there is a close similarity between the D. simulans and the D. melanogaster enzymes, but there are significant differences from the alcohol dehydrogenase of D. virilis. The N-terminal amino acid is blocked and the C-terminal amino acid is the same for all three alcohol dehydrogenases. The enzymes from the three species were carboxymethylated and digested with trypsin. The peptide 'maps' reveal, as expected, more homologies between the enzymes of D. simulans and D. melanogaster than with the enzyme of D. virilis.     

Neutral and Non-Neutral Evolution of Drosophila Mitochondrial DNA

To test hypotheses of neutral evolution of mitochondrial DNA (mtDNA), nucleotide sequences were determined for 1515 base pairs of the NADH dehydrogenase subunit 5 (ND5) gene in the mitochondrial DNA of 29 lines of Drosophila melanogaster and 9 lines of its sibling species Drosophila simulans. In contrast to the patterns for nuclear genes, where D. melanogaster generally exhibits much less nucleotide polymorphism, the number of segregating sites was slightly higher in a global sample of nine ND5 sequences in D. melanogaster (s = 8) than in the nine lines of D. simulans (s = 6). When compared to variation at nuclear loci, the mtDNA variation in D. melanogaster does not depart from neutral expectations. The ND5 sequences in D. simulans, however, show fewer than half the number of variable sites expected under neutrality when compared to sequences from the period locus. While this reduction in variation is not significant at the 5% level, HKA tests with published restriction data for mtDNA in D. simulans do show a significant reduction of variation suggesting a selective sweep of variation in the mtDNA in this species. Tests of neutral evolution based on the ratios of synonymous and replacement polymorphism and divergence are generally consistent with neutral expectations, although a significant excess of amino acid polymorphism within both species is localized in one region of the protein. The rate of mtDNA evolution has been faster in D. melanogaster than in D. simulans and the population structure of mtDNA is distinct in these species. The data reveal how different rates of mtDNA evolution between species and different histories of neutral and adaptive evolution within species can compromise historical inferences in population and evolutionary biology.     

Satellite Dnas in the Embryos of Various Species of the Genus Drosophila

A tentative evolutionary pattern has been found for two classes of the multiple satellite DNA's found in the genus Drosophila. The satellite DNA's from five Drosophila species (D. melanogaster, D. simulans, D. nasuta, D. virilis and D. hydei) were analyzed and found to fall into three arbitrary CsCl buoyant density classes: Class I, rho = 1.661-1.669 g cm(-3), DNA molecules composed of primarily dA and dT moieties; Class II, rho = 1.685 and rho = 1.692, DNA molecules of low GC content; and Class III, rho = 1.711, a DNA of high GC composition. The dAT satellite DNA's appear in all the species studied except D. hydei, the species of most recent evolutionary divergence, whereas the heavy satellite appears only in the two species of most recent divergence, D. virilis and D. hydei.     

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.     

Further study on selective transmission of mitochondrial DNA in heteroplasmic lines of Drosophila melanogaster.

The temperature-dependent transmission of mitochondrial DNA (mtDNA) was investigated in heteroplasmic lines of Drosophila melanogaster established by germ-plasm transplantation. Using D. melanogaster, D. simulans and D. mauritiana as germ-plasm donors, five recipient-donor combinations of heteroplasmy, differing from those previously examined (Matsuura et al., 1991), were constructed. For intraspecific reciprocal combinations, donor mtDNA in one combination was retained at 25 degrees C but was almost lost by the tenth generation at 19 degrees C. In the reciprocal, the proportion of the same type of recipient mtDNA decreased more quickly at 19 degrees C than 25 degrees C. Decreasing rates at 19 degrees C in the reciprocals differed from each other. For interspecific combinations, two species were used as germ-plasm donors. Donor mtDNA derived from D. simulans was lost at both temperatures and the rate of decrease was greater at 19 degrees C than 25 degrees C. The proportion of donor mtDNA derived from D. mauritiana increased at a greater rate at 25 degrees C than 19 degrees C when using two different strains of D. melanogaster as recipients. These results suggest that both the nuclear and two types of mitochondrial genomes are involved in the selective transmission of mtDNA.     

Regulative interactions between cells of wing discs from different dipteran species

The mechanism by which patterns are produced appears to be repeated in each segment of an animal, and it has been proposed that it may even have been conserved in evolution so that different species would have the same system of positional information. This idea has been tested by mixing cells of a defined fragment of the wing disc of Drosophila melanogaster with wing disc fragments of five other dipteran species to assay the ability of these disc fragments to stimulate intercalary regeneration of the D. melanogaster cells. The genetically marked (y; mwh) D. melanogaster fragment was mechanically mixed with wing discs or wing disc fragments of four drosophilids (D. melanogaster as a control, D. virilis, D. hydei, Zaprionus vittiger), of Musca domestica, and of Piophila casei. The mixed aggregates were cultured in vivo for 7 days, then metamorphosed in D. melanogaster larval hosts. The D. melanogaster fragments were only stimulated to regenerate when combined with complementary fragments from D. melanogaster or D. virilis wing discs. In the combination between D. melanogaster and D. hydei, the tissue formed integrated mosaic patterns, but no regeneration ensued. The one positive result (D. melanogaster mixed with D. virilis) shows that positional cues can be exchanged and correctly interpreted between cells of different species. The negative results do not prove that the mechanism for establishing patterns is different in the tested species, but may be due to incompatibilities that are not related to pattern formation.     

Analysis of nucleotide substitutions of mitochondrial DNAs in Drosophila melanogaster and its sibling species

To study the rate and pattern of nucleotide substitution in mitochondrial DNA (mtDNA), we cloned and sequenced a 975-bp segment of mtDNA from Drosophila melanogaster, D. simulans, and D. mauritiana containing the genes for three transfer RNAs and parts of two protein- coding genes, ND2 and COI. Statistical analysis of synonymous substitutions revealed a predominance of transitions over transversions among the three species, a finding differing from previous results obtained from a comparison of D. melanogaster and D. yakuba. The number of transitions observed was nearly the same for each species comparison, including D. yakuba, despite the differences in divergence times. However, transversions seemed to increase steadily with increasing divergence time. By contrast, nonsynonymous substitutions in the ND2 gene showed a predominance of transversions over transitions. Most transversions were between A and T and seemed to be due to some kind of mutational bias to which the A + T-rich mtDNA of Drosophila species may be subject. The overall rate of nucleotide substitution in Drosophila mtDNA appears to be slightly faster (approximately 1.4 times) than that of the Adh gene. This contrasts with the result obtained for mammals, in which the mtDNA evolves approximately 10 times faster than single-copy nuclear DNA. We have also shown that the start codon of the COI gene is GTGA in D. simulans and GTAA in D. mauritiana. These codons are different from that of D. melanogaster (ATAA).      

Mitochondrial DNA evolution in themelanogaster species subgroup ofDrosophila

Detailed restriction maps (40 cleavage sites on average) of mitochondrial DNAs (mtDNAs) from the eight species of the melanogaster species subgroup of Drosophila were established. Comparison of the cleavage sites allowed us to build a phylogenetic tree based on the matrix of nucleotide distances and to select the most parsimonious network. The two methods led to similar results, which were compared with those in the literature obtained from nuclear characters. The three chromosomally homosequential species D. simulans, D. mauritiana, and D. sechellia are mitochondrially very related, but exhibit complex phylogenetic relationships. D. melanogaster is their closest relative, and the four species form a monophyletic group (the D. melanogaster complex), which is confirmed by the shared unusual length of their mt genomes (18-19 kb). The other four species of the subgroup (D. yakuba, D. teissieri, D. erecta, and D. orena) are characterized by a much shorter mt genome (16-16.5 kb). The monophyletic character of the D. yakuba complex, however, is questionable. Two species of this complex, D. yakuba and D. teissieri, are mitochondrially indistinguishable (at the level of our investigation) in spite of their noticeable allozymic and chromosomal divergence. Finally, mtDNA distances were compared with the nuclear-DNA distances thus far established. These sequences seem to evolve at rather similar rates, the mtDNA rate being barely double that of nuclear DNA.     

Transgenic tools for members of the genus Drosophila with sequenced genomes

The sequencing of the genomes of 12 Drosophila species has created an opportunity for much in the way of comparative molecular analyses amongst these species. To aid that endeavor, we have made several transformation vectors based on the piggyBac transposon with 3xP3-EGFP and -ECFP transgenic markers that should be useful for mutagenesis and establishing the GAL4/UAS system in these species. We have tested the ability of mini-white to be used as a marker for insertional mutagenesis, and have observed mini-white derived pigmentation of the testes sheath in a subset of lines from D. pseudoobscura and D. virilis. We have incorporated a source of piggyBac transposase into nine Drosophila species, and have demonstrated the functionality of these transposase lines for mobilization of marked inserts in vivo. Additionally, we tested the ability of a D. melanogaster nanos enhancer element to drive expression of GAL4 in D. melanogaster, D. simulans, D. erecta, D. yakuba, D. pseudoobscura, and D. virilis. The efficacy of the nos-Gal4 transgene was determined by measuring the response of UAS-EGFPtub in all six species. Our results show that D. melanogaster nos-Gal4 drives expression in other species, to varying degrees, in similar spatiotemporal domains in the ovaries, testes, and embryos as seen in D. melanogaster. However, expression levels are variable, demonstrating the possible need to use species-specific promoters in some cases. In summary, we hope to provide a set of guidelines and basic tools, based upon this work, for both insertional mutagenesis and GAL4/UAS system-based experiments in multiple species of Drosophila.     

Distribution of transposable elements in neotropical species of Drosophila

The phylogenetic distribution of transposable families, P, gypsy, hobo, I, and mariner has been analyzed in 33 species of 11 groups of neotropical Drosophila and a Drosophilidae species Zygotrica vittimaculosa, using squash blot and dot blot. Genomic DNA of almost all neotropical species tested hybridized with gypsy probe and some species showed a particularly strong hybridization signal, as D. gaucha, D. virilis, and species of flavopilosa group. The hobo element was restricted to melanogaster group and some strains of D. willistoni. Only D. simulans DNA showed hybridization to mariner probe in all species tested and D. simulans and D. melanogaster showed hybridization with I element probe. P element homologous sequence was present in D. melanogaster and all species and strains of the willistoni and saltans groups tested. The presence of at least one P-homologous sequence was detected in Drosophila mediopunctata. This one was the only P-bearing species of all six tested from the tripunctata group. Four different pairs of primers homologous to segments of the canonical sequence of D. melanogaster's P were used to amplify specific sequences from D. mediopunctata DNA, showing the occurrence of seemingly well-conserved P-homologous sequences. This revised version was published online in July 2006 with corrections to the Cover Date.     

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.     

A comprehensive study of genic variation in natural populations of Drosophila melanogaster. VII. Varying rates of genic divergence as revealed by two-dimensional electrophoresis.

Four sibling species from the melanogaster subgroup (Drosophila melanogaster, D. simulans, D. sechellia, and D. mauritiana) were studied for genetic divergence, by high-resolution two-dimensional protein electrophoresis (2DE) coupled with ultrasensitive silver staining. A total of eight tissues from larval and adult developmental stages representing both gonadal (germ-line) and nongonadal (somatic) tissues were analyzed for protein divergence between species. Close to 400 polypeptides (protein spots) were scored from each tissue and species, and protein divergence was measured on the basis of qualitative differences (presence/absence) of protein spots in pairwise species comparisons. The observed levels of genic divergence varied among tissues and among species. When larval hemolymph proteins (which are known to be highly polymorphic) were excluded, there was no evidence to suggest that either the larval or adult-stage proteins, as a whole, are more diverged than the other; variation between different tissues rather than between developmental stages appears to be the most significant factor affecting genetic divergence between species. The reproductive tissue (testis and accessory gland) showed more divergence than did the nonreproductive tissue; D. melanogaster testis (from both larvae and adult males) showed the highest level of divergence. In view of the previous observation that D. simulans, D. mauritiana, and D. sechellia show similar but significantly less reproductive isolation from each other than from D. melanogaster, the present results suggest a correlation between the levels of reproductive-tract-protein divergence and the degree of reproductive isolation in these species.     

Unraveling Selection in the Mitochondrial Genome of Drosophila

We examine mitochondrial DNA variation at the cytochrome b locus within and between three species of Drosophila to determine whether patterns of variation conform to the predictions of neutral molecular evolution. The entire 1137-bp cytochrome b locus was sequenced in 16 lines of Drosophila melanogaster, 18 lines of Drosophila simulans and 13 lines of Drosophila yakuba. Patterns of variation depart from neutrality by several test criteria. Analysis of the evolutionary clock hypothesis shows unequal rates of change along D. simulans lineages. A comparison within and between species of the ratio of amino acid replacement change to synonymous change reveals a relative excess of amino acid replacement polymorphism compared to the neutral prediction, suggestive of slightly deleterious or diversifying selection. There is evidence for excess homozygosity in our world wide sample of D. melanogaster and D. simulans alleles, as well as a reduction in the number of segregating sites in D. simulans, indicative of selective sweeps. Furthermore, a test of neutrality for codon usage shows the direction of mutations at third positions differs among different topological regions of the gene tree. The analyses indicate that molecular variation and evolution of mtDNA are governed by many of the same selective forces that have been shown to govern nuclear genome evolution and suggest caution be taken in the use of mtDNA as a ``neutral' molecular marker.     

Nucleotide variation at the no-on-transient A gene in Drosophila littoralis

The no-on-transient A (nonA) gene encodes a putative RNA-binding protein, and mutations in this gene are known to affect vision, male courtship song and viability in Drosophila melanogaster. Here we have sequenced the coding region of the nonA gene of Drosophila littoralis and compared it with those of Drosophila virilis and D. melanogaster. All portions of nonA appeared to be conserved between D. littoralis and D. virilis, while the 5' region of the gene of these two species showed high divergence from that of a more distantly-related species, D. melanogaster. The same was true for the glycine repeat regions. No significant deviation from neutrality was observed in the analysis of intraspecific nucleotide variation in 5' or 3' region of the nonA gene in D. littoralis population. Also, comparison of D. littoralis sequences with homologous sequence of D. virilis suggests that the gene is evolving neutrally in D. virilis group. Divergence of the 5' regions between D. virilis group species and D. melanogaster could be a result of positive selection, but this finding is obscured by the long divergence time of the species groups.     

Historicity and the Population Genetics of Drosophila Melanogaster and D. Simulans

We summarize data showing that there is population structure in African populations of Drosophila from the melanogaster-simulans complex. In D. melanogaster, population structuring is found at individual loci, but is obscured by population structuring for large inversions that simultaneously affect several loci. In D. simulans, molecular polymorphism at the X-linked vermilion locus suggests that different groups of populations have been geographically isolated for some time. Invading populations are probably derived from different areas in Africa. European populations originate from an east African population that was probably not at a demographic equilibrium. The origin of the Antilles population is apparently different and is as yet unknown. In south-western France, populations from these two species undergo different population structuring at the scale of a few kilometres: D. melanogaster makes up a large panmictic population, whereas D. simulans forms a metapopulation that is divided into smaller demes.