Diversity of human tRNA genes from the 1000-genomes project

Pathways to disease from natural variations in human cytoplasmic tRNAs

Perfectly accurate translation of mRNA into protein is not a prerequisite for life. Resulting from errors in protein synthesis, mistranslation occurs in all cells, including human cells. The human genome encodes >600 tRNA genes, providing both the raw material for genetic variation and a buffer to ensure that resulting translation errors occur at tolerable levels. On the basis of data from the 1000 Genomes Project, we highlight the unanticipated prevalence of mistranslating tRNA variants in the human population and review studies on synthetic and natural tRNA mutations that cause mistranslation or de-regulate protein synthesis. Although mitochondrial tRNA variants are well known to drive human diseases, including developmental disorders, few studies have revealed a role for human cytoplasmic tRNA mutants in disease. In the context of the unexpectedly large number of tRNA variants in the human population, the emerging literature suggests that human diseases may be affected by natural tRNA variants that cause mistranslation or de-regulate tRNA expression and nucleotide modification. This review highlights examples relevant to genetic disorders, cancer, and neurodegeneration in which cytoplasmic tRNA variants directly cause or exacerbate disease and disease-linked phenotypes in cells, animal models, and humans. In the near future, tRNAs may be recognized as useful genetic markers to predict the onset or severity of human disease.

Translation regulation in skin cancer from a tRNA point of view

Protein synthesis is a central and dynamic process, frequently deregulated in cancer through aberrant activation or expression of translation initiation factors and tRNAs. The discovery of tRNA-derived fragments, a new class of abundant and, in some cases stress-induced, small Noncoding RNAs has perplexed the epigenomics landscape and highlights the emerging regulatory role of tRNAs in translation and beyond. Skin is the biggest organ in human body, which maintains homeostasis of its multilayers through regulatory networks that induce translational reprogramming, and modulate tRNA transcription, modification and fragmentation, in response to various stress signals, like UV irradiation. In this review, we summarize recent knowledge on the role of translation regulation and tRNA biology in the alarming prevalence of skin cancer.

Modifications and functional genomics of human transfer RNA

Transfer RNA (tRNA) is present at tens of millions of transcripts in a human cell and is the most abundant RNA in moles among all cellular RNAs. tRNA is also the most extensively modified RNA with, on an average, 13 modifications per molecule. The primary function of tRNA as the adaptor of amino acids and the genetic code in protein synthesis is well known. tRNA modifications play multi-faceted roles in decoding and other cellular processes. The abundance, modification, and aminoacylation (charging) levels of tRNAs contribute to mRNA decoding in ways that reflect the cell type and its environment; however, how these factors work together to maximize translation efficiency remains to be understood. tRNAs also interact with many proteins not involved in translation and this may coordinate translation activity and other processes in the cell. This review focuses on the modifications and the functional genomics of human tRNA and discusses future perspectives on the explorations of human tRNA biology.

Pharmacogenetics: An Important Part of Drug Development with A Focus on Its Application

Since the human genome project in 2003, the view of personalized medicine to improve diagnosis and cure diseases at the molecular level became more real. Sequencing the human genome brought some benefits in medicine such as early detection of diseases with a genetic predisposition, treating patients with rare diseases, the design of gene therapy and the understanding of pharmacogenetics in the metabolism of drugs. This review explains the concepts of pharmacogenetics, polymorphisms, mutations, variations, and alleles, and how this information has helped us better understand the metabolism of drugs. Multiple resources are presented to promote reducing the gap between scientists, physicians, and patients in understanding the use and benefits of pharmacogenetics. Some of the most common clinical examples of genetic variants and how pharmacogenetics was used to determine treatment options for patients having these variants were discussed. Finally, we evaluated some of the challenges of implementing pharmacogenetics in a clinical setting and proposed actions to be taken to make pharmacogenetics a standard diagnostic tool in personalized medicine.

Visualizing tRNA-dependent mistranslation in human cells

