The complete mitochondrial sequence of Tarsius bancanus: evidence for an extensive nucleotide compositional plasticity of primate mitochondrial DNA

Inconsistencies between phylogenetic interpretations obtained from independent sources of molecular data occasionally hamper the recovery of the true evolutionary history of certain taxa. One prominent example concerns the primate infraordinal relationships. Phylogenetic analyses based on nuclear DNA sequences traditionally represent Tarsius as a sister group to anthropoids. In contrast, mitochondrial DNA (mtDNA) data only marginally support this affiliation or even exclude Tarsius from primates. Two possible scenarios might cause this conflict: a period of adaptive molecular evolution or a shift in the nucleotide composition of higher primate mtDNAs through directional mutation pressure. To test these options, the entire mt genome of Tarsius bancanus was sequenced and compared with mtDNA of representatives of all major primate groups and mammals. Phylogenetic reconstructions at both the amino acid (AA) and DNA level of the protein-coding genes led to faulty tree topologies depending on the algorithms used for reconstruction. We propose that these artifactual affiliations rather reflect the nucleotide compositional similarity than phylogenetic relatedness and favor the directional mutation pressure hypothesis because: (1) the overall nucleotide composition changes dramatically on the lineage leading to higher primates at both silent and nonsilent sites, and (2) a highly significant correlation exists between codon usage and the nucleotide composition at the third, silent codon position. Comparisons of mt genes with mt pseudogenes that presumably transferred to the nucleus before the directional mutation pressure took place indicate that the ancestral DNA composition is retained in the relatively fossilized mtDNA-like sequences, and that the directed acceleration of the substitution rate in higher primates is restricted to mtDNA.     

Nucleotide sequence of a Xenopus laevis mitochondrial DNA fragment containing the D-loop, flanking tRNA genes and the apocytochrome b gene

Extensive corrections of the nucleotide sequence of the Xenopus laevis mitochondrial (mt) displacement (D) loop and surrounding genes [Wong et al., Nucl. Acids Res. 11 (1983) 4977-4995] are reported, including addition of two stretches of nucleotides and 60 scattered modifications. The additional sequences presented here correspond to the apocytochrome b gene, the tRNAGlu gene and part of URF6. This allows us to propose a conformational model for the X. laevis apocytochrome b protein and also permits comparisons with mammalian mtDNA. The D-loop sequence is poorly conserved except for sequences involved in the regulation of the mt genome (conserved sequence blocks and the DNA polymerase stop sequences). On the other hand, all genes show marked conservation both of their nucleotide sequence and their respective location on the mt genome. Organization of the genetic information described for mammalian mtDNA also holds for the X. laevis mtDNA. This result strongly suggests that all animal vertebrate mtDNAs have followed the same evolutionary pathway.     

Maternal inheritance of mitochondrial DNA polymorphisms in cultured human fibroblasts

We have isolated the total cellular DNA from the cultured diploid fibroblasts of a six-member, three-generation human family. Using a specific radioactive probe for mitochondrial (mt) sequences we have identified new polymorphic variants in this family for the Hhal restriction endonuclease cleavage pattern of the mtDNA. The inheritance of these cleavage patterns verifies the maternal inheritance of mtDNA through all three generations.     

Nucleotide sequence identity of mitochondrial DNA from different human tissues

Recombinant DNA techniques have been used to search for mitochondrial (mt) nucleotide (nt) sequence differences between human tissues within an individual. mtDNA isolated from brain, heart, liver, kidney, and skeletal muscle of two different individuals was cleaved with SacI and XbaI, and then cloned in bacteriophage M13. Partial nt sequence determination of 121 independently isolated recombinant M13 clones containing either the cytochrome oxidase subunit III gene or the D-loop region of human mtDNA revealed base substitution differences between individuals, and between each individual and the published human mtDNA sequence. A majority of these base substitutions were transitions. No systematic nt sequence differences were identified between tissues within an individual, however. These results suggest that mtDNA sequence alterations do not accompany organogenesis, and that somatic mutations do not accumulate in the mtDNA of different human tissues to a level of greater than one nt substitution per molecule.     

Rates of nuclear and cytoplasmic mitochondrial DNA sequence divergence in mammals

