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.     

Complete nucleotide sequence of the A+T-rich region of Drosophila mauritiana mitochondrial DNA

We determined the complete nucleotide sequence of the A+T-rich region of the maII type of mtDNA in D. mauritiana. The nucleotide sequence was found to contain 3,206 bp. Three types of conserved element, i.e., type I element, type II element, and T-stretch, were included in this sequence, as reported for D. melanogaster. Comparison between the two species revealed that the type I elements were less conserved than the type II elements. However, each of these type I elements contained a G-stretch within a loop of a putative stem-loop-forming sequence, which has also been observed in D. melanogaster. Moreover, in both type I and type II repeat arrays, the elements closest to the T-stretch diverged the most, due to nucleotide substitution and/or the insertion of short repeats. Sequence comparison of the two complete sequences of the A+T-rich region of D. melanogaster and the maII type of D. mauritiana, as well as comparison of partial sequences in other types of mtDNA within the melanogaster complex, suggested that the A+T-rich region in this complex has been maintained by concerted evolution after the duplication of two types of element, i.e., type I and type II.     

Mitochondrial DNA in the Drosophila Melanogaster Complex

Mitochondrial DNA in the complex Drosophila melanogaster was among the first studied in metazoans. The variability of the molecule was extensively studied using restriction enzymes, gene sequence and recently sequence of the whole coding region. Within the complex, seven major haplotypes have been described, one (me) in D. melanogaster, three in D. simulans (siI, siII, siIII), two in D. mauritiana (maI, maII), and one in D. sechellia (se). The molecular distance between the haplotypes is comprised between 1 and 5%, except for siII and maI, which are virtually identical. The nucleotide diversity within each of these haplotypes is very low, varying from 0 to 0.0005. Most of the cytoplasms are infected by the bacterium Wolbachia and different bacterial strains infect cytoplasms harboring different mtDNA types. mtDNA polymorphism is discussed in relation with Wolbachia, nuclear polymorphism and speciation events.     

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.     

Analysis of copia sequence variation within and between Drosophila species

The sequences of the 5' long-terminal repeat (LTR) and adjacent leader regions of 27 full-length copia elements isolated from natural populations of Drosophila melanogaster, D. simulans, and D. mauritiana are presented. Phylogenetic analyses indicate that although D. melanogaster copia elements are distinct from those of D. simulans and D. mauritiana, the elements of these latter two species are not distinguishable from one another. LTRs and adjacent 5' leader regions of elements isolated from D. simulans and D. mauritiana are structurally similar to one another and carry substantial deletional variation mapping to regions previously identified as being of potential importance for copia expression.      

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).      

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.     

Paleo-demography of the Drosophila melanogaster subgroup: application of the maximum likelihood method

The species divergence times and demographic histories of Drosophila melanogaster and its three sibling species, D. mauritiana, D. simulans, and D. yakuba, were investigated using a maximum likelihood (ML) method. Thirty-nine orthologous loci for these four species were retrieved from DDBJ/EMBL/GenBank database. Both autosomal and X-linked loci were used in this study. A significant degree of rate heterogeneity across loci was observed for each pair of species. Most loci have the GC content greater than 50% at the third codon position. The codon usage bias in Drosophila loci is considered to result in the high GC content and the heterogenous rates across loci. The chi-square, G, and Fisher's exact tests indicated that data sets with 11, 23, and 9 pairs of DNA sequences for the comparison of D. melanogaster with D. mauritiana, D. simulans, and D. yakuba, respectively, retain homogeneous rates across loci. We applied the ML method to these data sets to estimate the DNA sequence divergences before and after speciation of each species pair along with their standard deviations. Using 1.6 x 10(-8) as the rate of nucleotide substitutions per silent site per year, our results indicate that the D. melanogaster lineage split from D. yakuba approximately 5.1 +/- 0.8 million years ago (mya), D. mauritiana 2.7 +/- 0.4 mya, and D. simulans 2.3 +/- 0.3 mya. It implies that D. melanogaster became distinct from D. mauritiana and D. simulans at approximately the same time and from D. yakuba no earlier than 10 mya. The effective ancestral population size of D. melanogaster appears to be stable over evolutionary time. Assuming 10 generations per year for Drosophila, the effective population size in the ancestral lineage immediately prior to the time of species divergence is approximately 3 x 10(6), which is close to that estimated for the extant D. melanogaster population. The D. melanogaster did not encounter any obvious bottleneck during the past 10 million years.     

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.     

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.      

The 5S ribosomal genes in the Drosophila melanogaster species subgroup. Nucleotide sequence of a 5S unit from Drosophila simulans and Drosophila teissieri