High-fidelity translation and a strictly accurate proteome were originally assumed as essential to life and cellular viability. Yet recent studies in bacteria and eukaryotic model organisms suggest that proteome-wide mistranslation can provide selective advantages and is tolerated in the cell at higher levels than previously thought (one error in 6.9 × 10^-4 in yeast) with a limited impact on phenotype. Previously, we selected a tRNA^Pro containing a single mutation that induces mistranslation with alanine at proline codons in yeast. Yeast tolerate the mistranslation by inducing a heat-shock response and through the action of the proteasome. Here we found a homologous human tRNA^Pro (G3:U70) mutant that is not aminoacylated with proline, but is an efficient alanine acceptor. In live human cells, we visualized mistranslation using a green fluorescent protein reporter that fluoresces in response to mistranslation at proline codons. In agreement with measurements in yeast, quantitation based on the GFP reporter suggested a mistranslation rate of up to 2-5% in HEK 293 cells. Our findings suggest a stress-dependent phenomenon where mistranslation levels increased during nutrient starvation. Human cells did not mount a detectable heat-shock response and tolerated this level of mistranslation without apparent impact on cell viability. Because humans encode ∼600 tRNA genes and the natural population has greater tRNA sequence diversity than previously appreciated, our data also demonstrate a cell-based screen with the potential to elucidate mutations in tRNAs that may contribute to or alleviate disease.

The La and related RNA-binding proteins (LARPs): structures, functions, and evolving perspectives

La was first identified as a polypeptide component of ribonucleic protein complexes targeted by antibodies in autoimmune patients and is now known to be a eukaryote cell-ubiquitous protein. Structure and function studies have shown that La binds to a common terminal motif, UUU-3'-OH, of nascent RNA polymerase III (RNAP III) transcripts and protects them from exonucleolytic decay. For precursor-tRNAs, the most diverse and abundant of these transcripts, La also functions as an RNA chaperone that helps to prevent their misfolding. Related to this, we review evidence that suggests that La and its link to RNAP III were significant in the great expansions of the tRNAomes that occurred in eukaryotes. Four families of La-related proteins (LARPs) emerged during eukaryotic evolution with specialized functions. We provide an overview of the high-resolution structural biology of La and LARPs. LARP7 family members most closely resemble La but function with a single RNAP III nuclear transcript, 7SK, or telomerase RNA. A cytoplasmic isoform of La protein as well as LARPs 6, 4, and 1 function in mRNA metabolism and translation in distinct but similar ways, sometimes with the poly(A)-binding protein, and in some cases by direct binding to poly(A)-RNA. New structures of LARP domains, some complexed with RNA, provide novel insights into the functional versatility of these proteins. We also consider LARPs in relation to ancestral La protein and potential retention of links to specific RNA-related pathways. One such link may be tRNA surveillance and codon usage by LARP-associated mRNAs. WIREs RNA 2017, 8:e1430. doi: 10.1002/wrna.1430 For further resources related to this article, please visit the WIREs website.

The emerging complexity of the tRNA world: mammalian tRNAs beyond protein synthesis

The discovery of the genetic code and tRNAs as decoders of the code transformed life science. However, after establishing the role of tRNAs in protein synthesis, the field moved to other parts of the RNA world. Now, tRNA research is blooming again, with demonstration of the involvement of tRNAs in various other pathways beyond translation and in adapting translation to environmental cues. These roles are linked to the presence of tRNA sequence variants known as isoacceptors and isodecoders, various tRNA base modifications, the versatility of protein binding partners and tRNA fragmentation events, all of which collectively create an incalculable complexity. This complexity provides a vast repertoire of tRNA species that can serve various functions in cellular homeostasis and in adaptation of cellular functions to changing environments, and it likely arose from the fundamental role of RNAs in early evolution.

Determination of tRNA aminoacylation levels by high-throughput sequencing

Transfer RNA (tRNA) decodes mRNA codons when aminoacylated (charged) with an amino acid at its 3' end. Charged tRNAs turn over rapidly in cells, and variations in charged tRNA fractions are known to be a useful parameter in cellular responses to stress. tRNA charging fractions can be measured for individual tRNA species using acid denaturing gels, or comparatively at the genome level using microarrays. These hybridization-based approaches cannot be used for high resolution analysis of mammalian tRNAs due to their large sequence diversity. Here we develop a high-throughput sequencing method that enables accurate determination of charged tRNA fractions at single-base resolution (Charged DM-tRNA-seq). Our method takes advantage of the recently developed DM-tRNA-seq method, but includes additional chemical steps that specifically remove the 3'A residue in uncharged tRNA. Charging fraction is obtained by counting the fraction of A-ending reads versus A+C-ending reads for each tRNA species in the same sequencing reaction. In HEK293T cells, most cytosolic tRNAs are charged at >80% levels, whereas tRNASer and tRNAThr are charged at lower levels. These low charging levels were validated using acid denaturing gels. Our method should be widely applicable for investigations of tRNA charging as a parameter in biological regulation.

