The non-monophyletic origin of the tRNA molecule and the origin of genes only after the evolutionary stage of the last universal common ancestor (LUCA)

A model has been proposed suggesting that the tRNA molecule must have originated by direct duplication of an RNA hairpin structure [Di Giulio, M., 1992. On the origin of the transfer RNA molecule. J. Theor. Biol. 159, 199-214]. A non-monophyletic origin of this molecule has also been theorized [Di Giulio, M., 1999. The non-monophyletic origin of tRNA molecule. J. Theor. Biol. 197, 403-414]. In other words, the tRNA genes evolved only after the evolutionary stage of the last universal common ancestor (LUCA) through the assembly of two minigenes codifying for different RNA hairpin structures, which is what the exon theory of genes suggests when it is applied to the model of tRNA origin. Recent observations strongly corroborate this theorization because it has been found that some tRNA genes are completely separate in two minigenes codifying for the 5' and 3' halves of this molecule [Randau, L., et al., 2005a. Nanoarchaeum equitans creates functional tRNAs from separate genes for their 5'- and 3'-halves. Nature 433, 537-541]. In this paper it is shown that these tRNA genes codifying for the 5' and 3' halves of this molecule are the ancestral form from which the tRNA genes continuously codifying for the complete tRNA molecule are thought to have evolved. This, together with the very existence of completely separate tRNA genes codifying for their 5' and 3' halves, proves a non-monophyletic origin for tRNA genes, as a monophyletic origin would exclude the existence of these genes which have, on the contrary, been observed. Here the polyphyletic origin of genes codifying for proteins is also suggested and discussed. Moreover, a hypothesis is advanced to suggest that the LUCA might have had a fragmented genome made up of RNA and the possibility that 'Paleokaryotes' may exist is outlined. Finally, the characteristic of the indivisibility of homology that these polyphyletic origins seem to remove at the sequence level is discussed.     

The origin of the tRNA molecule: Independent data favor a specific model of its evolution

The properties, historical and empirical observations of a model of the origin of the tRNA molecule are discussed. This model would predict that this molecule originated by means of the assembly of two hairpin-like structures of RNA. The conclusion is that the model possesses a relevant part of the truth on the origin of the tRNA molecule.     

Permuted tRNA genes of Cyanidioschyzon merolae, the origin of the tRNA molecule and the root of the Eukarya domain

An evolutionary analysis is conducted on the permuted tRNA genes of Cyanidioschyzon merolae, in which the 5′ half of the tRNA molecule is codified at the 3′ end of the gene and its 3′ half is codified at the 5′ end. This analysis has shown that permuted genes cannot be considered as derived traits but seem to possess characteristics that suggest they are ancestral traits, i.e. they originated when tRNA molecule genes originated for the first time. In particular, if the hypothesis that permuted genes are a derived trait were true, then we should not have been able to observe that the most frequent class of permuted genes is that of the anticodon loop type, for the simple reason that this class would derive by random permutation from a class of non-permuted tRNA genes, which instead is the rarest. This would not explain the high frequency with which permuted tRNA genes with perfectly separate 5′ and 3′ halves were observed. Clearly the mechanism that produced this class of permuted genes would envisage the existence, in an advanced stage of evolution, of minigenes codifying for the 5′ and 3′ halves of tRNAs which were assembled in a permuted way at the origin of the tRNA molecule, thus producing a high frequency of permuted genes of the class here referred. Therefore, this evidence supports the hypothesis that the genes of the tRNA molecule were assembled by minigenes codifying for hairpin-like RNA molecules, as suggested by one model for the origin of tRNA [Di Giulio, M., 1992. On the origin of the transfer RNA molecule. J. Theor. Biol. 159, 199–214; Di Giulio, M., 1999. The non-monophyletic origin of tRNA molecule. J. Theor. Biol. 197, 403–414]. Moreover, the late assembly of the permuted genes of C. merolae, as well as their ancestrality, strengthens the hypothesis of the polyphyletic origins of these genes. Finally, on the basis of the uniqueness and the ancestrality of these permuted genes, I suggest that the root of the Eukarya domain is in the super-ensemble of the Plantae and that the Rhodophyta to which C. merolae belongs are the first line of divergence.     

Specific interaction of proteins with 5 S RNA and tRNA in the 42 S storage particle of Xenopus oocytes

