The Role of tRNA-Centered Translational Regulatory Mechanisms in Cancer

Cancer is a leading cause of morbidity and mortality worldwide. While numerous factors have been identified as contributing to the development of malignancy, our understanding of the mechanisms involved remains limited. Early cancer detection and the development of effective treatments are therefore critical areas of research. One class of molecules that play a crucial role in the transmission of genetic information are transfer RNAs (tRNAs), which are the most abundant RNA molecules in the human transcriptome. Dysregulated synthesis of tRNAs directly results in translation disorders and diseases, including cancer. Moreover, various types of tRNA modifications and the enzymes responsible for these modifications have been implicated in tumor biology. Furthermore, alterations in tRNA modification can impact tRNA stability, and impaired stability can prompt the cleavage of tRNAs into smaller fragments known as tRNA fragments (tRFs). Initially believed to be random byproducts lacking any physiological function, tRFs have now been redefined as non-coding RNA molecules with distinct roles in regulating RNA stability, translation, target gene expression, and other biological processes. In this review, we present recent findings on translational regulatory models centered around tRNAs in tumors, providing a deeper understanding of tumorigenesis and suggesting new directions for cancer treatment.

The Evolutionary Dynamics of the Mitochondrial tRNA in the Cichlid Fish Family

Simple Summary

Cichlids are a unique example of fish diversity and species richness which have been explained by sympatric speciation at different freshwater sources in Africa. The mitochondria contribute to cell vitality by providing energy. It contains a circular genome with an established translation system that is spatially independent of the cytosolic counterpart. The current study aimed to investigate the evolutionary dynamics of the mitochondrial tRNA and its role in cichlids’ diversity. The available cichlid mitogenomes in the public database were filtered, in addition to newly sequenced accessions from a specific cichlid group known as the haplotilapiine lineage that is widely distributed in the Egyptian sector of the Nile River. Based on the comparative analysis of mitogenomic data, we identified 22 tRNA genes, in which a single gene was D-armless, while the cloverleaf secondary structure subdivided into stem-loop formations was predicted and used to define the levels of genetic divergence for the remained tRNAs. Peculiarly, in cichlids, the formation known as “T-arm” showed the lowest polymorphism levels among other structures in contrast to other organisms (e.g., scorpions). Comparing the whole family to the specific haplotilapiine lineage showed that the tryptophan tRNA was the most conserved tRNA, with signatures of possible purifying selection.

Abstract

The mitochondrial transfer RNA genes (tRNAs) attract more attention due to their highly dynamic and rapidly evolving nature. The current study aimed to detect and evaluate the dynamics, characteristic patterns, and variations of mitochondrial tRNAs. The study was conducted in two main parts: first, the published mitogenomic sequences of cichlids mt tRNAs have been filtered. Second, the filtered mitochondrial tRNA and additional new mitogenomes representing the most prevalent Egyptian tilapiine were compared and analyzed. Our results revealed that all 22 tRNAs of cichlids folded into a classical cloverleaf secondary structure with four domains, except for trnS^GCU, missing the D domain in all cichlids. When consensus tRNAs were compared, most of the mutations were observed in the trnP at nucleotide levels (substitutions and indels), in contrast to trnL^UAA. From a structural perspective, the anticodon loop and T-loop formations were the most conserved structures among all parts of the tRNA in contrast to the A-stem and D-loop formations. The trnW was the lowest polymorphic unneutral tRNA among all cichlids (both the family and the haplotilapiine lineage), in contrast with the neutral trnD that was extremely polymorphic among and within the haplotilapiine lineage species compared to other cichlids species. From a phylogenetic perspective, the trnC was extremely hypervariable and neutral tRNA in both haplotilapiine lineage and cichlids but was unable to report correct phylogenetic signal for the cichlids. In contrast to trnI and trnY, less variable neutral tRNAs that were able to cluster the haplotilapiine lineage and cichlids species as previously reported. By observing the DNA polymorphism in the coding DNA sequences (CDS), the highest affected amino acid by non-synonymous mutations was isoleucine and was equally mutated to valine and vice versa; no correlation between mutations in CDS and tRNAs was statistically found. The current study provides an insight into the mitochondrial tRNA evolution and its effect on the cichlid diversity and speciation model at the maternal level.

Transfer RNA-derived small RNA: an emerging small non-coding RNA with key roles in cancer

Transfer RNAs (tRNAs) promote protein translation by binding to the corresponding amino acids and transporting them to the ribosome, which is essential in protein translation. tRNA-derived small RNAs (tsRNAs) are derived fragments of tRNAs that are cleaved explicitly under certain conditions. An increasing amount of research has demonstrated that tsRNAs have biological functions rather than just being degradation products. tsRNAs can exert functions such as regulating gene expression to influence cancer progression. Their dysregulation is closely associated with various cancers and can serve as diagnostic and prognostic biomarkers for cancer. This review summarizes the generation, classification, and biological functions of tsRNAs, and highlights the roles of tsRNAs in different cancers and their applications as tumor markers.

Protein-coding tRNA sequences?

Transfer RNAs (tRNAs) are ancient molecules likely predating the translation machinery. These extremely conserved RNA molecules transfer amino acids to the ribosome for the synthesis of proteins encoded by mRNAs, but canonical tRNAs are not protein-coding RNAs. Surprisely, when virtually translated, I observed that peptides derived from tRNA sequences match thousands of protein entries in databases. The analysis of these sequences indicates that the vast majority of these tRNA-derived proteins are annotated as small hypothetical peptides, likely arising from sequencing, prediction and/or annotation errors. But life often surpasses fiction. Importantly, tRNA-encoded amino acid domains were also found embedded in large functional proteins. Phylogenetic analysis of representative tRNA-derived protein domains may provide new insights into the origin, plasticity, and evolution of protein-coding genes.

A tRNA processing enzyme is a key regulator of the mitochondrial unfolded protein response

