Polymorphic Sites and the Mechanism of Evolution in Human Mitochondrial DNA

Twelve restriction enzymes were used to screen for the presence or absence of cleavage sites at 441 locations in the mitochondrial DNA of 112 humans from four continents. Cleavage maps were constructed by comparison of DNA fragment sizes with those expected from the published sequence for one human mtDNA. One hundred and sixty-three of the sites were polymorphic, i.e., present in some individuals but absent from others, 278 sites being invariant. These polymorphisms probably result from single base substitutions and occur in all functional regions of the genome.--In 77 cases, it was possible to specify the exact nature and location (within a restriction site) of the mutation responsible for the absence of a restriction site in a known human mtDNA sequence and its presence in another human mtDNA. Fifty-two of these 77 gain mutations occur in genes coding for proteins, 34 being silent and 18 causing amino acid replacements; moreover, nine of the replacements are radical.--Notable also is the anomalous ratio of transitions to transversions required to account for these 77 restriction site differences between the known human mtDNA sequences and other human mtDNAs. This ratio is lower for most groups of restriction sites than has been reported from sequence comparisons of limited parts of the mtDNA genome in closely related mammals, perhaps indicating a special functional role or sensitivity to mutagenesis for palindromic regions containing high levels of guanine and cytosine.--From the genomic distribution of the 163 polymorphic sites, it is inferred that the level of point mutational variability in tRNA and rRNA genes is nearly as high as in protein-coding genes but lower than in noncoding mtDNA. Thus, the functional constraints operating on components of the protein-synthetic apparatus may be lower for mitochondria than for other systems. Furthermore, the mitochondrial genes for tRNAs that recognize four codons are more variable than those recognizing only two codons.--Among the more variable of the human mitochondrial genes coding for proteins is that for subunit 2 of cytochrome oxidase; this polypeptide appears to have been evolving about five times faster in primates than in other mammals. Cytochrome c, a nuclearly encoded protein that interacts directly with the oxidase 2 subunit in electron transport, has also evolved faster in primates than in rodents or ungulates. This example, along with that for the mitochondrial rRNA genes and the nuclear genes coding for mitochondrial ribosomal proteins, provides evidence for coevolution between specific nuclear and mitochondrial genes.     

Large, rapidly evolving intergenic spacers in the mitochondrial DNA of the salamander family Ambystomatidae (Amphibia: Caudata)

We report the presence, in the mitochondrial DNA (mtDNA) of all of the sexual species of the salamander family Ambystomatidae, of a shared 240- bp intergenic spacer between tRNAThr and tRNAPro. We place the intergenic spacer in context by presenting the sequence of 1,746 bp of mtDNA from Ambystoma tigrinum tigrinum, describe the nucleotide composition of the intergenic spacer in all of the species of Ambystomatidae, and compare it to other coding and noncoding regions of Ambystoma and several other vertebrate mtDNAs. The nucleotide substitution rate of the intergenic spacer is approximately three times faster than the substitution rate of the control region, as shown by comparisons among six Ambystoma macrodactylum sequences and eight members of the Ambystoma tigrinum complex. We also found additional inserts within the intergenic spacers of five species that varied from 87-444 bp in length. The presence of the intergenic spacer in all sexual species of Ambystomatidae suggests that it arose at least 20 MYA and has been a stable component of the ambystomatid mtDNA ever since. As such, it represents one of the few examples of a large and persistent intergenic spacer in the mtDNA of any vertebrate clade.      

Identification and characterization of new long conserved noncoding sequences in vertebrates