During early oogenesis in amphibia, most of the 5 S RNA and tRNA is stored in a ribonucleoprotein particle that sediments at 42 S. In Xenopus laevis the 42 S particle contains two major proteins: of M_r 48 000 (P48) and 43 000 (P43). It is shown that heterogeneity in composition of the 42 S particle reflects a changing situation whereby initially, both 5 S RNA and tRNA are complexed with P48 (1 molecule 5 S RNA: 1 molecule P48; 2 or 3 molecules tRNA: 1 molecule P48), but later, tRNA becomes increasingly associated with P43 (in a 1:1 ratio) although 5 S RNA remains complexed with a cleavage product of P48. These changes relate to the eventual utilization of the excess 5 S RNA and tRNA in ribosome assembly and protein synthesis.     

Molecular evolution of transfer RNA from two precursor hairpins: Implications for the origin of protein synthesis

In this paper we are going to present a model for the coevolution of major components of the protein synthesis machinery in a primordial RNA world. We propose that the essential prerequisites for RNA-based protein synthesis, i.e., tRNA-like molecules, ribozymic charging catalysts, small-subunit(SSU) rRNA, and large-subunit(LSU) rRNA, evolved from the same ancestral RNA molecule. Several arguments are considered which suggest that tRNA-like molecules were derived by tandem joining of template-flanking hairpin structures involved in replication control. It is further argued that the ancestors of contemporary group I tRNA introns catalyzed such hairpin joining reactions, themselves also giving rise to the ribosomal RNAs. Our model includes a general stereochemical principle for the interaction between ribozymes and hairpin-derived recognition structures, which can be applied to such seemingly different processes as RNA polymerization, aminoacylation, tRNA decoding, and peptidyl transfer, implicating a common origin for these fundamental functions. These and other considerations suggest that generation and evolution of tRNA were coupled to the evolution of synthetases, ribosomal RNAs, and introns from the beginning and have been a consequence arising from the original function of tRNA precursor hairpins as replication and recombination control elements. Correspondence to: T.P. Dick     

Intensity fluctuation spectroscopy and transfer RNA conformation. III. Influence of NaCl concentration on the size and shape of the initially salt-free tRNA in solution.

The influence of sodium ion concentration in solution on the initially salt-free conformation of bulk tRNA from baker's yeast has been investigated by means of photon correlation spectroscopy. From the measured values of translational (D^T) and rotational (D^R) diffusion coefficients, the semiaxes of an ellipsoid of revolution, which are hydrodynamically equivalent to the tRNA molecule, were calculated for tRNA solutions in pure H₂O as well as in 0.005, 0.1, 0.5M NaCl and 0.01M MgCl₂ solutions at pH 4.2 and 7.5. These data, combined with our previous studies, suggested a model which describes the formation of an ordered tRNA structure due to increasing NaCl concentrations. Furthermore, we have obtained information concerning intermolecular interactions between tRNA molecules in solution. In low-salt or salt-free tRNA solutions, we detected in the linewidth distribution function an extra-fast component which can be attributed as possibly due to charge fluctuations related to the reaction of ionization of organic bases. In our light-scattering linewidth measurements, we do not see fluctuations of charged and uncharged states directly as concentration fluctuations. Rather, we postulate a modulation of long-range intermolecular electrostatic interactions between the tRNA molecules due to such charge fluctuations. It is this modulation which is related to the fast component of the time correlation function at finite concentrations. A quantitative theory is needed to provide a more definitive explanation of the dynamical behavior of tRNA in salt-free or low-salt solutions.     

A method using 3-O-methyl-D-glucose and phloretin for the determination of intracellular water space of cells in monolayer culture.

A method is presented for determining the extent of methylation of tRNAs synthesized in mammalian and bacterial cell systems and is based upon determining the distribution of radioactivity associated with the guanine constituents of total cellular tRNA preparations previously labeled with [2-¹⁴C]guanosine and with [methyl]-³H or -¹⁴C]methionine. Whereas labeling with guanosine provides a means of assessing the extent of methylation of the [2-¹⁴C]guanine residues incorporated into tRNA, methionine labeling provides a measure of the percentage of [methyl-³H or -¹⁴C]methylated constituents that are methylated guanines. Analyses such as the above reveal that the tRNA of KB cells acquires approximately three times as many methyl groups as that of E. coli B tRNA. Coupled with the knowledge that both mammalian and bacterial tRNA preparations contain an average of 24 guanine residues per molecule, the above analyses further reveal that 7.2 and 2.4 methyl groups are incorporated into each tRNA molecule synthesized in exponentially growing KB- and E. coli B-cells, respectively. Additional information regarding the extent of formation of individual methylated constituents per tRNA molecule synthesized is presented.     