The mitochondrial unfolded protein response (UPR^mt) has emerged as a predominant mechanism that preserves mitochondrial function. Consequently, multiple pathways likely exist to modulate UPR^mt. We discovered that the tRNA processing enzyme, homolog of ELAC2 (HOE-1), is key to UPR^mt regulation in Caenorhabditis elegans. We find that nuclear HOE-1 is necessary and sufficient to robustly activate UPR^mt. We show that HOE-1 acts via transcription factors ATFS-1 and DVE-1 that are crucial for UPR^mt. Mechanistically, we show that HOE-1 likely mediates its effects via tRNAs, as blocking tRNA export prevents HOE-1-induced UPR^mt. Interestingly, we find that HOE-1 does not act via the integrated stress response, which can be activated by uncharged tRNAs, pointing toward its reliance on a new mechanism. Finally, we show that the subcellular localization of HOE-1 is responsive to mitochondrial stress and is subject to negative regulation via ATFS-1. Together, we have discovered a novel RNA-based cellular pathway that modulates UPR^mt.

Streptococcus pneumoniae: a Plethora of Temperate Bacteriophages With a Role in Host Genome Rearrangement

Bacteriophages (phages) are viruses that infect bacteria. They are the most abundant biological entity on Earth (current estimates suggest there to be perhaps 10³¹ particles) and are found nearly everywhere. Temperate phages can integrate into the chromosome of their host, and prophages have been found in abundance in sequenced bacterial genomes. Prophages may modulate the virulence of their host in different ways, e.g., by the secretion of phage-encoded toxins or by mediating bacterial infectivity. Some 70% of Streptococcus pneumoniae (the pneumococcus)—a frequent cause of otitis media, pneumonia, bacteremia and meningitis—isolates harbor one or more prophages. In the present study, over 4000 S. pneumoniae genomes were examined for the presence of prophages, and nearly 90% were found to contain at least one prophage, either defective (47%) or present in full (43%). More than 7000 complete putative integrases, either of the tyrosine (6243) or serine (957) families, and 1210 full-sized endolysins (among them 1180 enzymes corresponding to 318 amino acid-long N-acetylmuramoyl-L-alanine amidases [LytA_PPH]) were found. Based on their integration site, 26 different pneumococcal prophage groups were documented. Prophages coding for tRNAs, putative virulence factors and different methyltransferases were also detected. The members of one group of diverse prophages (PPH090) were found to integrate into the 3’ end of the host lytA_Spn gene encoding the major S. pneumoniae autolysin without disrupting it. The great similarity of the lytA_Spnand lytA_PPH genes (85–92% identity) allowed them to recombine, via an apparent integrase-independent mechanism, to produce different DNA rearrangements within the pneumococcal chromosome. This study provides a complete dataset that can be used to further analyze pneumococcal prophages, their evolutionary relationships, and their role in the pathogenesis of pneumococcal disease.

Emerging Functions for snoRNAs and snoRNA-Derived Fragments

The widespread implementation of mass sequencing has revealed a diverse landscape of small RNAs derived from larger precursors. Whilst many of these are likely to be byproducts of degradation, there are nevertheless metabolically stable fragments derived from tRNAs, rRNAs, snoRNAs, and other non-coding RNA, with a number of examples of the production of such fragments being conserved across species. Coupled with specific interactions to RNA-binding proteins and a growing number of experimentally reported examples suggesting function, a case is emerging whereby the biological significance of small non-coding RNAs extends far beyond miRNAs and piRNAs. Related to this, a similarly complex picture is emerging of non-canonical roles for the non-coding precursors, such as for snoRNAs that are also implicated in such areas as the silencing of gene expression and the regulation of alternative splicing. This is in addition to a body of literature describing snoRNAs as an additional source of miRNA-like regulators. This review seeks to highlight emerging roles for such non-coding RNA, focusing specifically on “new” roles for snoRNAs and the small fragments derived from them.

Human Codon Usage: The Genetic Basis of Pathogen Latency

Infectious diseases pose two main compelling issues. First, the identification of the molecular factors that allow chronic infections, that is, the often completely asymptomatic coexistence of infectious agents with the human host. Second, the definition of the mechanisms that allow the switch from pathogen dormancy to pathologic (re)activation. Furthering previous studies, the present study (1) analyzed the frequency of occurrence of synonymous codons in coding DNA, that is, codon usage, as a genetic tool that rules protein expression; (2) described how human codon usage can inhibit protein expression of infectious agents during latency, so that pathogen genes the codon usage of which does not conform to the human codon usage cannot be translated; and (3) framed human codon usage among the front-line instruments of the innate immunity against infections. In parallel, it was shown that, while genetics can account for the molecular basis of pathogen latency, the changes of the quantitative relationship between codon frequencies and isoaccepting tRNAs during cell proliferation offer a biochemical mechanism that explains the pathogen switching to (re)activation. Immunologically, this study warns that using codon optimization methodologies can (re)activate, potentiate, and immortalize otherwise quiescent, asymptomatic pathogens, thus leading to uncontrollable pandemics.

The complete mitochondrial genome of the Honduran white bat Ectophylla alba (Allen 1982) (Chiroptera: Phyllostomidae)