Differential rates of nucleotide substitution among different gene segments and between distinct evolutionary lineages is well documented among mitochondrial genes and is likely a consequence of locus-specific selective constraints that delimit mutational divergence over evolutionary time. We compared sequence variation of 18 homologous loci (15 coding genes and 3 parts of the control region) among 10 mammalian mitochondrial DNA genomes which allowed us to describe different mitochondrial evolutionary patterns and to produce an estimation of the relative order of gene divergence. The relative rates of divergence of mitochondrial DNA genes in the family Felidae were estimated by comparing their divergence from homologous counterpart genes included in nuclear mitochondrial DNA (Numt, pronounced "new might"), a genomic fossil that represents an ancient transfer of 7.9 kb of mitochondrial DNA to the nuclear genome of an ancestral species of the domestic cat (Felis catus). Phylogenetic analyses of mitochondrial (mtDNA) sequences with multiple outgroup species were conducted to date the ancestral node common to the Numt and the cytoplasmic (Cymt) mtDNA genes and to calibrate the rate of sequence divergence of mitochondrial genes relative to nuclear homologous counterparts. By setting the fastest substitution rate as strictly mutational, an empirical "selective retardation index" is computed to quantify the sum of all constraints, selective and otherwise, that limit sequence divergence of mitochondrial gene sequences over time.      

A Method for Calibrating Molecular Clocks and Its Application to Animal Mitochondrial DNA

A generalized least-squares procedure is introduced for the calibration of molecular clocks and applied to the complete mitochondrial DNA sequences of 13 animal species. The proposed technique accounts for both nonindependence and heteroscedasticity of molecular-distance data, problems that have not been taken into to account in such analyses in the past. When sequence-identity data are transformed to account for multiple substitutions/site, the molecular divergence scales linearly with time, but with substantially more variation in the substitution rate than expected under a Poisson model. Significant levels of divergence are predicted at zero divergence time for most loci, suggesting high levels of site-specific heterozygosity among mtDNA molecules establishing in sister taxa. For nearly all loci, the baseline heterozygosity is lower and the substitution rate is higher in mammals relative to other animals. There is considerable variation in the evolutionary rate among loci but no compelling evidence that the average rate of mtDNA evolution is elevated with respect to that of nuclear DNA. Using the observed patterns of interspecific divergence, empirical estimates are derived for the mean coalescence times of organelles colonizing sister taxa.     

Quantifying the Elevation of Mitochondrial DNA Evolutionary Substitution Rates Over Nuclear Rates in the Intertidal Copepod Tigriopus californicus

Mitochondrial DNA (mtDNA) genomes generally evolve rapidly in animals, but considerable variation in the rates of evolution of mtDNA occurs among taxa. Higher levels of mutation will tend to increase the amount of polymorphism, which should also scale with population size, but there are mixed signals from previous studies on the evolutionary outcomes of the interactions of these processes. The copepod Tigriopus californicus provides an interesting model in which to study the evolution of mtDNA because it has high levels of divergence among populations and there is the suggestion that this divergence could be involved in reproductive isolation. This species also appears to have an elevated mtDNA substitution rate, but previous studies did not provide an accurate measurement. This article examines the rate of mtDNA substitution versus nuclear substitution in T. californicus and finds that the mtDNA rate for synonymous sites averages 55-fold higher, a level that exceeds the rates found in most other invertebrates. Levels of polymorphism are also examined in both mtDNA and nuclear genes, and it is shown that the effective population size of mtDNA genes is much lower than that of nuclear genes. In addition, no correlation between polymorphism in mtDNA and nuclear genes is found across populations, which suggest factors other than demography may shape polymorphism in this species. The results from this study suggest that mtDNA is evolving at a very rapid rate in this copepod species, and this could increase the likelihood that mtDNA evolution is involved in the generation of reproductive isolation.     

Relative rates of evolution in the coding and control regions of African mtDNAs

Reduced median networks of African haplogroup L mitochondrial DNA (mtDNA) sequences were analyzed to determine the pattern of substitutions in both the noncoding control and coding regions. In particular, we attempted to determine the causes of the previously reported (Howell et al. 2004) violation of the molecular clock during the evolution of these sequences. In the coding region, there was a significantly higher rate of substitution at synonymous sites than at nonsynonymous sites as well as in the tRNA and rRNA genes. This is further evidence for the operation of purifying selection during human mtDNA evolution. For most sites in the control region, the relative rate of substitution was similar to the rate of neutral evolution (assumed to be most closely approximated by the substitution rate at 4-fold degenerate sites). However, there are a number of mutational hot spots in the control region, approximately 3% of the total sites, that have a rate of substitution greater than the neutral rate, at some sites by more than an order of magnitude. It is possible either that these sites are evolving under conditions of positive selection or that the substitution rate at some sites in the control region is strongly dependent upon sequence context. Finally, we obtained preliminary evidence for "nonideal" evolution in the control region, including haplogroup-specific substitution patterns and a decoupling between relative rates of substitution in the control and coding regions.     

