Complete Mitochondrial DNA Sequences of Six Snakes: Phylogenetic Relationships and Molecular Evolution of Genomic Features

Complete mitochondrial DNA (mtDNA) sequences were determined for representative species from six snake families: the acrochordid little file snake, the bold boa constrictor, the cylindrophiid red pipe snake, the viperid himehabu, the pythonid ball python, and the xenopeltid sunbeam snake. Thirteen protein-coding genes, 22 tRNA genes, 2 rRNA genes, and 2 control regions were identified in these mtDNAs. Duplication of the control region and translocation of the tRNA^Leu gene were two notable features of the snake mtDNAs. The duplicate control regions had nearly identical nucleotide sequences within species but they were divergent among species, suggesting concerted sequence evolution of the two control regions. In addition, the duplicate control regions appear to have facilitated an interchange of some flanking tRNA genes in the viperid lineage. Phylogenetic analyses were conducted using a large number of sites (9570 sites in total) derived from the complete mtDNA sequences. Our data strongly suggested a new phylogenetic relationship among the major families of snakes: ((((Viperidae, Colubridae), Acrochordidae), (((Pythonidae, Xenopeltidae), Cylindrophiidae), Boidae)), Leptotyphlopidae). This conclusion was distinct from a widely accepted view based on morphological characters in denying the sister-group relationship of boids and pythonids, as well as the basal divergence of nonmacrostomatan cylindrophiids. These results imply the significance to reconstruct the snake phylogeny with ample molecular data, such as those from complete mtDNA sequences.[Reviewing Editor: Dr. Bill Ballard]     

The primary structure of the β-subunit of protocatechuate 3,4-dioxygenase from Pseudomonas aeruginosa

The complete amino acid sequence of the β-subunit of protocatechuate 3,4-dioxygenase was determined. The β-subunit contained four methionine residues. Thus, five peptides were obtained after cleavage of the carboxymethylated β-subunit with cyanogen bromide, and were isolated on Sephadex G-75 column chromatography. The amino acid sequences of the cyanogen bromide peptides were established by characterization of the peptides obtained after digestion with trypsin, chymotrypsin, thermolysin, or Staphylococcus aureus protease. The major sequencing techniques used were automated and manual Edman degradations. The five cyanogen bromide peptides were aligned by means of the amino acid sequences of the peptides containing methionine purified from the tryptic hydrolysate of the carboxymethylated β-subunit. The amino acid sequence of all the 238 residues was as follows: ProAlaGlnAspAsnSerArgPheValIleArgAsp ArgAsnTrpHis ProLysAlaLeuThrPro-Asp — TyrLysThrSerIleAlaArg SerProArgGlnAla LeuValSerIleProGlnSer — IleSerGluThrThrGly ProAsnPheSerHisLeu GlyPheGlyAlaHisAsp-His — AspLeuLeuLeuAsnPheAsn AsnGlyGlyLeu ProIleGlyGluArgIle-Ile — ValAlaGlyArgValValAsp GlnTyrGlyLysPro ValProAsnThrLeuValGluMet — TrpGlnAlaAsnAla GlyGlyArgTyrArg HisLysAsnAspArgTyrLeuAlaPro — LeuAspProAsn PheGlyGlyValGly ArgCysLeuThrAspSerAspGlyTyrTyr — SerPheArg ThrIleLysProGlyPro TyrProTrpArgAsnGlyProAsnAsp — TrpArgProAla HisIleHisPheGlyIle SerGlyProSerIleAlaThr-Lys — LeuIleThrGlnLeuTyr PheGluGlyAspPro LeuIleProMetCysProIleVal — LysSerIleAlaAsn ProGluAlaValGlnGln LeuIleAlaLysLeuAspMetAsnAsn — AlaAsnProMet AsnCysLeuAlaTyr ArgPheAspIleValLeuArgGlyGlnArgLysThrHis PheGluAsnCys. The sequence published earlier in summary form (Iwaki et al., 1979, J. Biochem.86, 1159–1162) contained a few errors which are pointed out in this paper.     

The complete nucleotide sequence of a snake (Dinodon semicarinatus) mitochondrial genome with two identical control regions.

