Predicting transfer RNA gene activity from sequence and genome context

The central role of transfer RNAs in mistranslation

Transfer RNAs (tRNA) are essential small non-coding RNAs that enable the translation of genomic information into proteins in all life forms. The principal function of tRNAs is to bring amino acid building blocks to the ribosomes for protein synthesis. In the ribosome, tRNAs interact with messenger RNA (mRNA) to mediate the incorporation of amino acids into a growing polypeptide chain following the rules of the genetic code. Accurate interpretation of the genetic code requires tRNAs to carry amino acids matching their anticodon identity and decode the correct codon on mRNAs. Errors in these steps cause the translation of codons with the wrong amino acids (mistranslation), compromising the accurate flow of information from DNA to proteins. Accumulation of mutant proteins due to mistranslation jeopardizes proteostasis and cellular viability. However, the concept of mistranslation is evolving, with increasing evidence indicating that mistranslation can be used as a mechanism for survival and acclimatization to environmental conditions. In this review, we discuss the central role of tRNAs in modulating translational fidelity through their dynamic and complex interplay with translation factors. We summarize recent discoveries of mistranslating tRNAs and describe the underlying molecular mechanisms and the specific conditions and environments that enable and promote mistranslation.

Mechanisms and Delivery of tRNA Therapeutics

Transfer ribonucleic acid (tRNA) therapeutics will provide personalized and mutation specific medicines to treat human genetic diseases for which no cures currently exist. The tRNAs are a family of adaptor molecules that interpret the nucleic acid sequences in our genes into the amino acid sequences of proteins that dictate cell function. Humans encode more than 600 tRNA genes. Interestingly, even healthy individuals contain some mutant tRNAs that make mistakes. Missense suppressor tRNAs insert the wrong amino acid in proteins, and nonsense suppressor tRNAs read through premature stop signals to generate full length proteins. Mutations that underlie many human diseases, including neurodegenerative diseases, cancers, and diverse rare genetic disorders, result from missense or nonsense mutations. Thus, specific tRNA variants can be strategically deployed as therapeutic agents to correct genetic defects. We review the mechanisms of tRNA therapeutic activity, the nature of the therapeutic window for nonsense and missense suppression as well as wild-type tRNA supplementation. We discuss the challenges and promises of delivering tRNAs as synthetic RNAs or as gene therapies. Together, tRNA medicines will provide novel treatments for common and rare genetic diseases in humans.

Genome-Wide Computational Prediction and Analysis of Noncoding RNAs in Oleidesulfovibrio alaskensis G20

Noncoding RNAs (ncRNAs) play key roles in the regulation of important pathways, including cellular growth, stress management, signaling, and biofilm formation. Sulfate-reducing bacteria (SRB) contribute to huge economic losses causing microbial-induced corrosion through biofilms on metal surfaces. To effectively combat the challenges posed by SRB, it is essential to understand their molecular mechanisms of biofilm formation. This study aimed to identify ncRNAs in the genome of a model SRB, Oleidesulfovibrio alaskensis G20 (OA G20). Three in silico approaches revealed genome-wide distribution of 37 ncRNAs excluding tRNAs in the OA G20. These ncRNAs belonged to 18 different Rfam families. This study identified riboswitches, sRNAs, RNP, and SRP. The analysis revealed that these ncRNAs could play key roles in the regulation of several pathways of biosynthesis and transport involved in biofilm formation by OA G20. Three sRNAs, Pseudomonas P10, Hammerhead type II, and sX4, which were found in OA G20, are rare and their roles have not been determined in SRB. These results suggest that applying various computational methods could enrich the results and lead to the discovery of additional novel ncRNAs, which could lead to understanding the "rules of life of OA G20" during biofilm formation.

Large-scale duplication events underpin population-level flexibility in tRNA gene copy number in Pseudomonas fluorescens SBW25