The origin of the tRNA molecule: implications for the origin of protein synthesis

A model for tRNA molecule origin is discussed. The model postulates that this molecule originated simply by direct duplication (and subsequent evolution) of a gene coding for an RNA hairpin structure, which can thus be hypothesized as the evolutionary precursor of the tRNA molecule. The main properties are defined for these hairpin structures and it is suggested that these structures might have housed, near their 3' end, anticodons that were transferred to the loop of the tRNA anticodon during duplication of the hairpin structures. Moreover, the main characteristics are given for the evolutionary intermediary formed by direct duplication of the hairpin structure, i.e. the double hairpin. The evolutionary stages envisaged by this model for tRNA origin seem to naturally imply some evolutionary transitions through which the origin of protein synthesis passed. Finally, some strong historical evidence is provided to corroborate the model.     

Mitochondrial genome sequence of <Emphasis Type="Italic">Unionicola parkeri</Emphasis> (Acari: Trombidiformes: Unionicolidae): molecular synapomorphies between closely-related <Emphasis Type="Italic">Unionicola</Emphasis> gill mites

The mitochondrial genome of Unionicola parkeri is a 14,734 bp circular DNA molecule. The sequence and annotation revealed a unique gene order, related to but distinct from the gene order in the closely related species U. foili. Mitochondrial tRNA sequences annotated in this genome predict non-canonical secondary structures for these molecules. The continuing patterns of unique gene orders and unusual tRNA structures in the Trombidiformes in general and Unionicola in particular support the use of phylogenetic approaches that use these types of molecular features as shared, derived character states. Further progress in using these molecular character states to reconstruct phylogeny will depend on careful annotation, especially of tRNA genes.     

Multimers of a suppressor transfer RNA: Supporting evidence for alternate conformations of the anticodon loop region

Multimeric aggregates of tyrosine suppressor tRNA have been identified. The formation of multimers requires a high concentration of tRNA, as is formed in an electrophoretic moving boundary, and a specific anticodon sequence.A simple model is presented where the anticodon loop, in an hf configuration (Crick et al., 1976) interacts with the CCA stem of another tyrosine suppressor tRNA molecule.     

Studies on Isoleucyl-tRNA Synthetase of an Extreme Thermophile,Thermus thermophilus HB8

Isoleucyl-tRNA synthetase (IRS) was partially purified from an extreme thermophile, T. thermophilus HB 8. The molecular weight (11.5 ×10⁴) and some kinetic constants were obtained and compared with IRS from other sources.

The present IRS catalyzed both isoleucine dependent and valine dependent ATP-PPi exchange reactions (optimum at around 80°C) but not valyl-tRNA formation. The optimum temperature for isoleucyl-tRNA formation was 62°C with E. coli tRNA and 75°C with T. thermophilus tRNA.

The enzyme showed a remarkable thermostability. The addition of E. coli or T. thermophilus tRNA enhanced the thermostability of the enzyme, which was shown to be fully active up to 77°C. When E. coli tRNA was used, the loading activity decreased in parallel to the unfolding of the substrate tRNA molecule. From these results the relation is discussed between tRNA conformation and function.     

A model for the interaction of the G3‐subdomain of Geobacillus stearothermophilus IF2 with the 30S ribosomal subunit

Bacterial translation initiation factor IF2 complexed with GTP binds to the 30S ribosomal subunit, promotes ribosomal binding of fMet‐tRNA, and favors the joining of the small and large ribosomal subunits yielding a 70S initiation complex ready to enter the translation elongation phase. Within the IF2 molecule subdomain G3, which is believed to play an important role in the IF2‐30S interaction, is positioned between the GTP‐binding G2 and the fMet‐tRNA binding C‐terminal subdomains. In this study the solution structure of subdomain G3 of Geobacillus stearothermophilus IF2 has been elucidated. G3 forms a core structure consisting of two β‐sheets with each four anti‐parallel strands, followed by a C‐terminal α‐helix. In line with its role as linker between G3 and subdomain C1, this helix has no well‐defined orientation but is endowed with a dynamic nature. The structure of the G3 core is that of a typical OB‐fold module, similar to that of the corresponding subdomain of Thermus thermophilus IF2, and to that of other known RNA‐binding modules such as IF2‐C2, IF1 and subdomains II of elongation factors EF‐Tu and EF‐G. Structural comparisons have resulted in a model that describes the interaction between IF2‐G3 and the 30S ribosomal subunit.     