The 17,191-bp mitochondrial DNA (mtDNA) of a Japanese colubrid snake, akamata (Dinodon semicarinatus), was cloned and sequenced. The snake mtDNA has some peculiar features that were found in our previous study using polymerase chain reaction: duplicate control regions that have completely identical sequences over 1 kbp, translocation of tRNALeu(UUR) gene, shortened TpsiC arm for most tRNA genes, and a pseudogene for tRNAPro. Phylogenetic analysis of amino acid sequences of protein genes suggested an unusually high rate of molecular evolution in the snake compared to other vertebrates. Southern hybridization experiments using mtDNAs purified from multiple akamata individuals showed that the duplicate state of the control region is not a transient or unstable feature found in a particular individual, but that it stably occurs in mitochondrial genomes of the species. This may, therefore, be regarded as an unprecedented example of stable functional redundancy in animal mtDNA. However, some of the examined individuals contain a rather scanty proportion of heteroplasmic mtDNAs with an organization of genes distinct from that of the major mtDNA. The gene organization of the minor mtDNA is in agreement with one of models that we present to account for the concerted evolution of duplicate control regions.     

Complete DNA sequence of the mitochondrial genome of the ascidian Halocynthia roretzi (Chordata, Urochordata)

The complete nucleotide sequence of the 14,771-bp-long mitochondrial (mt) DNA of a urochordate (Chordata)-the ascidian Halocynthia roretzi-was determined. All the Halocynthia mt-genes were found to be located on a single strand, which is rich in T and G rather than in A and C. Like nematode and Mytilus edulis mtDNAs, that of Halocynthia encodes no ATP synthetase subunit 8 gene. However, it does encode an additional tRNA gene for glycine (anticodon TCT) that enables Halocynthia mitochondria to use AGA and AGG codons for glycine. The mtDNA carries an unusual tRNA(Met) gene with a TAT anticodon instead of the usual tRNA(Met)(CAT) gene. As in other metazoan mtDNAs, there is not any long noncoding region. The gene order of Halocynthia mtDNA is completely different from that of vertebrate mtDNAs except for tRNA(His)-tRNA(Ser)(GCU), suggesting that evolutionary change in the mt-gene structure is much accelerated in the urochordate line compared with that in vertebrates. The amino acid sequences of Halocynthia mt-proteins deduced from their gene sequences are quite different from those in other metazoans, indicating that the substitution rate in Halocynthia mt-protein genes is also accelerated.     

The amino acid sequence of human chorionic gonadotropin. The alpha subunit and beta subunit.

The amino acid sequences of both the alpha and beta subunits of human chorionic gonadotropin have been determined. The amino acid sequence of the alpha subunit is: Ala - Asp - Val - Gln - Asp - Cys - Pro - Glu - Cys-10 - Thr - Leu - Gln - Asp - Pro - Phe - Ser - Gln-20 - Pro - Gly - Ala - Pro - Ile - Leu - Gln - Cys - Met - Gly-30 - Cys - Cys - Phe - Ser - Arg - Ala - Tyr - Pro - Thr - Pro-40 - Leu - Arg - Ser - Lys - Lys - Thr - Met - Leu - Val - Gln-50 - Lys - Asn - Val - Thr - Ser - Glu - Ser - Thr - Cys - Cys-60 - Val - Ala - Lys - Ser - Thr - Asn - Arg - Val - Thr - Val-70 - Met - Gly - Gly - Phe - Lys - Val - Glu - Asn - His - Thr-80 - Ala - Cys - His - Cys - Ser - Thr - Cys - Tyr - Tyr - His-90 - Lys - Ser. Oligosaccharide side chains are attached at residues 52 and 78. In the preparations studied approximately 10 and 30% of the chains lack the initial 2 and 3 NH2-terminal residues, respectively. This sequence is almost identical with that of human luteinizing hormone (Sairam, M. R., Papkoff, H., and Li, C. H. (1972) Biochem. Biophys. Res. Commun. 48, 530-537). The amino acid sequence of the beta subunit is: Ser - Lys - Glu - Pro - Leu - Arg - Pro - Arg - Cys - Arg-10 - Pro - Ile - Asn - Ala - Thr - Leu - Ala - Val - Glu - Lys-20 - Glu - Gly - Cys - Pro - Val - Cys - Ile - Thr - Val - Asn-30 - Thr - Thr - Ile - Cys - Ala - Gly - Tyr - Cys - Pro - Thr-40 - Met - Thr - Arg - Val - Leu - Gln - Gly - Val - Leu - Pro-50 - Ala - Leu - Pro - Gin - Val - Val - Cys - Asn - Tyr - Arg-60 - Asp - Val - Arg - Phe - Glu - Ser - Ile - Arg - Leu - Pro-70 - Gly - Cys - Pro - Arg - Gly - Val - Asn - Pro - Val - Val-80 - Ser - Tyr - Ala - Val - Ala - Leu - Ser - Cys - Gln - Cys-90 - Ala - Leu - Cys - Arg - Arg - Ser - Thr - Thr - Asp - Cys-100 - Gly - Gly - Pro - Lys - Asp - His - Pro - Leu - Thr - Cys-110 - Asp - Asp - Pro - Arg - Phe - Gln - Asp - Ser - Ser - Ser - Ser - Lys - Ala - Pro - Pro - Pro - Ser - Leu - Pro - Ser-130 - Pro - Ser - Arg - Leu - Pro - Gly - Pro - Ser - Asp - Thr-140 - Pro - Ile - Leu - Pro - Gln. Oligosaccharide side chains are found at residues 13, 30, 121, 127, 132, and 138. The proteolytic enzyme, thrombin, which appears to cleave a limited number of arginyl bonds, proved helpful in the determination of the beta sequence.     

