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

Drosophila Melanogaster Mitochondrial DNA: Gene Organization and Evolutionary Considerations

The sequence of a 8351-nucleotide mitochondrial DNA (mtDNA) fragment has been obtained extending the knowledge of the Drosophila melanogaster mitochondrial genome to 90% of its coding region. The sequence encodes seven polypeptides, 12 tRNAs and the 3' end of the 16S rRNA and CO III genes. The gene organization is strictly conserved with respect to the Drosophila yakuba mitochondrial genome, and different from that found in mammals and Xenopus. The high A + T content of D. melanogaster mitochondrial DNA is reflected in a reiterative codon usage, with more than 90% of the codons ending in T or A, G + C rich codons being practically absent. The average level of homology between the D. melanogaster and D. yakuba sequences is very high (roughly 94%), although insertion and deletions have been detected in protein, tRNA and large ribosomal genes. The analysis of nucleotide changes reveals a similar frequency for transitions and transversions, and reflects a strong bias against G + C on both strands. The predominant type of transition is strand specific.     

Nucleotide sequence of a segment of Drosophila mitochondrial DNA that contains the genes for cytochrome c oxidase subunits II and III and ATPase subunit 6.

The nucleotide sequence of a segment of the mtDNA molecule of Drosophila yakuba has been determined, within which have been identified the genes for tRNAleuUUR, cytochrome c oxidase subunit II (COII), tRNAlys, tRNAasp, URFA6L, ATPase subunit 6 (ATPase6), cytochrome c oxidase subunit III (COIII) and tRNAgly. The genes are arranged in the order given and all are transcribed from the same strand of the molecule in a direction opposite to that in which replication proceeds around the molecule. The tRNAlys gene is unusual among mitochondrial tRNAlys genes in that it contains a CTT anticodon. The triplet AGA is used to specify an amino acid in all of the COII, COIII, ATPase6, and URFA6L genes. However, the AGA codons found in these four polypeptide genes correspond in position to codons which specify nine different amino acids, but never arginine, in the equivalent polypeptide gene which have been sequenced from mtDNAs of mouse, yeast and Zea mays.     

A cluster of six tRNA genes in Drosophila mitochondrial DNA that includes a gene for an unusual tRNAserAGY.

Genes for URF3, tRNAala, tRNAarg, tRNAasn, tRNAserAGY, tRNAglu, tRNAphe, and the carboxyl terminal segment of the URF5 gene have been identified within a sequenced segment of the mtDNA molecule of Drosophila yakuba. The genes occur in the order given. The URF5 and tRNAphe genes are transcribed in the same direction as replication while the URF3 and remaining five tRNA genes are transcribed in the opposite direction. Considerable differences exist in the relative arrangement of these genes in D. yakuba and mammalian mtDNA molecules. In the tRNAserAGY gene an eleven nucleotide loop, within which secondary structure formation seems unlikely, replaces the dihydrouridine arm, and both the variable loop (six nucleotides) and the T phi C loop (nine nucleotides) are larger than in any other D. yakuba tRNA gene. As available evidence is consistent with AGA codons specifying serine rather than arginine in the Drosophila mitochondrial genetic code, the possibility is considered that the 5'GCU anticodon of the D. yakuba tRNAserAGY gene can recognize AGA as well as AGY codons.     

Nuclear-mitochondrial epistasis and drosophila aging: introgression of Drosophila simulans mtDNA modifies longevity in D. melanogaster nuclear backgrounds