Man's place in hominoidea revealed by mitochondrial DNA genealogy

Summary Molecular biology has resurrected C. Darwin and T.H. Huxley's question about the origin of humans, but the precise branching pattern and dating remain controversial. To settle this issue, a large amount of sequence information is required. We determined mitochondrial (mt) DNA sequences for five hominoids; pygmy and common chimpanzees, gorilla, orangutan, and siamang. The common region compared with the known human sequence is 4759 by long, encompassing genes for 11 transfer RNAs and 6 proteins. Because of the high substitution rates in mammalian mtDNA and an unprecedentedly large region compared, the sequence differences clearly indicate that the closest relatives to human are chimpanzees rather than gorilla. For dating the divergences of human, chimpanzee, and gorilla, we used only unsaturated parts of sequence differences in which the mtDNA genealogy is not obscured by multiple substitutions. The result suggests that gorilla branched off 7.7 ± 0.7 million years (Myr) ago and human 4.7 ± 0.5 Myr ago; the time difference between these divergences being as long as 3 Myr.Offprint requests to: S. Horai     

Mitochondrial genome nucleotide substitution pattern between domesticated silkmoth, Bombyx mori, and its wild ancestors, Chinese Bombyx mandarina and Japanese Bombyx mandarina

Bombyx mori and Bombyx mandarina are morphologically and physiologically similar. In this study, we compared the nucleotide variations in the complete mitochondrial (mt) genomes between the domesticated silkmoth, B. mori, and its wild ancestors, Chinese B. mandarina (ChBm) and Japanese B. mandarina (JaBm). The sequence divergence and transition mutation ratio between B. mori and ChBm are significantly smaller than those observed between B. mori and JaBm. The preference of transition by DNA strands between B. mori and ChBm is consistent with that between B. mori and JaBm, however, the regional variation in nucleotide substitution rate shows a different feature. These results suggest that the ChBm mt genome is not undergoing the same evolutionary process as JaBm, providing evidence for selection on mtDNA. Moreover, investigation of the nucleotide sequence divergence in the A+T-rich region of Bombyx mt genomes also provides evidence for the assumption that the A+T-rich region might not be the fastest evolving region of the mtDNA of insects.     

Extrachromosomal DNA in the Apicomplexa.

Malaria and related apicomplexan parasites have two highly conserved organellar genomes: one is of plastid (pl) origin, and the other is mitochondrial (mt). The organization of both organellar DNA molecules from the human malaria parasite Plasmodium falciparum has been determined, and they have been shown to be tightly packed with genes. The 35-kb circular DNA is the smallest known vestigial plastid genome and is presumed to be functional. All but two of its recognized genes are involved with genetic expression: one of the two encodes a member of the clp family of molecular chaperones, and the other encodes a conserved protein of unknown function found both in algal plastids and in eubacterial genomes. The possible evolutionary source and intracellular location of the plDNA are discussed. The 6-kb tandemly repeated mt genome is the smallest known and codes for only three proteins (cytochrome b and two subunits of cytochrome oxidase) as well as two bizarrely fragmented rRNAs. The organization of the mt genome differs somewhat among genera. The mtDNA sequence provides information not otherwise available about the structure of apicomplexan cytochrome b as well as the unusually fragmented rRNAs. The malarial mtDNA has a phage-like replication mechanism and undergoes extensive recombination like the mtDNA of some other lower eukaryotes.     

Identification of three populations ofOphiostoma ulmi (aggressive subgroup) by mitochondrial DNA restriction-site mapping and nuclear DNA fingerprinting