The 5S genes of the eight species of the D. melanogaster subgroup have been mapped. The spacers, in contrast with coding regions, differ markedly between most species. One 5S gene unit has been sequenced for both D. simulans and D. teissieri. The mature 5S RNA region in these two species is identical to the corresponding region of D. melanogaster. Only 5 nucleotide variations occur between the D. melanogaster and D. simulans 5S gene spacers. The spacer in D. teissieri is very different. Only two segments, located one at each side of the coding region, are clearly homologous to corresponding sequences of D. melanogaster and D. simulans.     

Divergence and heterogeneity of the histone gene repeating units in the Drosophila melanogaster species subgroup

The repeating units of the histone gene cluster containing the H1, H2A, H2B and H4 genes were amplified by PCR from the Drosophila melanogaster species subgroup, i.e., D. yakuba, D. erecta, D. sechellia, D. mauritiana, D. teissieri and D. orena. The PCR products were cloned and their nucleotide sequences of about 4.6-4.8kbp were determined to elucidate the mechanism of molecular evolution of the histone gene family. The heterogeneity among the histone gene repeating units was 0.6% and 0.7% for D. yakuba and D. sechellia, respectively, indicating the same level of heterogeneity as in the H3 gene region of D. melanogaster. Divergence of the genes among species even in the most closely related ones was much greater than the heterogeneity among family members, indicating a concerted mode of evolution for the histone gene repeating units. Among the species in the D. melanogaster species subgroup, the histone gene regions as well as 3rd codon position of the coding region showed nearly the same GC contents. These results suggested that the previous conclusion on analysis of the H3 gene regions, the gene family evolution in a concerted fashion, holds true for the whole histone gene repeating unit.     

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.     

Identical satellite DNA sequences in sibling species of Drosophila

The evolution of simple satellite DNAs was examined by DNA-DNA hybridization of ten Drosophila melanogaster satellite sequences to DNAs of the sibling species, Drosophila simulans and Drosophila erecta. Seven of these repeat types are present in tandem arrays in D. simulans and each of the ten sequences is repeated in D. erecta. In thermal melts, six of the seven satellite sequences in D. simulans and seven of the ten sequences in D. erecta melted within 1 deg.C of the corresponding values in D. melanogaster. The remaining sequences melted within 3 deg.C of the homologous hybrids. Therefore, there is little or no alteration in those satellite sequences held in common, despite a period of about ten million years since the divergence of D. melanogaster and D. simulans from a common ancestor. Simple satellite sequences appear to be more highly conserved than coding regions of the genome, on a per nucleotide basis. Since multiple copies of three satellite sequences could not be detected in D. simulans yet are present in D. erecta, a species more distantly related to D. melanogaster than is D. simulans, these sequences show discontinuities in evolution. There were major quantitative variations between species, showing that satellite DNAs are prone to massive amplification or diminution events over timespans as short as those separating sibling species. In D. melanogaster, these sequences amount to 21% of the genome but only 5% in D. simulans and 0.4% in D. erecta. There was a general trend of lower abundance with evolutionary distance for most satellites, suggesting that the amounts of different satellite sequences do not vary independently during evolution.     

The molecular evolution of the alcohol dehydrogenase and alcohol dehydrogenase-related genes in the Drosophila melanogaster species subgroup

The DNA sequences of the Adh genes of three members of the Drosophila melanogaster species subgroup have been determined. This completes the Adh sequences of the eight species of this subgroup. Two species, D. yakuba and D. teissieri, possess processed Adh pseudogenes. In all of the species of the subgroup, a gene of unknown function, Adhr, is located about 300 bp 3' to Adh. Although this gene is experiencing a higher rate of synonymous substitution than Adh, it is more constrained at the amino acid level. Phylogenetic relationships between all eight members of the melanogaster subgroup have been analyzed using a variety of methods. All analyses suggested that the D. yakuba and D. teissieri pseudogenes have a single common ancestor, rather than evolving independently in each species, and that D. melanogaster is the sister species to D. simulans, D. sechellia, and D. mauritiana. The evolutionary relationships of the latter three species remain equivocal.      

DNA from Drosophila melanogaster region 61C7/C8 has high evolution rate

The results of a comparative study of cloned DNA fragments of Drosophila simulans, D. mauritiana, D. teissieri, and D. erecta are presented. The fragments were amplified in PCR with primers specified to the region of D. melanogaster interband 61C7/C8. The uniqueness of all cloned fragments in the genomes of these species was confirmed. A comparative analysis of nucleotide sequences revealed that the rate of evolution of DNA from D. melanogaster interband 61C7/C8 is close to the rate of neutral evolution in the genus Drosophila.     

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.