Under the mitochondrial theory of aging, physiological decline with age results from the accumulated cellular damage produced by reactive oxygen species generated during electron transport in the mitochondrion. A large body of literature has documented age-specific declines in mitochondrial function that are consistent with this theory, but relatively few studies have been able to distinguish cause from consequence in the association between mitochondrial function and aging. Since mitochondrial function is jointly encoded by mitochondrial (mtDNA) and nuclear genes, the mitochondrial genetics of aging should be controlled by variation in (1) mtDNA, (2) nuclear genes, or (3) nuclear-mtDNA interactions. The goal of this study was to assess the relative contributions of these factors in causing variation in Drosophila longevity. We compared strains of flies carrying mtDNAs with varying levels of divergence: two strains from Zimbabwe (<20 bp substitutions between mtDNAs), strains from Crete and the United States (approximately 20-40 bp substitutions between mtDNAs), and introgression strains of Drosophila melanogaster carrying mtDNA from Drosophila simulans in a D. melanogaster Oregon-R chromosomal background (>500 silent and 80 amino acid substitutions between these mtDNAs). Longevity was studied in reciprocal cross genotypes between pairs of these strains to test for cytoplasmic (mtDNA) factors affecting aging. The intrapopulation crosses between Zimbabwe strains show no difference in longevity between mtDNAs; the interpopulation crosses between Crete and the United States show subtle but significant differences in longevity; and the interspecific introgression lines showed very significant differences between mtDNAs. However, the genotypes carrying the D. simulans mtDNA were not consistently short-lived, as might be predicted from the disruption of nuclear-mitochondrial coadaptation. Rather, the interspecific mtDNA strains showed a wide range of variation that flanked the longevities seen between intraspecific mtDNAs, resulting in very significant nuclear x mtDNA epistatic interaction effects. These results suggest that even "defective" mtDNA haplotypes could extend longevity in different nuclear allelic backgrounds, which could account for the variable effects attributable to mtDNA haplogroups in human aging.     

The Mitochondrial Genome of the Honeybee Apis Mellifera: Complete Sequence and Genome Organization

The complete sequence of honeybee (Apis mellifera) mitochondrial DNA is reported being 16,343 bp long in the strain sequenced. Relative to their positions in the Drosophila map, 11 of the tRNA genes are in altered positions, but the other genes and regions are in the same relative positions. Comparisons of the predicted protein sequences indicate that the honeybee mitochondrial genetic code is the same as that for Drosophila; but the anticodons of two tRNAs differ between these two insects. The base composition shows extreme bias, being 84.9% AT (cf. 78.6% in Drosophila yakuba). In protein-encoding genes, the AT bias is strongest at the third codon positions (which in some cases lack guanines altogether), and least in second codon positions. Multiple stepwise regression analysis of the predicted products of the protein-encoding genes shows a significant association between the numbers of occurrences of amino acids and %T in codon family, but not with the number of codons per codon family or other parameters associated with codon family base composition. Differences in amino acid abundances are apparent between the predicted Apis and Drosophila proteins, with a relative abundance in the Apis proteins of lysine and a relative deficiency of alanine. Drosophila alanine residues are as often replaced by serine as conserved in Apis. The differences in abundances between Drosophila and Apis are associated with %AT in the codon families, and the degree of divergence in amino acid composition between proteins correlates with the divergence in %AT at the second codon positions. Overall, transversions are about twice as abundant as transitions when comparing Drosophila and Apis protein-encoding genes, but this ratio varies between codon positions. Marked excesses of transitions over chance expectation are seen for the third positions of protein-coding genes and for the gene for the small subunit of ribosomal RNA. For the third codon positions the excess of transitions is adequately explained as due to the restriction of observable substitutions to transitions for conserved amino acids with two-codon families; the excess of transitions over expectation for the small ribosomal subunit suggests that the conservation of nucleotide size is favored by selection.     

Drosophila mitochondrial DNA: a novel gene order

Part of the replication origin-containing A+T-rich region of the Drosophila yakuba mtDNA molecule and segments on either side of this region have been sequenced, and the genes within them identified. The data confirm that the small and large rRNA genes lie in tandem adjacent to that side of the A+T-rich region which is replicated first, and establish that a tRNAval gene lies between the two rRNA genes and that URF1 follows the large rRNA gene. The data further establish that the genes for tRNAile, tRNAgln, tRNAf-met and URF2 lie in the order given, on the opposite side of the A+T-rich region to the rRNA genes and, except for tRNAgln, are contained in the opposite strand to the rRNA, tRNAval and URF1 genes. This is in contrast to mammalian mtDNAs where all of these genes are located on the side of the replication origin which is replicated last, within the order tRNAphe, small (12S) rRNA, tRNAval, large (16S) rRNA, tRNAleu, URF1, tRNAile, tRNAgln, tRNAf-met and URF2, and, except tRNAgln, are all contained in the same (H) strand. In D. yakuba URF1 and URF2, the triplet AGA appears to specify an amino acid, which is again different from the situation found in mammalian mtDNAs, where AGA is used only as a rare termination codon.     

Nucleotide Sequence and Gene Organization of the Starfish Asterina Pectinifera Mitochondrial Genome