The Honduran white bat, Ectophylla alba (Allen 1982), is one of eight species belonging to the family Phyllostomidae that exclusively roosts in tents. Due to its restricted distribution, habitat specificity, and diet requirements, E. alba has been strongly affected by habitat loss and fragmentation during the last decade. In this study, we developed the first genomic resource for this species; we assembled and analyzed in detail the complete mitochondrial genome of E. alba. The mitogenome of E. alba is 16,664 bp in length and is comprised of 13 protein coding genes (PCGs), 2 ribosomal RNA genes, 22 transfer RNA genes (tRNAs), and a putative Control Region (CR) 1,232 bp in length. Gene arrangement in the mitochondrial chromosome of E. alba is identical to that reported before in other species of co-familiar bats. All PCGs are under purifying selection, with atp8 experiencing the least selective pressure. In all PCGs, codons ending with adenine are preferred over others ending in thymine and cytosine. Except tRNA-Serine 1, all tRNAs exhibit a cloverleaf secondary structure. The CR of E. alba exhibits three domains commonly described in other mammals, including bats; extended terminal associated sequences (ETAS), central, and conserved sequence block (CSB). A ML phylogenetic reconstruction of the family Phyllostomidae based on all 13 mitochondrial PCGs confirms the monophyletic status of the subfamily Sternodermatinae and indicates the close relationship between E. alba and the genus Artibeus. This is the first genomic resource developed for E. alba and represents the first step to improving our understanding of the genomic underpinnings involved in the evolution of specialization as well as acclimatization and adaptation to local and global change of specialist bats.

Circulating Non-Coding RNAs as a Signature of Autism Spectrum Disorder Symptomatology

Autism spectrum disorder (ASD) is a multifaced neurodevelopmental disorder that becomes apparent during early childhood development. The complexity of ASD makes clinically diagnosing the condition difficult. Consequently, by identifying the biomarkers associated with ASD severity and combining them with clinical diagnosis, one may better factionalize within the spectrum and devise more targeted therapeutic strategies. Currently, there are no reliable biomarkers that can be used for precise ASD diagnosis. Consequently, our pilot experimental cohort was subdivided into three groups: healthy controls, individuals those that express severe symptoms of ASD, and individuals that exhibit mild symptoms of ASD. Using next-generation sequencing, we were able to identify several circulating non-coding RNAs (cir-ncRNAs) in plasma. To the best of our knowledge, this study is the first to show that miRNAs, piRNAs, snoRNAs, Y-RNAs, tRNAs, and lncRNAs are stably expressed in plasma. Our data identify cir-ncRNAs that are specific to ASD. Furthermore, several of the identified cir-ncRNAs were explicitly associated with either the severe or mild groups. Hence, our findings suggest that cir-ncRNAs have the potential to be utilized as objective diagnostic biomarkers and clinical targets.

Plant Small RNA World Growing Bigger: tRNA-Derived Fragments, Longstanding Players in Regulatory Processes

In the past 2 decades, the discovery of a new class of small RNAs, known as tRNA-derived fragments (tRFs), shed light on a new layer of regulation implicated in many biological processes. tRFs originate from mature tRNAs and are classified according to the tRNA regions that they derive from, namely 3′tRF, 5′tRF, and tRF-halves. Additionally, another tRF subgroup deriving from tRNA precursors has been reported, the 3′U tRFs. tRF length ranges from 17 to 26 nt for the 3′and 5′tRFs, and from 30 to 40 nt for tRF-halves. tRF biogenesis is still not yet elucidated, although there is strong evidence that Dicer (and DICER-LIKE) proteins, as well as other RNases such as Angiogenin in mammal and RNS proteins family in plants, are responsible for processing specific tRFs. In plants, the abundance of those molecules varies among tissues, developmental stages, and environmental conditions. More recently, several studies have contributed to elucidate the role that these intriguing molecules may play in all organisms. Among the recent discoveries, tRFs were found to be involved in distinctive regulatory layers, such as transcription and translation regulation, RNA degradation, ribosome biogenesis, stress response, regulatory signaling in plant nodulation, and genome protection against transposable elements. Although tRF biology is still poorly understood, the field has blossomed in the past few years, and this review summarizes the most recent developments in the tRF field in plants.

Adaptation of codon and amino acid use for translational functions in highly expressed cricket genes

BACKGROUND

For multicellular organisms, much remains unknown about the dynamics of synonymous codon and amino acid use in highly expressed genes, including whether their use varies with expression in different tissue types and sexes. Moreover, specific codons and amino acids may have translational functions in highly transcribed genes, that largely depend on their relationships to tRNA gene copies in the genome. However, these relationships and putative functions are poorly understood, particularly in multicellular systems.

RESULTS

Here, we studied codon and amino acid use in highly expressed genes from reproductive and nervous system tissues (male and female gonad, somatic reproductive system, brain and ventral nerve cord, and male accessory glands) in the cricket Gryllus bimaculatus. We report an optimal codon, defined as the codon preferentially used in highly expressed genes, for each of the 18 amino acids with synonymous codons in this organism. The optimal codons were mostly shared among tissue types and both sexes. However, the frequency of optimal codons was highest in gonadal genes. Concordant with translational selection, a majority of the optimal codons had abundant matching tRNA gene copies in the genome, but sometimes obligately required wobble tRNAs. We suggest the latter may comprise a mechanism for slowing translation of abundant transcripts, particularly for cell-cycle genes. Non-optimal codons, defined as those least commonly used in highly transcribed genes, intriguingly often had abundant tRNAs, and had elevated use in a subset of genes with specialized functions (gametic and apoptosis genes), suggesting their use promotes the translational upregulation of particular mRNAs. In terms of amino acids, we found evidence suggesting that amino acid frequency, tRNA gene copy number, and amino acid biosynthetic costs (size/complexity) had all interdependently evolved in this insect model, potentially for translational optimization.

CONCLUSIONS

Collectively, the results suggest a model whereby codon use in highly expressed genes, including optimal, wobble, and non-optimal codons, and their tRNA abundances, as well as amino acid use, have been influenced by adaptation for various functional roles in translation within this cricket. The effects of expression in different tissue types and the two sexes are discussed.

m5U54 tRNA Hypomodification by Lack of TRMT2A Drives the Generation of tRNA-Derived Small RNAs