Comparative sequence analyses have identified highly conserved genomic DNA sequences, including noncoding sequences, between humans and other species. By performing whole-genome comparisons of human and mouse, we have identified 611 conserved noncoding sequences longer than 500 bp, with more than 95% identity between the species. These long conserved noncoding sequences (LCNS) include 473 new sequences that do not overlap with previously reported ultraconserved elements (UCE), which are defined as aligned sequences longer than 200 bp with 100% identity in human, mouse, and rat. The LCNS were distributed throughout the genome except for the Y chromosome and often occurred in clusters within regions with a low density of coding genes. Many of the LCNS were also highly conserved in other mammals, chickens, frogs, and fish; however, we were unable to find orthologous sequences in the genomes of invertebrate species. In order to examine whether these conserved sequences are functionally important or merely mutational cold spots, we directly measured the frequencies of ENU-induced germline mutations in the LCNS of the mouse. By screening about 40.7 Mb, we found 35 mutations, including mutations at nucleotides that were conserved between human and fish. The mutation frequencies were equivalent to those found in other genomic regions, including coding sequences and introns, suggesting that the LCNS are not mutational cold spots at all. Taken together, these results suggest that mutations occur with equal frequency in LCNS but are eliminated by natural selection during the course of evolution.     

Identification and characterization of new long conserved noncoding sequences in vertebrates

Comparative sequence analyses have identified highly conserved genomic DNA sequences, including noncoding sequences, between humans and other species. By performing whole-genome comparisons of human and mouse, we have identified 611 conserved noncoding sequences longer than 500 bp, with more than 95% identity between the species. These long conserved noncoding sequences (LCNS) include 473 new sequences that do not overlap with previously reported ultraconserved elements (UCE), which are defined as aligned sequences longer than 200 bp with 100% identity in human, mouse, and rat. The LCNS were distributed throughout the genome except for the Y chromosome and often occurred in clusters within regions with a low density of coding genes. Many of the LCNS were also highly conserved in other mammals, chickens, frogs, and fish; however, we were unable to find orthologous sequences in the genomes of invertebrate species. In order to examine whether these conserved sequences are functionally important or merely mutational cold spots, we directly measured the frequencies of ENU-induced germline mutations in the LCNS of the mouse. By screening about 40.7 Mb, we found 35 mutations, including mutations at nucleotides that were conserved between human and fish. The mutation frequencies were equivalent to those found in other genomic regions, including coding sequences and introns, suggesting that the LCNS are not mutational cold spots at all. Taken together, these results suggest that mutations occur with equal frequency in LCNS but are eliminated by natural selection during the course of evolution.

     

Multiple nuclear-gene phylogenies: application to pinnipeds and comparison with a mitochondrial DNA gene phylogeny

Phylogenetic analyses of closely related species should use information from multiple, independent genes with relatively high rates of sequence evolution. To investigate species for which there are few prior sequence data for single-copy nuclear (scnDNA) genes, primers for gene amplification can be designed to highly conserved regions of exons in order to amplify both coding (exons) and noncoding (introns) sequences. We have explored this approach in a phylogenetic analysis of six species of pinnipeds that, together with terrestrial carnivore outgroups, encompass divergence times < or = 40-50 Mya. We sequenced one intron from each of the aldolase A (ALD-A), aldolase C (ALD-C), and histone H2AF genes; one exon from the major-histocompatibility-complex DQA gene; a H2AF processed pseudogene (psi H2AF); and, for comparison with the nuclear genes, the 5' portion of the mitochondrial DNA (mtDNA) control region. The pinniped psi H2AF genes were found to be of limited use because they were paralogous with the gene in the outgroup. The rate of silent substitution in scnDNA (primarily introns) was 5-10-fold lower than that for mtDNA control region I, and scnDNA sequence divergence increased linearly with time < or = 40-50 Mya. Alleles at three polymorphic scnDNA loci (ALD-A, H2AF, and DQA) in the southern elephant seal were paraphyletic with respect to the allele from the closely related northern elephant seal, while the more numerous mtDNA alleles were monophyletic. This we attribute to the consequences of a higher mutation rate rather than to a lower effective population size of mtDNA compared with scnDNA. Within the short (i.e., < 500-bp) sequences of individual scnDNA sequences, phylogenetically informative variation was insufficient to obtain robust phylogenies. However, the combined scnDNA sequences produced a well-supported phylogeny congruent with that derived from mtDNA. This analysis illustrates the high resolution of mtDNA sequences compared with a similar length of scnDNA sequence, but it also demonstrates the utility of combining information from multiple short scnDNA sequences obtained using broadly applicable primers.      