The complement of tRNA genes within a genome is typically considered to be a (relatively) stable characteristic of an organism. Here, we demonstrate that bacterial tRNA gene set composition can be more flexible than previously appreciated, particularly regarding tRNA gene copy number. We report the high-rate occurrence of spontaneous, large-scale, tandem duplication events in laboratory populations of the bacterium Pseudomonas fluorescens SBW25. The identified duplications are up to ∼1 Mb in size (∼15% of the wildtype genome) and are predicted to change the copy number of up to 917 genes, including several tRNA genes. The observed duplications are inherently unstable: they occur, and are subsequently lost, at extremely high rates. We propose that this unusually plastic type of mutation provides a mechanism by which tRNA gene set diversity can be rapidly generated, while simultaneously preserving the underlying tRNA gene set in the absence of continued selection. That is, if a tRNA set variant provides no fitness advantage, then high-rate segregation of the duplication ensures the maintenance of the original tRNA gene set. However, if a tRNA gene set variant is beneficial, the underlying duplication fragment(s) may persist for longer and provide raw material for further, more stable, evolutionary change.

A POLR3B-variant reveals a Pol III transcriptome response dependent on La protein/SSB

RNA polymerase III (Pol III, POLR3) synthesizes tRNAs and other small non-coding RNAs. Human POLR3 pathogenic variants cause a range of developmental disorders, recapitulated in part by mouse models, yet some aspects of POLR3 deficiency have not been explored. We characterized a human POLR3B:c.1625A>G;p.(Asn542Ser) disease variant that was found to cause mis-splicing of POLR3B. Genome-edited POLR3B1625A>G HEK293 cells acquired the mis-splicing with decreases in multiple POLR3 subunits and TFIIIB, although display auto-upregulation of the Pol III termination-reinitiation subunit POLR3E. La protein was increased relative to its abundant pre-tRNA ligands which bind via their U(n)U-3'-termini. Assays for cellular transcription revealed greater deficiencies for tRNA genes bearing terminators comprised of 4Ts than of ≥5Ts. La-knockdown decreased Pol III ncRNA expression unlinked to RNA stability. Consistent with these effects, small-RNAseq showed that POLR3B1625A>G and patient fibroblasts express more tRNA fragments (tRFs) derived from pre-tRNA 3'-trailers (tRF-1) than from mature-tRFs, and higher levels of multiple miRNAs, relative to control cells. The data indicate that decreased levels of Pol III transcripts can lead to functional excess of La protein which reshapes small ncRNA profiles revealing new depth in the Pol III system. Finally, patient cell RNA analysis uncovered a strategy for tRF-1/tRF-3 as POLR3-deficiency biomarkers.

tRNA therapeutics for genetic diseases

Transfer RNAs (tRNAs) have a crucial role in protein synthesis, and in recent years, their therapeutic potential for the treatment of genetic diseases - primarily those associated with a mutation altering mRNA translation - has gained significant attention. Engineering tRNAs to readthrough nonsense mutation-associated premature termination of mRNA translation can restore protein synthesis and function. In addition, supplementation of natural tRNAs can counteract effects of missense mutations in proteins crucial for tRNA biogenesis and function in translation. This Review will present advances in the development of tRNA therapeutics with high activity and safety in vivo and discuss different formulation approaches for single or chronic treatment modalities. The field of tRNA therapeutics is still in its early stages, and a series of challenges related to tRNA efficacy and stability in vivo, delivery systems with tissue-specific tropism, and safe and efficient manufacturing need to be addressed.

Selective gene expression maintains human tRNA anticodon pools during differentiation

Transfer RNAs are essential for translating genetic information into proteins. The human genome contains hundreds of predicted tRNA genes, many in multiple copies. How their expression is regulated to control tRNA repertoires is unknown. Here we combined quantitative tRNA profiling and chromatin immunoprecipitation with sequencing to measure tRNA expression following the differentiation of human induced pluripotent stem cells into neuronal and cardiac cells. We find that tRNA transcript levels vary substantially, whereas tRNA anticodon pools, which govern decoding rates, are more stable among cell types. Mechanistically, RNA polymerase III transcribes a wide range of tRNA genes in human induced pluripotent stem cells but on differentiation becomes constrained to a subset we define as housekeeping tRNAs. This shift is mediated by decreased mTORC1 signalling, which activates the RNA polymerase III repressor MAF1. Our data explain how tRNA anticodon pools are buffered to maintain decoding speed across cell types and reveal that mTORC1 drives selective tRNA expression during differentiation.

