A new estimate of sequence divergence of mitochondrial DNA using restriction endonuclease mappings

Summary A new estimate of the sequence divergence of mitochondrial DNA in related species using restriction enzyme maps is constructed. The estimate is derived assuming a simple Posisson-like model for the evolutionary process and is chosen to maximize an expression which is a reasonable approximation to the true likelihood of the restriction map data. Using this estimate, four sets of mitochondrial DNA data are analyzed and discussed.     

Complementary DNA sequence of a human cytoplasmic actin. Interspecies divergence of 3' non-coding regions

We have isolated and sequenced a cloned complementary DNA insert complementary to the messenger RNA of a cytoplasmic actin expressed in human epidermal cells. This provides the first cytoplasmic actin complementary DNA sequence for a vertebrate organism. The actin amino acid sequence predicted from this complementary DNA is identical to that of a bovine cytoplasmic actin and shows 98 and 85% homology with a Dictyostelium and a yeast actin, respectively. The complementary DNA sequence indicates that the 3' end of the mRNA contains an unusually long (greater than 400 nucleotides) 3' non-translated region. A comparison of this 3' non-coding region with those of recently determined actin complementary DNA sequences from other species reveals little or no homology among these sequences. Thus, these results indicate that although the actin amino acid sequences are extremely conserved, the non-coding regions of the mRNAs diverge rapidly.     

Nucleotide sequence divergence in the A+T-rich region of mitochondrial DNA in Drosophila simulans and Drosophila mauritiana

We have determined the nucleotide sequences of two regions within the A+T-rich region of mitochondrial DNA (mtDNA) in the siIII type of Drosophila simulans and the maI type of D. mauritiana. The sequences of the two regions in siIII and maI are almost identical. The sequences include elements corresponding to the type I and type II repeats elements and the T-stretches as reported in D. melanogaster; an approximately 340-bp region (A region) adjacent to the tRNA(Ile) gene includes a part of the type II repeat element, and an approximately 440- bp region (B region) includes a central portion of the A+T-rich region between the type I and type II repeat arrays. Each sequence of the two species was compared with those of D. melanogaster and D. yakuba. The sequences of the A region are relatively well conserved among the four species. The alignment of the two sequences of the B region with those of D. melanogaster and D. yakuba requires numerous insertions/deletions. For both regions, nucleotide differences between D. simulans or D. mauritiana and D. melanogaster are similar to those between the two and D. yakuba. The tendency is obvious in a subregion within the type II repeat element in the A region. These findings suggest that the rate of nucleotide substitution in the subregion is accelerated in the lineage leading to D. melanogaster. Loss of functional constraint in the stem-loop-forming sequence is proposed for this acceleration.      

DNA sequence divergence in the Drosophila virilis group

DNA sequence divergence was analyzed in some sibling species of the Drosophila virilis group. Clones comprising about 0.1% of the genome DNA were selected at random from a D. virilis library for a comparative study on DNA from D. lummei, D. novamexicana, D. borealis, and D. lacicola. Blot hybridization experiments indicated that about 70% of DNA from D. lummei and D. novamexicana and less than 50% of DNA from D. borealis and D. lacicola share sequences that are homologous to DNA in D. virilis. This finding is in excellent agreement with the genealogical tree based on cytological studies (Throckmorton 1982). - Four plasmids with inserts which are present in one or a few copies per genome were hybridized in situ to polytene chromosomes. These experiments demonstrate that (1) homologous "unique" DNA sequences are localized exclusively in homologous bands and (2) homologous bands that appear to be identical in different species may contain different DNA sequences.     

Population subdivision in Europe's great bustard inferred from mitochondrial and nuclear DNA sequence variation

A continent-wide survey of sequence variation in mitochondrial (mt) and nuclear (n) DNA of the endangered great bustard (Otis tarda) was conducted to assess the extent of phylogeographic structure in a morphologically monotypic bird. DNA sequence variation in a combined 809 bp segment of the mtDNA genome from 66 individuals from the last six breeding regions showed relatively low levels of intraspecific sequence diversity (n = 0.32%) but significant differences in the regional distribution of 11 haplotypes (phiST = 0.49). Despite their exceptional potential for dispersal, a complete and long-term historical separation between the populations from the Iberian Peninsula (Spain) and mainland Europe (Hungary, Slovakia, Germany, and Russia) was demonstrated. Divergence between populations based on a 3-bp insertion-deletion polymorphism within the intron region of the nuclear CHD-Z gene was geographically concordant with the primary subdivision identified within the mtDNA sequences. Inferred aspects of phylogeography were used to formulate conservation recommendations for this endangered species.     