Rates of DNA Duplication and Mitochondrial DNA Insertion in the Human Genome

The hundreds of mitochondrial pseudogenes in the human nuclear genome sequence (numts) constitute an excellent system for studying and dating DNA duplications and insertions. These pseudogenes are associated with many complete mitochondrial genome sequences and through those with a good fossil record. By comparing individual numts with primate and other mammalian mitochondrial genome sequences, we estimate that these numts arose continuously over the last 58 million years. Our pairwise comparisons between numts suggest that most human numts arose from different mitochondrial insertion events and not by DNA duplication within the nuclear genome. The nuclear genome appears to accumulate mtDNA insertions at a rate high enough to predict within-population polymorphism for the presence/absence of many recent mtDNA insertions. Pairwise analysis of numts and their flanking DNA produces an estimate for the DNA duplication rate in humans of 2.2 × 10^–9 per numt per year. Thus, a nucleotide site is about as likely to be involved in a duplication event as it is to change by point substitution. This estimate of the rate of DNA duplication of noncoding DNA is based on sequences that are not in duplication hotspots, and is close to the rate reported for functional genes in other species.     

Genomic selective constraints in murid noncoding DNA

Recent work has suggested that there are many more selectively constrained, functional noncoding than coding sites in mammalian genomes. However, little is known about how selective constraint varies amongst different classes of noncoding DNA. We estimated the magnitude of selective constraint on a large dataset of mouse-rat gene orthologs and their surrounding noncoding DNA. Our analysis indicates that there are more than three times as many selectively constrained, nonrepetitive sites within noncoding DNA as in coding DNA in murids. The majority of these constrained noncoding sites appear to be located within intergenic regions, at distances greater than 5 kilobases from known genes. Our study also shows that in murids, intron length and mean intronic selective constraint are negatively correlated with intron ordinal number. Our results therefore suggest that functional intronic sites tend to accumulate toward the 5' end of murid genes. Our analysis also reveals that mean number of selectively constrained noncoding sites varies substantially with the function of the adjacent gene. We find that, among others, developmental and neuronal genes are associated with the greatest numbers of putatively functional noncoding sites compared with genes involved in electron transport and a variety of metabolic processes. Combining our estimates of the total number of constrained coding and noncoding bases we calculate that over twice as many deleterious mutations have occurred in intergenic regions as in known genic sequence and that the total genomic deleterious point mutation rate is 0.91 per diploid genome, per generation. This estimated rate is over twice as large as a previous estimate in murids.     

How rapidly does the human mitochondrial genome evolve?

The results of an empirical nucleotide-sequencing approach indicate that the evolution of the human mitochondrial noncoding D-loop is both more rapid and more complex than is revealed by standard phylogenetic approaches. The nucleotide sequence of the D-loop region of the mitochondrial genome was determined for 45 members of a large matrilineal Leber hereditary optic neuropathy pedigree. Two germ-line mutations have arisen in members of one branch of the family, thereby leading to triplasmic descendants with three mitochondrial genotypes. Segregation toward the homoplasmic state can occur within a single generation in some of these descendants, a result that suggests rapid fixation of mitochondrial mutations as a result of developmental bottlenecking. However, slow segregation was observed in other offspring, and therefore no single or simple pattern of segregation can be generalized from the available data. Evidence for rare mtDNA recombination within the D-loop was obtained for one family member. In addition to these germ-line mutations, a somatic mutation was found in the D-loop of one family member. When this genealogical approach was applied to the nucleotide sequences of mitochondrial coding regions, the results again indicated a very rapid rate of evolution.     

The contribution of slippage-like processes to genome evolution

