Mitochondrial DNA polymorphism: evidence that variants detected by restriction enzymes differ in nucleotide sequence rather than in methylation

MITOCHONDRIAL-DNA RESTRICTION-SITE POLYMORPHISMS IN THE TELEOST FUNDULUS HETEROCLITUS SUPPORT SECONDARY INTERGRADATION

Fundulus heteroclitus is a highly polymorphic fish distributed along the Atlantic coast of North America. Several loci show directional changes in gene frequency with latitude (i.e., clines). Such directional changes have classically been described by two general models: primary and secondary intergradation. Previously, it has not been possible to distinguish between these models for Fundulus heteroclitus on the basis of allelic isozymes or morphological data. However, recent analysis of mitochondrial-DNA (mtDNA) restriction electromorphs helps resolve this issue. Mitochondrial-DNA samples from 48 individuals representing four populations were digested with 17 restriction endonucleases. After electrophoresis, the sizes of the mtDNA fragments were used to analyze the phylogenetic relatedness of fish collected over most of the species range. The analysis clearly identified two major races within the species: a northern and a southern form. The distribution of the mtDNA electromorphs, combined with zoogeographical changes in allelic isozymes and in eggs and adult morphologies (published elsewhere), makes the secondary-intergradation hypothesis most compelling.

Evolution of DNA sequence organization in mitochondrial genomes of Zea

Molecular mechanisms underlying the evolution of mtDNA in maize and three species of teosinte have been investigated. By using DNA transfer techniques and cloned fragments of maize mtDNA, changes in the position of homologous sequences in restriction digests were analyzed in related taxa. The resulting patterns indicate general conservation of sequence homology in all four species. One-third of the cloned fragments showed sequences conserved in homology and position of BamHI restriction fragments. Other fragments produced patterns indicating that extensive rearrangement of DNA sequences has taken place. In two cases, hybridization patterns revealed major changes in relative abundance of specific sequences that we believe are related to the molecular heterogeneity phenomenon of maize mtDNA. The evolution of mtDNA in these plants appears dramatically different from that in animal mtDNA, where base substitutions may account for most of the observed changes in restriction patterns.

Co-ordinate decrease in the expression of the mitochondrial genome and nuclear genes for mitochondrial proteins in the lactation-induced mitochondrial hypotrophy of rat brown fat

The relative abundance of the mitochondrial-encoded mRNAs for cytochrome c oxidase subunit II and NADH dehydrogenase subunit I was lower in brown adipose tissue (BAT) from lactating rats than in virgin controls. This decrease was in parallel with a significant decrease in mitochondrial 16 S rRNA levels and in the relative content of mitochondrial DNA in the tissue. BAT from lactating rats showed lowered mRNA expression of the nuclear-encoded genes for the mitochondrial uncoupling protein, subunit IV of cytochrome c oxidase and the adenine nucleotide translocase isoforms ANT1 and ANT2, whereas mRNA levels for the ATP synthase beta-subunit were unchanged. However, the relative content of this last protein was lower in BAT mitochondria from lactating rats than in virgin controls. It is concluded that lactation-induced mitochondrial hypotrophy in BAT is associated with a co-ordinate decrease in the expression of the mitochondrial genome and nuclear genes for mitochondrial proteins. This decrease is caused by regulatory events acting at different levels, including pre- and post-transcriptional regulation. BAT appears to be a useful model with which to investigate the molecular mechanisms involved in the co-ordination of the expression of the mitochondrial and nuclear genomes during mitochondrial biogenesis.

Effects of cold environment on mitochondrial genome expression in the rat: evidence for a tissue-specific increase in the liver, independent of changes in mitochondrial gene abundance

The abundance of the mitochondrially encoded mRNA for subunit II of cytochrome c oxidase (COII mRNA) increases in the liver of rats exposed to environmental cold stress (4 degrees C ambient temperature). Only transient increases or no changes in COII mRNA levels were observed in brown fat and soleus muscle respectively. The time course of the liver COII mRNA increase was compared with the effects of cold stress on mitochondrial 16S ribosomal RNA expression and indicated that cold induces a rapid (few hours) increase in the liver mitochondrial mRNA levels and high levels of both messenger and ribosomal RNA mitochondrial transcripts are present after a few days of cold exposure. No changes in mitochondrial DNA abundance relative to total cellular DNA were observed in the liver of rats at any time during cold stress. It is concluded that mitochondrial genome expression is specifically increased in the liver of cold-exposed rats through different mechanisms, independent of changes in mitochondrial genome abundance.