Design, construction, and functional characterization of a tRNA neochromosome in yeast

Here, we report the design, construction, and characterization of a tRNA neochromosome, a designer chromosome that functions as an additional, de novo counterpart to the native complement of Saccharomyces cerevisiae. Intending to address one of the central design principles of the Sc2.0 project, the ∼190-kb tRNA neochromosome houses all 275 relocated nuclear tRNA genes. To maximize stability, the design incorporates orthogonal genetic elements from non-S. cerevisiae yeast species. Furthermore, the presence of 283 rox recombination sites enables an orthogonal tRNA SCRaMbLE system. Following construction in yeast, we obtained evidence of a potent selective force, manifesting as a spontaneous doubling in cell ploidy. Furthermore, tRNA sequencing, transcriptomics, proteomics, nucleosome mapping, replication profiling, FISH, and Hi-C were undertaken to investigate questions of tRNA neochromosome behavior and function. Its construction demonstrates the remarkable tractability of the yeast model and opens up opportunities to directly test hypotheses surrounding these essential non-coding RNAs.

Effect of mRNA/tRNA mutations on translation speed: Implications for human diseases

Recent discoveries establish tRNAs as central regulators of mRNA translation dynamics, and therefore cotranslational folding and function of the encoded protein. The tRNA pool, whose composition and abundance change in a cell- and tissue-dependent manner, is the main factor which determines mRNA translation velocity. In this review, we discuss a group of pathogenic mutations, in the coding sequences of either protein-coding genes or in tRNA genes, that alter mRNA translation dynamics. We also summarize advances in tRNA biology that have uncovered how variations in tRNA levels on account of genetic mutations affect protein folding and function, and thereby contribute to phenotypic diversity in clinical manifestations.

Selection of tRNA Genes in Human Breast Tumours Varies Substantially between Individuals

Abnormally elevated expression of tRNA is a common feature of breast tumours. Rather than a uniform increase in all tRNAs, some are deregulated more strongly than others. Elevation of particular tRNAs has been associated with poor prognosis for patients, and experimental models have demonstrated the ability of some tRNAs to promote proliferation or metastasis. Each tRNA isoacceptor is encoded redundantly by multiple genes, which are commonly dispersed across several chromosomes. An unanswered question is whether the consistently high expression of a tRNA in a cancer type reflects the consistent activation of the same members of a gene family, or whether different family members are activated from one patient to the next. To address this question, we interrogated ChIP-seq data to determine which tRNA genes were active in individual breast tumours. This revealed that distinct sets of tRNA genes become activated in individual cancers, whereas there is much less variation in the expression patterns of families. Several pathways have been described that are likely to contribute to increases in tRNA gene transcription in breast tumours, but none of these can adequately explain the observed variation in the choice of genes between tumours. Current models may therefore lack at least one level of regulation.

Translation and mRNA Stability Control

Messenger RNA (mRNA) stability and translational efficiency are two crucial aspects of the post-transcriptional process that profoundly impact protein production in a cell. While it is widely known that ribosomes produce proteins, studies during the past decade have surprisingly revealed that ribosomes also control mRNA stability in a codon-dependent manner, a process referred to as codon optimality. Therefore, codons, the three-nucleotide words read by the ribosome, have a potent effect on mRNA stability and provide cis-regulatory information that extends beyond the amino acids they encode. While the codon optimality molecular mechanism is still unclear, the translation elongation rate appears to trigger mRNA decay. Thus, transfer RNAs emerge as potential master gene regulators affecting mRNA stability. Furthermore, while few factors related to codon optimality have been identified in yeast, the orthologous genes in vertebrates do not necessary share the same functions. Here, we discuss codon optimality findings and gene regulation layers related to codon composition in different eukaryotic species.

Perseverance of protein homeostasis despite mistranslation of glycine codons with alanine