Simple sequences present in long (>30 kb) sequences representative of the single-copy genome of five species (Homo sapiens, Caenorhabditis elegans Saccharomyces cerevisiae, E. coli, and Mycobacterium leprae) have been analyzed. A close relationship was observed between genome size and the overall level of sequence repetition. This suggested that the incorporation of simple sequences had accompanied increases of genome size during evolution. Densities of simple sequence motifs were higher in noncoding regions than in coding regions in eukaryotes but not in eubacteria. All five genomes showed very biased frequency distributions of simple sequence motifs in all species, particularly in eukaryotes where AAA and TTT predominated. Interspecific comparisons showed that noncoding sequences in eukaryotes showed highly significantly similar frequency distributions of simple sequence motifs but this was not true of coding sequences. ANOVA of the frequency distributions of simple sequence motifs indicated strong contributions from motif base composition and repeat unit length, but much of the variation remained unexplained by these parameters. The sequence composition of simple sequences therefore appears to reflect both underlying sequence biases in slippage-like processes and the action of selection. Frequency distributions of simple sequence motifs in coding sequences correlated weakly or not at all with those in noncoding sequences. Selection on coding sequences to eliminate undesirable sequences may therefore have been strong, particularly in the human lineage.     

Markov chain analysis finds a significant influence of neighboring bases on the occurrence of a base in eucaryotic nuclear DNA sequences both protein-coding and noncoding

Summary Sixty-four eucaryotic nuclear DNA sequences, half of them coding and half noncoding, have been examined as expressions of first-, second-, or third-order Markov chains. Standard statistical tests found that most of the sequences required at least second-order Markov chains for their representation, and some required chains of third order. For all 64 sequences the observed one-step second-order transition count matrices were effective in predicting the two-step transition count matrices, and 56 of 64 were effective in predicting the three-step transition count matrices. The departure from random expectation of the observed first- and second-order transition count matrices meant that a considerable sample of eucaryotic nuclear DNA sequences, both protein coding and noncoding, have significant local structure over subsequences of three to five contiguous bases, and that this structure occurs throughout the total length of the sequence. These results suggested that present DNA sequences may have arisen from the duplication, concatenation, and gradual modification of very early short sequences.     

Markov chain analysis finds a significant influence of neighboring bases on the occurrence of a base in eucaryotic nuclear DNA sequences both protein-coding and noncoding

Sixty-four eucaryotic nuclear DNA sequences, half of them coding and half noncoding, have been examined as expressions of first-, second-, or third-order Markov chains. Standard statistical tests found that most of the sequences required at least second-order Markov chains for their representation, and some required chains of third order. For all 64 sequences the observed one-step second-order transition count matrices were effective in predicting the two-step transition count matrices, and 56 of 64 were effective in predicting the three-step transition count matrices. The departure from random expectation of the observed first- and second-order transition count matrices meant that a considerable sample of eucaryotic nuclear DNA sequences, both protein coding and noncoding, have significant local structure over subsequences of three to five contiguous bases, and that this structure occurs throughout the total length of the sequence. These results suggested that present DNA sequences may have arisen from the duplication, concatenation, and gradual modification of very early short sequences.     

Molecular analyses of mtDNA deletion mutations in microdissected skeletal muscle fibers from aged rhesus monkeys

Mitochondrial DNA (mtDNA) deletion mutations co-localize with electron transport system (ETS) abnormalities in rhesus monkey skeletal muscle fibers. Using laser capture microdissection in conjunction with PCR and DNA sequence analysis, mitochondrial genomes from single sections of ETS abnormal fibers were characterized. All ETS abnormal fibers contained mtDNA deletion mutations. Deletions were large, removing 20-78% of the genome, with some to nearly all of the functional genes lost. In one-third of the deleted genomes, the light strand origin was deleted, whereas the heavy strand origin of replication was conserved in all fibers. A majority (27/39) of the deletion mutations had direct repeat sequences at their breakpoints and most (36/39) had one breakpoint within or in close proximity to the cytochrome b gene. Several pieces of evidence support the clonality of the mtDNA deletion mutation within an ETS abnormal region of a fiber: (a) only single, smaller than wild-type, PCR products were obtained from each ETS abnormal region; (b) the amplification of mtDNA from two regions of the same ETS abnormal fiber identified identical deletion mutations, and (c) a polymorphism was observed at nucleotide position 16103 (A and G) in the wild-type mtDNA of one animal (sequence analysis of an ETS abnormal region revealed that mtDNA deletion mutations contained only A or G at this position). Species-specific differences in the regions of the genomes lost as well as the presence of direct repeat sequences at the breakpoints suggest mechanistic differences in deletion mutation formation between rodents and primates.     