Mitochondrial gene order is not conserved in arthropods: prostriate and metastriate tick mitochondrial genomes

The entire mitochondrial genome was sequenced in a prostriate tick, Ixodes hexagonus, and a metastriate tick, Rhipicephalus sanguineus. Both genomes encode 22 tRNAs, 13 proteins, and two ribosomal RNAs. Prostriate ticks are basal members of Ixodidae and have the same gene order as Limulus polyphemus. In contrast, in R. sanguineus, a block of genes encoding NADH dehydrogenase subunit 1 (ND1), tRNA(Leu)(UUR), tRNA(Leu)(CUN), 16S rDNA, tRNA(Val), 12S rDNA, the control region, and the tRNA(Ile) and tRNA(Gln) have translocated to a position between the tRNA(Glu) and tRNA(Phe) genes. The tRNA(Cys) gene has translocated between the control region and the tRNA(Met) gene, and the tRNA(Leu)(CUN) gene has translocated between the tRNA(Ser)(UCN) gene and the control region. Furthermore, the control region is duplicated, and both copies undergo concerted evolution. Primers that flank these rearrangements confirm that this gene order is conserved in all metastriate ticks examined. Correspondence analysis of amino acid and codon use in the two ticks and in nine other arthropod mitochondrial genomes indicate a strong bias in R. sanguineus towards amino acids encoded by AT-rich codons.      

Sequence variation and structural conservation in the D-loop region and flanking genes of mitochondrial DNA from Japanese pond frogs

The nucleotide sequences of the D-loop region and its flanking genes of the mitochondrial DNA (mtDNA) from Japanese pond frogs were determined by the methods of PCR, cloning, and sequencing. The frogs belonged to two species, one subspecies, and one local race. The gene arrangements adjacent to the D-loop region were analyzed. The frogs shared a unique mitochondrial gene order that was found in Rana catesbeiana; i.e., cyt b--D-loop region--tRNA(Leu(CUN))--tRNA(Thr)--tRNA(Pro)--tRNA(Phe)--12S rRNA. The arrangements of the three tRNA genes of these frogs were different from those of X. laevis, a species which has the same overall structure as in mammals. Highly repetitive sequences with repeat units (16-bp or 17-bp sequence specific for each taxon) were found in the D-loop region. The length of repetitive sequences varied from 0.6 kbp to 1.2 kbp, and caused the extensive size variation in mtDNA. Several short sequence elements such as putative TAS, OH, CSB-1, and CSB-2 were found in the D-loop region of these frogs. The sequences of these short regulatory elements were conserved in R. catesbeiana, X. laevis, and also in human. The comparison of sequence divergences of the D-loop region and its adjacent genes among various taxa revealed that the rates of nucleotide substitutions depend on genes. The nucleotide sequences of the 3'-side segment of the D-loop region were the most variable among taxa, whereas those of the tRNA and 12S rRNA genes were the most conservative.     

Variations in mitochondrial tRNA gene organization of reptiles as phylogenetic markers

