Metabolic rate and directional nucleotide substitution in animal mitochondrial DNA

There is marked heterogeneity of nucleotide composition in mitochondrial DNA across divergent animals. Differences in nucleotide composition presumably reflect differences in directional nucleotide substitution for A+T or G+C nucleotides. In mitochondrial DNA, there is A+T directional nucleotide substitution in most (if not all) animals surveyed, and the magnitude of directional A+T nucleotide substitution differs greatly within and among groups. Differences in directional nucleotide substitution among lineages of mammals can be explained by changes in metabolic physiology. This relationship is thought to be mediated by the effect of oxygen radicals because these toxic compounds are by-products of aerobic metabolism and are known mutagens. Association between metabolism and nucleotide composition provides additional evidence in favor of the hypothesis that rates and patterns of nucleotide substitution in mitochondrial DNA can be influenced by factors that impinge on rates of endogenous DNA damage.      

On the estimation of the rate of nucleotide substitution for the control region of human mitochondrial DNA

To apply molecular clock for studying human evolution, the pattern of nucleotide substitution for the control region of human mtDNA was analyzed in detail. It is well known that the rate of nucleotide substitution for the control region is much higher than that for any other part of mtDNA. In this study, the higher substitution rate was attributed to the higher rate of transition-type substitution between pyrimidines within the D-loop part, whereas the rates of other types of substitution were essentially the same over the entire mtDNA molecule. Even within the control region, the rate and pattern of nucleotide substitution were different between the D-loop part and the rest. The rate and pattern for the non-D-loop part were very similar to those for fourfold-degenerate sites in the protein-coding region. In contrast, the D-loop and non-D-loop parts showed similarities in the base composition, whereas the base composition of fourfold-degenerate sites slightly different from that of the both parts of the control region. It is concluded, therefore, that the nucleotide frequencies of the control region should be used to estimate the number of substitutions (d) between the control region sequences. However, a method to verify the accuracy of the estimation of d by means of the transition/transversion (s/v) ratio was theoretically studied. It was suggested that the s/v ratio becomes constant over a wide range of d values only when the estimation of d is unbiased. On the basis of this result, the estimates of d previously obtained between human sequences were evaluated.     

The influence of adjacent nucleotides on the pattern of nucleotide substitution in mitochondrial introns of angiosperms

It has been known that in noncoding regions of the chloroplast genome, the pattern of nucleotide substitution is influenced by the two nucleotides flanking the substitution site. In a GC-rich environment, a bias toward transition was observed, whereas in an AT-rich environment, a bias toward transversion was observed. In this study, the influence of the two adjacent neighbors on the substitution pattern was observed in the first intron of the mitochondrial nad4 gene, although the AT content of this intron is only 48%. The proportion of transversions increases from 0.32 to 0.75 as the A + T content (number of A's + T's) of the two nearest neighbors increases from 0 to 2. This trend was also observed in another mitochondrial group I intron with an AT content of 64%. In addition, a similar, though weaker, effect was observed in vertebrate pseudogenes. So this effect is present in all three types of genomes. Furthermore, in contrast to the situation in the noncoding regions of chloroplast DNA, where most nucleotide substitutions occurred in the categories with an A + T content of either 1 or 2, nucleotide substitutions in the mitochondrial first nad4 intron occurred more evenly in three categories of different A + T contents. This might be due largely to the difference in the AT content (0.48 vs. 0.72) between the mitochondrial first nad4 intron and the chloroplast DNA regions studied.     

Complex mitochondrial DNA in Drosophila.

The larval mtDNA isolated from D. virilis, D. simulans and D. melanogaster exists in complex molecular forms in addition to the simple monomeric circular form. The frequency of circular dimers and oligomers is highly elevated in apparently normal larval tissues. These complex forms of mtDNA are separable on agarose gels. Hind III restriction endonuclease and electron microscopic analyses used in the present study have revealed that circular dimers are simply the circular concatemers of two monomeric circles which are arranged in a head-to-tail structure with no detectable heterologous regions such as insertions or deletions. The electrophoretic patterns of Hind III digested mtDNAs of D. simulans and D. melanogaster (sibling species) are identical and distinguishable from that of distantly related species, D. virilis.     

