tRNA biology in the omics era: Stress signalling dynamics and cancer progression

Advances in methods for tRNA sequencing and quantification

In the past decade tRNA sequencing (tRNA-seq) has attracted considerable attention as an important tool for the development of novel approaches to quantify highly modified tRNA species and to propel tRNA research aimed at understanding the cellular physiology and disease and development of tRNA-based therapeutics. Many methods are available to quantify tRNA abundance while accounting for modifications and tRNA charging/acylation. Advances in both library preparation methods and bioinformatic workflows have enabled developments in next-generation sequencing (NGS) workflows. Other approaches forgo NGS applications in favor of hybridization-based approaches. In this review we provide a brief comparative overview of various tRNA quantification approaches, focusing on the advantages and disadvantages of these methods, which together facilitate reliable tRNA quantification.

PtRNAdb: a web resource of plant tRNA genes from a wide range of plant species

tRNA, as well as their derived products such as short interspersed nuclear elements (SINEs), pseudogenes, and transfer RNA (tRNA)-derived fragments (tRFs), have now been shown to be vital for cellular life, functioning, and adaptation during different stress conditions in all diverse life forms. In this study, we have developed PtRNAdb (www.nipgr.ac.in/PtRNAdb), a plant-exclusive tRNA database containing 113,849 tRNA gene sequences from phylogenetically diverse plant species. We have analyzed a total of 106 nuclear, 89 plastidial, and 38 mitochondrial genomes of plants by the tRNAscan-SE software package, and after careful curation of the output data, we integrated the data and developed this database. The information about the tRNA gene sequences obtained was further enriched with a consensus sequence-based study of tRNA genes based on their isoacceptors and isodecoders. We have also built covariance models based on the isoacceptors and isodecoders of all the tRNA sequences using the infernal tool. The user can also perform BLAST not only against PtRNAdb entries but also against all the tRNA sequences stored in the PlantRNA database and annotated tRNA genes across the plant kingdom available at NCBI. This resource is believed to be of high utility for plant researchers as well as molecular biologists to carry out further exploration of the plant tRNAome on a wider spectrum, as well as for performing comparative and evolutionary studies related to tRNAs, and their derivatives across all domains of life.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03255-7.

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.

A dual role of human tRNA methyltransferase hTrmt13 in regulating translation and transcription

Since numerous RNAs and RBPs prevalently localize to active chromatin regions, many RNA-binding proteins (RBPs) may be potential transcriptional regulators. RBPs are generally thought to regulate transcription via noncoding RNAs. Here, we describe a distinct, dual mechanism of transcriptional regulation by the previously uncharacterized tRNA-modifying enzyme, hTrmt13. On one hand, hTrmt13 acts in the cytoplasm to catalyze 2'-O-methylation of tRNAs, thus regulating translation in a manner depending on its tRNA-modification activity. On the other hand, nucleus-localized hTrmt13 directly binds DNA as a transcriptional co-activator of key epithelial-mesenchymal transition factors, thereby promoting cell migration independent of tRNA-modification activity. These dual functions of hTrmt13 are mutually exclusive, as it can bind either DNA or tRNA through its CHHC zinc finger domain. Finally, we find that hTrmt13 expression is tightly associated with poor prognosis and survival in diverse cancer patients. Our discovery of the noncatalytic roles of an RNA-modifying enzyme provides a new perspective for understanding epitranscriptomic regulation.

Increased expression of tryptophan and tyrosine tRNAs elevates stop codon readthrough of reporter systems in human cell lines