Amplification and sequencing of mitochondrial DNA regions corresponding to three major clusters of transfer RNA genes from a variety of species representing major groups of birds and reptiles revealed some new variations in tRNA gene organization. First, a gene rearrangement from tRNA(His)-tRNA(Ser)(AGY)-tRNA(Leu)(CUN) to tRNA(Ser)(AGY)- tRNA(His)tRNA(Leu)(CUN) occurs in all three crocodilians examined (alligator, caiman, and crocodile). In addition an exceptionally long spacer region between the genes for NADH dehydrogenase subunit 4 and tRNA(Ser)(AGY) is found in caiman. Second, in congruence with a recent finding by Seutin et al., a characteristic stem-and-loop structure for the putative light-strand replication origin located between tRNA(Asn) and tRNA(Cys) genes is absent for all the birds and crocodilians. This stem-and-loop structure is absent in an additional species, the Texas blind snake, whereas the stem-and-loop structure is present in other snakes, lizards, turtles, mammals, and a frog. The disappearance of the stem-and-loop structure in the blind snake most likely occurred independently of that on the lineage leading to birds and crocodilians. Finally, the blind snake has a novel type of tRNA gene arrangement in which the tRNA(Gln) gene moved from one tRNA cluster to another. Sequence substitution rates for the tRNA genes appeared to be somewhat higher in crocodialians than in birds and mammals. As regards the controversial phylogenetic relationship among the Aves, Crocodilia, and Mammalia, a sister group relationship of birds and crocodilians relative to mammals, as suggested from the common loss of the stem-and- loop structure, was supported with statistical significance by molecular phylogenetic analyses using the tRNA gene sequence data.      

Identification of coding polymorphisms in human circadian rhythm genes PER1, PER2, PER3, CLOCK, ARNTL, CRY1, CRY2 and TIMELESS in a multi-ethnic screening panel.

STUDY OBJECTIVE: In this study, the exonic regions of the circadian rhythm genes PER1, PER2, PER3, CLOCK, ARNTL, CRY1, CRY2 and TIMELESS were re-sequenced and coding changes identified in a panel of 95 individuals varying in ethnicity. STUDY PARTICIPANTS: DNA screening panel consisting of 95 DNA samples (17 American Caucasians, 17 African Americans, 8 Ashkenazi Jews, 8 Chinese, 8 Japanese, 5 Mexican Indians, 8 Mexicans, 8 Northern Europeans, 8 Puerto Ricans, and 8 South Americans) selected from the Coriell Institute Human Variation Panel. RESULTS: In addition to coding changes already identified in the database dbSNP, novel coding changes were identified, including PER1: Pro37Ser, Pro351Ser, Gln988Pro, Ala998Thr; PER2: Leu83Arg, Leu157Leu, Thre174Ile, Phe400Phe, Pro822Pro, Ala828Thr, Ala861Val, Phe876Leu, Val883Met, Val903Ile, Ala923Pro; PER3: Pro67Pro, Val90Ile, His638His, Ala820Ala, Leu929Leu; ARNTL: Arg166Gln, Ser459Phe; CLOCK: Ala34Ala, Ser208Cys, Phe233Phe, Ser632Thr, Ser816Ser; TIMELESS: Met870Val and CRY2: His35His. No coding polymorphisms were identified in CRY1. CONCLUSIONS: Considerable genetic variation occurs within the coding region of the genes regulating circadian rhythm. Many of the non-synonymous coding polymorphisms could affect protein structure/function with the potential to affect molecular regulation of the sleep/wake cycle. Many of the potential functional effects could be ethnic group specific.     

Complete DNA Sequence of the Mitochondrial Genome of the Black Chiton, Katharina Tunicata

The DNA sequence of the 15,532-base pair (bp) mitochondrial DNA (mtDNA) of the chiton Katharina tunicata has been determined. The 37 genes typical of metazoan mtDNA are present: 13 for protein subunits involved in oxidative phosphorylation, 2 for rRNAs and 22 for tRNAs. The gene arrangement resembles those of arthropods much more than that of another mollusc, the bivalve Mytilus edulis. Most genes abut directly or overlap, and abbreviated stop codons are inferred for four genes. Four junctions between adjacent pairs of protein genes lack intervening tRNA genes; however, at each of these junctions there is a sequence immediately adjacent to the start codon of the downstream gene that is capable of forming a stem-and-loop structure. Analysis of the tRNA gene sequences suggests that the D arm is unpaired in tRNA(ser(AGN)), which is typical of metazoan mtDNAs, and also in tRNA(ser(UCN)), a condition found previously only in nematode mtDNAs. There are two additional sequences in Katharina mtDNA that can be folded into structures resembling tRNAs; whether these are functional genes is unknown. All possible codons except the stop codons TAA and TAG are used in the protein-encoding genes, and Katharina mtDNA appears to use the same variation of the mitochondrial genetic code that is used in Drosophila and Mytilus. Translation initiates at the codons ATG, ATA and GTG. A + T richness appears to have affected codon usage patterns and, perhaps, the amino acid composition of the encoded proteins. A 142-bp non-coding region between tRNA(glu) and CO3 contains a 72-bp tract of alternating A and T.     