Genetic selection for mistranslation rescues a defective co-chaperone in yeast

Despite the general requirement for translation fidelity, mistranslation can be an adaptive response. We selected spontaneous second site mutations that suppress the stress sensitivity caused by a Saccharomyces cerevisiae tti2 allele with a Leu to Pro mutation at residue 187, identifying a single nucleotide mutation at the same position (C70U) in four tRNA^Pro_UGG genes. Linkage analysis and suppression by SUF9_G3:U70 expressed from a centromeric plasmid confirmed the causative nature of the suppressor mutation. Since the mutation incorporates the G3:U70 identity element for alanyl-tRNA synthetase into tRNA^Pro, we hypothesized that suppression results from mistranslation of Pro187 in Tti2_L187P as Ala. A strain expressing Tti2_L187A was not stress sensitive. In vitro, tRNA^Pro_UGG (C70U) was mis-aminoacylated with alanine by alanyl–tRNA synthetase, but was not a substrate for prolyl–tRNA synthetase. Mass spectrometry from protein expressed in vivo and a novel GFP reporter for mistranslation confirmed substitution of alanine for proline at a rate of ∼6%. Mistranslating cells expressing SUF9_G3:U70 induce a partial heat shock response but grow nearly identically to wild-type. Introducing the same G3:U70 mutation in SUF2 (tRNA^Pro_AGG) suppressed a second tti2 allele (tti2_L50P). We have thus identified a strategy that allows mistranslation to suppress deleterious missense Pro mutations in Tti2.

Role of the tRNA-Derived Small RNAs in Cancer: New Potential Biomarkers and Target for Therapy

Noncoding RNAs are untranslated RNA molecules that can be divided into two main types: infrastructural, including transfer RNAs (tRNAs) and ribosomal RNAs (rRNAs), and regulatory, including long ncRNAs (lncRNAs) and small ncRNAs (sRNA). Among small ncRNA, the role of microRNAs (miRNAs) and Piwi-interacting RNAs (piRNAs) in cancer is well documented. Recently, other small ncRNAs have been described. In particular, tRNA-derived small RNAs (tsRNA) have been found to be frequently dysregulated in cancer. Since tsRNAs can be considered unique sequences and are able to bind both Argonaute proteins (like miRNAs) and Piwi proteins (like piRNAs), their dysregulation could play a critical role in cancer by interfering with gene expression regulation at different levels. Like microRNAs, ts-53 (previously known as miR-3676) interacts with the 3'UTR of TCL1, therefore supporting a role for tsRNAs on the posttranscriptional regulation of gene expression. Like piRNAs, tsRNAs are produced as single-stranded molecules and can interact with DNA and histone methylation machinery, suggesting a role in the pretranscriptional regulation of gene expression. Herein, we describe the most recent findings about the role of tsRNAs in cancer.

Evolving Mistranslating tRNAs Through a Phenotypically Ambivalent Intermediate in Saccharomyces cerevisiae

The genetic code converts information from nucleic acid into protein. The genetic code was thought to be immutable, yet many examples in nature indicate that variations to the code provide a selective advantage. We used a sensitive selection system involving suppression of a deleterious allele (tti2-L187P) in Saccharomyces cerevisiae to detect mistranslation and identify mechanisms that allow genetic code evolution. Though tRNA^Ser containing a proline anticodon (UGG) is toxic, using our selection system we identified four tRNA^Ser_UGG variants, each with a single mutation, that mistranslate at a tolerable level. Mistranslating tRNA^Leu_UGG variants were also obtained, demonstrating the generality of the approach. We characterized two of the tRNA^Ser_UGG variants. One contained a G26A mutation, which reduced cell growth to 70% of the wild-type rate, induced a heat shock response, and was lost in the absence of selection. The reduced toxicity of tRNA^Ser_UGG-G26A is likely through increased turnover of the tRNA, as lack of methylation at G26 leads to degradation via the rapid tRNA decay pathway. The second tRNA^Ser_UGG variant, with a G9A mutation, had minimal effect on cell growth, was relatively stable in cells, and gave rise to less of a heat shock response. In vitro, the G9A mutation decreases aminoacylation and affects folding of the tRNA. Notably, the G26A and G9A mutations were phenotypically neutral in the context of an otherwise wild-type tRNA^Ser These experiments reveal a model for genetic code evolution in which tRNA anticodon mutations and mistranslation evolve through phenotypically ambivalent intermediates that reduce tRNA function.