Regulation of translation via stop codon readthrough (SC-RT) expands not only tissue-specific but also viral proteomes in humans and, therefore, represents an important subject of study. Understanding this mechanism and all involved players is critical also from a point of view of prospective medical therapies of hereditary diseases caused by a premature termination codon. tRNAs were considered for a long time to be just passive players delivering amino acid residues according to the genetic code to ribosomes without any active regulatory roles. In contrast, our recent yeast work identified several endogenous tRNAs implicated in the regulation of SC-RT. Swiftly emerging studies of human tRNA-ome also advocate that tRNAs have unprecedented regulatory potential. Here, we developed a universal U6 promotor-based system expressing various human endogenous tRNA iso-decoders to study consequences of their increased dosage on SC-RT employing various reporter systems in vivo. This system combined with siRNA-mediated downregulations of selected aminoacyl-tRNA synthetases demonstrated that changing levels of human tryptophan and tyrosine tRNAs do modulate efficiency of SC-RT. Overall, our results suggest that tissue-to-tissue specific levels of selected near-cognate tRNAs may have a vital potential to fine-tune the final landscape of the human proteome, as well as that of its viral pathogens.

On the Track of the Missing tRNA Genes: A Source of Non-Canonical Functions?

Cellular tRNAs appear today as a diverse population of informative macromolecules with conserved general elements ensuring essential common functions and different and distinctive features securing specific interactions and activities. Their differential expression and the variety of post-transcriptional modifications they are subject to, lead to the existence of complex repertoires of tRNA populations adjusted to defined cellular states. Despite the tRNA-coding genes redundancy in prokaryote and eukaryote genomes, it is surprising to note the absence of genes coding specific translational-active isoacceptors throughout the phylogeny. Through the analysis of different releases of tRNA databases, this review aims to provide a general summary about those “missing tRNA genes.” This absence refers to both tRNAs that are not encoded in the genome, as well as others that show critical sequence variations that would prevent their activity as canonical translation adaptor molecules. Notably, while a group of genes are universally missing, others are absent in particular kingdoms. Functional information available allows to hypothesize that the exclusion of isodecoding molecules would be linked to: 1) reduce ambiguities of signals that define the specificity of the interactions in which the tRNAs are involved; 2) ensure the adaptation of the translational apparatus to the cellular state; 3) divert particular tRNA variants from ribosomal protein synthesis to other cellular functions. This leads to consider the “missing tRNA genes” as a source of putative non-canonical tRNA functions and to broaden the concept of adapter molecules in ribosomal-dependent protein synthesis.

tRNA Biology in the Pathogenesis of Diabetes: Role of Genetic and Environmental Factors

The global rise in type 2 diabetes results from a combination of genetic predisposition with environmental assaults that negatively affect insulin action in peripheral tissues and impair pancreatic β-cell function and survival. Nongenetic heritability of metabolic traits may be an important contributor to the diabetes epidemic. Transfer RNAs (tRNAs) are noncoding RNA molecules that play a crucial role in protein synthesis. tRNAs also have noncanonical functions through which they control a variety of biological processes. Genetic and environmental effects on tRNAs have emerged as novel contributors to the pathogenesis of diabetes. Indeed, altered tRNA aminoacylation, modification, and fragmentation are associated with β-cell failure, obesity, and insulin resistance. Moreover, diet-induced tRNA fragments have been linked with intergenerational inheritance of metabolic traits. Here, we provide a comprehensive review of how perturbations in tRNA biology play a role in the pathogenesis of monogenic and type 2 diabetes.

RNA polymerase III transcription as a disease factor

In this review, Yeganeh et al. summarize different human diseases that have been linked to defects in the Pol III transcription apparatus or to Pol III products imbalance and discuss the possible underlying mechanisms.

RNA polymerase (Pol) III is responsible for transcription of different noncoding genes in eukaryotic cells, whose RNA products have well-defined functions in translation and other biological processes for some, and functions that remain to be defined for others. For all of them, however, new functions are being described. For example, Pol III products have been reported to regulate certain proteins such as protein kinase R (PKR) by direct association, to constitute the source of very short RNAs with regulatory roles in gene expression, or to control microRNA levels by sequestration. Consistent with these many functions, deregulation of Pol III transcribed genes is associated with a large variety of human disorders. Here we review different human diseases that have been linked to defects in the Pol III transcription apparatus or to Pol III products imbalance and discuss the possible underlying mechanisms.

Role of the RNA-binding protein La in cancer pathobiology