Transfer RNA (tRNA) molecules contain various post-transcriptional modifications that are crucial for tRNA stability, translation efficiency, and fidelity. Besides their canonical roles in translation, tRNAs also originate tRNA-derived small RNAs (tsRNAs), a class of small non-coding RNAs with regulatory functions ranging from translation regulation to gene expression control and cellular stress response. Recent evidence indicates that tsRNAs are also modified, however, the impact of tRNA epitranscriptome deregulation on tsRNAs generation is only now beginning to be uncovered. The 5-methyluridine (m⁵U) modification at position 54 of cytosolic tRNAs is one of the most common and conserved tRNA modifications among species. The tRNA methyltransferase TRMT2A catalyzes this modification, but its biological role remains mostly unexplored. Here, we show that TRMT2A knockdown in human cells induces m⁵U54 tRNA hypomodification and tsRNA formation. More specifically, m⁵U54 hypomodification is followed by overexpression of the ribonuclease angiogenin (ANG) that cleaves tRNAs near the anticodon, resulting in accumulation of 5′tRNA-derived stress-induced RNAs (5′tiRNAs), namely 5′tiRNA-Gly^GCC and 5′tiRNA-Glu^CTC, among others. Additionally, transcriptomic analysis confirms that down-regulation of TRMT2A and consequently m⁵U54 hypomodification impacts the cellular stress response and RNA stability, which is often correlated with tiRNA generation. Accordingly, exposure to oxidative stress conditions induces TRMT2A down-regulation and tiRNA formation in mammalian cells. These results establish a link between tRNA hypomethylation and ANG-dependent tsRNAs formation and unravel m⁵U54 as a tRNA cleavage protective mark.

Differential expression of microRNAs and tRNA fragments mediate the adaptation of the liver fluke Fasciola gigantica to its intermediate snail and definitive mammalian hosts

The tropical liver fluke Fasciola gigantica affects livestock and humans in many Asian countries, large parts of Africa, and parts of Europe. Despite the public health and economic impacts of F. gigantica, understanding of F. gigantica biology and how the complex lifecycle of this liver fluke is transcriptionally regulated remain unknown. Here, we tested the hypothesis that the regulatory small non-coding RNAs (sncRNAs), microRNAs (miRNAs) and tRNA-derived fragments (tRFs) play roles in the adaptation of F. gigantica to its intermediate and definitive hosts. We sequenced sncRNAs of eight lifecycle stages of F. gigantica. In total, 56 miRNAs from 33 conserved families and four Fasciola-specific miRNAs were identified. Expression analysis of miRNAs suggested clear stage-related patterns. By leveraging the existing transcriptomic data, we predicted a miRNA-based regulation of metabolism, transport, growth and developmental processes. Also, by comparing miRNA complement of F. gigantica with that of Fasciola hepatica, we detected a high level of conservation and identified differences in some miRNAs, which can be used to distinguish the two species. Moreover, we found that tRFs at each lifecycle stage were predominantly derived by tRNA^-Lys and tRNA^-Gly at 5' half sites, but relatively high expression was related to the buffalo-infecting stages. Taken together, we provided a comprehensive overview of the dynamic transcriptional changes of small RNAs that occur during the developmental stages of F. gigantica. This global analysis of F. gigantica lifecycle stages revealed new roles of miRNAs and tRFs in parasite development and will facilitate future research into understanding of fasciolosis pathobiology.

High-Throughput Detection of mtDNA Mutations Leading to tRNA Processing Errors

Some mutations in the tRNA genes of mitochondrial DNA (mtDNA) have been demonstrated to affect the processing of the mitochondrial transcriptome in human patients with mitochondrial disease. A recent analysis of mtDNA mutations in 527 human tumors revealed that approximately a quarter of the somatic mt-tRNA gene mutations lead to aberrant processing of the mitochondrial transcriptome in these tumors. Here, we describe a method, based on mtDNA mutations induced by the mtDNA mutator mouse, to map the sites that lead to transcript processing abnormalities. Mutations in the mtDNA are identified and quantified by amplicon-based mtDNA sequencing, and compared to the allelic ratios observed in matched RNASeq data. Strong deviation in the variant allele frequencies between the amplicon and RNASeq data suggests that such mutations lead to disruptions in mitochondrial transcript processing.

The role of Transfer RNA-Derived Small RNAs (tsRNAs) in Digestive System Tumors

Transfer RNA-derived small RNA(tsRNA) is a type of non-coding tRNA undergoing cleavage by specific nucleases such as Dicer. TsRNAs comprise of tRNA-derived fragments (tRFs) and tRNA halves (tiRNAs). Based on the splicing site within the tRNA, tRFs can be classified into tRF-1, tRF-2, tRF-3, tRF-5, and i-tRF. TiRNAs can be classified into 5′-tiRNA and 3′-tiRNA. Both tRFs and tiRNAs have important roles in carcinogenesis, especially cancer of digestive system. TRFs and tiRNAs can promote cell proliferation and cell cycle progression by regulating the expression of oncogenes, combining with RNA binding proteins such as Y-box binding protein 1 (YBX1) to prevent transcription. Despite many reviews on the basic biological function of tRFs and tiRNAs, few have described their correlation with tumors especially gastrointestinal tumor. This review focused on the relationship of tRFs and tiRNAs with the biological behavior, clinicopathological characteristics, diagnosis, treatment and prognosis of digestive system tumors, and would provide novel insights for the early detection and treatment of digestive system tumors.

Contribution of miRNAs, tRNAs and tRFs to Aberrant Signaling and Translation Deregulation in Lung Cancer

Simple Summary

The profiles of miRNAs, tRNA-derived fragments and tRNAs from lung cancer biopsy specimens indicate involvement of gene networks that modulate signaling and translation initiation. The current study highlights the important role of several regulatory small non-coding RNAs in aberrant signaling and translation deregulation in lung cancer.

Abstract