Hypervariable noncoding sequences in Saccharomyces cerevisiae

Compared to protein-coding sequences, the evolution of noncoding sequences and the selective constraints placed on these sequences is not well characterized. To compare the evolution of coding and noncoding sequences, we have conducted a survey for DNA polymorphism at five randomly chosen loci among a diverse collection of 81 strains of Saccharomyces cerevisiae. Average rates of both polymorphism and divergence are 40% lower at noncoding sites and 90% lower at nonsynonymous sites in comparison to synonymous sites. Although noncoding and coding sequences show substantial variability in ratios of polymorphism to divergence, two of the loci, MLS1 and PDR10, show a higher rate of polymorphism at noncoding compared to synonymous sites. The high rate of polymorphism is not accompanied by a high rate of divergence and is limited to a few small regions. These hypervariable regions include sites with three segregating bases at a single site and adjacent polymorphic sites. We show that this clustering of polymorphic sites is significantly greater than one would expect on the basis of the spacing between polymorphic fourfold degenerate sites. Although hypervariable noncoding sequences could result from selection on regulatory mutations, they could also result from transient mutational hotspots.     

EVOLUTION AND MITOCHONDRIAL DNA IN NEUROSPORA CRASSA

We studied mitochondrial DNA variability in 19 natural Neurospora crassa isolates and one wild-type isolate to examine evolution of these fungi and their mitochondrial DNA (mtDNA). We combined restriction endonuclease analysis of natural isolate mtDNA with DNA-DNA hybridization to cloned EcoR I fragments of a wild-type genome to discriminate between length mutations and site changes due to nucleotide substitution. Most variability was due to length mutations (insertions and deletions); genome size could vary 25% between pairs of isolates. Length-mutation distribution was not random, nor simply explained by the presence of coding versus noncoding regions. Restriction-site changes were few; the estimated amount of nucleotide substitution per nucleotide between the most divergent pair of isolates was 0.78%. Evolutionary relationships among isolates based on both types of mutations were compatible, and suggest that geographically distinct populations of mitochondrial DNA exist in the biological species, N. crassa. In contrast, no such correlation was shown by the previously determined distribution of nuclear heterokaryon incompatibility genes in the same isolates (Mylyk, 1975, 1976).     

Lineage-specific selection in human mtDNA: lack of polymorphisms in a segment of MTND5 gene in haplogroup J

Human mitochondrial DNA (mtDNA) is a nonrecombining genome that codes for 13 subunits of the mitochondrial oxidative phosphorylation system, 2 rRNAs, and 22 tRNAs. Mutations have accumulated sequentially in mtDNA lineages that diverged tens of thousands of years ago. The genes in mtDNA are subject to different functional constraints and are therefore expected to evolve at different rates, but the rank order of these rates should be the same in all lineages of a phylogeny. Previous studies have indicated, however, that specific regions of mtDNA may have experienced different histories of selection in different lineages, possibly because of lineage-specific interactions or environmental factors such as climate. We report here on a survey for lineage-specific patterns of nucleotide polymorphism in human mtDNA. We calculated molecular polymorphism indices and neutrality tests for classes of functional sites and genes in 837 human mtDNA sequences, compared the results between continent-specific mtDNA lineages, and used two sliding window methods to identify differences in the patterns of polymorphism between haplogroups. A general correlation between nucleotide position and the level of nucleotide polymorphism was identified in the coding region of the mitochondrial genome. Nucleotide diversity in the protein-coding sequence of mtDNA was generally not much higher than that found for many genes in nuclear DNA. A comparison of nonsynonymous/synonymous rate ratios in the 13 protein-coding genes suggested differences in the relative levels of selection between haplogroups, including the European haplogroup clusters. Interestingly, a segment of the MTND5 gene was found to be almost void of segregating sites and nonsynonymous mutations in haplogroup J, which has been associated with susceptibility to certain complex diseases. Our results suggest that there are haplogroup-specific differences in the intensity of selection against particular regions of the mitochondrial genome, indicating that some mutations may be non-neutral within specific phylogenetic lineages but neutral within others.     