RNA-binding proteins are important regulators of RNA metabolism and are of critical importance in all steps of the gene expression cascade. The role of aberrantly expressed RBPs in human disease is an exciting research field and the potential application of RBPs as a therapeutic target or a diagnostic marker represents a fast-growing area of research.Aberrant overexpression of the human RNA-binding protein La has been found in various cancer entities including lung, cervical, head and neck, and chronic myelogenous leukaemia. Cancer-associated La protein supports tumour-promoting processes such as proliferation, mobility, invasiveness and tumour growth. Moreover, the La protein maintains the survival of cancer cells by supporting an anti-apoptotic state that may cause resistance to chemotherapeutic therapy.The human La protein represents a multifunctional post-translationally modified RNA-binding protein with RNA chaperone activity that promotes processing of non-coding precursor RNAs but also stimulates the translation of selective messenger RNAs encoding tumour-promoting and anti-apoptotic factors. In our model, La facilitates the expression of those factors and helps cancer cells to cope with cellular stress. In contrast to oncogenes, able to initiate tumorigenesis, we postulate that the aberrantly elevated expression of the human La protein contributes to the non-oncogenic addiction of cancer cells. In this review, we summarize the current understanding about the implications of the RNA-binding protein La in cancer progression and therapeutic resistance. The concept of exploiting the RBP La as a cancer drug target will be discussed.

Robust landscapes of ribosome dwell times and aminoacyl-tRNAs in response to nutrient stress in liver

Translation depends on messenger RNA (mRNA)-specific initiation, elongation, and termination rates. While translation elongation is well studied in bacteria and yeast, less is known in higher eukaryotes. Here we combined ribosome and transfer RNA (tRNA) profiling to investigate the relations between translation elongation rates, (aminoacyl-) tRNA levels, and codon usage in mammals. We modeled codon-specific ribosome dwell times from ribosome profiling, considering codon pair interactions between ribosome sites. In mouse liver, the model revealed site- and codon-specific dwell times that differed from those in yeast, as well as pairs of adjacent codons in the P and A site that markedly slow down or speed up elongation. While translation efficiencies vary across diurnal time and feeding regimen, codon dwell times were highly stable and conserved in human. Measured tRNA levels correlated with codon usage and several tRNAs showed reduced aminoacylation, which was conserved in fasted mice. Finally, we uncovered that the longest codon dwell times could be explained by aminoacylation levels or high codon usage relative to tRNA abundance.

LOTTE-seq (Long hairpin oligonucleotide based tRNA high-throughput sequencing): specific selection of tRNAs with 3'-CCA end for high-throughput sequencing

Transfer RNAs belong to the most abundant type of ribonucleic acid in the cell, and detailed investigations revealed correlations between alterations in the tRNA pool composition and certain diseases like breast cancer. However, currently available methods do not sample the entire tRNA pool or lack specificity for tRNAs. A specific disadvantage of such methods is that only full-length tRNAs are analysed, while tRNA fragments or incomplete cDNAs due to RT stops at modified nucleosides are lost. Another drawback in certain approaches is that the tRNA fraction has to be isolated and separated from high molecular weight RNA, resulting in considerable labour costs and loss of material. Based on a hairpin-shaped adapter oligonucleotide selective for tRNA transcripts, we developed a highly specific protocol for efficient and comprehensive high-throughput analysis of tRNAs that combines the benefits of existing methods and eliminates their disadvantages. Due to a 3'-TGG overhang, the adapter is specifically ligated to the tRNA 3'-CCA end. Reverse transcription prior to the ligation of a second adapter allows to include prematurely terminated cDNA products, increasing the number of tRNA reads. This strategy renders this approach a powerful and universal tool to analyse the tRNA pool of cells and organisms under different conditions in health and disease.

Archaeal NSUN6 catalyzes m5C72 modification on a wide-range of specific tRNAs