By linking amino acids to their codon assignments, transfer RNAs (tRNAs) are essential for protein synthesis and translation fidelity. Some human tRNA variants cause amino acid mis-incorporation at a codon or set of codons. We recently found that a naturally occurring tRNASer variant decodes phenylalanine codons with serine and inhibits protein synthesis. Here, we hypothesized that human tRNA variants that misread glycine (Gly) codons with alanine (Ala) will also disrupt protein homeostasis. The A3G mutation occurs naturally in tRNAGly variants (tRNAGlyCCC, tRNAGlyGCC) and creates an alanyl-tRNA synthetase (AlaRS) identity element (G3 : U70). Because AlaRS does not recognize the anticodon, the human tRNAAlaAGC G35C (tRNAAlaACC) variant may function similarly to mis-incorporate Ala at Gly codons. The tRNAGly and tRNAAla variants had no effect on protein synthesis in mammalian cells under normal growth conditions; however, tRNAGlyGCC A3G depressed protein synthesis in the context of proteasome inhibition. Mass spectrometry confirmed Ala mistranslation at multiple Gly codons caused by the tRNAGlyGCC A3G and tRNAAlaAGC G35C mutants, and in some cases, we observed multiple mistranslation events in the same peptide. The data reveal mistranslation of Ala at Gly codons and defects in protein homeostasis generated by natural human tRNA variants that are tolerated under normal conditions. This article is part of the theme issue 'Reactivity and mechanism in chemical and synthetic biology'.

TFIIIC-based chromatin insulators through eukaryotic evolution

Eukaryotic chromosomes are divided into domains with distinct structural and functional properties, such as differing levels of chromatin compaction and gene transcription. Domains of relatively compact chromatin and minimal transcription are termed heterochromatic, whereas euchromatin is more open and actively transcribed. Insulators separate these domains and maintain their distinct features. Disruption of insulators can cause diseases such as cancer. Many insulators contain tRNA genes (tDNAs), examples of which have been shown to block the spread of activating or silencing activities. This characteristic of specific tDNAs is conserved through evolution, such that human tDNAs can serve as barriers to the spread of silencing in fission yeast. Here we demonstrate that tDNAs from the methylotrophic fungus Pichia pastoris can function effectively as insulators in distantly-related budding yeast. Key to the function of tDNAs as insulators is TFIIIC, a transcription factor that is also required for their expression. TFIIIC binds additional loci besides tDNAs, some of which have insulator activity. Although the mechanistic basis of TFIIIC-based insulation has been studied extensively in yeast, it is largely uncharacterized in metazoa. Utilising publicly-available genome-wide ChIP-seq data, we consider the extent to which mechanisms conserved from yeast to man may suffice to allow efficient insulation by TFIIIC in the more challenging chromatin environments of metazoa and suggest features that may have been acquired during evolution to cope with new challenges. We demonstrate the widespread presence at human tDNAs of USF1, a transcription factor with well-established barrier activity in vertebrates. We predict that tDNA-based insulators in higher organisms have evolved through incorporation of modules, such as binding sites for factors like USF1 and CTCF that are absent from yeasts, thereby strengthening function and providing opportunities for regulation between cell types.

Eukaryotic tRNA sequences present conserved and amino acid-specific structural signatures

Metazoan organisms have many tRNA genes responsible for decoding amino acids. The set of all tRNA genes can be grouped in sets of common amino acids and isoacceptor tRNAs that are aminoacylated by corresponding aminoacyl-tRNA synthetases. Analysis of tRNA alignments shows that, despite the high number of tRNA genes, specific tRNA sequence motifs are highly conserved across multicellular eukaryotes. The conservation often extends throughout the isoacceptors and isodecoders with, in some cases, two sets of conserved isodecoders. This study is focused on non-Watson–Crick base pairs in the helical stems, especially GoU pairs. Each of the four helical stems may contain one or more conserved GoU pairs. Some are amino acid specific and could represent identity elements for the cognate aminoacyl tRNA synthetases. Other GoU pairs are found in more than a single amino acid and could be critical for native folding of the tRNAs. Interestingly, some GoU pairs are anticodon-specific, and others are found in phylogenetically-specific clades. Although the distribution of conservation likely reflects a balance between accommodating isotype-specific functions as well as those shared by all tRNAs essential for ribosomal translation, such conservations may indicate the existence of specialized tRNAs for specific translation targets, cellular conditions, or alternative functions.