草鱼(Ctenopharyngodon idellus)线粒体DNA控制区及其两侧tRNA基因的克隆与结构分析

动物线粒体ＤＮＡ控制区是线粒体基因组复制与基因表达的最主要的调控区．采用杂交和测序的方法对草鱼线粒体ＤＮＡ控制区进行定位、克隆并测定了控制区及其旁侧的ｔＲＮＡＰｈｅ、ｒＲＮＡＰｒｏ和ｒＲＮＡＴｈｒ三个基因的序列,与多种脊椎动物的相应序列进行了比较,并进行了结构分析．草鱼线粒体控制区全长９２７ｂｐ,含有与酵母和爪蟾线粒体启动子相似的序列,其ＣＳＢⅠ、ＣＳＢⅡ和ＣＳＢⅢ序列与其他几种动物的ＣＳＢ比较相当保守,ＴＡＳ与其回文基序可形成稳定的茎环结构,成为Ｈ－链复制的终止信号．草鱼线粒体ｔＲＮＡＰｈｅ、ｔＲＮＡＰｒｏ和ｔＲＮＡＴｈｒ可折叠成三叶草形二级结构,其基因具有许多不同于细胞质ｔＲＮＡ基因的结构特点,可能反映了线粒体ｔＲＮＡ与线粒体核糖体具有不寻常的作用方式     

Nucleotide Sequence and Gene Organization of the Starfish Asterina Pectinifera Mitochondrial Genome

The 16,260-bp mitochondrial DNA (mtDNA) from the starfish Asterina pectinifera has been sequenced. The genes for 13 proteins, two rRNAs and 22 tRNAs are organized in an extremely economical fashion, similar to those of other animal mtDNAs, with some of the genes overlapping each other. The gene organization is the same as that for another echinoderm, sea urchin, except for the inversion of a 4.6-kb segment that contains genes for two proteins, 13 tRNAs and the 16S rRNA. Judging from the organization of the protein coding genes, mammalian mtDNAs resemble the sea urchin mtDNA more than that of the starfish. The region around the 3' end of the 12S rRNA gene of the starfish shows a high similarity with those for vertebrates. This region encodes a possible stem and loop structure; similar potential structures occur in this region of vertebrate mtDNAs and also in nonmitochondrial small subunit rRNA. A similar stem and loop structure is also found at the 3' end of the 16S rRNA genes in A. pectinifera, in another starfish Pisaster ochraceus, in vertebrates and in Drosophila, but not in sea urchins. The full sequence data confirm the presumption that AGA/AGG, AUA and AAA codons, respectively, code for serine, isoleucine, and asparagine in the starfish mitochondria, and that AGA/AGG codons are read by tRNA(GCU)(Ser), which possesses a truncated dihydrouridine arm, that was previously suggested from a partial mtDNA sequence. The structural characteristics of tRNAs and possible mechanisms for the change in the mitochondrial genetic code are also discussed.     

Sequence evidence in the archaeal genomes that tRNAs emerged through the combination of ancestral genes as 5' and 3' tRNA halves

The discovery of separate 5' and 3' halves of transfer RNA (tRNA) molecules-so-called split tRNA-in the archaeal parasite Nanoarchaeum equitans made us wonder whether ancestral tRNA was encoded on 1 or 2 genes. We performed a comprehensive phylogenetic analysis of tRNAs in 45 archaeal species to explore the relationship between the three types of tRNAs (nonintronic, intronic and split). We classified 1953 mature tRNA sequences into 22 clusters. All split tRNAs have shown phylogenetic relationships with other tRNAs possessing the same anticodon. We also mimicked split tRNA by artificially separating the tRNA sequences of 7 primitive archaeal species at the anticodon and analyzed the sequence similarity and diversity of the 5' and 3' tRNA halves. Network analysis revealed specific characteristics of and topological differences between the 5' and 3' tRNA halves: the 5' half sequences were categorized into 6 distinct groups with a sequence similarity of >80%, while the 3' half sequences were categorized into 9 groups with a higher sequence similarity of >88%, suggesting different evolutionary backgrounds of the 2 halves. Furthermore, the combinations of 5' and 3' halves corresponded with the variation of amino acids in the codon table. We found not only universally conserved combinations of 5'-3' tRNA halves in tRNA(iMet), tRNA(Thr), tRNA(Ile), tRNA(Gly), tRNA(Gln), tRNA(Glu), tRNA(Asp), tRNA(Lys), tRNA(Arg) and tRNA(Leu) but also phylum-specific combinations in tRNA(Pro), tRNA(Ala), and tRNA(Trp). Our results support the idea that tRNA emerged through the combination of separate genes and explain the sequence diversity that arose during archaeal tRNA evolution.     