Evolutionary divergence of mitochondrial DNA from Paramecium aurelia

Summary Mitochondrial (mt) DNA from four sibling species within the Paramecium aurelia complex, including stocks of different geographic origin and mutants, were analyzed using four 6-bp recognition site and one 4-bp recognition site endonucleases and the sequence divergence was estimated using Upholt's (1977) statistical procedure. All four species were readily distinguishable regardless of the restriction endonuclease employed. With intraspecies comparisons, no differences were observed which could be accounted for on the basis of geographic origin. Except for species 4, each stock (and mutant) gave a species-specific fragment pattern. For species 4, while the patterns were distinct from the other species, two species-specific type of patterns were found, designated A and B. The sequence divergence between these was estimated to be between 1 and 2 percent. With interspecies comparisons, the sequence divergence ranged from 3.9 to 10.3% with the greatest divergence being between species 1 and 4, and the least between species 1 and 5. The similarity between species 1 and 5 is in accord with other criteria for interspecies comparisons. The degree of sequence divergence measured here in Paramecium mt DNA is well within the range reported for rodents and primates. All four species mt DNA were cleaved to many DNA fragments by DPN II, an enzyme which recognizes non-methylated sites, and not by DPNI, the methyl-site specific counterpart of DPN II, suggesting that mt DNA from Paramecium aurelia is not appreciably methylated, if at all.     

Oxidative DNA damage processing in nuclear and mitochondrial DNA

Living organisms are constantly exposed to oxidative stress from environmental agents and from endogenous metabolic processes. The resulting oxidative modifications occur in proteins, lipids and DNA. Since proteins and lipids are readily degraded and resynthesized, the most significant consequence of the oxidative stress is thought to be the DNA modifications, which can become permanent via the formation of mutations and other types of genomic instability. Many different DNA base changes have been seen following some form of oxidative stress, and these lesions are widely considered as instigators for the development of cancer and are also implicated in the process of aging. Several studies have documented that oxidative DNA lesions accumulate with aging, and it appears that the major site of this accumulation is mitochondrial DNA rather than nuclear DNA. The DNA repair mechanisms involved in the removal of oxidative DNA lesions are much more complex than previously considered. They involve base excision repair (BER) pathways and nucleotide excision repair (NER) pathways, and there is currently a great deal of interest in clarification of the pathways and their interactions. We have used a number of different approaches to explore the mechanism of the repair processes, to examine the repair of different types of oxidative lesions and to measure different steps of the repair processes. Furthermore, we can measure the DNA damage processing in the nuclear DNA and separately, in the mitochondrial DNA. Contrary to widely held notions, mitochondria have efficient DNA repair of oxidative DNA damage.     

Catarrhine primate divergence dates estimated from complete mitochondrial genomes: concordance with fossil and nuclear DNA evidence

Accurate divergence date estimates improve scenarios of primate evolutionary history and aid in interpretation of the natural history of disease-causing agents. While molecule-based estimates of divergence dates of taxa within the superfamily Hominoidea (apes and humans) are common in the literature, few such estimates are available for the Cercopithecoidea (Old World monkeys), the sister taxon of the hominoids in the primate infraorder Catarrhini. To help fill this gap, we have sequenced the entire mitochondrial DNA (mtDNA) genomes from a representative of three cercopithecoid tribes, Cercopithecini (Chlorocebus aethiops), Colobini (Colobus guereza), and Presbytini (Trachypithecus obscurus), and analyzed these new data together with other catarrhine mtDNA genomes available in public databases. Molecular divergence date estimates are dependent on calibration points gleaned from the paleontological record. We defined criteria for the selection of good calibration points and identified three points meeting these criteria: Homo-Pan, 6.0 Ma; Pongo-hominines, 14.0 Ma; hominoid/cercopithecoid, 23.0 Ma. Because a uniform molecular clock does not fit the catarrhine mtDNA data, we estimated divergence dates using a penalized likelihood and a Bayesian method, both of which take into account the effects of rate differences on lineages, phylogenetic tree structure, and multiple calibration points. The penalized likelihood method applied to the coding regions of the mtDNA genome yielded the following divergence date estimates, with approximate 95% confidence intervals: cercopithecine-colobine, 16.2 (14.4-17.9) Ma; colobin-presbytin, 10.9 (9.6-12.3) Ma; cercopithecin-papionin, 11.6 (10.3-12.9) Ma; and Macaca-Papio, 9.8 (8.6-10.9) Ma. Within the hominoids, the following dates were inferred: hylobatid-hominid, 16.8 (15.0-18.5) Ma; Gorilla-Homo+Pan, 8.1 (7.1-9.0) Ma; Pongo pygmaeus pygmaeus-P. p. abelii, 4.1 (3.5-4.7) Ma; and Pan troglodytes-P. paniscus, 2.4 (2.0-2.7) Ma. These dates were similar to those found using penalized likelihood on other subsets of the data, but slightly younger than several of the Bayesian estimates.     