Epigenetic regulation of human non-coding RNA gene transcription

Recent investigations on the non-protein-coding transcriptome of human cells have revealed previously hidden layers of gene regulation relying on regulatory non-protein-coding (nc) RNAs, including the widespread ncRNA-dependent regulation of epigenetic chromatin states and of mRNA translation and stability. However, despite its centrality, the epigenetic regulation of ncRNA genes has received relatively little attention. In this mini-review, we attempt to provide a synthetic account of recent literature suggesting an unexpected complexity in chromatin-dependent regulation of ncRNA gene transcription by the three human nuclear RNA polymerases. Emerging common features, like the heterogeneity of chromatin states within ncRNA multigene families and their influence on 3D genome organization, point to unexplored issues whose investigation could lead to a better understanding of the whole human epigenomic network.

Dynamic changes in tRNA modifications and abundance during T cell activation

The tRNA pool determines the efficiency, throughput, and accuracy of translation. Previous studies have identified dynamic changes in the tRNA (transfer RNA) supply and mRNA (messenger RNA) demand during cancerous proliferation. Yet dynamic changes may also occur during physiologically normal proliferation, and these are less well characterized. We examined the tRNA and mRNA pools of T cells during their vigorous proliferation and differentiation upon triggering their antigen receptor. We observed a global signature of switch in demand for codons at the early proliferation phase of the response, accompanied by corresponding changes in tRNA expression levels. In the later phase, upon differentiation, the response of the tRNA pool relaxed back to the basal level, potentially restraining excessive proliferation. Sequencing of tRNAs allowed us to evaluate their diverse base-modifications. We found that two types of tRNA modifications, wybutosine and ms²t6A, are reduced dramatically during T cell activation. These modifications occur in the anticodon loops of two tRNAs that decode "slippery codons," which are prone to ribosomal frameshifting. Attenuation of these frameshift-protective modifications is expected to increase the potential for proteome-wide frameshifting during T cell proliferation. Indeed, human cell lines deleted of a wybutosine writer showed increased ribosomal frameshifting, as detected with an HIV gag-pol frameshifting site reporter. These results may explain HIV's specific tropism toward proliferating T cells since it requires ribosomal frameshift exactly on the corresponding codon for infection. The changes in tRNA expression and modifications uncover a layer of translation regulation during T cell proliferation and expose a potential tradeoff between cellular growth and translation fidelity.

Endogenous tRNA-derived small RNA (tRF3-Thr-AGT) inhibits ZBP1/NLRP3 pathway-mediated cell pyroptosis to attenuate acute pancreatitis (AP)

Endogenous transfer RNA‐derived small RNAs (tsRNAs) are newly identified RNAs that are closely associated with the pathogenesis of multiple diseases, but the involvement of tsRNAs in regulating acute pancreatitis (AP) development has not been reported. In this study, we screened out a novel tsRNA, tRF3‐Thr‐AGT, that was aberrantly downregulated in the acinar cell line AR42J treated with sodium taurocholate (STC) and the pancreatic tissues of STC‐induced AP rat models. In addition, STC treatment suppressed cell viability, induced pyroptotic cell death and cellular inflammation in AP models in vitro and in vivo. Overexpression of tRF3‐Thr‐AGT partially reversed STC‐induced detrimental effects on the AR42J cells. Next, Z‐DNA‐binding protein 1 (ZBP1) was identified as the downstream target of tRF3‐Thr‐AGT. Interestingly, upregulation of tRF3‐Thr‐AGT suppressed NOD‐like receptor protein 3 (NLRP3)‐mediated pyroptotic cell death in STC‐treated AR42J cells via degrading ZBP1. Moreover, the effects of tRF3‐Thr‐AGT overexpression on cell viability and inflammation in AR42J cells were abrogated by upregulating ZBP1 and NLRP3. Collectively, our data indicated that tRF3‐Thr‐AGT suppressed ZBP1 expressions to restrain NLRP3‐mediated pyroptotic cell death and inflammation in AP models. This study, for the first time, identified the role and potential underlying mechanisms by which tRF3‐Thr‐AGT regulated AP pathogenesis.

tRNAscan-SE 2.0: improved detection and functional classification of transfer RNA genes