Three genetically distinct populations of the Dutch elm pathogenOphistoma ulmi within the aggressive subgroup were defined by the hybridization of a human minisatellite DNA sequence (HVR 33.6) to polymorphic dispersed DNA sequences within theO. ulmi nuclear genomes. For the 10 isolates examined there was a close correlation between nuclear DNA fingerprints and mitochondrial (mt) DNA restriction patterns. A restriction-site map was constructed for the mitochondrial genomes for each of these populations. The three mt DNA maps corresponded to genome sizes of 49.1 (Type I), 49.9 (Type II), and 53.9 (Type III) kilobase pairs (kbp) of DNA. The Type I and Type II mt genomes differed from the Type III mt genome by discrete length mutations of 4.8 and 4.0 kbp, respectively. It is unknown whether these length mutations resulted from insertions into or deletions from a progenitor mitochondrial genome. There was no correlation between the mitochondrial or nuclear genotypes and the geographical source of the isolates.     

The Use of Hybrid Somatic Cells as an Approach to Mitochondrial Genetics in Animals

We have studied the fate of parental mitochondrial DNA (mtDNA) in hybrid somatic cells derived by Sendai virus-induced fusion of human cells and mouse or rat cells. Many hybrid cell strains were obtained which contained sequences from both human and rodent mtDNA after 40 to 60 population doublings. Some strains were subcloned and cultured further for up to 150 doublings; a large fraction of these strains contained both parental mtDNA sequences at that time.The relation between human and rodent mtDNA sequences was tested in some of the hybrid cell strains. In a high fraction of strains tested the human and rodent mtDNA sequences were linked to each other by what are most likely covalent bonds. This linkage may be described as "recombination" of mtDNA sequences from two different animals.     

Evolution of duplicate control regions in the mitochondrial genomes of metazoa: a case study with Australasian Ixodes ticks

To investigate the evolution pattern and phylogenetic utility of duplicate control regions (CRs) in mitochondrial (mt) genomes, we sequenced the entire mt genomes of three Ixodes species and part of the mt genomes of another 11 species. All the species from the Australasian lineage have duplicate CRs, whereas the other species have one CR. Sequence analyses indicate that the two CRs of the Australasian Ixodes ticks have evolved in concert in each species. In addition to the Australasian Ixodes ticks, species from seven other lineages of metazoa also have mt genomes with duplicate CRs. Accumulated mtDNA sequence data from these metazoans and two recent experiments on replication of mt genomes in human cell lines with duplicate CRs allowed us to re-examine four intriguing questions about the presence of duplicate CRs in the mt genomes of metazoa: (1) Why do some mt genomes, but not others, have duplicate CRs? (2) How did mt genomes with duplicate CRs evolve? (3) How could the nucleotide sequences of duplicate CRs remain identical or very similar over evolutionary time? (4) Are duplicate CRs phylogenetic markers? It appears that mt genomes with duplicate CRs have a selective advantage in replication over mt genomes with one CR. Tandem duplication followed by deletion of genes is the most plausible mechanism for the generation of mt genomes with duplicate CRs. Once duplicate CRs occur in an mt genome, they tend to evolve in concert, probably by gene conversion. However, there are lineages where gene conversion may not always occur, and, thus, the two CRs may evolve independently in these lineages. Duplicate CRs have much potential as phylogenetic markers at low taxonomic levels, such as within genera, within families, or among families, but not at high taxonomic levels, such as among orders.     

Evolutionary distance estimation under heterogeneous substitution pattern among lineages

Most of the sophisticated methods to estimate evolutionary divergence between DNA sequences assume that the two sequences have evolved with the same pattern of nucleotide substitution after their divergence from their most recent common ancestor (homogeneity assumption). If this assumption is violated, the evolutionary distance estimated will be biased, which may result in biased estimates of divergence times and substitution rates, and may lead to erroneous branching patterns in the inferred phylogenies. Here we present a simple modification for existing distance estimation methods to relax the assumption of the substitution pattern homogeneity among lineages when analyzing DNA and protein sequences. Results from computer simulations and empirical data analyses for human and mouse genes are presented to demonstrate that the proposed modification reduces the estimation bias considerably and that the modified method performs much better than the LogDet methods, which do not require the homogeneity assumption in estimating the number of substitutions per site. We also discuss the relationship of the substitution and mutation rate estimates when the substitution pattern is not the same in the lineages leading to the two sequences compared.     

Evidence for maternal inheritance of mitochondrial DNA in allotetraploid.