Transcriptomics profiles of miRNAs, tRNAs or tRFs are used as biomarkers, after separate examination of several cancer cell lines, blood samples or biopsies. However, the possible contribution of all three profiles on oncogenic signaling and translation as a net regulatory effect, is under investigation. The present analysis of miRNAs and tRFs from lung cancer biopsies indicated putative targets, which belong to gene networks involved in cell proliferation, transcription and translation regulation. In addition, we observed differential expression of specific tRNAs along with several tRNA-related genes with possible involvement in carcinogenesis. Transfection of lung adenocarcinoma cells with two identified tRFs and subsequent NGS analysis indicated gene targets that mediate signaling and translation regulation. Broader analysis of all major signaling and translation factors in several biopsy specimens revealed a crosstalk between the PI3K/AKT and MAPK pathways and downstream activation of eIF4E and eEF2. Subsequent polysome profile analysis and 48S pre-initiation reconstitution experiments showed increased global translation rates and indicated that aberrant expression patterns of translation initiation factors could contribute to elevated protein synthesis. Overall, our results outline the modulatory effects that possibly correlate the expression of important regulatory non-coding RNAs with aberrant signaling and translation deregulation in lung cancer.

Diversity of tRNA Clusters in the Chloroviruses

Viruses rely on their host’s translation machinery for the synthesis of their own proteins. Problems belie viral translation when the host has a codon usage bias (CUB) that is different from an infecting virus due to differences in the GC content between the host and virus genomes. Here, we examine the hypothesis that chloroviruses adapted to host CUB by acquisition and selection of tRNAs that at least partially favor their own CUB. The genomes of 41 chloroviruses comprising three clades, each infecting a different algal host, have been sequenced, assembled and annotated. All 41 viruses not only encode tRNAs, but their tRNA genes are located in clusters. While differences were observed between clades and even within clades, seven tRNA genes were common to all three clades of chloroviruses, including the tRNA^Arg gene, which was found in all 41 chloroviruses. By comparing the codon usage of one chlorovirus algal host, in which the genome has been sequenced and annotated (67% GC content), to that of two of its viruses (40% GC content), we found that the viruses were able to at least partially overcome the host’s CUB by encoding tRNAs that recognize AU-rich codons. Evidence presented herein supports the hypothesis that a chlorovirus tRNA cluster was present in the most recent common ancestor (MRCA) prior to divergence into three clades. In addition, the MRCA encoded a putative isoleucine lysidine synthase (TilS) that remains in 39/41 chloroviruses examined herein, suggesting a strong evolutionary pressure to retain the gene. TilS alters the anticodon of tRNA^Met that normally recognizes AUG to then recognize AUA, a codon for isoleucine. This is advantageous to the chloroviruses because the AUA codon is 12–13 times more common in the chloroviruses than their host, further helping the chloroviruses to overcome CUB. Among large DNA viruses infecting eukaryotes, the presence of tRNA genes and tRNA clusters appear to be most common in the Phycodnaviridae and, to a lesser extent, in the Mimiviridae.

Arabidopsis tRNA-derived fragments as potential modulators of translation

Transfer RNA-derived fragments (tRFs) exist in all branches of life. They are involved in RNA degradation, regulation of gene expression, ribosome biogenesis. In archaebacteria, kinetoplastid, yeast, and human cells, they were also shown to regulate translation. In Arabidopsis, the tRFs population fluctuates under developmental or environmental conditions but their functions are yet poorly understood. Here, we show that populations of long (30-35 nt) or short (19-25 nt) tRFs produced from Arabidopsis tRNAs can inhibit in vitro translation of a reporter gene. Analysing a series of oligoribonucleotides mimicking natural tRFs, we demonstrate that only a limited set of tRFs possess the ability to affect protein synthesis. Out of a dozen of tRFs, only two deriving from tRNA^Ala(AGC) and tRNA^Asn(GUU) strongly attenuate translation in vitro. Contrary to human tRF(Ala), the 4 Gs present at the 5' extremity of Arabidopsis tRF(Ala) are not implicated in this inhibition while the G₁₈ and G₁₉ residues are essential. Protein synthesis inhibition by tRFs does not require complementarity with the translated mRNA but, having the capability to be associated with polyribosomes, tRFs likely act as general modulation factors of the translation process in plants.

Human cells adapt to translational errors by modulating protein synthesis rate and protein turnover

Deregulation of tRNAs, aminoacyl-tRNA synthetases (aaRS) or tRNA modifying enzymes, increase the level of protein synthesis errors (PSE) and are associated with several diseases, but the cause-effect mechanisms of these pathologies remain elusive. To clarify the role of PSE in human biology, we have engineered a HEK293 cell line to overexpress a wild type (Wt) tRNASer and two tRNASer mutants that misincorporate serine at non-cognate codon sites. Then, we followed long-term adaptation to PSE of such recombinant cells by analysing cell viability, protein synthesis rate and activation of protein quality control mechanisms (PQC). Engineered cells showed higher level of misfolded and aggregated proteins; activated the ubiquitin-proteasome system (UPS) and the unfolded protein response (UPR), indicative of proteotoxic stress. Adaptation to PSE involved increased protein turnover, UPR up-regulation and altered protein synthesis rate. Gene expression analysis showed that engineered cells presented recurrent alterations in the endoplasmic reticulum, cell adhesion and calcium homeostasis. Herein, we unveil new phenotypic consequences of protein synthesis errors in human cells and identify the protein quality control processes that are necessary for long-term adaptation to PSE and proteotoxic stress. Our data provide important insight on how chronic proteotoxic stress may cause disease and highlight potential biological pathways that support the association of PSE with disease.

tRNA functional signatures classify plastids as late-branching cyanobacteria

BACKGROUND

Eukaryotes acquired the trait of oxygenic photosynthesis through endosymbiosis of the cyanobacterial progenitor of plastid organelles. Despite recent advances in the phylogenomics of Cyanobacteria, the phylogenetic root of plastids remains controversial. Although a single origin of plastids by endosymbiosis is broadly supported, recent phylogenomic studies are contradictory on whether plastids branch early or late within Cyanobacteria. One underlying cause may be poor fit of evolutionary models to complex phylogenomic data.

RESULTS