Stabilisation of some of the protein synthesis components in the thermophilic fungus,Humicola lanuginosa

The thermal stabilities of tRNA from the thermophilic fungus,Humicola lanuginose were compared with that from the mesophilic yeast,Candida utilis, by measuring the increase in the optical density with temperature. tRNAs from both the species were stable in the presence of millimolar quantities of magnesium chloride upto 50°C, the optimum growth temperature of the fungus. Aminoacyl tRNA synthetases were maximally active at 40°C under thein vitro assay conditions. They were fractionated and one species of valine tRNA synthetase was purified to homogeneity. The purified enzyme was protected against inactivation to varying degrees when preincubated with the substrates valine, tRNA and ATP as well as spermine. Protein turnover studies showed that the rate of turnover was higher at higher temperatures. It was concluded from these results that the protein synthesizing machinery of this fungus has no intrinsic stability but it is stabilised by intracellular factors. Higher rate of protein turnover also plays a role for growth at higher temperature.     

A critical role of water in the specific cleavage of the anticodon loop of some eukaryotic methionine initiator tRNAs

We have noticed that during a long storage and handling, the plant methionine initiator tRNA is spontaneously hydrolyzed within the anticodon loop at the C34-A35 phosphodiester bond. A literature search indicated that there is also the case for human initiator tRNA^Met but not for yeast tRNA^Met _i or E. coli tRNA^Met _f. All these tRNAs have an identical nucleotide sequence of the anticodon stems and loops with only one difference at position 33 within the loop. It means that cytosine 33 (C33) makes the anticodon loop of plant and human tRNA^Met _i susceptible to the specific cleavage reaction. Using crystallographic data of tRNA^Met _f of E. coli with U33, we modeled the anticodon loop of this tRNA with C33. We found that C33 within the anticodon loop creates a pocket that can accomodate a hydrogen bonded water molecule that acts as a general base and catalyzes a hydrolysis of C-A bond. We conclude that a single nucleotide change in the primary structure of tRNA^Met _i made changes in hydration pattern and readjustment in hydrogen bonding which lead to a cleavage of the phosphodiester bond.     

Extension of Watson's Model for the Elongation Cycle of Protein Biosynthesis

Abstract

The scheme for the elongation cycle of protein biosynthesis is proposed based on modern quantitative data on the interactions of mRNA and different functional forms of tRNA with 70S ribosomes and their 30S and 50S subunits. This scheme takes into account recently discovered third ribosomal (E) site with presumable exit function. The E site is introduced into 70S ribosome by its 50S subunit, the codon-anticodon interaction does not take place at the E site, and the affinity of tRNA for the E site is considerably lower than that for the P site. On the other hand, the P and A sites are located mainly on a 30S subunit, the codon-anticodon interactions being realized on both these sites. An mRNA molecule is placed exclusively on a 30S subunit where it makes U-turn. The proposed scheme does not contradict to any data but includes all main postulates of the initial Watson's model (J. D. Watson, Bull Soc. Chim. Biol. 46, 1399 (1964), and is considered as a natural extension of the later according to modern experimental data.     

Complete DNA sequence of the mitochondrial genome ofCepaea nemoralis (Gastropoda: Pulmonata)