The effect of nucleotide substitution on DNA denaturation profiles.

The melting profiles were obtained for DNA restriction fragments of approx. 1150 bp with deletion of one, five or six base pairs making them different from each other. In all cases the deletions caused a shift of one melting peak without affecting the positions of the other three peaks. The effect amounted to 0.28 +/- 0.03C upon the deletion of one GC pair. The melting of DNA fragments was also studied by electrophoresis in denaturing gradient gels. The deletion of one GC pair was shown to cause an appreciable shift of the electrophoretic denaturation profile.     

Estimating the pattern of nucleotide substitution

Knowledge of the pattern of nucleotide substitution is important both to our understanding of molecular sequence evolution and to reliable estimation of phylogenetic relationships. The method of parsimony analysis, which has been used to estimate substitution patterns in real sequences, has serious drawbacks and leads to results difficult to interpret. In this paper a model-based maximum likelihood approach is proposed for estimating substitution patterns in real sequences. Nucleotide substitution is assumed to follow a homogeneous Markov process, and the general reversible process model (REV) and the unrestricted model without the reversibility assumption are used. These models are also applied to examine the adequacy of the model of Hasegawa et al. (J. Mol. Evol. 1985;22:160–174) (HKY85). Two data sets are analyzed. For the -globin pseudogenes of six primate species, the REV model fits the data much better than HKY85, while, for a segment of mtDNA sequences from nine primates, REV cannot provide a significantly better fit than HKY85 when rate variation over sites is taken into account in the models. It is concluded that the use of the REV model in phylogenetic analysis can be recommended, especially for large data sets or for sequences with extreme substitution patterns, while HKY85 may be expected to provide a good approximation. The use of the unrestricted model does not appear to be worthwhile.     

Direct estimation of the mitochondrial DNA mutation rate in Drosophila melanogaster

Mitochondrial DNA (mtDNA) variants are widely used in evolutionary genetics as markers for population history and to estimate divergence times among taxa. Inferences of species history are generally based on phylogenetic comparisons, which assume that molecular evolution is clock-like. Between-species comparisons have also been used to estimate the mutation rate, using sites that are thought to evolve neutrally. We directly estimated the mtDNA mutation rate by scanning the mitochondrial genome of Drosophila melanogaster lines that had undergone approximately 200 generations of spontaneous mutation accumulation (MA). We detected a total of 28 point mutations and eight insertion-deletion (indel) mutations, yielding an estimate for the single-nucleotide mutation rate of 6.2 × 10⁻⁸ per site per fly generation. Most mutations were heteroplasmic within a line, and their frequency distribution suggests that the effective number of mitochondrial genomes transmitted per female per generation is about 30. We observed repeated occurrences of some indel mutations, suggesting that indel mutational hotspots are common. Among the point mutations, there is a large excess of G→A mutations on the major strand (the sense strand for the majority of mitochondrial genes). These mutations tend to occur at nonsynonymous sites of protein-coding genes, and they are expected to be deleterious, so do not become fixed between species. The overall mtDNA mutation rate per base pair per fly generation in Drosophila is estimated to be about 10× higher than the nuclear mutation rate, but the mitochondrial major strand G→A mutation rate is about 70× higher than the nuclear rate. Silent sites are substantially more strongly biased towards A and T than nonsynonymous sites, consistent with the extreme mutation bias towards A+T. Strand-asymmetric mutation bias, coupled with selection to maintain specific nonsynonymous bases, therefore provides an explanation for the extreme base composition of the mitochondrial genome of Drosophila.     

The ribosomal RNA genes of Drosophila mitochondrial DNA.