Using Posterior Predictive Analysis, we show that recently applied evolutionary models poorly fit three phylogenomic datasets curated from cyanobacteria and plastid genomes because of heterogeneities in both substitution processes across sites and of compositions across lineages. To circumvent these sources of bias, we developed CYANO-MLP, a machine learning algorithm that consistently and accurately phylogenetically classifies ("phyloclassifies") cyanobacterial genomes to their clade of origin based on bioinformatically predicted function-informative features in tRNA gene complements. Classification of cyanobacterial genomes with CYANO-MLP is accurate and robust to deletion of clades, unbalanced sampling, and compositional heterogeneity in input tRNA data. CYANO-MLP consistently classifies plastid genomes into a late-branching cyanobacterial sub-clade containing single-cell, starch-producing, nitrogen-fixing ecotypes, consistent with metabolic and gene transfer data.

CONCLUSIONS

Phylogenomic data of cyanobacteria and plastids exhibit both site-process heterogeneities and compositional heterogeneities across lineages. These aspects of the data require careful modeling to avoid bias in phylogenomic estimation. Furthermore, we show that amino acid recoding strategies may be insufficient to mitigate bias from compositional heterogeneities. However, the combination of our novel tRNA-specific strategy with machine learning in CYANO-MLP appears robust to these sources of bias with high accuracy in phyloclassification of cyanobacterial genomes. CYANO-MLP consistently classifies plastids as late-branching Cyanobacteria, consistent with independent evidence from signature-based approaches and some previous phylogenetic studies.

Purification, characterization and cytotoxic activities of individual tRNAs from Escherichia coli

Transfer RNAs (tRNAs) are the most abundant class in small non-coding RNAs which have been proved to be pharmacologically active. In the present study, we evaluated the potential anticancer activities of tRNAs from Escherichia coli MRE 600 to investigate the relationship between non-pathogenic Escherichia coli strain and colorectal cancer. To purify individual tRNAs, we firstly developed a two-dimensional liquid chromatography (2D-LC) and successfully obtained two pure tRNAs. Nuclease mediated base-specific digestions coupled with UHPLC-MS/MS techniques led to an identification of these two tRNAs as tRNA-Val(UAC) and tRNA-Leu(CAG) with typical cloverleaf-like secondary structure. MTT assay demonstrated that both tRNA-1 and tRNA-2 exhibit strong cytotoxicity with IC₅₀ of 113.0 nM and 124.8 nM on HCT-8 cells in a dose-dependent manner. Further clonogenic assay revealed that the purified tRNAs exhibit significant inhibition in colony formation with survival percentage of 79.0 ± 1.6 and 71.2 ± 2.2 at the concentration of 100 nM. These findings provided evidences of anticancer activities of tRNAs from non-pathogenic Escherichia coli strain, indicating that the pharmacological effects of these neglected biomacromolecules from microorganisms should be emphasized. This study put new insights into the therapeutic effects of intestinal microorganism on human diseases, therefore broadened our knowledge of the biological functions of gut microbiota.

tRNA Deregulation and Its Consequences in Cancer

The expression of transfer RNAs (tRNAs) is deregulated in cancer cells but the mechanisms and functional meaning of such deregulation are poorly understood. The proteome of cancer cells is not fully encoded by their transcriptome, however, the contribution of mRNA translation to such diversity remains to be elucidated. We review data supporting the hypothesis that tRNA expression deregulation and translational error rate is an important contributor to proteome diversity and cell population heterogeneity, genome instability, and drug resistance in tumors. This hypothesis is aligned with recent data in various model organisms, showing unanticipated adaptive roles of translational errors (adaptive mistranslation), expression control of specific gene subsets by tRNAs, and proteome diversification by elevation of translational error rates in tumors.

Gordon H. Dixon's trace in my personal career and the quantic jump experienced in regulatory information

Even before Rosalin Franklin had discovered the DNA double helix, in her impressive X-ray diffraction image pattern, Erwin Schröedinger, described, in his excellent book, What is Life, how the finding of aperiodic crystals in biological systems surprised him (an aperiodic crystal, which, in my opinion is the material carrier of life). In the 21st century and still far from being able to define life, we are attending to a quick acceleration of knowledge on regulatory information. With the discovery of new codes and punctuation marks, we will greatly increase our understanding in front of an impressive avalanche of genomic sequences. Trifonov et al. defined a genetic code as a widespread DNA sequence pattern that carries a message with an impact on biology. These patterns are largely captured in transcribed messages that give meaning and identity to the particular cells. In this review, I will go through my personal career in and after my years of work in the laboratory of Gordon H. Dixon, extending toward the impressive acquisition of new knowledge on regulatory information and genetic codes provided by remarkable scientists in the field. Abbreviations: CA II: carbonic anhydridase II (chicken); Car2: carbonic anhydridase 2 (mouse); CpG islands: short (>0.5 kb) stretches of DNA with a G+C content ≥55%; DNMT1: DNA methyltransferases 1; DNMT3b: DNA methyltransferases 3B; DSB: double-strand DNA breaks; ERT: endogenous retrotransposon; ERV: endogenous retroviruses; ES cells: embryonic stem cells; GAPDH: glyceraldehide phosphate dehydrogenase; H1: histone H1; HATs: histone acetyltransferases; HDACs: histone deacetylases; H3K4me3: histone 3 trimethylated at lys 4; H3K79me2: histone 3 dimethylated at lys 79; HMG: high mobility group proteins; HMT: histone methyltransferase; HP1: heterochromatin protein 1; HR: homologous recombination; HSE: heat-shock element; ICRs: imprinted control regions; IRF: interferon regulatory factor; LDH-A/-B: lactate dehydrogenase A/B; LTR: long terminal repeats; MeCP2: methyl CpG binding protein 2; OCT4: octamer-binding transcription factor 4; PAF1: RNA Polymerase II associated factor 1; piRNA: PIWI-interacting RNA; poly(A) tails: poly-adenine tails; PRC2: polycomb repressive complex 2; PTMs: post-translational modifications; SIRT 1: sirtuin 1, silent information regulator; STAT3: signal transducer and activator of transcription; tRNAs: transfer RNA; tRFs: tRNA-derived fragments; TSS: transcription start site; TE: transposable elements; UB I: polyubiquitin I; UB II: polyubiquitin II; UBE 2N: ubiquitin conjugating enzyme E2N; 5'-UTR: 5'-untranslated sequences; 3'-UTR: 3'-untranslated sequences.