Human NOL1/NOP2/Sun RNA methyltransferase family member 6 (hNSun6) generates 5-methylcytosine (m5C) at C72 of four specific tRNAs, and its homologs are present only in higher eukaryotes and hyperthermophilic archaea. Archaeal NSun6 homologs possess conserved catalytic residues, but have distinct differences in their RNA recognition motifs from eukaryotic NSun6s. Until now, the biochemical properties and functions of archaeal NSun6 homologs were unknown. In archaeon Pyrococcus horikoshii OT3, the gene encoding the NSun6 homolog is PH1991. We demonstrated that the PH1991 protein could catalyze m5C72 formation on some specific PhtRNAs in vitro and was thus named as PhNSun6. Remarkably, PhNSun6 has a much wider range of tRNA substrates than hNSun6, which was attributed to its tRNA substrate specificity. The mechanism was further elucidated using biochemical and crystallographic experiments. Structurally, the binding pocket for nucleotide 73 in PhNSun6 is specific to accommodate U73 or G73-containing PhtRNAs. Furthermore, PhNSun6 lacks the eukaryotic NSun6-specific Lys-rich loop, resulting in the non-recognition of D-stem region by PhNSun6. Functionally, the m5C72 modification could slightly promote the thermal stability of PhtRNAs, but did not affect the amino acid accepting activity of PhtRNAs.

Drosophila tsRNAs preferentially suppress general translation machinery via antisense pairing and participate in cellular starvation response

Transfer RNA-derived small RNAs (tsRNAs) are an emerging class of small RNAs, yet their regulatory roles have not been well understood. Here we studied the molecular mechanisms and consequences of tsRNA-mediated regulation in Drosophila. By analyzing 495 public small RNA libraries, we demonstrate that most tsRNAs are conserved, prevalent and abundant in Drosophila. By carrying out mRNA sequencing and ribosome profiling of S2 cells transfected with single-stranded tsRNA mimics and mocks, we show that tsRNAs recognize target mRNAs through conserved complementary sequence matching and suppress target genes by translational inhibition. The target prediction suggests that tsRNAs preferentially suppress translation of the key components of the general translation machinery, which explains how tsRNAs inhibit the global mRNA translation. Serum starvation experiments confirm tsRNAs participate in cellular starvation responses by preferential targeting the ribosomal proteins and translational initiation or elongation factors. Knock-down of AGO2 in S2 cells under normal and starved conditions reveals a dependence of the tsRNA-mediated regulation on AGO2. We also validated the repressive effects of representative tsRNAs on cellular global translation and specific targets with luciferase reporter assays. Our study suggests the tsRNA-mediated regulation might be crucial for the energy homeostasis and the metabolic adaptation in the cellular systems.

Cytosine-5 RNA methylation links protein synthesis to cell metabolism

Posttranscriptional modifications in transfer RNA (tRNA) are often critical for normal development because they adapt protein synthesis rates to a dynamically changing microenvironment. However, the precise cellular mechanisms linking the extrinsic stimulus to the intrinsic RNA modification pathways remain largely unclear. Here, we identified the cytosine-5 RNA methyltransferase NSUN2 as a sensor for external stress stimuli. Exposure to oxidative stress efficiently repressed NSUN2, causing a reduction of methylation at specific tRNA sites. Using metabolic profiling, we showed that loss of tRNA methylation captured cells in a distinct catabolic state. Mechanistically, loss of NSUN2 altered the biogenesis of tRNA-derived noncoding fragments (tRFs) in response to stress, leading to impaired regulation of protein synthesis. The intracellular accumulation of a specific subset of tRFs correlated with the dynamic repression of global protein synthesis. Finally, NSUN2-driven RNA methylation was functionally required to adapt cell cycle progression to the early stress response. In summary, we revealed that changes in tRNA methylation profiles were sufficient to specify cellular metabolic states and efficiently adapt protein synthesis rates to cell stress.

The tRNA Epitranscriptome and Diabetes: Emergence of tRNA Hypomodifications as a Cause of Pancreatic β-Cell Failure

