Evolution of Pulmonate Gastropod Mitochondrial Genomes: Comparisons of Gene Organizations of Euhadra, Cepaea and Albinaria and Implications of Unusual Trna Secondary Structures

Complete gene organizations of the mitochondrial genomes of three pulmonate gastropods, Euhadra herklotsi, Cepaea nemoralis and Albinaria coerulea, permit comparisons of their gene organizations. Euhadra and Cepaea are classified in the same superfamily, Helicoidea, yet they show several differences in the order of tRNA and protein coding genes. Albinaria is distantly related to the other two genera but shares the same gene order in one part of its mitochondrial genome with Euhadra and in another part with Cepaea. Despite their small size (14.1-14.5 kbp), these snail mtDNAs encode 13 protein genes, two rRNA genes and at least 22 tRNA genes. These genomes exhibit several unusual or unique features compared to other published metazoan mitochondrial genomes, including those of other molluscs. Several tRNAs predicted from the DNA sequences possess bizarre structures lacking either the T stem or the D stem, similar to the situation seen in nematode mt-tRNAs. The acceptor stems of many tRNAs show a considerable number of mismatched basepairs, indicating that the RNA editing process recently demonstrated in Euhadra is widespread in the pulmonate gastropods. Strong selection acting on mitochondrial genomes of these animals would have resulted in frequent occurrence of the mismatched basepairs in regions of overlapping genes.     

Three divergent mitochondrial genomes from California populations of the copepod Tigriopus californicus

Previous work on the harpacticoid copepod Tigriopus californicus has focused on the extensive population differentiation in three mtDNA protein coding genes (COXI, COXII, Cytb). In order to get a more complete understanding of mtDNA evolution in this species, we sequenced three complete mitochondrial genomes (one from each of three California populations) and compared them to two published mtDNA genomes from an Asian congener, Tigriopus japonicus. Several features of the mtDNA genome appear to be conserved within the genus: 1) the unique order of the protein coding genes, rRNA genes and most of the tRNA genes, 2) the genome is compact, varying between 14.3 and 14.6 kb, and 3) all genes are encoded on the same strand of the mtDNA. Within T. californicus, extremely high levels of nucleotide divergence (>20%) are observed across much of the mitochondrial genome. Inferred amino acid sequences of the proteins encoded in the mtDNAs also show high levels of divergence; at the extreme, the three ND3 variants in T. californicus showed >25% amino acid substitutions, compared with <3% amino acid divergence at the previously studied COXI locus. Unusual secondary structures make functional assignments of some tRNAs difficult. The only apparent tRNA(trp) in these genomes completely overlaps the 5' end of the 16S rRNA in all three T. californicus mtDNAs. Although not previously noted, this feature is also conserved in T. japonicus mtDNAs; whether this sequence is processed into a functional tRNA has not been determined. The putative control region contains a duplicated segment of different length (from 88 to 155 bp) in each of the T. californicus sequences. In each case, the duplicated segments are not tandem repeats; despite their different lengths, the distance between the start of the first and the start of the second repeat is conserved (520 bp). The functional significance, if any, of this repeat structure remains unknown.     

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.     

Parallel evolution of truncated transfer RNA genes in arachnid mitochondrial genomes

The cloverleaf secondary structure of transfer RNA (tRNA) is highly conserved across all forms of life. Here, we provide sequence data and inferred secondary structures for all tRNA genes from 8 new arachnid mitochondrial genomes, including representatives from 6 orders. These data show remarkable reductions in tRNA gene sequences, indicating that T-arms are missing from many of the 22 tRNAs in the genomes of 4 out of 7 orders of arachnids. Additionally, all opisthothele spiders possess some tRNA genes that lack sequences that could form well-paired aminoacyl acceptor stems. We trace the evolution of T-arm loss onto phylogenies of arachnids and show that a genome-wide propensity to lose sequences that encode canonical cloverleaf structures likely evolved multiple times within arachnids. Mapping of structural characters also shows that certain tRNA genes appear more evolutionarily prone to lose the sequence coding for the T-arm and that once a T-arm is lost, it is not regained. We use tRNA structural data to construct a phylogeny of arachnids and find high bootstrap support for a clade that is not supported in phylogenies that are based on more traditional morphological characters. Together, our data demonstrate variability in structural evolution among different tRNAs as well as evidence for parallel evolution of the loss of sequence coding for tRNA arms within an ancient and diverse group of animals.     