Codon misreading tRNAs promote tumor growth in mice

Deregulation of tRNAs, aminoacyl-tRNA synthetases and tRNA modifying enzymes are common in cancer, raising the hypothesis that protein synthesis efficiency and accuracy (mistranslation) are compromised in tumors. We show here that human colon tumors and xenograft tumors produced in mice by two epithelial cancer cell lines mistranslate 2- to 4-fold more frequently than normal tissue. To clarify if protein mistranslation plays a role in tumor biology, we expressed mutant Ser-tRNAs that misincorporate Ser-at-Ala (frequent error) and Ser-at-Leu (infrequent error) in NIH3T3 cells and investigated how they responded to the proteome instability generated by the amino acid misincorporations. There was high tolerance to both misreading tRNAs, but the Ser-to-Ala misreading tRNA was a more potent inducer of cell transformation, stimulated angiogenesis and produced faster growing tumors in mice than the Ser-to-Leu misincorporating tRNA. Upregulation of the Akt pathway and the UPR were also observed. Most surprisingly, the relative expression of both misreading tRNAs increased during tumor growth, suggesting that protein mistranslation is advantageous in cancer contexts. These data highlight new features of protein synthesis deregulation in tumor biology.

tRNA-derived small RNAs: New players in genome protection against retrotransposons

Previously believed to be mere random degradation products, tRNA-derived small RNAs have been lately connected to a series of functions that include, surprisingly, genome protection against retrotransposons. tRNAs have been known for a long time to be involved in the replication cycle of retroviruses, pararetroviruses and retrotransposons as primers of their reverse transcription. tRNA-derived small RNAs, as functional small RNAs or as mere tRNA degradation products, have emerged as important players in the regulation of genic transcription. Nevertheless, the involvement of functional sRNAs derived from tRNA transcripts in transposon posttranscriptional control is a regulatory layer that remained elusive until now. Here I review the recent discoveries in the field that connect tRNA-derived small RNAs and retrotransposon control.

tRNAs as primers and inhibitors of retrotransposons

The functional relationship between tRNAs and retrotransposons have been known for more than 35 years. tRNAs are used as primer molecules to guide the reverse transcription of retrotransposons. Recently, tRNAs have also emerge as important players in the postranscriptional regulation of retrotransposons by means of tRNA-derived small RNAs. This surprisingly new layer of regulation indicates that tRNAs are used both in the promotion and the suppression of the reverse transcription of retrotransposons indicating their primary role in the life cycle of LTR retrotransposons. This adds another level of translational control to tRNAs. Here we review the different known levels of interactions of tRNAs and retrotransposons and highlight the unknown parts of this interaction.

SMORE: Synteny Modulator of Repetitive Elements

Several families of multicopy genes, such as transfer ribonucleic acids (tRNAs) and ribosomal RNAs (rRNAs), are subject to concerted evolution, an effect that keeps sequences of paralogous genes effectively identical. Under these circumstances, it is impossible to distinguish orthologs from paralogs on the basis of sequence similarity alone. Synteny, the preservation of relative genomic locations, however, also remains informative for the disambiguation of evolutionary relationships in this situation. In this contribution, we describe an automatic pipeline for the evolutionary analysis of such cases that use genome-wide alignments as a starting point to assign orthology relationships determined by synteny. The evolution of tRNAs in primates as well as the history of the Y RNA family in vertebrates and nematodes are used to showcase the method. The pipeline is freely available.

Origin and Evolution of RNA-Dependent RNA Polymerase

RNA-dependent RNA polymerases (RdRp) are very ancient enzymes and are essential for all viruses with RNA genomes. We reconstruct the origin and evolution of this polymerase since the initial stages of the origin of life. The origin of the RdRp was traced back from tRNA ancestors. At the origin of the RdRp the most ancient part of the protein is the cofactor-binding site that had the capacity of binding to simple molecules as magnesium, calcium, and ribonucleotides. Our results suggest that RdRp originated from junctions of proto-tRNAs that worked as the first genes at the emergence of the primitive translation system, where the RNA was the informational molecule. The initial domain, worked as a building block for the emergence of the fingers and thumb domains. From the ancestral RdRp, we could establish the evolutionary stages of viral evolution from a rooted ancestor to modern viruses. It was observed that the selective pressure under the RdRp was the organization and functioning of the genome, where RNA double-stranded and RNA single-stranded virus formed a separate group. We propose an evolutionary route to the polymerases and the results suggest an ancient scenario for the origin of RNA viruses.

Involvement of tRNAs in replication of human mitochondrial DNA and modifying effects of telomerase

Overexpression of telomerase has been shown to significantly increase the lifespan of mice. When mechanistically attributed to repair of critically short telomeres, the lifespan extending action of telomerase cannot be reconciled with the observation that telomerase-null mice do not exhibit shortening of lifespan for at least two generations. We hypothesized that telomerase may interfere with replication of mitochondrial DNA (mtDNA) in a way that reduces formation of deletions - the primary cause of age-dependent cell attrition in non-renewable cells such as myocytes and neurons. Here we show that several tRNA genes may function as alternative origins of replication (ORIs). We also show that telomerase interacts with canonical light strand ORI (ORI_L) and tRNAs and modifies their activities. Our results suggest that replication of mitochondrial DNA (mtDNA) proceeds through a variety of mechanisms resulting in a mixture of classic strand-displacement mode, and coupled replication of heavy and light strands. Our results also suggest that effects of telomerase may arise from binding ORI_L and thus limiting contribution of the deletion-prone strand displacement mode to mtDNA synthesis. These findings imply that it may be possible to uncouple detrimental and beneficial effects of telomerase, and thereby to improve telomerase-based strategies to extend lifespan.