Discrepancy in divergence of the mitochondrial and nuclear genomes ofDrosophila teissieri andDrosophila yakuba

Summary Restriction sites were compared in the mitochondrial DNA (mtDNA) molecules from representatives of two closely related species of fruit flies: nine strains ofDrosophila teissieri and eight strains ofDrosophila yakuba. Nucleotide diversities amongD. teissieri strains and amongD. yakuba strains were 0.07% and 0.03%, respectively, and the nucleotide distance between the species was 0.22%. Also determined was the nucleotide sequence of a 2305-nucleotide pari (ntp) segment of the mtDNA molecule ofD. teissieri that contains the noncoding adenine+thymine (A+T)-rich region (1091 ntp) as well as the genes for the mitochondrial small-subunit rRNA, tRNA^f-met, tRNA^gln, and tRNA^ile, and portions of the ND2 and tRNA^val genes. This sequence differs from the corresponding segment of theD. yakuba mtDNA by base substitutions at 0.1% and 0.8% of the positions in the coding and noncoding regions, respectively. The higher divergence due to base substitutions in the A+T-rich region is accompanied by a greater number of insertions/deletions than in the coding regions. From alignment of theD. teissieri A+T-rich sequence with those ofD. yakuba andDrosophila virilis, it appears that the 40% of this sequence that lies adjacent to the tRNA^ile gene has been highly conserved. Divergence between the entireD. teissieri andD. yakuba mtDNA molecules, estimated from the sequences, was 0.3%; this value is close to the value (0.22%) obtained from the restriction analysis, but 10 times lower than the value estimated from published DNA hybridization results. From consideration of the relationships of mitochondrial nucleotide distance and allozyme genetic distance found among seven species of theDrosophila melanogaster subgroup, the mitochondrial nucleotide distance observed forD. teissieri andD. yakuba is anomalously low in relation to the nuclear genetic distance.     

Reduction in the rate of DNA reassociation by sequence divergence

An estimate is made of the effect of imperfectly complementary sequences on the rate of reassociation of DNA. Rate measurements are reported for the reassociation of deaminated DNA and for the pairing of DNAs from related bacteria. A method is presented for separating the effect of the incubation temperature on the rate from the effect of sequence divergence. After correction to the optimum incubation temperature, the rate of DNA reassociation appears to be reduced by a factor of two for each 10 deg. C reduction in melting temperature due to sequence divergence. For most typical cases this effect is modest. However it can be quite important for measurements of the relation between the DNAs of different species.     

Phylogenetic systematics and evolution of primate-derived Pneumocystis based on mitochondrial or nuclear DNA sequence comparison

Previous studies have demonstrated that the agent of Pneumocystis pneumonia (PcP), Pneumocystis carinii, is actually a complex of eukaryotic organisms, and cophylogeny could explain the distribution of the hosts and parasites. In the present work, we tested the hypothesis of cophylogeny between the primate-derived Pneumocystis group and their hosts. Specific strains isolated from 20 primate species, including humans, were used to produce a phylogeny of the parasites. Aligned sequences corresponding to DNA sequences of three genes (DHPS, mtSSU-rRNA, and mtLSU-rRNA) were separately analyzed and then combined in a single data set. The resulting parasite phylogeny was compared with different controversial phylogenies for the hosts. This comparison demonstrated that, depending upon which topology is accepted for the hosts, at least 61% and perhaps 77% of the homologous nodes of the respective cladograms of the hosts and parasites may be interpreted as resulting from codivergence events. This finding and the high specificity of these parasites suggests that cophylogeny may be considered the dominant pattern of evolution for Pneumocystis organisms, representing a new example of parallel evolution between primates and their specific parasites. Because the phylogeny of Pneumocystis followed very closely the differentiation of their hosts at the species level, the study of the parasites could provide valuable information on the phylogeny of their hosts. We used this information to explore controversial hypotheses of the phylogeny of the Platyrrhini by comparison with the phylogeny of their specific Pneumocystis parasites. If these organisms were closely associated as lung parasites with primates through the ages, the hypothesis of the Pneumocystis spp. being new pathogenic agents could be refuted. However, these organisms are opportunistic symbionts, becoming pathogenic whenever the immunological defences of their hosts decline. This study also provides support for the hypothesis that the different Pneumocystis species are genetically independent organisms, helping to clarify their taxonomic status.     