Sequence evolution of mitochondrial tRNA genes and deep-branch animal phylogenetics

Mitochondrial DNA sequences are often used to construct molecular phylogenetic trees among closely related animals. In order to examine the usefulness of mtDNA sequences for deep-branch phylogenetics, genes in previously reported mtDNA sequences were analyzed among several animals that diverged 20–600 million years ago. Unambiguous alignment was achieved for stem-forming regions of mitochondrial tRNA genes by virtue of their conservative secondary structures. Sequences derived from stem parts of the mitochondrial tRNA genes appeared to accumulate much variation linearly for a long period of time: nearly 100 Myr for transition differences and more than 350 Myr for transversion differences. This characteristic could be attributed, in part, to the structural variability of mitochondrial tRNAs, which have fewer restrictions on their tertiary structure than do nonmitochondrial tRNAs. The tRNA sequence data served to reconstruct a well-established phylogeny of the animals with 100% bootstrap probabilities by both maximum parsimony and neighbor joining methods. By contrast, mitochondrial protein genes coding for cytochrome b and cytochrome oxidase subunit I did not reconstruct the established phylogeny or did so only weakly, although a variety of fractions of the protein gene sequences were subjected to tree-building. This discouraging phylogenetic performance of mitochondrial protein genes, especially with respect to branches originating over 300 Myr ago, was not simply due to high randomness in the data. It may have been due to the relative susceptibility of the protein genes to natural selection as compared with the stem parts of mitochondrial tRNA genes. On the basis of these results, it is proposed that mitochondrial tRNA genes may be useful in resolving deep branches in animal phylogenies with divergences that occurred some hundreds of Myr ago. For this purpose, we designed a set of primers with which mtDNA fragments encompassing clustered tRNA genes were successfully amplified from various vertebrates by the polymerase chain reaction.Abbreviations AA stem amino acid-acceptor stem - AC stem anticodon stem - COI cytochrome oxidase subunit I - cytb cytochrome b - D stem dihydrouridine stem - MP maximum parsimony - mtDNA mitochondrial DNA - Myr million years - NJ neighbor joining - PCR polymerase chain reaction - Ti transition - T stem tC stem - Tv transversion Correspondence to: Y. Kumazawa     

Nucleotides of tRNA governing the specificity of Escherichia coli methionyl-tRNA(fMet) formyltransferase.

In Escherichia coli, the free amino group of the aminoacyl moiety of methionyl-tRNA(fMet) is specifically modified by a transformylation reaction. To identify the nucleotides governing the recognition of the tRNA substrate by the formylase, initiator tRNA(fMet) was changed into an elongator tRNA with the help of an in vivo selection method. All the mutations isolated were in the tRNA acceptor arm, at positions 72 and 73. The major role of the acceptor arm was further established by the demonstration of the full formylability of a chimaeric tRNA(Met) containing the acceptor stem of tRNA(fMet) and the remaining of the structure of tRNA(mMet). In addition, more than 30 variants of the genes encoding tRNA(mMet) or tRNA(fMet) have been constructed, the corresponding mutant tRNA products purified and the parameters of the formylation reaction measured. tRNA(mMet) became formylatable by the only change of the G1.C72 base-pair into C1-A72. It was possible to render tRNA(mMet) as good a substrate as tRNA(fMet) for the formylase by the introduction of a limited number of additional changes in the acceptor stem. In conclusion, A73, G2.C71, C3.G70 and G4.C69 are positive determinants for the specific processing of methionyl-tRNA(fMet) by the formylase while the occurrence of a G.C or C.G base-pair between positions 1 and 72 acts as a major negative determinant. This pattern appears to account fully for the specificity of the formylase and the lack of formylation of any aminoacylated tRNA, excepting the methionyl-tRNA(fMet).     