tRNAs are crucial noncoding RNA molecules that serve as amino acid carriers during protein synthesis. The transcription of tRNA genes is a highly regulated process. The tRNA pool is tissue and cell specific, it varies during development, and it is modulated by the environment. tRNAs are highly posttranscriptionally modified by specific tRNA-modifying enzymes. The tRNA modification signature of a cell determines the tRNA epitranscriptome. Perturbations in the tRNA epitranscriptome, as a consequence of mutations in tRNAs and tRNA-modifying enzymes or environmental exposure, have been associated with human disease, including diabetes. tRNA fragmentation induced by impaired tRNA modifications or dietary factors has been linked to pancreatic β-cell demise and paternal inheritance of metabolic traits. Herein, we review recent findings that associate tRNA epitranscriptome perturbations with diabetes.

Translational Control under Stress: Reshaping the Translatome

Adequate reprogramming of cellular metabolism in response to stresses or suboptimal growth conditions involves a myriad of coordinated changes that serve to promote cell survival. As protein synthesis is an energetically expensive process, its regulation under stress is of critical importance. Reprogramming of messenger RNA (mRNA) translation involves well-understood stress-activated kinases that target components of translation initiation machinery, resulting in the robust inhibition of general translation and promotion of the translation of stress-responsive proteins. Translational arrest of mRNAs also results in the accumulation of transcripts in cytoplasmic foci called stress granules. Recent studies focus on the key roles of transfer RNA (tRNA) in stress-induced translational reprogramming. These include stress-specific regulation of tRNA pools, codon-biased translation influenced by tRNA modifications, tRNA miscoding, and tRNA cleavage. In combination, signal transduction pathways and tRNA metabolism changes regulate translation during stress, resulting in adaptation and cell survival. This review examines molecular mechanisms that regulate protein synthesis in response to stress.

The cellular landscape of mid-size noncoding RNA

Noncoding RNA plays an important role in all aspects of the cellular life cycle, from the very basic process of protein synthesis to specialized roles in cell development and differentiation. However, many noncoding RNAs remain uncharacterized and the function of most of them remains unknown. Mid-size noncoding RNAs (mncRNAs), which range in length from 50 to 400 nucleotides, have diverse regulatory functions but share many fundamental characteristics. Most mncRNAs are produced from independent promoters although others are produced from the introns of other genes. Many are found in multiple copies in genomes. mncRNAs are highly structured and carry many posttranscriptional modifications. Both of these facets dictate their RNA-binding protein partners and ultimately their function. mncRNAs have already been implicated in translation, catalysis, as guides for RNA modification, as spliceosome components and regulatory RNA. However, recent studies are adding new mncRNA functions including regulation of gene expression and alternative splicing. In this review, we describe the different classes, characteristics and emerging functions of mncRNAs and their relative expression patterns. Finally, we provide a portrait of the challenges facing their detection and annotation in databases. This article is categorized under: Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs RNA Structure and Dynamics > RNA Structure, Dynamics, and Chemistry RNA Structure and Dynamics > Influence of RNA Structure in Biological Systems RNA Evolution and Genomics > RNA and Ribonucleoprotein Evolution.

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.

tRNA-Derived Small RNAs: Biogenesis, Modification, Function and Potential Impact on Human Disease Development

Transfer RNAs (tRNAs) are abundant small non-coding RNAs that are crucially important for decoding genetic information. Besides fulfilling canonical roles as adaptor molecules during protein synthesis, tRNAs are also the source of a heterogeneous class of small RNAs, tRNA-derived small RNAs (tsRNAs). Occurrence and the relatively high abundance of tsRNAs has been noted in many high-throughput sequencing data sets, leading to largely correlative assumptions about their potential as biologically active entities. tRNAs are also the most modified RNAs in any cell type. Mutations in tRNA biogenesis factors including tRNA modification enzymes correlate with a variety of human disease syndromes. However, whether it is the lack of tRNAs or the activity of functionally relevant tsRNAs that are causative for human disease development remains to be elucidated. Here, we review the current knowledge in regard to tsRNAs biogenesis, including the impact of RNA modifications on tRNA stability and discuss the existing experimental evidence in support for the seemingly large functional spectrum being proposed for tsRNAs. We also argue that improved methodology allowing exact quantification and specific manipulation of tsRNAs will be necessary before developing these small RNAs into diagnostic biomarkers and when aiming to harness them for therapeutic purposes.