Factors That Shape Eukaryotic tRNAomes: Processing, Modification and Anticodon-Codon Use

Transfer RNAs (tRNAs) contain sequence diversity beyond their anticodons and the large variety of nucleotide modifications found in all kingdoms of life. Some modifications stabilize structure and fit in the ribosome whereas those to the anticodon loop modulate messenger RNA (mRNA) decoding activity more directly. The identities of tRNAs with some universal anticodon loop modifications vary among distant and parallel species, likely to accommodate fine tuning for their translation systems. This plasticity in positions 34 (wobble) and 37 is reflected in codon use bias. Here, we review convergent evidence that suggest that expansion of the eukaryotic tRNAome was supported by its dedicated RNA polymerase III transcription system and coupling to the precursor-tRNA chaperone, La protein. We also review aspects of eukaryotic tRNAome evolution involving G34/A34 anticodon-sparing, relation to A34 modification to inosine, biased codon use and regulatory information in the redundancy (synonymous) component of the genetic code. We then review interdependent anticodon loop modifications involving position 37 in eukaryotes. This includes the eukaryote-specific tRNA modification, 3-methylcytidine-32 (m³C₃₂) and the responsible gene, TRM140 and homologs which were duplicated and subspecialized for isoacceptor-specific substrates and dependence on i⁶A₃₇ or t⁶A₃₇. The genetics of tRNA function is relevant to health directly and as disease modifiers.

An Enumerative Combinatorics Model for Fragmentation Patterns in RNA Sequencing Provides Insights into Nonuniformity of the Expected Fragment Starting-Point and Coverage Profile

RNA sequencing (RNA-seq) has emerged as the method of choice for measuring the expression of RNAs in a given cell population. In most RNA-seq technologies, sequencing the full length of RNA molecules requires fragmentation into smaller pieces. Unfortunately, the issue of nonuniform sequencing coverage across a genomic feature has been a concern in RNA-seq and is attributed to biases for certain fragments in RNA-seq library preparation and sequencing. To investigate the expected coverage obtained from fragmentation, we develop a simple fragmentation model that is independent of bias from the experimental method and is not specific to the transcript sequence. Essentially, we enumerate all configurations for maximal placement of a given fragment length, F, on transcript length, T, to represent every possible fragmentation pattern, from which we compute the expected coverage profile across a transcript. We extend this model to incorporate general empirical attributes such as read length, fragment length distribution, and number of molecules of the transcript. We further introduce the fragment starting-point, fragment coverage, and read coverage profiles. We find that the expected profiles are not uniform and that factors such as fragment length to transcript length ratio, read length to fragment length ratio, fragment length distribution, and number of molecules influence the variability of coverage across a transcript. Finally, we explore a potential application of the model where, with simulations, we show that it is possible to correctly estimate the transcript copy number for any transcript in the RNA-seq experiment.

Mistranslation: from adaptations to applications

BACKGROUND

The conservation of the genetic code indicates that there was a single origin, but like all genetic material, the cell's interpretation of the code is subject to evolutionary pressure. Single nucleotide variations in tRNA sequences can modulate codon assignments by altering codon-anticodon pairing or tRNA charging. Either can increase translation errors and even change the code. The frozen accident hypothesis argued that changes to the code would destabilize the proteome and reduce fitness. In studies of model organisms, mistranslation often acts as an adaptive response. These studies reveal evolutionary conserved mechanisms to maintain proteostasis even during high rates of mistranslation.

SCOPE OF REVIEW

This review discusses the evolutionary basis of altered genetic codes, how mistranslation is identified, and how deviations to the genetic code are exploited. We revisit early discoveries of genetic code deviations and provide examples of adaptive mistranslation events in nature. Lastly, we highlight innovations in synthetic biology to expand the genetic code.

MAJOR CONCLUSIONS

The genetic code is still evolving. Mistranslation increases proteomic diversity that enables cells to survive stress conditions or suppress a deleterious allele. Genetic code variants have been identified by genome and metagenome sequence analyses, suppressor genetics, and biochemical characterization.

GENERAL SIGNIFICANCE