tRNAscan-SE has been widely used for transfer RNA (tRNA) gene prediction for over twenty years, developed just as the first genomes were decoded. With the massive increase in quantity and phylogenetic diversity of genomes, the accurate detection and functional prediction of tRNAs has become more challenging. Utilizing a vastly larger training set, we created nearly one hundred specialized isotype- and clade-specific models, greatly improving tRNAscan-SE’s ability to identify and classify both typical and atypical tRNAs. We employ a new comparative multi-model strategy where predicted tRNAs are scored against a full set of isotype-specific covariance models, allowing functional prediction based on both the anticodon and the highest-scoring isotype model. Comparative model scoring has also enhanced the program's ability to detect tRNA-derived SINEs and other likely pseudogenes. For the first time, tRNAscan-SE also includes fast and highly accurate detection of mitochondrial tRNAs using newly developed models. Overall, tRNA detection sensitivity and specificity is improved for all isotypes, particularly those utilizing specialized models for selenocysteine and the three subtypes of tRNA genes encoding a CAU anticodon. These enhancements will provide researchers with more accurate and detailed tRNA annotation for a wider variety of tRNAs, and may direct attention to tRNAs with novel traits.

The genomic loci of specific human tRNA genes exhibit ageing-related DNA hypermethylation

The epigenome has been shown to deteriorate with age, potentially impacting on ageing-related disease. tRNA, while arising from only ˜46 kb (<0.002% genome), is the second most abundant cellular transcript. tRNAs also control metabolic processes known to affect ageing, through core translational and additional regulatory roles. Here, we interrogate the DNA methylation state of the genomic loci of human tRNA. We identify a genomic enrichment for age-related DNA hypermethylation at tRNA loci. Analysis in 4,350 MeDIP-seq peripheral-blood DNA methylomes (16-82 years), identifies 44 and 21 hypermethylating specific tRNAs at study-and genome-wide significance, respectively, contrasting with none hypomethylating. Validation and replication (450k array and independent targeted Bisuphite-sequencing) supported the hypermethylation of this functional unit. Tissue-specificity is a significant driver, although the strongest consistent signals, also independent of major cell-type change, occur in tRNA-iMet-CAT-1-4 and tRNA-Ser-AGA-2-6. This study presents a comprehensive evaluation of the genomic DNA methylation state of human tRNA genes and reveals a discreet hypermethylation with advancing age.

Social Networking of Quasi-Species Consortia drive Virolution via Persistence

The emergence of cooperative quasi-species consortia (QS-C) thinking from the more accepted quasispecies equations of Manfred Eigen, provides a conceptual foundation from which concerted action of RNA agents can now be understood. As group membership becomes a basic criteria for the emergence of living systems, we also start to understand why the history and context of social RNA networks become crucial for survival and function. History and context of social RNA networks also lead to the emergence of a natural genetic code. Indeed, this QS-C thinking can also provide us with a transition point between the chemical world of RNA replicators and the living world of RNA agents that actively differentiate self from non-self and generate group identity with membership roles. Importantly the social force of a consortia to solve complex, multilevel problems also depend on using opposing and minority functions. The consortial action of social networks of RNA stem-loops subsequently lead to the evolution of cellular organisms representing a tree of life.

Tissue-specific reprogramming of host tRNA transcriptome by the microbiome

Transfer RNAs (tRNAs) are essential for translation, and tRNA expression and modifications are regulated by many factors. However, the interplay between the microbiome and host tRNA profiles through host-microbiome interactions has not been explored. In this study, we investigated host-microbiome interactions via the tRNA profiling of four tissue types from germ-free and specific pathogen-free mice. Our analyses reveal that cytosolic and mitochondrial tRNA expression and tRNA modifications in the host are reprogrammed in a tissue-specific and microbiome-dependent manner. In terms of tRNA expression, the intestines and brains are more sensitive to the influence of the microbiome than the livers and kidneys. In terms of tRNA modifications, cytosolic tRNAs show more obvious changes in the livers and kidneys in the presence of the microbiome. Our findings reveal a previously unexplored relationship among the microbiome, tRNA abundance, and epitranscriptome in a mammalian host.