Mitochondrial COII Introgression into the Nuclear Genome of Gorilla gorilla

Numts are nonfunctional mitochondrial sequences that have translocated into nuclear DNA, where they evolve independently from the original mitochondrial DNA (mtDNA) sequence. Numts can be unintentionally amplified in addition to authentic mtDNA, complicating both the analysis and interpretation of mtDNA-based studies. Amplification of numts creates particular issues for studies on the noncoding, hypervariable 1 mtDNA region of gorillas. We provide data on putative numt sequences of the coding mitochondrial gene cytochrome oxidase subunit II (COII). Via polymerase chain reaction (PCR) and cloning, we obtained COII sequences for gorilla, orangutan, and human high-quality DNA and also from a gorilla fecal DNA sample. Both gorilla and orangutan samples yielded putative numt sequences. Phylogenetically more anciently transferred numts were amplified with a greater incidence from the gorilla fecal DNA sample than from the high-quality gorilla sample. Data on phylogenetically more recently transferred numts are equivocal. We further demonstrate the need for additional investigations into the use of mtDNA markers for noninvasively collected samples from gorillas and other primates.     

PREDICTING NUCLEAR GENE COALESCENCE FROM MITOCHONDRIAL DATA: THE THREE-TIMES RULE

Abstract.— Coalescence theory predicts when genetic drift at nuclear loci will result in fixation of sequence differences to produce monophyletic gene trees. However, the theory is difficult to apply to particular taxa because it hinges on genetically effective population size, which is generally unknown. Neutral theory also predicts that evolution of monophyly will be four times slower in nuclear than in mitochondrial genes primarily because genetic drift is slower at nuclear loci. Variation in mitochondrial DNA (mtDNA) within and between species has been studied extensively, but can these mtDNA data be used to predict coalescence in nuclear loci? Comparison of neutral theories of coalescence of mitochondrial and nuclear loci suggests a simple rule of thumb. The “three‐times rule” states that, on average, most nuclear loci will be monophyletic when the branch length leading to the mtDNA sequences of a species is three times longer than the average mtDNA sequence diversity observed within that species. A test using mitochondrial and nuclear intron data from seven species of whales and dolphins suggests general agreement with predictions of the three‐times rule. We define the coalescence ratio as the mitochondrial branch length for a species divided by intraspecific mtDNA diversity. We show that species with high coalescence ratios show nuclear monophyly, whereas species with low ratios have polyphyletic nuclear gene trees. As expected, species with intermediate coalescence ratios show a variety of patterns. Especially at very high or low coalescence ratios, the three‐times rule predicts nuclear gene patterns that can help detect the action of selection. The three‐times rule may be useful as an empirical benchmark for evaluating evolutionary processes occurring at multiple loci.     

Analysis of heteroplasmy in the major noncoding region of mitochondrial DNA in the blood and atherosclerotic plaques of carotid arteries

For identification of somatic mitochondrial DNA (mtDNA) mutations, the mtDNA major noncoding region (D-loop) sequence in blood samples and carotid atherosclerosis plaques from patients with atherosclerosis was analyzed. Five point heteroplasmic positions were observed in 4 of 23 individuals (17%). Only in two cases could heteroplasmy have resulted from somatic mutation, whereas three heteroplasmic positions were found in both vascular tissue and blood. In addition, length heteroplasmy in a polycytosine stretches was registered at nucleotide positions 303–315 in 16 individuals, and also in the 16184–16193 region in four patients. The results suggest that somatic mtDNA mutations can occur during atherosclerosis, but some heteroplasmic mutations may appear in all tissues, possibly being inherited.