Understanding the mechanisms of translation and genetic code deviations enables the design of new codes to produce novel proteins. Engineering the translation machinery and expanding the genetic code to incorporate non-canonical amino acids are valuable tools in synthetic biology that are impacting biomedical research. This article is part of a Special Issue entitled "Biochemistry of Synthetic Biology - Recent Developments" Guest Editor: Dr. Ilka Heinemann and Dr. Patrick O'Donoghue.

Genome-wide identification and characterization of tRNA-derived RNA fragments in land plants

The manuscript by Alves et al. entitled "Genome-wide identification and characterization of tRNA-derived RNA fragments in land plants" describes the identification and characterization of tRNAderived sRNA fragments in plants. By combining bioinformatic analysis and genetic and molecular approaches, we show that tRF biogenesis does not rely on canonical microRNA/siRNA processing machinery (i.e., independent of DICER-LIKE proteins). Moreover, we provide evidences that the Arabidopsis S-like Ribonuclease 1 (RNS1) might be involved in the biogenesis of tRFs. Detailed analyses showed that plant tRFs are sorted into different types of ARGONAUTE proteins and that they have potential target candidate genes. Our work advances the understanding of the tRF biology in plants by providing evidences that plant and animal tRFs shared common features and raising the hypothesis that an interplay between tRFs and other sRNAs might be important to fine-tune gene expression and protein biosynthesis in plant cells. Small RNA (sRNA) fragments derived from tRNAs (3'-loop, 5'-loop, anti-codon loop), named tRFs, have been reported in several organisms, including humans and plants. Although they may interfere with gene expression, their biogenesis and biological functions in plants remain poorly understood. Here, we capitalized on small RNA sequencing data from distinct species such as Arabidopsis thaliana, Oryza sativa, and Physcomitrella patens to examine the diversity of plant tRFs and provide insight into their properties. In silico analyzes of 19 to 25-nt tRFs derived from 5' (tRF-5s) and 3'CCA (tRF-3s) tRNA loops in these three evolutionary distant species showed that they are conserved and their abundance did not correlate with the number of genomic copies of the parental tRNAs. Moreover, tRF-5 is the most abundant variant in all three species. In silico and in vivo expression analyses unraveled differential accumulation of tRFs in Arabidopsis tissues/organs, suggesting that they are not byproducts of tRNA degradation. We also verified that the biogenesis of most Arabidopsis 19-25 nt tRF-5s and tRF-3s is not primarily dependent on DICER-LIKE proteins, though they seem to be associated with ARGONAUTE proteins and have few potential targets. Finally, we provide evidence that Arabidopsis ribonuclease RNS1 might be involved in the processing and/or degradation of tRFs. Our data support the notion that an interplay between tRFs and other sRNAs might be important to fine tune gene expression and protein biosynthesis in plant cells.

Interaction of tRNA with MEK2 in pancreatic cancer cells

Although the translational function of tRNA has long been established, extra translational functions of tRNA are still being discovered. We previously developed a computational method to systematically predict new tRNA-protein complexes and experimentally validated six candidate proteins, including the mitogen-activated protein kinase kinase 2 (MEK2), that interact with tRNA in HEK293T cells. However, consequences of the interaction between tRNA and these proteins remain to be elucidated. Here we tested the consequence of the interaction between tRNA and MEK2 in pancreatic cancer cell lines. We also generated disease and drug resistance-derived MEK2 mutants (Q60P, P128Q, S154F, E207K) to evaluate the function of the tRNA-MEK2 interaction. Our results demonstrate that tRNA interacts with the wild-type and mutant MEK2 in pancreatic cancer cells; furthermore, the MEK2 inhibitor U0126 significantly reduces the tRNA-MEK2 interaction. In addition, tRNA affects the catalytic activity of the wild type and mutant MEK2 proteins in different ways. Overall, our findings demonstrate the interaction of tRNA with MEK2 in pancreatic cancer cells and suggest that tRNA may impact MEK2 activity in cancer cells.

RNA Polymerase III Advances: Structural and tRNA Functional Views