The nucleotide sequence of a segment of the mtDNA molecule of Drosophila yakuba which contains the A+T-rich region and the small and large rRNA genes separated by the tRNAval gene has been determined. The 5' end of the small rRNA gene was located by S1 protection analysis. In contrast to mammalian mtDNA, a tRNA gene was not found at the 5' end of the D. yakuba small rRNA gene. The small and large rRNA genes are 20.7% and 16.7% G+C and contain only 789 and 1326 nucleotides. The 5' regions of the small rRNA gene (371 nucleotides) and of the large rRNA gene (643 nucleotides) are extremely low in G+C (14.6% and 9.5%, respectively) and convincing sequence homologies between these regions and the corresponding regions of mouse mt-rRNA genes were found only for a few short segments. Nevertheless, the entire lengths of both of the D. yakuba mt-rRNA genes can be folded into secondary structures which are remarkably similar to secondary structures proposed for the rRNAs of mouse mtDNA. The replication origin-containing, A+T-rich region (1077 nucleotides; 92.8% A+T), which lies between the tRNAile gene and the small rRNA gene, lacks open reading frames greater than 123 nucleotides.     

Intralineage variation in the pattern ofrbcL nucleotide substitution

Variation in chloroplastrbcL sequences was studied in representative species of four different lineages: the tribeRubieae (Rubiaceae), and the generaDrosera (Droseraceae),Nothofagus (Nothofagaceae) andIlex (Aquifoliaceae). Each lineage has its particular non-overlapping set ofrbcL polymorphic sites, indicating that common unconstrainedrbcL sites are not shared. Large differences in the rate and pattern of nucleotide substitution are observed among the four lineages. The genusIlex has the lowest rate of substitution, the lowest transition/transversion ratio, the lowest synonymous/replacement ratio and the lowest number of substitutions at the third codon position. An apparent relationship of these measures to the age of the lineages is observed. The A + T content and codon use among the four lineages are very similar and, apparently, cannot account for the observed differences in patterns of nucleotide substitution. However, the A + T content of the two bases immediately flanking the polymorphic sites is higher inIlex than in the other lineages. This could be correlated with the transversion/transition bias observed inIlex. The particularly low synonymous/replacement ratio found inIlex could also be explained by the small population sizes of species in this genus.     

The genomic rate of adaptive amino acid substitution in Drosophila

The proportion of amino acid substitutions driven by adaptive evolution can potentially be estimated from polymorphism and divergence data by an extension of the McDonald-Kreitman test. We have developed a maximum-likelihood method to do this and have applied our method to several data sets from three Drosophila species: D. melanogaster, D. simulans, and D. yakuba. The estimated number of adaptive substitutions per codon is not uniformly distributed among genes, but follows a leptokurtic distribution. However, the proportion of amino acid substitutions fixed by adaptive evolution seems to be remarkably constant across the genome (i.e., the proportion of amino acid substitutions that are adaptive appears to be the same in fast-evolving and slow-evolving genes; fast-evolving genes have higher numbers of both adaptive and neutral substitutions). Our estimates do not seem to be significantly biased by selection on synonymous codon use or by the assumption of independence among sites. Nevertheless, an accurate estimate is hampered by the existence of slightly deleterious mutations and variations in effective population size. The analysis of several Drosophila data sets suggests that approximately 25% +/- 20% of amino acid substitutions were driven by positive selection in the divergence between D. simulans and D. yakuba.     

Site-by-site estimation of the rate of substitution and the correlation of rates in mitochondrial DNA

A maximum likelihood method for independently estimating the relative rate of substitution at different nucleotide sites is presented. With this method, the evolution of DNA sequences can be analyzed without assuming a specific distribution of rates among sites. To investigate the pattern of correlation of rates among sites, the method was applied to a data set consisting of the protein-coding regions of the mitochondrial genome from 10 vertebrate species. Rates appear to be strongly correlated at distances up to 40 codons apart. Furthermore, there appears to be some higher order correlation of sites approximately 75 codons apart. The method of site-by-site estimation of the rate of substitution may also be applied to examine other aspects of rate variation along a DNA sequence and to assess the difference in the support of a tree along the sequence.     

Methylation pattern of mouse mitochondrial DNA.

The pattern of methylation of mouse mitochondrial DNA (mtDNA) was studied using several techniques. By employing a sensitive analytical procedure it was possible to show that this DNA contains the modified base 5-methylcytosine (m5Cyt). This residue occurred exclusively at the dinucleotide sequence CpG at a frequency of 3 to 5%. The pattern of methylation was further investigated by determining the state of methylation of several MspI (HpaII) sites. Different sites were found to be methylated to a different extent, implying that methylation of mtDNA is nonrandom. Based on the known base composition and nucleotide sequence of mouse mtDNA, the dinucleotide sequence CpG was found to be underrepresented in this DNA. The features of mtDNA methylation (CpG methylation, partial methylation of specific sites and CpG underrepresentation) are also characteristic of vertebrate nuclear DNA. This resemblance may reflect functional relationship between the mitochondrial and nuclear genomes.     

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.     

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.     