Accurate Profiling and Quantification of tRNA Fragments from RNA-Seq Data: A Vade Mecum for MINTmap

There is an increasing interest within the scientific community in identifying tRNA-derived fragments (tRFs) and elucidating the roles they play in the cell. Such endeavors can be greatly facilitated by mining the numerous datasets from many cellular contexts that exist publicly. However, the standard mapping tools cannot be used for the purpose. Several factors complicate this endeavor including: the presence of multiple identical or nearly identical isodecoders at various genomic locations; the presence of identical sequence segments that are shared by isodecoders of the same or even different anticodons; the existence of numerous partial tRNA sequences across the genome; the existence of hundreds of "lookalike" sequences that resemble true tRNAs; and others. This is generating a need for specialized tools that can mine deep sequencing data to identify and quantify tRFs. We discuss the various complicating factors and their ramifications, and how to use and run MINTmap, a tool that addresses these considerations.

LARP4 mRNA codon-tRNA match contributes to LARP4 activity for ribosomal protein mRNA poly(A) tail length protection

Messenger RNA function is controlled by the 3' poly(A) tail (PAT) and poly(A)-binding protein (PABP). La-related protein-4 (LARP4) binds poly(A) and PABP. LARP4 mRNA contains a translation-dependent, coding region determinant (CRD) of instability that limits its expression. Although the CRD comprises <10% of LARP4 codons, the mRNA levels vary >20 fold with synonymous CRD substitutions that accommodate tRNA dynamics. Separately, overexpression of the most limiting tRNA increases LARP4 levels and reveals its functional activity, net lengthening of the PATs of heterologous mRNAs with concomitant stabilization, including ribosomal protein (RP) mRNAs. Genetic deletion of cellular LARP4 decreases PAT length and RPmRNA stability. This LARP4 activity requires its PABP-interaction domain and the RNA-binding module which we show is sensitive to poly(A) 3'-termini, consistent with protection from deadenylation. The results indicate that LARP4 is a posttranscriptional regulator of ribosomal protein production in mammalian cells and suggest that this activity can be controlled by tRNA levels.

Absence of neurological abnormalities in mice homozygous for the Polr3a G672E hypomyelinating leukodystrophy mutation

Recessive mutations in the ubiquitously expressed POLR3A gene cause one of the most frequent forms of childhood-onset hypomyelinating leukodystrophy (HLD): POLR3-HLD. POLR3A encodes the largest subunit of RNA Polymerase III (Pol III), which is responsible for the transcription of transfer RNAs (tRNAs) and a large array of other small non-coding RNAs. In order to study the central nervous system pathophysiology of the disease, we introduced the French Canadian founder Polr3a mutation c.2015G > A (p.G672E) in mice, generating homozygous knock-in (KI/KI) as well as compound heterozygous mice for one Polr3a KI and one null allele (KI/KO). Both KI/KI and KI/KO mice are viable and are able to reproduce. To establish if they manifest a motor phenotype, WT, KI/KI and KI/KO mice were submitted to a battery of behavioral tests over one year. The KI/KI and KI/KO mice have overall normal balance, muscle strength and general locomotion. Cerebral and cerebellar Luxol Fast Blue staining and measurement of levels of myelin proteins showed no significant differences between the three groups, suggesting that myelination is not overtly impaired in Polr3a KI/KI and KI/KO mice. Finally, expression levels of several Pol III transcripts in the brain showed no statistically significant differences. We conclude that the first transgenic mice with a leukodystrophy-causing Polr3a mutation do not recapitulate the childhood-onset HLD observed in the majority of human patients with POLR3A mutations, and provide essential information to guide selection of Polr3a mutations for developing future mouse models of the disease.