RNA synthesis in eukaryotes is divided among three RNA polymerases (RNAPs). RNAP III transcribes hundreds of tRNA genes and fewer additional short RNA genes. We survey recent work on transcription by RNAP III including an atomic structure, mechanisms of action, interactions with chromatin and retroposons, and a conserved link between its activity and a tRNA modification that enhances mRNA decoding. Other new work suggests important mechanistic connections to oxidative stress, autoimmunity and cancer, embryonic stem cell pluripotency, and tissue-specific developmental effects. We consider that, for some of its complex functions, variation in RNAP III activity levels lead to nonuniform changes in tRNAs that can shift the translation profiles of key codon-biased mRNAs with resultant phenotypes or disease states.

A dual fluorescent reporter for the investigation of methionine mistranslation in live cells

In mammalian cells under oxidative stress, the methionyl-tRNA synthetase (MetRS) misacylates noncognate tRNAs at frequencies as high as 10% distributed among up to 28 tRNA species. Instead of being detrimental for the cell, misincorporation of methionine residues in the proteome reduces the risk of oxidative damage to proteins, which aids the oxidative stress response. tRNA microarrays have been essential for the detection of the full pattern of misacylated tRNAs, but have limited capacity to investigate the misacylation and mistranslation mechanisms in live cells. Here we develop a dual-fluorescence reporter to specifically measure methionine misincorporation at glutamic acid codons GAA and GAG via tRNA(Glu) mismethionylation in human cells. Our method relies on mutating a specific Met codon in the active site of the fluorescent protein mCherry to a Glu codon that renders mCherry nonfluorescent when translation follows the genetic code. Mistranslation utilizing mismethionylated tRNA(Glu) restores fluorescence in proportion to the amount of misacylated tRNA(Glu). This cellular approach works well for both transient transfection and established stable HEK293 lines. It is rapid, straightforward, and well suited for high-throughput activity analysis under a wide range of physiological conditions. As a proof of concept, we apply this method to characterize the effect of human tRNA(Glu) isodecoders on mistranslation and discuss the implications of our findings.

tRNA-Derived Small Non-Coding RNAs in Response to Ischemia Inhibit Angiogenesis

Ischemic injuries will lead to necrotic tissue damage, and post-ischemia angiogenesis plays critical roles in blood flow restoration and tissue recovery. Recently, several types of small RNAs have been reported to be involved in this process. In this study, we first generated a rat brain ischemic model to investigate the involvement of new types of small RNAs in ischemia. We utilized deep sequencing and bioinformatics analyses to demonstrate that the level of small RNA fragments derived from tRNAs strikingly increased in the ischemic rat brain. Among these sequences, tRNA^Val- and tRNA^Gly-derived small RNAs account for the most abundant segments. The up-regulation of tRNA^Val- and tRNA^Gly-derived fragments was verified through northern blot and quantitative PCR analyses. The levels of these two fragments also increased in a mouse hindlimb ischemia model and cellular hypoxia model. Importantly, up-regulation of the tRNA^Val- and tRNA^Gly-derived fragments in endothelial cells inhibited cell proliferation, migration and tube formation. Furthermore, we showed that these small RNAs are generated by angiogenin cleavage. Our results indicate that tRNA-derived fragments are involved in tissue ischemia, and we demonstrate for the first time that tRNA^Val- and tRNA^Gly-derived fragments inhibit angiogenesis by modulating the function of endothelial cells.

Small Non-coding Transfer RNA-Derived RNA Fragments (tRFs): Their Biogenesis, Function and Implication in Human Diseases

tRNA-derived RNA fragments (tRFs) are an emerging class of non-coding RNAs (ncRNAs). A growing number of reports have shown that tRFs are not random degradation products but are functional ncRNAs made of specific tRNA cleavage. They play regulatory roles in several biological contexts such as cancer, innate immunity, stress responses, and neurological disorders. In this review, we summarize the biogenesis and functions of tRFs.

Stabilization of eukaryotic ribosomal termination complexes by deacylated tRNA

Stabilization of the ribosomal complexes plays an important role in translational control. Mechanisms of ribosome stabilization have been studied in detail for initiation and elongation of eukaryotic translation, but almost nothing is known about stabilization of eukaryotic termination ribosomal complexes. Here, we present one of the mechanisms of fine-tuning of the translation termination process in eukaryotes. We show that certain deacylated tRNAs, remaining in the E site of the ribosome at the end of the elongation cycle, increase the stability of the termination and posttermination complexes. Moreover, only the part of eRF1 recognizing the stop codon is stabilized in the A site of the ribosome, and the stabilization is not dependent on the hydrolysis of peptidyl-tRNA. The determinants, defining this property of the tRNA, reside in the acceptor stem. It was demonstrated by site-directed mutagenesis of tRNA^Val and construction of a mini-helix structure identical to the acceptor stem of tRNA. The mechanism of this stabilization is different from the fixation of the unrotated state of the ribosome by CCA end of tRNA or by cycloheximide in the E site. Our data allow to reveal the possible functions of the isodecoder tRNAs in eukaryotes.