Aminoacyl RNA domain of turnip yellow mosaic virus Val-RNA interacting with elongation factor Tu

Turnip yellow mosaic virus (TYMV) Val-RNA forms a complex with the peptide elongation factor Tu (EF-Tu) in the presence of GTP: the Val-RNA is protected by EF-Tu·GTP from non-enzymatic deacylation and nuclease digestion. The determination of the length of the shortest TYMV Val-RNA fragment that binds EF-Tu·GTP leads us to conclude that the valylated aminoacyl RNA domain equivalent in tRNAs to the continuous helix formed by the acceptor stem and the T arm is sufficient for complex formation. Since the aminoacyl RNA domain is also sufficient for adenylation by the ATP(CTP):tRNA nucleotidyltransferase, an analogy can be drawn between these two tRNA-specific proteins.     

RNA editing in trans-splicing intron sequences of nad2 mRNAs in Oenothera mitochondria.

The complete open reading frame of subunit 2 of the NADH dehydrogenase in Oenothera mitochondria is split into five exons. The first two and the last three exons are encoded in distant genomic locations and are transcribed separately. Three tRNA genes coding for tRNA(Cys), tRNA(Asn), and tRNA(Tyr) are located upstream of the terminal three exons c, d, and e. The genomic distance, the interspersed tRNA genes, and the group II intron sequences flanking the two separated exons suggest trans-splicing to be required to connect exons b and c. Maturation of the mRNA includes RNA editing at 36 sites in the open reading frame. Three RNA editing events are observed in the split group II intron sequences. Two of these events allow after editing additional base pairings in the secondary structure, one in the stem of domain I, the other in the putative trans-pairing region of domain IV. These RNA editings may thus be involved in the trans-splicing reaction.     

A unique serine-specific elongation factor Tu found in nematode mitochondria

The translation elongation factor Tu (EF-Tu) delivers aminoacyl-tRNAs to ribosomes by recognizing the tRNA acceptor and T stems. However, the unusual truncation observed in some animal mitochondrial tRNAs seems to prevent recognition by a canonical EF-Tu. For instance, nematode mitochondria contain tRNAs lacking a T or D arm. We recently found an atypical EF-Tu (EF-Tu1) specific for nematode mitochondrial tRNAs that lack the T arm. We have now discovered a second factor, EF-Tu2, which binds only to tRNAs that lack a D arm. EF-Tu2 seems unique in its amino acid specificity because it recognizes the aminoacyl moiety of seryl-tRNAs and the tRNA structure itself. Such EF-Tu evolution might explain tRNA structural divergence in animal mitochondria.     

Seventeen new complete mtDNA sequences reveal extensive mitochondrial genome evolution within the Demospongiae