Homologous recombination and the pattern of nucleotide substitution in Ehrlichia ruminantium

Patterns of nucleotide substitution at orthologous loci were examined between three genomes of Ehrlichia ruminantium, the causative agent of heartwater disease of ruminants. The most recent common ancestor of two genomes (Erwe and Erwo) belonging to the Welgevonden strain was estimated to have occurred 26,500-57,000 years ago, while the most recent common ancestor of these two genomes and the Erga genome (Gardel strain) was estimated to have occurred 2.1-4.7 million years ago. The search for genes showing extremely high values of the number of synonymous substitutions per site was used to identify genes involved in past homologous recombination. The most striking case involved the map1 gene, encoding major antigenic protein-1; evidence for homologous recombination is consistent with previous phylogenetic analysis of map1 alleles. At this and certain other loci, homologous recombination may have contributed to the evolution of host-pathogen interactions. In addition, comparison of the patterns of synonymous and nonsynonymous substitution provided evidence for positive selection favoring a high level of amino acid change between the Welgevonden and Gardel strains at a locus of unknown function (designated Erum4340 in the Erwo genome).     

Baculovirus expression reconstitutes Drosophila mitochondrial DNA polymerase.

Drosophila mitochondrial DNA polymerase has been reconstituted and purified from baculovirus-infected insect cells. Baculoviruses encoding full-length and mature forms of the catalytic and accessory subunits were generated and used in single and co-infection studies. Recombinant heterodimeric holoenzyme was reconstituted in both the mitochondria and cytoplasm of Sf9 cells and required the mitochondrial presequences in both subunits. The recombinant holoenzyme contains DNA polymerase and 3'-5' exonuclease that are stimulated substantially by both salt and mitochondrial single-stranded DNA-binding protein. Thus, the recombinant enzyme exhibits biochemical properties indistinguishable from those of the native enzyme from Drosophila embryos. Production of the catalytic subunit alone yielded soluble protein with the chromatographic properties of the heterodimeric holoenzyme. However, the purified catalytic core has a 50-fold lower specific activity. This provides evidence of a critical role for the accessory subunit in the catalytic efficiency of Drosophila mitochondrial DNA polymerase.     

3'-->5' exonuclease in Drosophila mitochondrial DNA polymerase. Substrate specificity and functional coordination of nucleotide polymerization and mispair hydrolysis.

A mispair-specific 3'-->5' exonuclease copurifies quantitatively with the near-homogeneous Drosophila gamma polymerase (Kaguni, L.S., and Olson, M.W. (1989) Proc. Natl. Acad. Sci. U.S.A. 86, 6469-6473). The exonuclease and polymerase exhibit similar reaction requirements and optima, suggesting functional coordination of their activities. Under nonpolymerization conditions, the 3'-->5' exonuclease hydrolyzes 3'-terminal mispairs approximately 15-fold more efficiently than 3'-terminal base pairs on primed single-stranded DNA substrates, whereas it does not discriminate between any of three specific mispairs (dAMP:dAMP;dGMP:dGMP; dGMP:dAMP). Under polymerization conditions, gamma polymerase does not extend a 3'-terminal mispair from the "stationary" state, even in the presence of a large excess of the next correct nucleotide. Instead, 3'-terminal mispairs are hydrolyzed quantitatively by the 3'-->5' exonuclease over the reaction time course. During DNA synthesis by gamma polymerase in the "polymerization" mode, limited misincorporation and subsequent mispair extension do occur. Here, it appears that misincorporation and not mispair extension is rate-limiting. Template-primer challenge experiments suggest that the mechanism of template-primer transfer from the 3'-->5' exonuclease active site to the DNA polymerase active site is intermolecular; transfer from the exonuclease to polymerase mode appears to require dissociation and reassociation of mitochondrial DNA polymerase.