Real-time DNA quantification of nuclear and mitochondrial DNA in forensic analysis

The rapid development of molecular genetic analysis tools has made it possible to analyze most biological materialfound at the scene of a crime. Evidence materials containing DNA quantities too low to be analyzed using nuclear markers can be analyzed using the highly abundant mtDNA. However, there is a shortage of sensitive nDNA and mtDNA quantification assays. In this study, an assay for the quantification of very small amounts of DNA, based on the real-time Taq-Man assay, has been developed. This analysis will provide an estimate of the total number of nDNA copies and the total number of mtDNA molecules in a particular evidence material. The quantification is easy to perform, fast, and requires a minimum of the valuable DNA extracted from the evidence materiaL The results will aid in the evaluation of whether the specific sample is suitable for nDNA or mtDNA analysis. Furthermore, the optimal amount of DNA to be used in further analysis can be estimated ensuring that the analysis is successful and that the DNA is retained for future independent analysis. This assay has significant advantages over existing techniques because of its high sensitivity, accuracy, and the combined analysis of nDNA and mtDNA. Moreover, it has the potential to provide additional information about the presence of inhibitors in forensic samples. Subsequent mitochondrial and nuclear analysis of quantified samples illustrated the potential to predict the number of DNA copies required for a successful analysis in a certain typing assay.     

Sequence divergence of the mitochondrial DNA of Pacific Ocean salmon

Restriction assay of mtDNA has been made in 6 salmon species form the genus Oncorhynchus and one species from the related genus, i.e. Salvelinus malma. The size of the mitochondrial genome was found to be identical and equal to 16.7 kilobases. The digestion patterns of mtDNA cleaved by 5 restriction endonucleases (Eco RI, Bgl I, Bgl II, Hind III, and Pst I) were used for analysis of the levels of interspecific variation and for estimation the matrix of mtDNA sequence differentiation. It was found that the level of nucleotide sequence divergence (p) in the genus Oncorhynchus varies within 1.7-6.7%. Minimum p value was observed in a pair O. keta--O. gorbuscha, maximum one--between O. masu and other species. With respect to similarity in their mtDNA, three groups may be distinguished: 1) O. gorbuscha--O. keta; 2) O. nerka--O. kisutch, O. tschawytscha; 3) O. masu. Mean value of intergeneric level of sequence divergence between Oncorhynchus and Salvelinus was found to be equal to 8%. On the basis of mtDNA analysis, the dendrogram of similarity of the species was plotted which is consistent in principle with current viewpoints on phylogenetic relations among the Pacific salmon.     

Revisiting the mouse mitochondrial DNA sequence

The existence of reliable mtDNA reference sequences for each species is of great relevance in a variety of fields, from phylogenetic and population genetics studies to pathogenetic determination of mtDNA variants in humans or in animal models of mtDNA-linked diseases. We present compelling evidence for the existence of sequencing errors on the current mouse mtDNA reference sequence. This includes the deletion of a full codon in two genes, the substitution of one amino acid on five occasions and also the involvement of tRNA and rRNA genes. The conclusions are supported by: (i) the re-sequencing of the original cell line used by Bibb and Clayton, the LA9 cell line, (ii) the sequencing of a second L-derivative clone (L929), and (iii) the comparison with 12 other mtDNA sequences from live mice, 10 of them maternally related with the mouse from which the L cells were generated. Two of the latest sequences are reported for the first time in this study (Balb/cJ and C57BL/6J). In addition, we found that both the LA9 and L929 mtDNAs also contain private clone polymorphic variants that, at least in the case of L929, promote functional impairment of the oxidative phosphorylation system. Conse quently, the mtDNA of the strain used for the mouse genome project (C57BL/6J) is proposed as the new standard for the mouse mtDNA sequence.