The 16,260-bp mitochondrial DNA (mtDNA) from the starfish Asterina pectinifera has been sequenced. The genes for 13 proteins, two rRNAs and 22 tRNAs are organized in an extremely economical fashion, similar to those of other animal mtDNAs, with some of the genes overlapping each other. The gene organization is the same as that for another echinoderm, sea urchin, except for the inversion of a 4.6-kb segment that contains genes for two proteins, 13 tRNAs and the 16S rRNA. Judging from the organization of the protein coding genes, mammalian mtDNAs resemble the sea urchin mtDNA more than that of the starfish. The region around the 3' end of the 12S rRNA gene of the starfish shows a high similarity with those for vertebrates. This region encodes a possible stem and loop structure; similar potential structures occur in this region of vertebrate mtDNAs and also in nonmitochondrial small subunit rRNA. A similar stem and loop structure is also found at the 3' end of the 16S rRNA genes in A. pectinifera, in another starfish Pisaster ochraceus, in vertebrates and in Drosophila, but not in sea urchins. The full sequence data confirm the presumption that AGA/AGG, AUA and AAA codons, respectively, code for serine, isoleucine, and asparagine in the starfish mitochondria, and that AGA/AGG codons are read by tRNA(GCU)(Ser), which possesses a truncated dihydrouridine arm, that was previously suggested from a partial mtDNA sequence. The structural characteristics of tRNAs and possible mechanisms for the change in the mitochondrial genetic code are also discussed.     

Transitions,transversions, and the molecular evolutionary clock

Summary Nucleotide substitutions in the form of transitions (purine-purine or pyrimidine-pyrimidine interchanges) and transversions (purine-pyrimidine interchanges) occur during evolution and may be complied by aligning the sequences of homologous genes. Referring to the genetic code tables, silent transitions take place in third positions of codons in family boxes and two-codon sets. Silent transversions in third positions occur only in family boxes, except for AC transversions between AGR and CGR arginine codons (R=A or G). Comparisons of several protein genes have been made, and various subclasses of transitional and transversional nucleotide substitutions have been compiled. Considerable variations occur among the relative proportions of transitions and transversions. Such variations could possibly be caused by mutator genes, favoring either transitions or, conversely, transversions, during DNA replication. At earlier stages of evolutionary divergence, transitions are usually more frequent, but there are exceptions. No indication was found that transversions usually originate from multiple substitutions in transitions.     

Evoluation of Drosophila mitochondrial DNAs. Analysis of heteroduplex molecules.

We have mapped the single block of non-homologous sequences and measured the extent and distribution of base-pair substitutions within the homologous sequences in Drosophila melanogaster: Drosophila virilis heteroduplex mitochondrial DNAs (mtDNAs). Of the 4.8 kilobases long, unusually (A + T)-rich region in D. melanogaster mtDNA, only 0.5 kilobases can react with related, but not identical sequences in D. virilis mtDNA, while the rest (4.3 kilobases in the long arm of a heteroduplex loop) is replaced by a shorter, non-homologous region (1.0 kilobases in the short arm of the loop). No additional heterologous regions are evident. Homologous sequences have accumulated on the average 15.5% base-pair changes. Regionally, these substitutions are relatively uniformly distributed (14.5--16.5%) except for a single, more conserved region (10--13%), which presumably represents the ribosomal cistrons. The lack of general sequence stability suggests that the invariant topographic organization of the nucleotide sequence, previously recognized among Drosophila mtDNAs, is under more stringent selection than the sequence per se.     

Complete DNA Sequence of the Mitochondrial Genome of the Black Chiton, Katharina Tunicata

The DNA sequence of the 15,532-base pair (bp) mitochondrial DNA (mtDNA) of the chiton Katharina tunicata has been determined. The 37 genes typical of metazoan mtDNA are present: 13 for protein subunits involved in oxidative phosphorylation, 2 for rRNAs and 22 for tRNAs. The gene arrangement resembles those of arthropods much more than that of another mollusc, the bivalve Mytilus edulis. Most genes abut directly or overlap, and abbreviated stop codons are inferred for four genes. Four junctions between adjacent pairs of protein genes lack intervening tRNA genes; however, at each of these junctions there is a sequence immediately adjacent to the start codon of the downstream gene that is capable of forming a stem-and-loop structure. Analysis of the tRNA gene sequences suggests that the D arm is unpaired in tRNA(ser(AGN)), which is typical of metazoan mtDNAs, and also in tRNA(ser(UCN)), a condition found previously only in nematode mtDNAs. There are two additional sequences in Katharina mtDNA that can be folded into structures resembling tRNAs; whether these are functional genes is unknown. All possible codons except the stop codons TAA and TAG are used in the protein-encoding genes, and Katharina mtDNA appears to use the same variation of the mitochondrial genetic code that is used in Drosophila and Mytilus. Translation initiates at the codons ATG, ATA and GTG. A + T richness appears to have affected codon usage patterns and, perhaps, the amino acid composition of the encoded proteins. A 142-bp non-coding region between tRNA(glu) and CO3 contains a 72-bp tract of alternating A and T.     