The mutational spectrum of single base-pair substitutions causing human genetic disease: patterns and predictions

Reports of single base-pair substitutions that cause human genetic disease and that have been located and characterized in an unbiased fashion were collated; 32% of point mutations were CG----TG or CG----CA transitions consistent with a chemical model of mutation via methylation-mediated deamination. This represents a 12-fold higher frequency than that predicted from random expectation, confirming that CG dinucleotides are indeed hotspots of mutation causing human genetic disease. However, since CG also appears hypermutable irrespective of methylation-mediated deamination, a second mechanism may also be involved in generating CG mutations. The spectrum of point mutations occurring outwith CG dinucleotides is also non-random, at both the mono- and dinucleotide, levels. An intrinsic bias in clinical detection was excluded since frequencies of specific amino acid substitutions did not correlate with the 'chemical difference' between the amino acids exchanged. Instead, a strong correlation was observed with the mutational spectrum predicted from the experimentally measured mispairing frequencies of vertebrate DNA polymerases alpha and beta in vitro. This correlation appears to be independent of any difference in the efficiency of enzymatic proofreading/mismatch-repair mechanisms but is consistent with a physical model of mutation through nucleotide misincorporation as a result of transient misalignment of bases at the replication fork. This model is further supported by an observed correlation between dinucleotide mutability and stability, possibly because transient misalignment must be stabilized long enough for misincorporation to occur. Since point mutations in human genes causing genetic disease neither arise by random error nor are independent of their local sequence environment, predictive models may be considered. We present a computer model (MUTPRED) based upon empirical data; it is designed to predict the location of point mutations within gene coding regions causing human genetic disease. The mutational spectrum predicted for the human factor IX gene was shown to resemble closely the observed spectrum of point mutations causing haemophilia B. Further, the model was able to predict successfully the rank order of disease prevalence and/or mutation rates associated with various human autosomal dominant and sex-linked recessive conditions. Although still imperfect, this model nevertheless represents an initial attempt to relate the variable prevalence of human genetic disease to the mutability inherent in the nucleotide sequences of the underlying genes.

RNAs containing mitochondrial ND6 and COI sequences present an abnormal structure in chemically induced rat hepatomas

We have constructed a cDNA library prepared from an hepatoma cell line (HTC cells) and isolated a clone, pHT 13, which corresponds to mRNAs present at a much higher level in rat hepatomas than in normal hepatocytes. The sequence of the pHT 13 insert has been previously published (Nucleic Acids Res. 1988, 16,10935). This clone contains mitochondrial DNA sequences with an abnormal organization, since it includes part of the NADH dehydrogenase subunit 6 (ND6) and of the cytochrome oxidase subunit I (COI) genes separated by 230 bases instead of 9 kb in mitochondrial genome from normal hepatocytes. In this work we show (1) that RNAs homologous to this clone are present in hepatoma cells but not in normal hepatocytes, (2) that a 3 kb fragment of tumor mitochondrial DNA contains both the ND6 and the COI sequences. The abnormal structure of the DNA is confirmed by Southern blot analysis which shows that distinct types of mitochondrial DNAs are present in hepatoma cells.

ATP inhibits nuclear and mitochondrial type I topoisomerases from human leukemia cells

Type I topoisomerases have been purified from nuclei and mitochondria of human acute lymphoblastic leukemia cells. Both of these ATP-independent enzymes are actually found to be inhibited by ATP at physiologically significant concentrations. Other adenine nucleotides showed varying effects: ADP inhibited only at high concentrations; AMP had no effect on either topoisomerase. Both enzymes were also inhibited by dATP. The importance of the adenine ring structure was confirmed by the lack of an inhibitory effect observed with equivalent levels of GTP, UTP, CTP, or their deoxy counterparts. Assays performed in the presence of nonhydrolyzable analogs of ATP suggest that hydrolysis of ATP does not accompany this enzyme inhibition. This was supported by direct determination of the ATPase activity of the purified enzymes. Type I topoisomerase from calf thymus and HeLa cells were also found to be sensitive to ATP. These results suggest that mammalian type I topoisomerases in general may possess a nucleotide-binding site that may be involved in regulation of enzyme activity.