tRNA Derived smallRNAs: smallRNAs Repertoire Has Yet to Be Decoded in Plants

Among several smallRNAs classes, microRNAs play an important role in controlling the post-transcriptional events. Next generation sequencing has played a major role in extending the landscape of miRNAs and revealing their spatio-temporal roles in development and abiotic stress. Lateral evolution of these smallRNAs classes have widely been seen with the recently emerging knowledge on tRNA derived smallRNAs. In the present perspective, we discussed classification, identification and roles of tRNA derived smallRNAs across plants and their potential involvement in abiotic and biotic stresses.

Identification of urinary exosomal noncoding RNAs as novel biomarkers in chronic kidney disease

In chronic kidney disease (CKD), the decline in the glomerular filtration rate is associated with increased morbidity and mortality and thus poses a major challenge for healthcare systems. While the contribution of tissue-derived miRNAs and mRNAs to CKD progression has been extensively studied, little is known about the role of urinary exosomes and their association with CKD. Exosomes are small, membrane-derived endocytic vesicles that contribute to cell-to-cell communication and are present in various body fluids, such as blood or urine. Next-generation sequencing approaches have revealed that exosomes are enriched in noncoding RNAs and thus exhibit great potential for sensitive nucleic acid biomarkers in various human diseases. Therefore, in this study we aimed to identify urinary exosomal ncRNAs as novel biomarkers for diagnosis of CKD. Since up to now most approaches have focused on the class of miRNAs, we extended our analysis to several other noncoding RNA classes, such as tRNAs, tRNA fragments (tRFs), mitochondrial tRNAs, or lincRNAs. For their computational identification from RNA-seq data, we developed a novel computational pipeline, designated as ncRNASeqScan. By these analyses, in CKD patients we identified 30 differentially expressed ncRNAs, derived from urinary exosomes, as suitable biomarkers for early diagnosis. Thereby, miRNA-181a appeared as the most robust and stable potential biomarker, being significantly decreased by about 200-fold in exosomes of CKD patients compared to healthy controls. Using a cell culture system for CKD indicated that urinary exosomes might indeed originate from renal proximal tubular epithelial cells.

tRFs: miRNAs in disguise

tRFs and tiRNAs are two new classes of regulatory non-coding small RNAs that are derived from the cleavage of pre-existing tRNAs. tRFs are 18-22 nt long and are classified into the tRF-5, tRF-3, and tRF-1 series. Here, we discuss in detail the regulatory roles of tRFs in translation, viral infections, and carcinogenesis. Moreover, we have reviewed the association of tRFs with Argonaute proteins, including their potential to function as miRNAs. Interestingly, few miRNAs are generated from pre-existing tRNAs. Hence, tRNAs generate similar-sized tRFs and miRNAs, leading to misannotations due to cross mapping of tRFs and tRNA-derived miRNAs during deep sequencing data analysis. Therefore, it is important to catalogue the overlapping sequences between tRNA-derived miRNAs and tRFs. We have catalogued the miRNAs that overlap with tRFs sequences in humans using miRBase. We identified 20 tRNA-derived miRNAs that share sequences with tRFs. Of the 20 miRNAs, 5 miRNAs (miR-3182, miR-4521, miR-1260a, miR-1260b, and miR-7977) showed significant prediction scores. Furthermore, we have identified a lysine degradation pathway as a common regulatory pathway for miR-1260a, miR-1260b, and miR-3182 by using DIANA-mirPath.

Exploring GpG bases next to anticodon in tRNA subsets

Transfer RNA (tRNA) structure, modifications and functions are evolutionary and established in bacteria, archaea and eukaryotes. Typically the tRNA modifications are indispensable for its stability and are required for decoding the mRNA into amino acids for protein synthesis. A conserved methylation has been located on the anticodon loop specifically at the 37^th position and it is next to the anticodon bases. This modification is called as m1G37 and it is catalyzed by tRNA (m¹G37) methyltransferase (TrmD). It is deciphered that G37 positions occur on few additional amino acids specific tRNA subsets in bacteria. Furthermore, Archaea and Eukaryotes have more number of tRNA subsets which contains G37 position next to the anticodon and the G residue are located at different positions such as G36, G37, G38, 39, and G40. In eight bacterial species, G (guanosine) residues are presents at the 37^th and 38^th position except three tRNA subsets having G residues at 36^th and 39^th positions. Therefore we propose that m1G37 modification may be feasible at 36^th, 37^th, 38^th, 39^th and 40^th positions next to the anticodon of tRNAs. Collectively, methylation at G residues close to the anticodon may be possible at different positions and without restriction of anticodon 3^rd base A, C, U or G.

Fragile genomic sites are associated with origins of replication

Genome rearrangements are mediators of evolution and disease. Such rearrangements are frequently bounded by transfer RNAs (tRNAs), transposable elements, and other repeated elements, suggesting a functional role for these elements in creating or repairing breakpoints. Though not well explored, there is evidence that origins of replication also colocalize with breakpoints. To investigate a potential correlation between breakpoints and origins, we analyzed evolutionary breakpoints defined between Saccharomyces cerevisiae and Kluyveromyces waltii and S. cerevisiae and a hypothetical ancestor of both yeasts, as well as breakpoints reported in the experimental literature. We find that origins correlate strongly with both evolutionary breakpoints and those described in the literature. Specifically, we find that origins firing earlier in S phase are more strongly correlated with breakpoints than are later-firing origins. Despite origins being located in genomic regions also bearing tRNAs and Ty elements, the correlation we observe between origins and breakpoints appears to be independent of these genomic features. This study lays the groundwork for understanding the mechanisms by which origins of replication may impact genome architecture and disease.