tRNA fragments in human health and disease

Transfer RNA (tRNA) is traditionally considered to be an adaptor molecule that helps ribosomes to decode messenger RNA (mRNA) and synthesize protein. Recent studies have demonstrated that tRNAs also serve as a major source of small non-coding RNAs that possess distinct and varied functions. These tRNA fragments are heterogeneous in size, nucleotide composition, biogenesis and function. Here we describe multiple roles that tRNA fragments play in cell physiology and discuss their relevance to human health and disease.

Nuclear and mitochondrial tRNA-lookalikes in the human genome

We are interested in identifying and characterizing loci of the human genome that harbor sequences resembling known mitochondrial and nuclear tRNAs. To this end, we used the known nuclear and mitochondrial tRNA genes (the “tRNA-Reference” set) to search for “tRNA-lookalikes” and found many such loci at different levels of sequence conservation. We find that the large majority of these tRNA-lookalikes resemble mitochondrial tRNAs and exhibit a skewed over-representation in favor of some mitochondrial anticodons. Our analysis shows that the tRNA-lookalikes have infiltrated specific chromosomes and are preferentially located in close proximity to known nuclear tRNAs (z-score ≤ −2.54, P-value ≤ 0.00394). Examination of the transcriptional potential of these tRNA-lookalike loci using public transcript annotations revealed that more than 20% of the lookalikes are transcribed as part of either known protein-coding pre-mRNAs, known lncRNAs, or known non-protein-coding RNAs, while public RNA-seq data perfectly agreed with the endpoints of tRNA-lookalikes. Interestingly, we found that tRNA-lookalikes are significantly depleted in known genetic variations associated with human health and disease whereas the known tRNAs are enriched in such variations. Lastly, a manual comparative analysis of the cloverleaf structure of several of the transcribed tRNA-lookalikes revealed no disruptive mutations suggesting the possibility that these loci give rise to functioning tRNA molecules.

Different types of secondary information in the genetic code

Whole-genome and functional analyses suggest a wealth of secondary or auxiliary genetic information (AGI) within the redundancy component of the genetic code. Although there are multiple aspects of biased codon use, we focus on two types of auxiliary information: codon-specific translational pauses that can be used by particular proteins toward their unique folding and biased codon patterns shared by groups of functionally related mRNAs with coordinate regulation. AGI is important to genetics in general and to human disease; here, we consider influences of its three major components, biased codon use itself, variations in the tRNAome, and anticodon modifications that distinguish synonymous decoding. AGI is plastic and can be used by different species to different extents, with tissue-specificity and in stress responses. Because AGI is species-specific, it is important to consider codon-sensitive experiments when using heterologous systems; for this we focus on the tRNA anticodon loop modification enzyme, CDKAL1, and its link to type 2 diabetes. Newly uncovered tRNAome variability among humans suggests roles in penetrance and as a genetic modifier and disease modifier. Development of experimental and bioinformatics methods are needed to uncover additional means of auxiliary genetic information.

A novel tRNA variable number tandem repeat at human chromosome 1q23.3 is implicated as a boundary element based on conservation of a CTCF motif in mouse

The human genome contains numerous large tandem repeats, many of which remain poorly characterized. Here we report a novel transfer RNA (tRNA) tandem repeat on human chromosome 1q23.3 that shows extensive copy number variation with 9-43 repeat units per allele and displays evidence of meiotic and mitotic instability. Each repeat unit consists of a 7.3 kb GC-rich sequence that binds the insulator protein CTCF and bears the chromatin hallmarks of a bivalent domain in human embryonic stem cells. A tRNA containing tandem repeat composed of at least three 7.6-kb GC-rich repeat units reside within a syntenic region of mouse chromosome 1. However, DNA sequence analysis reveals that, with the exception of the tRNA genes that account for less than 6% of a repeat unit, the remaining 7.2 kb is not conserved with the notable exception of a 24 base pair sequence corresponding to the CTCF binding site, suggesting an important role for this protein at the locus.