Mitochondrial DNA and nuclear DNA from normal rat liver have a common sequence

Although Pst I does not cut the circular mitochondrial genome of the rat, BamHI generates from this genome two unequal fragments of DNA. Each of these fragments was cloned in pBR322. Nuclear DNA was digested from rat liver singly or doubly with Pst I and BamHI, and it was demonstrated that nuclear DNA shared a common sequence with the larger mitochondrial DNA BamHI fragment. The cloned larger mitochondrial DNA fragment was further subdivided with HindIII into four pieces that were labeled and then used to probe the double-digested nuclear DNA. The hybridization data showed that the common sequence is less than 3 kilobase pairs long and lies within the part of the mitochondrial genome containing the D-loop and a portion of the rRNA genes. It therefore appears that, as in lower eukaryotes, there are shared sequences between the nuclear and mitochondrial genomes in mammals.

Rat mitochondrial DNA polymorphism: sequence analysis of a hypervariable site for insertions/deletions

Three closely related variants of rat (Rattus norvegicus) mtDNA have been shown to differ in the number of T residues found in a run of Ts (light strand) which spans the junction between the tRNACys and tRNATyr genes. The number of Ts in the repeat varies from 6 to 8 in these DNAs. Another, less closely related, R. norvegicus variant has a run of 5 Ts at this site and in the related species, Rattus rattus, a run of 4 Ts is found. In R. norvegicus mtDNA runs of 5 As and 5 Gs are found just to the 3' side of the variable T repeat, and it is suggested that the three runs of repeated nucleotides may stabilize heteroduplexes which result from strand slippage and which give rise to the insertions and/or deletions. Among 17 mtDNA clones derived from an individual with the 8T repeat, one clone was found which possessed a 9T repeat. This variant may represent an additional DNA type originally present within the individual.

Eukaryotic DNA methylation

Eukaryotic genomes contain 5-methylcytosine (5mC) as a rare base.5mC arises by postsynthetic modification of cytosine and occurs, at least in animals, predominantly in the dinucleotide CpG. The base is not distributed randomly in these genomes but conforms to a pattern. This pattern varies between taxa but appears to be inherited in a semi-conservative fashion. At the level of the genome, gross changes in the level of DNA methylation have been noted. This has encouraged speculation that the modification may play a role in cellular differentiation. Tissue-specific patterns of DNA methylation, predicted by various models of differentiation, have been found for most vertebrate genes so far examined. A correlation has emerged between the undermethylation of these regions and their transcription, but this is not always the case. While data for eukaryotic viral sequences are less equivocal, studies of this kind cannot in isolation distinguish between undermethylation being a cause or a consequence of gene activity. If it were a cause, it is probable that the demethylation of specific CpG sites would be a necessary yet not a sufficient condition for transcription to occur. The introduction of artificially methylated DNA sequences into individual eukaryotic cells by microinjection or transformation may provide the means to elucidate these questions in the future. In the meantime, the study of eukaryotic DNA methylation promises to contribute much to our understanding of the regulation of gene expression in these organisms.

Radiation of human mitochondria DNA types analyzed by restriction endonuclease cleavage patterns

Human mitochondrial DNA (mtDNA) restriction endonuclease fragment patterns were analyzed using total blood cell DNA isolated from 200 individuals representing five different populations. Thirty-two fragment patterns (morphs) were observed with the enzymes Hpa I, Bam HI, Hae II, Msp I and Ava II yielding thirty-five different combinations of fragment patterns (mt DNA types). The major ethnic groups exhibit quantitative as well as qualitative differences in their mtDNA types, all of which are related to each other by a tree in which the closely related mtDNA types cluster according to geographic origin. Three mtDNA types are postulated to be 'central' to ethnic radiations due to their high frequencies, their appearance in more than one ethnic group, or their presence in other primate species. Genetic distances among populations were computed and employed in construction of an average linkage tree. If one of the three central mtDNA types is the root of the tree, differences in evolutionary rates among the branches become apparent. In particular, the Bushmen appear to have a higher evolutionary rate for mtDNA than the other four populations. Comparisons with nuclear gene frequencies suggest that this higher evolutionary rate may be the product of an elevated mutation rate or fixation of mutations in mtDNA.