Two major transitions in animal evolution--the origins of multicellularity and bilaterality--correlate with major changes in mitochondrial DNA (mtDNA) organization. Demosponges, the largest class in the phylum Porifera, underwent only the first of these transitions and their mitochondrial genomes display a peculiar combination of ancestral and animal-specific features. To get an insight into the evolution of mitochondrial genomes within the Demospongiae, we determined 17 new mtDNA sequences from this group and analyzing them with five previously published sequences. Our analysis revealed that all demosponge mtDNAs are 16- to 25-kbp circular molecules, containing 13-15 protein genes, 2 rRNA genes, and 2-27 tRNA genes. All but four pairs of sampled genomes had unique gene orders, with the number of shared gene boundaries ranging from 1 to 41. Although most demosponge species displayed low rates of mitochondrial sequence evolution, a significant acceleration in evolutionary rates occurred in the G1 group (orders Dendroceratida, Dictyoceratida, and Verticillitida). Large variation in mtDNA organization was also observed within the G0 group (order Homosclerophorida) including gene rearrangements, loss of tRNA genes, and the presence of two introns in Plakortis angulospiculatus. While introns are rare in modern-day demosponge mtDNA, we inferred that at least one intron was present in cox1 of the common ancestor of all demosponges. Our study uncovered an extensive mitochondrial genomic diversity within the Demospongiae. Although all sampled mitochondrial genomes retained some ancestral features, including a minimally modified genetic code, conserved structures of tRNA genes, and presence of multiple non-coding regions, they vary considerably in their size, gene content, gene order, and the rates of sequence evolution. Some of the changes in demosponge mtDNA, such as the loss of tRNA genes and the appearance of hairpin-containing repetitive elements, occurred in parallel in several lineages and suggest general trends in demosponge mtDNA evolution.     

Aminoacylation of RNA Minihelices: Implications for tRNA Synthetase Structural Design and Evolution

Abstract

The genetic code is based on the aminoacylation of tRNA with amino acids catalyzed by the aminoacyl-tRNA synthetases. The synthetases are constructed from discrete domains and all synthetases possess a core catalytic domain that catalyzes amino acid activation, binds the acceptor stem of tRNA, and transfers the amino acid to tRNA. Fused to the core domain are additional domains that mediate RNA interactions distal to the acceptor stem. Several synthetases catalyze the aminoacylation of RNA oligonucleotide substrates that recreate only the tRNA acceptor stems. In one case, a relatively small catalytic domain catalyzes the aminoacylation of these substrates independent of the rest of the protein. Thus, the active site domain may represent a primordial synthetase in which polypeptide insertions that mediate RNA acceptor stem interactions are tightly integrated with determinants for aminoacyl adenylate synthesis. The relationship between nucleotide sequences in small RNA oligonucleotides and the specific amino acids that are attached to these oligonucleotides could constitute a second genetic code.     

The complete DNA sequence of the mitochondrial genome of Physarum polycephalum

The complete sequence of the mitochondrial DNA (mtDNA) of the true slime mold Physarun polycephalum has been determined. The mtDNA is a circular 62,862-bp molecule with an A+T content of 74.1%. A search with the program BLAST X identified the protein-coding regions. The mitochondrial genome of P. polycephalum was predicted to contain genes coding for 12 known proteins [for three cytochrome c oxidase subunits, apocytochrome b, two F1Fo-ATPase subunits, five NADH dehydrogenase (nad) subunits, and one ribosomal protein], two rRNA genes, and five tRNA genes. However, the predicted ORFs are not all in the same frame, because mitochondrial RNA in P. polycephalum undergoes RNA editing to produce functional RNAs. The nucleotide sequence of an nad7 cDNA showed that 51 nucleotides were inserted at 46 sites in the mRNA. No guide RNA-like sequences were observed in the mtDNA of P. polycephalum. Comparison with reported Physarum mtDNA sequences suggested that sites of RNA editing vary among strains. In the Physarum mtDNA, 20 ORFs of over 300 nucleotides were found and ORFs 14 19 are transcribed.     

Crystal structure of an Escherichia coli tRNA(Gly) microhelix at 2.0 A resolution

tRNA identity elements determine the correct aminoacylation by the cognate aminoacyl-tRNA synthetase. In class II aminoacyl tRNA synthetase systems, tRNA specificity is assured by rather few and simple recognition elements, mostly located in the acceptor stem of the tRNA. Here we present the crystal structure of an Escherichia coli tRNA(Gly) aminoacyl stem microhelix at 2.0 A resolution. The tRNA(Gly) microhelix crystallizes in the space group P3(2)21 with the cell constants a=b=35.35 A, c=130.82 A, gamma=120 degrees . The helical parameters, solvent molecules and a potential magnesium binding site are discussed.