The nucleotide sequence of a mitochondrial genome of the pulmonate gastropod molluscCepaea nemoralis has been determined. Contained within the 14,100 basepairs (bp) are the two ribosomal RNA genes and 13 protein coding genes typical of metazoan mitochondrial genomes. TheCepaea mtDNA does contain a gene for ATPase subunit 8, like the clausiliid pulmonate,Albinaria, and the chiton,Katharina, but unlike the bivalve mollusc,Mytilus. The mitochondrial genetic code ofCepaea is proposed to be the same as that ofMytilus, Katharina, andDrosophila. Only 14 putative tRNA genes are presented, although there is sufficient unassigned sequence to encode the remainder of the expected total of 22 tRNA genes. These 14 tRNA genes are a mixture of standard cloverleaf structures and nonstandard structures containing TV replacement loops as seen in nematode and mosquito mitochondrial genomes. If the eight unidentified tRNA genes are indeed present, very little unassigned sequence would remain to serve as a control region. Genes are transcribed from both strands of the molecule. Base composition is the least biased for any reported animal mitochondrial genome and is also very little skewed between strands using measures independent of base composition. TheCepaea mitochondrial gene order is quite unlike that of any other reported metazoan mtDNA, with the exception of the recently reported partial sequences ofAlbinaria. No gene bound-aries are shared among all the reported molluscan taxa, demonstrating a complete lack of conservation of mitochondrial gene order across the phylum Mollusca.     

Yeast tRNAAsp-Aspartyl-tRNA Synthetase Complex: Low Resolution Crystal Structure

Abstract

Yeast aspartyl-tRNA synthetase, a dimer of molecular weight 125000, and two molecules of its cognate tRNA (M_r = 24160) cocrystallize in the cubic space group 1432 (a = 354 Å). The crystal structure was solved to low resolution using neutron and X-ray diffraction data. Neutron single crystal diffraction data were collected in five solvents differing by their D₂O content in order to use the contrast variation method to distinguish between the protein and tRNA The synthetase was first located at 40 Å resolution using the 65% D₂O neutron data (tRNA matched). tRNA molecules were found at 20 Å resolution using both neutron and X-ray data. The resulting model was refined against 10 Å resolution X-ray data, using density modification and least-squares refinement of the tRNA positions. The crystal structure, solved without a priori phase knowledge, was confirmed later by isomorphous replacement. The molecular model of the complex is in good agreement with results obtained in solution by probing the protected part of the tRNA by chemical reagents.     

Increased tRNA modification and gene-specific codon usage regulate cell cycle progression during the DNA damage response

S-phase and DNA damage promote increased ribonucleotide reductase (RNR) activity. Translation of RNR1 has been linked to the wobble uridine modifying enzyme tRNA methyltransferase 9 (Trm9). We predicted that changes in tRNA modification would translationally regulate RNR1 after DNA damage to promote cell cycle progression. In support, we demonstrate that the Trm9-dependent tRNA modification 5-methoxycarbonylmethyluridine (mcm?U) is increased in hydroxyurea (HU)-induced S-phase cells, relative to G? and G?, and that mcm?U is one of 16 tRNA modifications whose levels oscillate during the cell cycle. Codon-reporter data matches the mcm?U increase to Trm9 and the efficient translation of AGA codons and RNR1. Further, we show that in trm9Δ cells reduced Rnr1 protein levels cause delayed transition into S-phase after damage. Codon re-engineering of RNR1 increased the number of trm9Δ cells that have transitioned into S-phase 1 h after DNA damage and that have increased Rnr1 protein levels, similar to that of wild-type cells expressing native RNR1. Our data supports a model in which codon usage and tRNA modification are regulatory components of the DNA damage response, with both playing vital roles in cell cycle progression.     

Sequential folding of transfer RNA: A nuclear magnetic resonance study of successively longer tRNA fragments with a common 5′ end

Most folding studies on proteins and nucleic acids have been addressed to the transition between the folded and unfolded states of an intact molecule, where an entire residue sequence is present during the folding event. However, since these polymers are synthesized sequentially from one terminus to the other in vivo, their folding pathways may be influenced greatly by the sequential appearance of the residues as a function of time.The three-dimensional structure of yeast tRNA^Phe in the crystalline state is correlated with 360 MHz proton nuclear magnetic resonances from three fragments plus an intact molecule of the tRNA that share a common 5′ end and are in a solution condition similar to that of the crystal structure. This has allowed identification of folded structures present in the fragments and presumably present in the growing tRNA molecule as it is being synthesized from the 5′ end. The experiments show that only the correct stems are formed in the fragments; no additional or competing helical region is produced. This suggests that in the biosynthesis of this tRNA, correct folding of helical stems occurs before the entire molecule is formed. Further, some of the tertiary interactions (hydrogen bonds) found in the crystal structure are also probably present before the synthesis is completed. These findings are generalized to consider the precursor of the tRNA as well as other tRNAs.