Amino acid change associated with the major polymorphic Hinc II site of Oriental and Caucasian mitochondrial DNAs

The mitochondrial DNAs (mtDNAs) from 116 Oriental and Caucasian blood samples were analyzed for their Hinc II restriction endonuclease cleavage patterns using Southern analysis and 32P human mtDNA probes. Seven distinct patterns were found, all of which could be interrelated by single nucleotide changes. The predominant pattern (mtHinc II-2) was found in 97% of the Caucasians and 73% of the Orientals. This mtDNA morph had one more Hinc II site than did the second most common morph (mtHinc II-1), which was found only in 20% of the Orientals. Three additional patterns were in a single Oriental sample, a fourth in a single Caucasian sample, and a fifth in one member of each population. The polymorphic site that differentiated mtHinc II-1 and mtHinc II-2 was cloned and sequenced. A single nucleotide change was found that created an Hinc II site and changed the amino acid sequence of the URF5 gene. Comparison of these sequences with those of other primates [15] revealed that the Asian mtHinc II-1 and mtHinc II-4 mtDNAs were identical in this region with those of chimpanzees and orangutans. These results suggest that the Asian mtHinc II-1 mtDNA may have been ancestral to other human mtDNAs.

Mitochondrial DNA sequences of primates: tempo and mode of evolution

We cloned and sequenced a segment of mitochondrial DNA from human, chimpanzee, gorilla, orangutan, and gibbon. This segment is 896 bp in length, contains the genes for three transfer RNAs and parts of two proteins, and is homologous in all 5 primates. The 5 sequences differ from one another by base substitutions at 283 positions and by a deletion of one base pair. The sequence differences range from 9 to 19% among species, in agreement with estimates from cleavage map comparisons, thus confirming that the rate of mtDNA evolution in primates is 5 to 10 times higher than in nuclear DNA. The most striking new finding to emerge from these comparisons is that transitions greatly outnumber transversions. Ninety-two percent of the differences among the most closely related species (human, chimpanzee, and gorilla) are transitions. For pairs of species with longer divergence times, the observed percentage of transitions falls until, in the case of comparisons between primates and non-primates, it reaches a value of 45. The time dependence is probably due to obliteration of the record of transitions by multiple substitutions at the same nucleotide site. This finding illustrates the importance of choosing closely related species for analysis of evolutionary process. The remarkable bias toward transitions in mtDNA evolution necessitates the revision of equations that correct for multiple substitutions at the same site. With revised equations, we calculated the incidence of silent and replacement substitutions in the two protein-coding genes. The silent substitution rate is 4 to 6 times higher than the replacement rate, indicating strong functional constraints at replacement sites. Moreover, the silent rate for these two genes is about 10% per million years, a value 10 times higher than the silent rate for the nuclear genes studied so far. In addition, the mean substitution rate in the three mitochondrial tRNA genes is at least 100 times higher than in nuclear tRNA genes. Finally, genealogical analysis of the sequence differences supports the view that the human lineage branched off only slightly before the gorilla and chimpanzee lineages diverged and strengthens the hypothesis that humans are more related to gorillas and chimpanzees than is the orangutan.

Novel features of animal mtDNA evolution as shown by sequences of two rat cytochrome oxidase subunit II genes

The sequence of the region of the mitochondrial genome that encodes cytochrome oxidase subunit II (COII) has been determined for each of two closely related rat species, Rattus norvegicus and R. rattus. Comparison of the two sequences shows that 94.4% of the nucleotide substitutions are silent. The occurrence of this high proportion of silent substitutions leads us to propose that the rapid evolution of mtDNA relative to nuclear DNA is due only to silent changes and that amino acid-altering substitutions accumulate in nuclear and mtDNA at comparable rates. Other novel features of the nucleotide substitution pattern in the rat COII gene are a high transition/transversion ratio (8.0:1) and a strong bias toward C in equilibrium T transitions in the light strand. Comparison of the R. norvegicus COII sequence with the bovine and human sequences show that there may be selective constraints on some silent positions within the gene and that its rate of evolution may be different in different mammalian lineages.

Nicotiana chloroplast genome III. Chloroplast DNA evolution