Platyhelminth mitochondrial DNA: Evidence for early evolutionary origin of a tRNAserAGN that contains a dihydrouridine arm replacement loop,and of serine-specifying AGA and AGG codons

Summary The nucleotide sequence of a segment of the mitochondrial DNA (mtDNA) molecule of the liver flukeFasciola hepatica (phylum Platyhelminthes, class Trematoda) has been determined, within which have been identified the genes for tRNA^ala, tRNA^asp, respiratory chain NADH dehydrogenase subunit I (ND1), tRNA^asn, tRNA^pro, tRNA^ile, tRNA^lys, ND3, tRNA^serAGN, tRNA^trp, and cytochromec oxidase subunit I (COI). The 11 genes are arranged in the order given and are all transcribed from the same strand of the molecule. The overall order of theF. hepatica mitochondrial genes differs from what is found in other metazoan mtDNAs. All of the sequenced tRNA genes except the one for tRNA^serAGN can be folded into a secondary structure with four arms resembling most other metazoan mitochondrial tRNAs, rather than the tRNAs that contain a TψC arm replacement loop, found in nematode mtDNAs. TheF. hepatica mitochondrial tRNA^serAGN gene contains a dihydrouridine arm replacement loop, as is the case in all other metazoan mtDNAs examined to date. AGA and AGG are found in theF. hepatica mitochondrial protein genes and both codons appear to specify serine. These findings concerningF. hepatica mtDNA indicate that both a dihydrouridine arm replacement loop-containing tRNA^serAGN gene and the use of AGA and AGG codons to specify serine must first have occurred very early in, or before, the evolution of metazoa.     

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.     

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.     

The complete mitochondrial DNA sequence of the crustacean Artemia franciscana

The complete mitochondrial DNA (mtDNA) sequence of the brine shrimp Artemia franciscana has been determined. It extends the present knowledge of mitochondrial genomes to the crustacean class and supplies molecular markers for future comparative studies in this large branch of the arthropod phylum. Artemia mtDNA is 15,822 nucleotides long, and when compared with its Drosophila counterpart, it shows very few gene rearrangements, merely affecting two tRNAs placed 3 downstream of the ND 2 gene. In this position a stem-loop secondary structure with characteristics similar to the vertebrate mtDNA L-strand origin of replication is found. This suggests that, associated with tRNA changes, the diversification of the mitochondrial genome from an ancestor common to crustacea and insects could be explained by errors in the mtDNA replication process. Although the gene content is the same as in most animal mtDNAs, the sizes of the protein coding genes are in some cases considerably smaller. Artemia mtDNA uses the same genetic code as found in insects, ATN and GTG are used as initiation codons, and several genes end in incomplete T or TA codons.Correspondence to: R. Garesse     

Characteristics of Nucleotide Substitution in the Hepatitis C Virus Genome: Constraints on Sequence Change in Coding Regions at Both Ends of the Genome

Comparison of complete genome sequences for different variants of hepatitis C virus (HCV) reveals several different constraints on sequence change. Synonymous changes are suppressed in coding regions at both 5′ and 3′ ends of the genome. No evidence was found for the existence of alternative reading frames or for a lower mutation frequency in these regions. Instead, suppression may be due to constraints imposed by RNA secondary structures identified within the core and NS5b genes. Nonsynonymous substitutions are less frequent than synonymous ones except in the hypervariable region of E2 and, to a lesser extent, in E1, NS2, and NS5b. Transitions are more frequent than transversions, particularly at the third position of codons where the bias is 16:1. In addition, nucleotide substitutions may not occur symmetrically since there is a bias toward G or C at the third position of codons, while T ↔ C transitions were twice as frequent as A ↔ G transitions. These different biases do not affect the phylogenetic analysis of HCV variants but need to be taken into account in interpreting sequence change in longitudinal studies. Received: 9 September 1996 / Accepted: 20 April 1997