The complete mitochondrial DNA (mtDNA) sequences of the allotetraploid and red crucian carp were determined in this paper. We compared the complete mtDNA sequences between the allotetraploid and its female parent red crucian carp, and between the allotetraploid and its male parent common carp. The results indicated that the complete mtDNA nucleotide identity (99.7%) between the allotetraploid and its female parent red crucian carp was higher than that (89.0%) between the allotetraploid and its male parent common carp. Moreover, the analysis on the start and stop codons, overlaps and spacers, and phylogeny of the mt genomes indicated the genetic relationship between the allotetraploid and its female parent red crucian carp was closer than that between the allotetraploid and its male parent common carp. Our results indicated that the allotetraploid mt genome was strictly maternally inherited. Through maternal inheritance, the mt genome in the F(11) allotetraploid displayed extremely high similarity to that in the female parent red crucian carp after 11 generations (from F(1) to F(11) hybrids). Such results indicated that the F(11) allotetraploid possessed the stable inheritance characteristic. Thus the tetraploid stocks possessed the good base to form a new tetraploid species in the future. Since the establishment of the new tetraploid stocks has the great significance in analyzing evolutionary theory of vertebrate and in improving aquaculture industry, analysis of the mt genome and the elucidation of the variation of the mt genome in the allotetraploid and its parents proved that it was a useful genetic marker to monitor the variations in the progeny of the crosses.     

Protein bound to mitochondrial DNA from tumor and normal tissues of humans and animals

Stable mitochondrial (mt) DNA-protein complexes have been reported repeatedly in the last decade. We have found that the amount of proteins bound to mt DNA was increased in human HL-60 promyelocytic and chronic myelocytic leukemia cells. Mt of human tumor cells and bovine and rate liver cells were isolated by differential centrifugation. The DNA was purified by SDS-NaCl precipitation of protein, alcohol precipitation, and sucrose density gradient centrifugation. Mt DNA-bound protein obtained by enzymatic digestion of the DNA and purified by gel exclusion and reverse-phase HPLC chomatographies and SDS-polyacrylamide gel electrophoresis, showed a major protein band at 70 kD and a minor band at 35 kD. Hydrolysis of the mtDNA-protein complex in formic acid yielded DNA bases and peptides by HPLC on a C18 reverse-phase column. Tumor cell mtDNA contained 5-10 times more bound protein than mtDNA from normal tissue.     

Evidence for local DNA influences on patterns of substitutions in the human alpha-interferon gene family

The evolutionary history of genes can be used to examine patterns of spontaneous mutation if the sequences are sufficiently extensive to provide reliable data. Many human alpha-interferon genes have been sequenced and they form a large multigene family including several pseudogenes. A phylogenetic history for 15 human interferon sequences was reconstructed and their ancestral sequences inferred using a maximum parsimony method. This evolutionary history provided a record of more than 738 spontaneous mutations that have occurred in man's recent evolution. Of these mutations, more than 267 base substitution and deletion-insertion events were analyzed to determine the possible effects of nearby DNA sequences. Many substitutions occur at the end of long runs of identical bases and some dinucleotide pairs may mutate more often than others. Because templating by local DNA sequences has been implicated in prokaryotic mutation, the sequences were also examined for nearby repeats that include the substituted nucleotide and hence are potentially capable of templating the substitution. The majority of sequence alterations examined have either a similar direct repeat or palindrome nearby. Often such templates can account for simultaneous multiple mutations. These results suggest that sequence-directed events may occur occasionally in eukaryotes and that neighbouring DNA sequences can influence both the occurrence and types of mutations in several different ways.     

DNA Polymorphism Detectable by Restriction Endonucleases

Data on DNA polymorphisms detected by restriction endonucleases are rapidly accumulating. With the aim of analyzing these data, several different measures of nucleon (DNA segment) diversity within and between populations are proposed, and statistical methods for estimating these quantities are developed. These statistical methods are applicable to both nuclear and nonnuclear DNAs. When evolutionary change of nucleons occurs mainly by mutation and genetic drift, all the measures can be expressed in terms of the product of mutation rate per nucleon and effective population size. A method for estimating nucleotide diversity from nucleon diversity is also presented under certain assumptions. It is shown that DNA divergence between two populations can be studied either by the average number of restriction site differences or by the average number of nucleotide differences. In either case, a large number of different restriction enzymes should be used for studying phylogenetic relationships among related organisms, since the effect of stochastic factors on these quantities is very large. The statistical methods developed have been applied to data of Shah and Langley on mitochondrial (mt)DNA from Drosophila melanogaster, simulans and virilis. This application has suggested that the evolutionary change of mtDNA in higher animals occurs mainly by nucleotide substitution rather than by deletion and insertion. The evolutionary distances among the three species have also been estimated.