Nicotiana chloroplast genomes exhibit a high degree of diversity and a general similarity as revealed by restriction enzyme analysis. This property can be measured accurately by restriction enzymes which generate over 20 fragments. However, the restriction enzymes which generate a small number (about 10) of fragments are extremely useful not only in constructing the restriction maps but also in establishing the sequence of ct-DNA evolution. By using a single enzyme, Sma I, a elimination and sequential gain of its recognition sites during the course of ct-DNA evolution is clearly demonstrated. Thus, a sequence of ct-DNA evolution for many Nicotiana species is formulated. The observed changes are all clustered in one region to form a "hot spot" in the circular molecule of ct-DNA. The mechanisms involved for such alterations are mostly point mutations but inversion and deficiency are also indicated. Since there is a close correlation between the ct-DNA evolution and speciation in Nicotiana a high degree of cooperation and coordination betwen organellar and nuclear genomes is evident.

Polymorphisms of mitochondrial DNAs in Norway rats (Rattus norvegicus): cleavage site variations and length polymorphism of restriction fragments

Extensive polymorphism was found in mitochondrial DNAs (mtDNAs) of Norway rats (Rattus norvegicus). The restriction endonuclease cleavage patterns of mtDNAs of laboratory rats, wild rats, tumor cells, and culture cells were compared. The polymorphism is defined by two criteria; one is cleavage site variation and the other is length polymorphism of restriction fragments. The cleavage site variation may be caused by point mutation, and the length polymorphism by sequence deletions or insertions. At least five types, types A-E, were identified by cleavage site variations, and two groups, a and b, were identified by length polymorphism of one HpaII fragment, Hpa5. All types except type C belonged to either group-a or group-b, whereas both groups were found in type C. Differentiation of polymorphic Norway rat mtDNA types and the experimental use of the polymorphism are discussed.

Ethnic variation in Hpa 1 endonuclease cleavage patterns of human mitochondrial DNA

The mtDNAs of 235 individuals from five ethnic groups were analyzed for restriction site variation by digestion with restriction endonuclease Hpa I, Southern transfer, and hybridization with 32P-labeled human mtDNA. Six different cleavage patterns (morphs) were found, all of which could be related to each other by single nucleotide substitutions. Differences were found in the frequency of these morphs among the populations. The largest difference observed was in the frequency of the morph most common in Caucasians and Orientals compared to the frequency of that found in Africans. This difference apparently originated by the sequence change G-T-C-A-A-C to G-T-T-A-A-C. This alteration permitted recognition by Hpa I but did not alter the amino acid sequence. Two other observed differences were due to separate substitutions occurring in the ribosomal RNA genes. Comparison with primate data shows that the morph with two fragments, found in 12.5% of Oriental and 4% of Bantu samples, might be the ancestral type common to all hominoids. These two conserved sites were localized in tRNA genes in the anticodon loop. Assuming that the two-fragment morph is ancestral, this finding is consistent with previous data suggesting that Asia is genetically central to the radiations that are thought to have given rise to the human ethnic groups.

Intra- and interspecific variation of the mitochondrial genome in Rattus norvegicus and Rattus rattus: restriction enzyme analysis of variant mitochondrial DNA molecules and their evolutionary relationships

Restriction endonuclease analysis has revealed extensive mtDNA polymorphism in two species of rats, Rattus rattus and Rattus norvegicus. Sequence divergence values for the eight detected R. norvegicus variants range from 0.2% to 1.8% and for the eight R. rattus variants, from 0.2% to 9.6%. Three of the most closely related R. norvegicus mtDNA's appear to differ by deletions/insertions of about 4 base pairs apiece. Restriction sites for seven enzymes have been mapped for 11 of these variants. The 31 intraspecific and 41 interspecific variant sites appear to be evenly distributed on the mtDNA molecule outside of the rRNA cistrons. The location of sites present in all the DNAs suggests that the rRNA genes and possibly the light strand origin of replication may be more highly evolutionarily conserved than other parts of the molecule. The sequence divergences among the mtDNAs of animals whose geographic origins are separated by major barriers, such as oceans, were significantly greater than those among animals found within large land masses, such as the continental United States. Dendrograms (phenograms), which have been constructed to depict the relationships among the various DNAs, indicate that East Asian members of the R. rattus species are more closely related to American rats of this species than to other Asian R. rattus animals from Sri Lanka. Moreover, it appears that R. norvegicus comprises a group taxonomically distinct from any of the R. rattus subspecies.