The CO-I and CO-II region of honeybee mitochondrial DNA: evidence for variation in insect mitochondrial evolutionary rates

The sequence of a region of honeybee (Apis mellifera ligustica) mitochondrial DNA, which contains the genes for cytochrome c oxidase subunits I and II (CO-I and CO-II) and inferred genes for tRNA(Asp), tRNA(Leu)UUR, tRNA(Lys), and tRNA(Trp), is presented. The region includes the segment previously identified as incurring a length increase in some other bee strains, including Africanized bees. The sequence information of this study and of that by Vlasak et al. shows that several shifts of tRNA genes have occurred between Apis and Drosophila, but shifts of other kinds of genes have yet to be demonstrated. The CO-I and CO-II gene sequences are both more A+T rich than are the corresponding Drosophila genes. Parsimony analyses using the mouse and Xenopus sequences as outgroups show significantly more amino acid substitutions on the branch to Apis (120) than on that to Drosophila (44), indicating a difference in the long-term evolutionary rates of hymenopteran and dipteran mtDNA.     

Identification and location of nine T5 bacteriophage tRNA genes by DNA sequence analysis.

Sequence analysis of two DNA fragments generated from bacteriophage T5 DNA by restriction with Hpa I and Hae III has resulted in the detection and localization of nine tRNA genes (His, two Ser genes, Leu, Val, Lys, fMet, Pro, and Ile). The genes which code for tRNAs His and Leu are partials, whereas the remaining genes are complete. A majority of the tRNA genes are located in close proximity to one another. A unique feature of the Pro and Ile genes is that their DNA sequence overlap.     

Transfer RNA genes of Euglena gracilis chloroplast DNA

Summary Transfer RNA genes have been mapped to at least nine different loci on the physical map of the Euglena gracilis chloroplast genome. One of these loci in the ribosomal RNA operons is present three times per genome. The DNA sequences of six of the nine different loci, containing 21 different tRNA genes, have been determined. Genes corresponding to the amino acids Ala, Arg, Asn, Cys, Gln, Gly (2), Glu, His, Ile, Leu (2), Met (2), Phe, Ser, Thr, Trp, Tyr, Val, and one unassigned species have been identified. All genes except one are found in clusters of 2–6 genes. None of the known genes contains introns, nor codes for the 3-CCA terminus. In addition to these genes, two pseudo tRNA genes are present in the rDNA leader region.     

Editing corrects mispairing in the acceptor stem of bean and potato mitochondrial phenylalanine transfer RNAs.

Editing is a general event in plant mitochondrial messenger RNAs, but has never been detected in a plant mitochondrial transfer RNA (tRNA). We demonstrate here the occurrence of a tRNA editing event in higher plant mitochondria: in both bean and potato, the C encoded at position 4 in the mitochondrial tRNA(Phe)(GAA) gene is converted into a U in the mature tRNA. This nucleotide change corrects the mismatched C4-A69 base-pair which appears when folding the gene sequence into the cloverleaf structure and it is consistent with the fact that C to U transitions constitute the common editing events affecting plant mitochondrial messenger RNAs. The tRNA(Phe)(GAA) gene is located upstream of the single copy tRNA(Pro)(UGG) gene in both the potato and the bean mitochondrial DNAs. The sequences of potato and bean tRNA(Pro)(UGG) genes are colinear with the sequence of the mature bean mitochondrial tRNA(Pro)(UGG), demonstrating that this tRNA is not edited. A single copy tRNA(Ser)(GCU) gene was found upstream of the tRNA(Phe) gene in the potato mitochondrial DNA. A U6-U67 mismatched base-pair appears in the cloverleaf folding of this gene and is maintained in the mature potato mitochondrial tRNA(Ser)(GCU), which argues in favour of the hypothesis that the editing system of plant mitochondria can only perform C to U or occasionally U to C changes.