tRNA biology charges to the front

MncR: Late Integration Machine Learning Model for Classification of ncRNA Classes Using Sequence and Structural Encoding

Non-coding RNA (ncRNA) classes take over important housekeeping and regulatory functions and are quite heterogeneous in terms of length, sequence conservation and secondary structure. High-throughput sequencing reveals that the expressed novel ncRNAs and their classification are important to understand cell regulation and identify potential diagnostic and therapeutic biomarkers. To improve the classification of ncRNAs, we investigated different approaches of utilizing primary sequences and secondary structures as well as the late integration of both using machine learning models, including different neural network architectures. As input, we used the newest version of RNAcentral, focusing on six ncRNA classes, including lncRNA, rRNA, tRNA, miRNA, snRNA and snoRNA. The late integration of graph-encoded structural features and primary sequences in our MncR classifier achieved an overall accuracy of >97%, which could not be increased by more fine-grained subclassification. In comparison to the actual best-performing tool ncRDense, we had a minimal increase of 0.5% in all four overlapping ncRNA classes on a similar test set of sequences. In summary, MncR is not only more accurate than current ncRNA prediction tools but also allows the prediction of long ncRNA classes (lncRNAs, certain rRNAs) up to 12.000 nts and is trained on a more diverse ncRNA dataset retrieved from RNAcentral.

XRN2 suppresses aberrant entry of tRNA trailers into argonaute in humans and Arabidopsis

MicroRNAs (miRNAs) are a well-characterized class of small RNAs (sRNAs) that regulate gene expression post-transcriptionally. miRNAs function within a complex milieu of other sRNAs of similar size and abundance, with the best characterized being tRNA fragments or tRFs. The mechanism by which the RNA-induced silencing complex (RISC) selects for specific sRNAs over others is not entirely understood in human cells. Several highly expressed tRNA trailers (tRF-1s) are strikingly similar to microRNAs in length but are generally excluded from the microRNA effector pathway. This exclusion provides a paradigm for identifying mechanisms of RISC selectivity. Here, we show that 5' to 3' exoribonuclease XRN2 contributes to human RISC selectivity. Although highly abundant, tRF-1s are highly unstable and degraded by XRN2 which blocks tRF-1 accumulation in RISC. We also find that XRN mediated degradation of tRF-1s and subsequent exclusion from RISC is conserved in plants. Our findings reveal a conserved mechanism that prevents aberrant entry of a class of highly produced sRNAs into Ago2.

Emerging roles of tRNA in cancer

Transfer RNAs (tRNAs) play pivotal roles in the transmission of genetic information, and abnormality of tRNAs directly leads to translation disorders and causes diseases, including cancer. The complex modifications enable tRNA to execute its delicate biological function. Alteration of appropriate modifications may affect the stability of tRNA, impair its ability to carry amino acids, and disrupt the pairing between anticodons and codons. Studies confirmed that dysregulation of tRNA modifications plays an important role in carcinogenesis. Furthermore, when the stability of tRNA is impaired, tRNAs are cleaved into small tRNA fragments (tRFs) by specific RNases. Though tRFs have been found to play vital regulatory roles in tumorigenesis, its formation process is far from clear. Understanding improper tRNA modifications and abnormal formation of tRFs in cancer is conducive to uncovering the role of metabolic process of tRNA under pathological conditions, which may open up new avenues for cancer prevention and treatment.

PANDORA-Seq unveils the hidden small noncoding RNA landscape in atherosclerosis of LDL receptor-deficient mice

Small noncoding RNAs (sncRNAs) play diverse roles in numerous biological processes. While the widely used RNA sequencing (RNA-Seq) method has advanced sncRNA discovery, RNA modifications can interfere with the complementary DNA library construction process, preventing the discovery of highly modified sncRNAs including transfer RNA-derived small RNAs (tsRNAs) and ribosomal RNA-derived small RNAs (rsRNAs) that may have important functions in disease development. To address this technical obstacle, we recently developed a novel PANDORA-Seq (Panoramic RNA Display by Overcoming RNA Modification Aborted Sequencing) method to overcome RNA modification-elicited sequence interferences. To identify novel sncRNAs associated with atherosclerosis development, LDL receptor-deficient (LDLR-/-) mice were fed a low-cholesterol diet or high-cholesterol diet (HCD) for 9 weeks. Total RNAs isolated from the intima were subjected to PANDORA-Seq and traditional RNA-Seq. By overcoming RNA modification-elicited limitations, PANDORA-Seq unveiled an rsRNA/tsRNA-enriched sncRNA landscape in the atherosclerotic intima of LDLR-/- mice, which was strikingly different from that detected by traditional RNA-Seq. While microRNAs were the dominant sncRNAs detected by traditional RNA-Seq, PANDORA-Seq substantially increased the reads of rsRNAs and tsRNAs. PANDORA-Seq also detected 1,383 differentially expressed sncRNAs induced by HCD feeding, including 1,160 rsRNAs and 195 tsRNAs. One of HCD-induced intimal tsRNAs, tsRNA-Arg-CCG, may contribute to atherosclerosis development by regulating the proatherogenic gene expression in endothelial cells. Overall, PANDORA-Seq revealed a hidden rsRNA and tsRNA population associated with atherosclerosis development. These understudied tsRNAs and rsRNAs, which are much more abundant than microRNAs in the atherosclerotic intima of LDLR-/- mice, warrant further investigations.

RNA methylation and cellular response to oxidative stress-promoting anticancer agents

Disruption of the complex network that regulates redox homeostasis often underlies resistant phenotypes, which hinder effective and long-lasting cancer eradication. In addition, the RNA methylome-dependent control of gene expression also critically affects traits of cellular resistance to anti-cancer agents. However, few investigations aimed at establishing whether the epitranscriptome-directed adaptations underlying acquired and/or innate resistance traits in cancer could be implemented through the involvement of redox-dependent or -responsive signaling pathways. This is unexpected mainly because: i) the effectiveness of many anti-cancer approaches relies on their capacity to promote oxidative stress (OS); ii) altered redox milieu and reprogramming of mitochondrial function have been acknowledged as critical mediators of the RNA methylome-mediated response to OS. Here we summarize the current state of understanding on this topic, as well as we offer new perspectives that might lead to original approaches and strategies to delay or prevent the problem of refractory cancer and tumor recurrence.

The importance of pseudouridylation: human disorders related to the fifth nucleoside

Pseudouridylation is one of the most abundant RNA modifications in eukaryotes, making pseudouridine known as the "fifth nucleoside." This highly conserved alteration affects all non-coding and coding RNA types. Its role and importance have been increasingly widely researched, especially considering that its absence or damage leads to serious hereditary diseases. Here, we summarize the human genetic disorders described to date that are related to the participants of the pseudouridylation process.

Excision of Integrated Human Herpesvirus 6A Genomes Using CRISPR/Cas9 Technology

Human herpesviruses 6A and 6B are betaherpesviruses that can integrate their genomes into the telomeres of latently infected cells. Integration can also occur in germ cells, resulting in individuals who harbor the integrated virus in every cell of their body and can pass it on to their offspring. This condition is termed inherited chromosomally integrated HHV-6 (iciHHV-6) and affects about 1% of the human population. The integrated HHV-6A/B genome can reactivate in iciHHV-6 patients and in rare cases can also cause severe diseases including encephalitis and graft-versus-host disease. Until now, it has remained impossible to prevent virus reactivation or remove the integrated virus genome. Therefore, we developed a system that allows the removal of HHV-6A from the host telomeres using the CRISPR/Cas9 system. We used specific guide RNAs (gRNAs) targeting the direct repeat region at the ends of the viral genome to remove the virus from latently infected cells generated in vitro and iciHHV-6A patient cells. Fluorescence-activated cell sorting (FACS), quantitative PCR (qPCR), and fluorescence in situ hybridization (FISH) analyses revealed that the virus genome was efficiently excised and lost in most cells. Efficient excision was achieved with both constitutive and transient expression of Cas9. In addition, reverse transcription-qPCR (RT-qPCR) revealed that the virus genome did not reactivate upon excision. Taken together, our data show that our CRISPR/Cas9 approach allows efficient removal of the integrated virus genome from host telomeres. IMPORTANCE Human herpesvirus 6 (HHV-6) infects almost all humans and integrates into the telomeres of latently infected cells to persist in the host for life. In addition, HHV-6 can also integrate into the telomeres of germ cells, which results in about 80 million individuals worldwide who carry the virus in every cell of their body and can pass it on to their offspring. In this study, we develop the first system that allows excision of the integrated HHV-6 genome from host telomeres using CRISPR/Cas9 technology. Our data revealed that the integrated HHV-6 genome can be efficiently removed from the telomeres of latently infected cells and cells of patients harboring the virus in their germ line. Virus removal could be achieved with both stable and transient Cas9 expression, without inducing viral reactivation.

Elucidation of transfer RNAs as stress regulating agents and the experimental strategies to conceive the functional role of tRNA-derived fragments in plants

In plants, the transfer RNAs (tRNAs) exhibit their profound influence in orchestrating diverse physiological activities like cell growth, development, and response to several surrounding stimuli. The tRNAs, which were known to restrict their function solely in deciphering the codons, are now emerging as frontline defenders in stress biology. The plants that are constantly confronted with a huge panoply of stresses rely on tRNA-mediated stress regulation by altering the tRNA abundance, curbing the transport of tRNAs, fragmenting the mature tRNAs during stress. Among them, the studies on the generation of transfer RNA-derived fragments (tRFs) and their biological implication in stress response have attained huge interest. In plants, the tRFs hold stable expression patterns and regulate biological functions under diverse environmental conditions. In this review, we discuss the fate of plant tRNAs upon stress and thereafter how the tRFs are metamorphosed into sharp ammunition to wrestle with stress. We also address the various methods developed to date for uncovering the role of tRFs and their function in plants.

Systemic inflammation as a biomarker of seizure propensity and a target for treatment to reduce seizure propensity

People with diabetes can wear a device that measures blood glucose and delivers just the amount of insulin needed to return the glucose level to within bounds. Currently, people with epilepsy do not have access to an equivalent wearable device that measures a systemic indicator of an impending seizure and delivers a rapidly acting medication or other intervention (e.g., an electrical stimulus) to terminate or prevent a seizure. Given that seizure susceptibility is reliably increased in systemic inflammatory states, we propose a novel closed-loop device where release of a fast-acting therapy is governed by sensors that quantify the magnitude of systemic inflammation. Here, we review the evidence that patients with epilepsy have raised levels of systemic indicators of inflammation than controls, and that some anti-inflammatory drugs have reduced seizure occurrence in animals and humans. We then consider the options of what might be incorporated into a responsive anti-seizure system.

Transfer RNA-derived small RNAs in tumor microenvironment

Transfer RNAs (tRNAs) are a class of non-coding RNAs responsible for amino acid translocation during protein synthesis and are ubiquitously found in organisms. With certain modifications and under specific conditions, tRNAs can be sheared and fragmented into small non-coding RNAs, also known as tRNA-derived small RNAs (tDRs). With the development of high-throughput sequencing technologies and bioinformatic strategies, more and more tDRs have been identified and their functions in organisms have been characterized. tRNA and it derived tDRs, have been shown to be essential not only for transcription and translation, but also for regulating cell proliferation, apoptosis, metastasis, and immunity. Aberrant expression of tDRs is associated with a wide range of human diseases, especially with tumorigenesis and tumor progression. The tumor microenvironment (TME) is a complex ecosystem consisting of various cellular and cell-free components that are mutually compatible with the tumor. It has been shown that tDRs regulate the TME by regulating cancer stem cells, immunity, energy metabolism, epithelial mesenchymal transition, and extracellular matrix remodeling, playing a pro-tumor or tumor suppressor role. In this review, the biogenesis, classification, and function of tDRs, as well as their effects on the TME and the clinical application prospects will be summarized and discussed based on up to date available knowledge.

Emerging roles of tRNA-derived fragments in cancer

tRNA-derived fragments (tRFs) are an emerging category of small non-coding RNAs that are generated from cleavage of mature tRNAs or tRNA precursors. The advance in high-throughput sequencing has contributed to the identification of increasing number of tRFs with critical functions in distinct physiological and pathophysiological processes. tRFs can regulate cell viability, differentiation, and homeostasis through multiple mechanisms and are thus considered as critical regulators of human diseases including cancer. In addition, increasing evidence suggest the extracellular tRFs may be utilized as promising diagnostic and prognostic biomarkers for cancer liquid biopsy. In this review, we focus on the biogenesis, classification and modification of tRFs, and summarize the multifaceted functions of tRFs with an emphasis on the current research status and perspectives of tRFs in cancer.

Paternal phthalate exposure-elicited offspring metabolic disorders are associated with altered sperm small RNAs in mice

Exposure to ubiquitous plastic-associated endocrine disrupting chemicals (EDCs) is associated with the increased risk of many chronic diseases. For example, phthalate exposure is associated with cardiometabolic mortality in humans, with societal costs ∼ $39 billion/year or more. We recently demonstrated that several widely used plastic-associated EDCs increase cardiometabolic disease in appropriate mouse models. In addition to affecting adult health, parental exposure to EDCs has also been shown to cause metabolic disorders, including obesity and diabetes, in the offspring. While most studies have focused on the impact of maternal EDC exposure on the offspring's health, little is known about the effects of paternal EDC exposure. In the current study, we investigated the adverse impact of paternal exposure to a ubiquitous but understudied phthalate, dicyclohexyl phthalate (DCHP) on the metabolic health of F1 and F2 offspring in mice. Paternal DCHP exposure led to exacerbated insulin resistance and impaired insulin signaling in F1 offspring without affecting diet-induced obesity. We previously showed that sperm small non-coding RNAs including tRNA-derived small RNAs (tsRNAs) and rRNA-derived small RNAs (rsRNAs) contribute to the intergenerational transmission of paternally acquired metabolic disorders. Using a novel PANDORA-seq, we revealed that DCHP exposure can lead to sperm tsRNA/rsRNA landscape changes that were undetected by traditional RNA-seq, which may contribute to DCHP-elicited adverse effects. Lastly, we found that paternal DCHP can also cause sex-specific transgenerational adverse effects in F2 offspring and elicited glucose intolerance in female F2 descendants. Our results suggest that exposure to endocrine disrupting phthalates may have intergenerational and transgenerational adverse effects on the metabolic health of their offspring. These findings increase our understanding of the etiology of chronic human diseases originating from chemical-elicited intergenerational and transgenerational effects.

tRNA-derived small RNAs in plant response to biotic and abiotic stresses

tRNA-derived small RNAs (tsRNAs) represent a novel category of small non-coding RNAs and serve as a new regulator of gene expression at both transcriptional and post-transcriptional levels. Growing evidence indicates that tsRNAs can be induced by diverse stimuli and regulate stress-responsive target genes, allowing plants to adapt to unfavorable environments. Here, we discuss the latest developments about the biogenesis and classification of tsRNAs and highlight the expression regulation and potential function of tsRNAs in plant biotic and abiotic stress responses. Of note, we also collect useful bioinformatics tools and resources for tsRNAs study in plants. Finally, we propose current limitations and future directions for plant tsRNAs research. These recent discoveries have refined our understanding of whether and how tsRNAs enhance plant stress tolerance.

tRNA-Derived RNA Fragments Are Novel Biomarkers for Diagnosis, Prognosis, and Tumor Subtypes in Prostate Cancer

BACKGROUND

tRNA-derived RNA fragments (tRFs) are a novel class of small ncRNA that are derived from precursor or mature tRNAs. Recently, the general relevance of their roles and clinical values in tumorigenesis, metastasis, and recurrence have been increasingly highlighted. However, there has been no specific systematic study to elucidate any potential clinical significance for these tRFs in prostate adenocarcinoma (PRAD), one of the most common and malignant cancers that threatens male health worldwide. Here, we investigate the clinical value of 5'-tRFs in PRAD.

METHODS

Small RNA sequencing data were analyzed to discover new 5'-tRFs biomarkers for PRAD. Machine learning algorithms were used to identify 5'-tRF classifiers to distinguish PRAD tumors from normal tissues. LASSO and Cox regression analyses were used to construct 5'-tRF prognostic predictive models. NMF and consensus clustering analyses were performed on 5'-tRF profiles to identify molecular subtypes of PRAD.

RESULTS

The overall levels of 5'-tRFs were significantly upregulated in the PRAD tumor samples compared to their adjacent normal samples. tRF classifiers composed of 13 5'-tRFs achieved AUC values as high as 0.963, showing high sensitivity and specificity in distinguishing PRAD tumors from normal samples. Multiple 5'-tRFs were identified as being associated with the PRAD prognosis. The tRF score, defined by a set of eight 5'-tRFs, was highly predictive of survival in PRAD patients. The combination of tRF and Gleason scores showed a significantly better performance than the Gleason score alone, suggesting that 5'-tRFs can offer PRAD patients additional and improved prognostic information. Four molecular subtypes of the PRAD tumor were identified based on their 5'-tRF expression profiles. Genetically, these 5'-tRFs PRAD tumor subtypes exhibited distinct genomic landscapes in tumor cells. Clinically, they showed marked differences in survival and clinicopathological features.

CONCLUSIONS

5'-tRFs are potential clinical biomarkers for the diagnosis, prognosis, and classification of tumor subtypes on a molecular level. These can help clinicians formulate personalized treatment plans for PRAD patients and may have similar potential applications for other disease types.

Synthetic biology tools to promote the folding and function of RNA aptamers in mammalian cells

RNA aptamers are structured RNAs that can bind to diverse ligands, including proteins, metabolites, and other small molecules. RNA aptamers are widely used as in vitro affinity reagents. However, RNA aptamers have not been highly successful as bioactive intracellular molecules that can bind target molecules and influence cellular processes. We describe how poor RNA aptamer expression and especially poor RNA aptamer folding have limited the use of RNA aptamers in RNA synthetic biology applications. We discuss innovative new approaches that promote RNA aptamer folding in living cells and how these approaches have improved the function of aptamers in mammalian cells. These new approaches are making RNA aptamer-based synthetic biology and RNA aptamer therapeutic applications much more achievable.

Noncanonical Amino Acid Incorporation in Mice

Genetic code expansion enables in cellulo biosynthesis of curative proteins with enhanced specificity, improved stability, and even novel functions, due to the incorporation of artificial, designed, noncanonical amino acids (ncAAs). In addition, this orthogonal system also holds great potential for in vivo suppressing nonsense mutations during protein translation, providing an alternative strategy for alleviating inherited diseases caused by premature termination codons (PTCs). Here we describe the approach to explore the therapeutic efficacy and long-term safety of this strategy in transgenic mdx mice with stably expanded genetic codes. Theoretically, this method is applicable to about 11% of monogenic diseases involving nonsense mutations.

High-quality genome of the basidiomycete yeast Dioszegia hungarica PDD-24b-2 isolated from cloud water

The genome of the basidiomycete yeast Dioszegia hungarica strain PDD-24b-2 isolated from cloud water at the summit of puy de Dôme (France) was sequenced using a hybrid PacBio and Illumina sequencing strategy. The obtained assembled genome of 20.98 Mb and a GC content of 57% is structured in 16 large-scale contigs ranging from 90 kb to 5.56 Mb, and another 27.2 kb contig representing the complete circular mitochondrial genome. In total, 8,234 proteins were predicted from the genome sequence. The mitochondrial genome shows 16.2% cgu codon usage for arginine but has no canonical cognate tRNA to translate this codon. Detected transposable element (TE)-related sequences account for about 0.63% of the assembled genome. A dataset of 2,068 hand-picked public environmental metagenomes, representing over 20 Tbp of raw reads, was probed for D. hungarica related ITS sequences, and revealed worldwide distribution of this species, particularly in aerial habitats. Growth experiments suggested a psychrophilic phenotype and the ability to disperse by producing ballistospores. The high-quality assembled genome obtained for this D. hungarica strain will help investigate the behavior and ecological functions of this species in the environment.

tModBase: deciphering the landscape of tRNA modifications and their dynamic changes from epitranscriptome data

tRNA molecules contain dense, abundant modifications that affect tRNA structure, stability, mRNA decoding and tsRNA formation. tRNA modifications and related enzymes are responsive to environmental cues and are associated with a range of physiological and pathological processes. However, there is a lack of resources that can be used to mine and analyse these dynamically changing tRNA modifications. In this study, we established tModBase (https://www.tmodbase.com/) for deciphering the landscape of tRNA modification profiles from epitranscriptome data. We analysed 103 datasets generated with second- and third-generation sequencing technologies and illustrated the misincorporation and termination signals of tRNA modification sites in ten species. We thus systematically demonstrate the modification profiles across different tissues/cell lines and summarize the characteristics of tRNA-associated human diseases. By integrating transcriptome data from 32 cancers, we developed novel tools for analysing the relationships between tRNA modifications and RNA modification enzymes, the expression of 1442 tRNA-derived small RNAs (tsRNAs), and 654 DNA variations. Our database will provide new insights into the features of tRNA modifications and the biological pathways in which they participate.

Conservation and Diversification of tRNA t6A-Modifying Enzymes across the Three Domains of Life

The universal N6-threonylcarbamoyladenosine (t6A) modification occurs at position 37 of tRNAs that decipher codons starting with adenosine. Mechanistically, t6A stabilizes structural configurations of the anticodon stem loop, promotes anticodon-codon pairing and safeguards the translational fidelity. The biosynthesis of tRNA t6A is co-catalyzed by two universally conserved protein families of TsaC/Sua5 (COG0009) and TsaD/Kae1/Qri7 (COG0533). Enzymatically, TsaC/Sua5 protein utilizes the substrates of L-threonine, HCO3-/CO2 and ATP to synthesize an intermediate L-threonylcarbamoyladenylate, of which the threonylcarbamoyl-moiety is subsequently transferred onto the A37 of substrate tRNAs by the TsaD-TsaB -TsaE complex in bacteria or by the KEOPS complex in archaea and eukaryotic cytoplasm, whereas Qri7/OSGEPL1 protein functions on its own in mitochondria. Depletion of tRNA t6A interferes with protein homeostasis and gravely affects the life of unicellular organisms and the fitness of higher eukaryotes. Pathogenic mutations of YRDC, OSGEPL1 and KEOPS are implicated in a number of human mitochondrial and neurological diseases, including autosomal recessive Galloway-Mowat syndrome. The molecular mechanisms underscoring both the biosynthesis and cellular roles of tRNA t6A are presently not well elucidated. This review summarizes current mechanistic understandings of the catalysis, regulation and disease implications of tRNA t6A-biosynthetic machineries of three kingdoms of life, with a special focus on delineating the structure-function relationship from perspectives of conservation and diversity.

Emerging functions of tRNA modifications in mRNA translation and diseases

tRNAs are essential modulators that recognize mRNA codons and bridge amino acids for mRNA translation. The tRNAs are heavily modified, which is essential for forming a complex secondary structure that facilitates codon recognition and mRNA translation. In recent years, studies have identified the regulatory roles of tRNA modifications in mRNA translation networks. Misregulation of tRNA modifications is closely related to the progression of developmental diseases and cancers. In this review, we summarize the tRNA biogenesis process and then discuss the effects and mechanisms of tRNA modifications on tRNA processing and mRNA translation. Finally, we provide a comprehensive overview of tRNA modifications' physiological and pathological functions, focusing on diseases including cancers.

Efficient CRISPR/Cas9 system based on autonomously replicating plasmid with an AMA1 sequence and precisely targeted gene deletion in the edible fungus, Cordyceps militaris

Cordyceps militaris is a popular edible fungus with important economic value worldwide. In this study, an efficient CRISPR/Cas9 genome-editing system based on an autonomously replicating plasmid with an AMA1 sequence was constructed. Further, a precisely targeted gene deletion via homology-directed repair was effectively introduced in C. militaris. Gene editing was successful, with efficiencies of 55.1% and 89% for Cmwc-1 and Cmvvd, respectively. Precisely targeted gene deletion was achieved at an efficiency of 73.9% by a single guide RNA supplementation with donor DNAs. Double genes, Cmwc-1 and Cmvvd, were edited simultaneously with an efficiency of 10%. Plasmid loss was observed under non-selective culture conditions, which could permit recycling of the selectable marker and avoid the adverse effects of the CRISPR/Cas9 system on the fungus, which is beneficial for the generation of new cultivars. RNA Pol III promoters, endogenous tRNAPro of C. militaris, and chimeric AfU6-tRNAGly can be used to improve the efficiency. Polyethylene glycol-mediated protoplast transformation was markedly more efficient than Agrobacterium tumefaciens-mediated transformation of C. militaris. To our knowledge, this is the first description of genome editing and precisely targeted gene deletion in mushrooms based on AMA1 plasmids. Our findings will enable the modification of multiple genes in both functional genomics research and strain breeding.

Tied up in knots: Untangling substrate recognition by the SPOUT methyltransferases

The SpoU-TrmD (SPOUT) methyltransferase superfamily was designated when structural similarity was identified between the transfer RNA-modifying enzymes TrmH (SpoU) and TrmD. SPOUT methyltransferases are found in all domains of life and predominantly modify transfer RNA or ribosomal RNA substrates, though one instance of an enzyme with a protein substrate has been reported. Modifications placed by SPOUT methyltransferases play diverse roles in regulating cellular processes such as ensuring translational fidelity, altering RNA stability, and conferring bacterial resistance to antibiotics. This large collection of S-adenosyl-L-methionine-dependent methyltransferases is defined by a unique α/β fold with a deep trefoil knot in their catalytic (SPOUT) domain. Herein, we describe current knowledge of SPOUT enzyme structure, domain architecture, and key elements of catalytic function, including S-adenosyl-L-methionine co-substrate binding, beginning with a new sequence alignment that divides the SPOUT methyltransferase superfamily into four major clades. Finally, a major focus of this review will be on our growing understanding of how these diverse enzymes accomplish the molecular feat of specific substrate recognition and modification, as highlighted by recent advances in our knowledge of protein-RNA complex structures and the discovery of the dependence of one SPOUT methyltransferase on metal ion binding for catalysis. Considering the broad biological roles of RNA modifications, developing a deeper understanding of the process of substrate recognition by the SPOUT enzymes will be critical for defining many facets of fundamental RNA biology with implications for human disease.

A Test and Refinement of Folding Free Energy Nearest Neighbor Parameters for RNA Including N6-Methyladenosine

RNA folding free energy change parameters are widely used to predict RNA secondary structure and to design RNA sequences. These parameters include terms for the folding free energies of helices and loops. Although the full set of parameters has only been traditionally available for the four common bases and backbone, it is well known that covalent modifications of nucleotides are widespread in natural RNAs. Covalent modifications are also widely used in engineered sequences. We recently derived a full set of nearest neighbor terms for RNA that includes N6-methyladenosine (m6A). In this work, we test the model using 98 optical melting experiments, matching duplexes with or without N6-methylation of A. Most experiments place RRACH, the consensus site of N6-methylation, in a variety of contexts, including helices, bulge loops, internal loops, dangling ends, and terminal mismatches. For matched sets of experiments that include either A or m6A in the same context, we find that the parameters for m6A are as accurate as those for A. Across all experiments, the root mean squared deviation between estimated and experimental free energy changes is 0.67 kcal/mol. We used the new experimental data to refine the set of nearest neighbor parameter terms for m6A. These parameters enable prediction of RNA secondary structures including m6A, which can be used to model how N6-methylation of A affects RNA structure.

Transfer RNA-derived small RNAs and their potential roles in the therapeutic heterogeneity of sacubitril/valsartan in heart failure patients after acute myocardial infarction

Background

It has been reported that sacubitril/valsartan can improve cardiac function in acute myocardial infarction (AMI) patients complicated by heart failure (HF). However, a number of patients cannot be treated successfully; this phenomenon is called sacubitril/valsartan resistance (SVR), and the mechanisms remain unclear.

Methods

In our present research, the expression profiles of transfer RNA (tRNA)-derived small RNAs (tsRNAs) in SVR along with no sacubitril/valsartan resistance (NSVR) patients were determined by RNA sequencing. Through bioinformatics, quantitative real-time PCR (qRT-PCR), and cell-based experiments, we identified SVR-related tsRNAs and confirmed their diagnostic value, predicted their targeted genes, and explored the enriched signal pathways as well as regulatory roles of tsRNAs in SVR.

Results

Our research indicated that 36 tsRNAs were upregulated and that 21 tsRNAs were downregulated in SVR. Among these tsRNAs, the expression of tRF-59:76-Tyr-GTA-2-M3 and tRF-60:76-Val-AAC-1-M5 was upregulated, while the expression of tRF-1:29-Gly-GCC-1 was downregulated in the group of SVR. Receiver operating characteristic (ROC) curve analysis demonstrated that these three tsRNAs were potential biomarkers of the therapeutic heterogeneity of sacubitril/valsartan. Moreover, tRF-60:76-Val-AAC-1-M5 might target Tnfrsf10b and Bcl2l1 to influence the observed therapeutic heterogeneity through the lipid and atherosclerosis signaling pathways.

Conclusion

Hence, tsRNA might play a vital role in SVR. These discoveries provide new insights for the mechanistic investigation of responsiveness to sacubitril/valsartan.

Three tRNA nuclear exporters in S. cerevisiae: parallel pathways, preferences, and precision

tRNAs that are transcribed in the nucleus are exported to the cytoplasm to perform their iterative essential function in translation. However, the complex set of tRNA post-transcriptional processing and subcellular trafficking steps are not completely understood. In particular, proteins involved in tRNA nuclear export remain unknown since the canonical tRNA nuclear exportin, Los1/Exportin-t, is unessential in all tested organisms. We previously reported that budding yeast Mex67-Mtr2, a mRNA nuclear exporter, co-functions with Los1 in tRNA nuclear export. Here we employed in vivo co-purification of tRNAs with endogenously expressed nuclear exporters to document that Crm1 also is a bona fide tRNA nuclear exporter. We document that Los1, Mex67-Mtr2 and Crm1 possess individual tRNA preferences for forming nuclear export complexes with members of the 10 families of intron-containing pre-tRNAs. Remarkably, Mex67-Mtr2, but not Los1 or Crm1, is error-prone, delivering tRNAs to the cytoplasm prior to 5′ leader removal. tRNA retrograde nuclear import functions to monitor the aberrant leader-containing spliced tRNAs, returning them to the nucleus where they are degraded by 3′ to 5′ exonucleases. Overall, our work identifies a new tRNA nuclear exporter, uncovers exporter preferences for specific tRNA families, and documents contribution of tRNA nuclear import to tRNA quality control.

Recent insights into the structure, function, and regulation of the eukaryotic transfer RNA splicing endonuclease complex

The splicing of transfer RNA (tRNA) introns is a critical step of tRNA maturation, for intron-containing tRNAs. In eukaryotes, tRNA splicing is a multi-step process that relies on several RNA processing enzymes to facilitate intron removal and exon ligation. Splicing is initiated by the tRNA splicing endonuclease (TSEN) complex which catalyzes the excision of the intron through its two nuclease subunits. Mutations in all four subunits of the TSEN complex are linked to a family of neurodegenerative and neurodevelopmental diseases known as pontocerebellar hypoplasia (PCH). Recent studies provide molecular insights into the structure, function, and regulation of the eukaryotic TSEN complex and are beginning to illuminate how mutations in the TSEN complex lead to neurodegenerative disease. Using new advancements in the prediction of protein structure, we created a three-dimensional model of the human TSEN complex. We review functions of the TSEN complex beyond tRNA splicing by highlighting recently identified substrates of the eukaryotic TSEN complex and discuss mechanisms for the regulation of tRNA splicing, by enzymes that modify cleaved tRNA exons and introns. Finally, we review recent biochemical and animal models that have worked to address the mechanisms that drive PCH and synthesize these studies with previous studies to try to better understand PCH pathogenesis. This article is categorized under: RNA Processing > tRNA Processing RNA in Disease and Development > RNA in Disease RNA Interactions with Proteins and Other Molecules > Protein-RNA Recognition.

Human TRUB1 is a highly conserved pseudouridine synthase responsible for the formation of Ψ55 in mitochondrial tRNAAsn, tRNAGln, tRNAGlu and tRNAPro

Pseudouridine (Ψ) at position 55 in tRNAs plays an important role in their structure and function. This modification is catalyzed by TruB/Pus4/Cbf5 family of pseudouridine synthases in bacteria and yeast. However, the mechanism of TRUB family underlying the formation of Ψ55 in the mammalian tRNAs is largely unknown. In this report, the CMC/reverse transcription assays demonstrated the presence of Ψ55 in the human mitochondrial tRNAAsn, tRNAGln, tRNAGlu, tRNAPro, tRNAMet, tRNALeu(UUR) and tRNASer(UCN). TRUB1 knockout (KO) cell lines generated by CRISPR/Cas9 technology exhibited the loss of Ψ55 modification in mitochondrial tRNAAsn, tRNAGln, tRNAGlu and tRNAPro but did not affect other 18 mitochondrial tRNAs. An in vitro assay revealed that recombinant TRUB1 protein can catalyze the efficient formation of Ψ55 in tRNAAsn and tRNAGln, but not in tRNAMet and tRNAArg. Notably, the overexpression of TRUB1 cDNA reversed the deficient Ψ55 modifications in these tRNAs in TRUB1KO HeLa cells. TRUB1 deficiency affected the base-pairing (18A/G-Ψ55), conformation and stability but not aminoacylation capacity of these tRNAs. Furthermore, TRUB1 deficiency impacted mitochondrial translation and biogenesis of oxidative phosphorylation system. Our findings demonstrated that human TRUB1 is a highly conserved mitochondrial pseudouridine synthase responsible for the Ψ55 modification in the mitochondrial tRNAAsn, tRNAGln, tRNAGlu and tRNAPro.

Oligodendrocyte differentiation alters tRNA modifications and codon optimality-mediated mRNA decay

Oligodendrocytes are specialized cells that confer neuronal myelination in the central nervous system. Leukodystrophies associated with oligodendrocyte deficits and hypomyelination are known to result when a number of tRNA metabolism genes are mutated. Thus, for unknown reasons, oligodendrocytes may be hypersensitive to perturbations in tRNA biology. In this study, we survey the tRNA transcriptome in the murine oligodendrocyte cell lineage and find that specific tRNAs are hypomodified in oligodendrocytes within or near the anticodon compared to oligodendrocyte progenitor cells (OPCs). This hypomodified state may be the result of differential expression of key modification enzymes during oligodendrocyte differentiation. Moreover, we observe a concomitant relationship between tRNA hypomodification and tRNA decoding potential; observing oligodendrocyte specific alterations in codon optimality-mediated mRNA decay and ribosome transit. Our results reveal that oligodendrocytes naturally maintain a delicate, hypersensitized tRNA/mRNA axis. We suggest this axis is a potential mediator of pathology in leukodystrophies and white matter disease when further insult to tRNA metabolism is introduced.

A unifying network modeling approach for codon optimization

MOTIVATION

Synthesizing genes to be expressed in other organisms is an essential tool in biotechnology. While the many-to-one mapping from codons to amino acids makes the genetic code degenerate, codon usage in a particular organism is not random either. This bias in codon use may have a remarkable effect on the level of gene expression. A number of measures have been developed to quantify a given codon sequence's strength to express a gene in a host organism. Codon optimization aims to find a codon sequence that will optimize one or more of these measures. Efficient computational approaches are needed since the possible number of codon sequences grows exponentially as the number of amino acids increases.

RESULTS

We develop a unifying modeling approach for codon optimization. With our mathematical formulations based on graph/network representations of amino acid sequences, any combination of measures can be optimized in the same framework by finding a path satisfying additional limitations in an acyclic layered network. We tested our approach on bi-objectives commonly used in the literature, namely, Codon Pair Bias versus Codon Adaptation Index and Relative Codon Pair Bias versus Relative Codon Bias. However, our framework is general enough to handle any number of objectives concurrently with certain restrictions or preferences on the use of specific nucleotide sequences. We implemented our models using Python's Gurobi interface and showed the efficacy of our approach even for the largest proteins available. We also provided experimentation showing that highly expressed genes have objective values close to the optimized values in the bi-objective codon design problem.

AVAILABILITY AND IMPLEMENTATION

http://alpersen.bilkent.edu.tr/NetworkCodon.zip.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Cloning and Sequencing Eukaryotic Small RNAs

Small RNAs are ubiquitous regulators of gene expression that participate in nearly all aspects of physiology in a wide range of organisms. There are many different classes of eukaryotic small RNAs that play regulatory roles at every level of gene expression, including transcription, RNA stability, and translation. While eukaryotic small RNAs display diverse functions across and within classes, they are generally grouped functionally based on the machinery required for their biogenesis, the effector proteins they associate with, and their molecular characteristics. The development of techniques to clone and sequence small RNAs has been critical for their identification, yet the ligation-dependent addition of RNA adapters and the use of reverse transcriptase to generate cDNA in traditional library preparation protocols can be unsuitable to detect certain small RNA subtypes. In particular, 3' or 5' chemical modifications that are characteristic of specific types of small RNAs can impede the ligation-dependent addition of RNA adapters, while internal RNA modifications can interfere with accurate reverse transcription. The inability to clone certain small RNA subtypes with traditional protocols results in an inaccurate assessment of small RNA abundance and diversity, where some RNAs appear over-represented and others are not detected. This overview aims to guide users on how to design small RNA cloning workflows in eukaryotes to more accurately capture specific small RNAs of interest. Hence, we discuss the molecular biology underlying the identification and quantitation of small RNAs, explore the limitations of commonly used protocols, and detail the alternative approaches that can be used to enrich specific small RNA classes. © 2022 Wiley Periodicals LLC.

The tRNA regulome in neurodevelopmental and neuropsychiatric disease

Transfer (t)RNAs are 70-90 nucleotide small RNAs highly regulated by 43 different types of epitranscriptomic modifications and requiring aminoacylation ('charging') for mRNA decoding and protein synthesis. Smaller cleavage products of mature tRNAs, or tRNA fragments, have been linked to a broad variety of noncanonical functions, including translational inhibition and modulation of the immune response. Traditionally, knowledge about tRNA regulation in brain is derived from phenotypic exploration of monogenic neurodevelopmental and neurodegenerative diseases associated with rare mutations in tRNA modification genes. More recent studies point to the previously unrecognized potential of the tRNA regulome to affect memory, synaptic plasticity, and affective states. For example, in mature cortical neurons, cytosine methylation sensitivity of the glycine tRNA family (tRNA^Gly) is coupled to glycine biosynthesis and codon-specific alterations in ribosomal translation together with robust changes in cognition and depression-related behaviors. In this Review, we will discuss the emerging knowledge of the neuronal tRNA landscape, with a focus on epitranscriptomic tRNA modifications and downstream molecular pathways affected by alterations in tRNA expression, charging levels, and cleavage while mechanistically linking these pathways to neuropsychiatric disease and provide insight into future areas of study for this field.

Initiator tRNA lacking 1-methyladenosine is targeted by the rapid tRNA decay pathway in evolutionarily distant yeast species

All tRNAs have numerous modifications, lack of which often results in growth defects in the budding yeast Saccharomyces cerevisiae and neurological or other disorders in humans. In S. cerevisiae, lack of tRNA body modifications can lead to impaired tRNA stability and decay of a subset of the hypomodified tRNAs. Mutants lacking 7-methylguanosine at G46 (m7G46), N2,N2-dimethylguanosine (m2,2G26), or 4-acetylcytidine (ac4C12), in combination with other body modification mutants, target certain mature hypomodified tRNAs to the rapid tRNA decay (RTD) pathway, catalyzed by 5’-3’ exonucleases Xrn1 and Rat1, and regulated by Met22. The RTD pathway is conserved in the phylogenetically distant fission yeast Schizosaccharomyces pombe for mutants lacking m7G46. In contrast, S. cerevisiae trm6/gcd10 mutants with reduced 1-methyladenosine (m1A58) specifically target pre-tRNAiMet(CAU) to the nuclear surveillance pathway for 3’-5’ exonucleolytic decay by the TRAMP complex and nuclear exosome. We show here that the RTD pathway has an unexpected major role in the biology of m1A58 and tRNAiMet(CAU) in both S. pombe and S. cerevisiae. We find that S. pombe trm6Δ mutants lacking m1A58 are temperature sensitive due to decay of tRNAiMet(CAU) by the RTD pathway. Thus, trm6Δ mutants had reduced levels of tRNAiMet(CAU) and not of eight other tested tRNAs, overexpression of tRNAiMet(CAU) restored growth, and spontaneous suppressors that restored tRNAiMet(CAU) levels had mutations in dhp1/RAT1 or tol1/MET22. In addition, deletion of cid14/TRF4 in the nuclear surveillance pathway did not restore growth. Furthermore, re-examination of S. cerevisiae trm6 mutants revealed a major role of the RTD pathway in maintaining tRNAiMet(CAU) levels, in addition to the known role of the nuclear surveillance pathway. These findings provide evidence for the importance of m1A58 in the biology of tRNAiMet(CAU) throughout eukaryotes, and fuel speculation that the RTD pathway has a major role in quality control of body modification mutants throughout fungi and other eukaryotes.

Artificial Intelligence, Healthcare, Clinical Genomics, and Pharmacogenomics Approaches in Precision Medicine

Precision medicine has greatly aided in improving health outcomes using earlier diagnosis and better prognosis for chronic diseases. It makes use of clinical data associated with the patient as well as their multi-omics/genomic data to reach a conclusion regarding how a physician should proceed with a specific treatment. Compared to the symptom-driven approach in medicine, precision medicine considers the critical fact that all patients do not react to the same treatment or medication in the same way. When considering the intersection of traditionally distinct arenas of medicine, that is, artificial intelligence, healthcare, clinical genomics, and pharmacogenomics—what ties them together is their impact on the development of precision medicine as a field and how they each contribute to patient-specific, rather than symptom-specific patient outcomes. This study discusses the impact and integration of these different fields in the scope of precision medicine and how they can be used in preventing and predicting acute or chronic diseases. Additionally, this study also discusses the advantages as well as the current challenges associated with artificial intelligence, healthcare, clinical genomics, and pharmacogenomics.

Coexpression of MEIOTIC-TOPOISOMERASE VIB-dCas9 with guide RNAs specific to a recombination hotspot is insufficient to increase crossover frequency in Arabidopsis

During meiosis, homologous chromosomes pair and recombine, which can result in reciprocal crossovers that increase genetic diversity. Crossovers are unevenly distributed along eukaryote chromosomes and show repression in heterochromatin and the centromeres. Within the chromosome arms, crossovers are often concentrated in hotspots, which are typically in the kilobase range. The uneven distribution of crossovers along chromosomes, together with their low number per meiosis, creates a limitation during crop breeding, where recombination can be beneficial. Therefore, targeting crossovers to specific genome locations has the potential to accelerate crop improvement. In plants, meiotic crossovers are initiated by DNA double-strand breaks that are catalyzed by SPO11 complexes, which consist of 2 catalytic (SPO11-1 and SPO11-2) and 2 noncatalytic subunits (MTOPVIB). We used the model plant Arabidopsis thaliana to coexpress an MTOPVIB-dCas9 fusion protein with guide RNAs specific to the 3a crossover hotspot. We observed that this was insufficient to significantly change meiotic crossover frequency or pattern within 3a. We discuss the implications of our findings for targeting meiotic recombination within plant genomes.

Structural and mechanistic basis of RNA processing by protein-only ribonuclease P enzymes

Ribonuclease P (RNase P) enzymes are responsible for the 5' processing of tRNA precursors. In addition to the well-characterised ribozyme-based RNase P enzymes, an evolutionarily distinct group of protein-only RNase Ps exists. These proteinaceous RNase Ps (PRORPs) can be found in all three domains of life and can be divided into two structurally different types: eukaryotic and prokaryotic. Recent structural studies on members of both families reveal a surprising diversity of molecular architectures, but also highlight conceptual and mechanistic similarities. Here, we provide a comparison between the different types of PRORP enzymes and review how the combination of structural, biochemical, and biophysical studies has led to a molecular picture of protein-mediated tRNA processing.

Eukaryotic tRNA splicing - one goal, two strategies, many players

Transfer RNAs (tRNAs) are transcribed as precursor molecules that undergo several maturation steps before becoming functional for protein synthesis. One such processing mechanism is the enzyme-catalysed splicing of intron-containing pre-tRNAs. Eukaryotic tRNA splicing is an essential process since intron-containing tRNAs cannot fulfil their canonical function at the ribosome. Splicing of pre-tRNAs occurs in two steps: The introns are first excised by a tRNA-splicing endonuclease and the exons are subsequently sealed by an RNA ligase. An intriguing complexity has emerged from newly identified tRNA splicing factors and their interplay with other RNA processing pathways during the past few years. This review summarises our current understanding of eukaryotic tRNA splicing and the underlying enzyme machinery. We highlight recent structural advances and how they have shaped our mechanistic understanding of tRNA splicing in eukaryotic cells. A special focus lies on biochemically distinct strategies for exon-exon ligation in fungi versus metazoans.

How RNases Shape Mitochondrial Transcriptomes

Mitochondria are the power houses of eukaryote cells. These endosymbiotic organelles of prokaryote origin are considered as semi-autonomous since they have retained a genome and fully functional gene expression mechanisms. These pathways are particularly interesting because they combine features inherited from the bacterial ancestor of mitochondria with characteristics that appeared during eukaryote evolution. RNA biology is thus particularly diverse in mitochondria. It involves an unexpectedly vast array of factors, some of which being universal to all mitochondria and others being specific from specific eukaryote clades. Among them, ribonucleases are particularly prominent. They play pivotal functions such as the maturation of transcript ends, RNA degradation and surveillance functions that are required to attain the pool of mature RNAs required to synthesize essential mitochondrial proteins such as respiratory chain proteins. Beyond these functions, mitochondrial ribonucleases are also involved in the maintenance and replication of mitochondrial DNA, and even possibly in the biogenesis of mitochondrial ribosomes. The diversity of mitochondrial RNases is reviewed here, showing for instance how in some cases a bacterial-type enzyme was kept in some eukaryotes, while in other clades, eukaryote specific enzymes were recruited for the same function.

tiRNAs: Insights into Their Biogenesis, Functions, and Future Applications in Livestock Research

Transfer RNA (tRNA)-derived small RNAs (tsRNAs) belong to a group of transfer ribonucleic acid (tRNA)-derived fragments that have recently gained interest as molecules with specific biological functions. Their involvement in the regulation of physiological processes and pathological phenotypes suggests molecular roles similar to those of miRNAs. tsRNA biogenesis under specific physiological conditions will offer new perspectives in understanding diseases, and may provide new sources for biological marker design to determine and monitor the health status of farm animals. In this review, we focus on the latest discoveries about tsRNAs and give special attention to molecules initially thought to be mainly associated with tRNA-derived stress-induced RNAs (tiRNAs). We present an outline of their biological functions, offer a collection of useful databases, and discuss future research perspectives and applications in livestock basic and applied research.

PltRNAdb: Plant transfer RNA database

Transfer RNAs (tRNAs) are intermediate-sized non-coding RNAs found in all organisms that help translate messenger RNA into protein. Recently, the number of sequenced plant genomes has increased dramatically. The availability of this extensive data greatly accelerates the study of tRNAs on a large scale. Here, 8,768,261 scaffolds/chromosomes containing 229,093 giga-base pairs representing whole-genome sequences of 256 plant species were analyzed to identify tRNA genes. As a result, 331,242 nuclear, 3,216 chloroplast, and 1,467 mitochondrial tRNA genes were identified. The nuclear tRNA genes include 275,134 tRNAs decoding 20 standard amino acids, 1,325 suppressor tRNAs, 6,273 tRNAs with unknown isotypes, 48,475 predicted pseudogenes, and 37,873 tRNAs with introns. Efforts also extended to the creation of PltRNAdb (https://bioinformatics.um6p.ma/PltRNAdb/index.php), a data source for tRNA genes from 256 plant species. PltRNAdb website allows researchers to search, browse, visualize, BLAST, and download predicted tRNA genes. PltRNAdb will help improve our understanding of plant tRNAs and open the door to discovering the unknown regulatory roles of tRNAs in plant genomes.

A multiplex platform for small RNA sequencing elucidates multifaceted tRNA stress response and translational regulation

Small RNAs include tRNA, snRNA, micro-RNA, tRNA fragments and others that constitute > 90% of RNA copy numbers in a human cell and perform many essential functions. Popular small RNA-seq strategies limit the insights into coordinated small RNA response to cellular stress. Small RNA-seq also lacks multiplexing capabilities. Here, we report a multiplex small RNA-seq library preparation method (MSR-seq) to investigate cellular small RNA and mRNA response to heat shock, hydrogen peroxide, and arsenite stress. Comparing stress-induced changes of total cellular RNA and polysome-associated RNA, we identify a coordinated tRNA response that involves polysome-specific tRNA abundance and synergistic N3-methylcytosine (m³C) tRNA modification. Combining tRNA and mRNA response to stress we reveal a mechanism of stress-induced down-regulation in translational elongation. We also find that native tRNA molecules lacking several modifications are biased reservoirs for the biogenesis of tRNA fragments. Our results demonstrate the importance of simultaneous investigation of small RNAs and their modifications in response to varying biological conditions.

Reversible RNA phosphorylation stabilizes tRNA for cellular thermotolerance

Post-transcriptional modifications have critical roles in tRNA stability and function^1–4. In thermophiles, tRNAs are heavily modified to maintain their thermal stability under extreme growth temperatures^5,6. Here we identified 2′-phosphouridine (U^p) at position 47 of tRNAs from thermophilic archaea. U^p47 confers thermal stability and nuclease resistance to tRNAs. Atomic structures of native archaeal tRNA showed a unique metastable core structure stabilized by U^p47. The 2′-phosphate of U^p47 protrudes from the tRNA core and prevents backbone rotation during thermal denaturation. In addition, we identified the arkI gene, which encodes an archaeal RNA kinase responsible for U^p47 formation. Structural studies showed that ArkI has a non-canonical kinase motif surrounded by a positively charged patch for tRNA binding. A knockout strain of arkI grew slowly at high temperatures and exhibited a synthetic growth defect when a second tRNA-modifying enzyme was depleted. We also identified an archaeal homologue of KptA as an eraser that efficiently dephosphorylates U^p47 in vitro and in vivo. Taken together, our findings show that U^p47 is a reversible RNA modification mediated by ArkI and KptA that fine-tunes the structural rigidity of tRNAs under extreme environmental conditions.

Reversible internal RNA phosphrylation contributes to thermal stability and nuclease resistance of tRNA, and cellular thermotolerance of hyperthermophiles. 

Identification of a Trm732 Motif Required for 2'-O-methylation of the tRNA Anticodon Loop by Trm7

Posttranscriptional tRNA modifications are essential for proper gene expression, and defects in the enzymes that perform tRNA modifications are associated with numerous human disorders. Throughout eukaryotes, 2'-O-methylation of residues 32 and 34 of the anticodon loop of tRNA is important for proper translation, and in humans, a lack of these modifications results in non-syndromic X-linked intellectual disability. In yeast, the methyltransferase Trm7 forms a complex with Trm732 to 2'-O-methylate tRNA residue 32 and with Trm734 to 2'-O-methylate tRNA residue 34. Trm732 and Trm734 are required for the methylation activity of Trm7, but the role of these auxiliary proteins is not clear. Additionally, Trm732 and Trm734 homologs are implicated in biological processes not directly related to translation, suggesting that these proteins may have additional cellular functions. To identify critical amino acids in Trm732, we generated variants and tested their ability to function in yeast cells. We identified a conserved RRSAGLP motif in the conserved DUF2428 domain of Trm732 that is required for tRNA modification activity by both yeast Trm732 and its human homolog, THADA. The identification of Trm732 variants that lack tRNA modification activity will help to determine if other biological functions ascribed to Trm732 and THADA are directly due to tRNA modification or to secondary effects due to other functions of these proteins.

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.

Synthetic riboswitches for the analysis of tRNA processing by eukaryotic RNase P enzymes

Removal of the 5'-leader region is an essential step in the maturation of tRNA molecules in all domains of life. This reaction is catalyzed by various RNase P activities, ranging from ribonucleoproteins with ribozyme activity to protein-only forms. In Escherichia coli, the efficiency of RNase P-mediated cleavage can be controlled by computationally designed riboswitch elements in a ligand-dependent way, where the 5'-leader sequence of a tRNA precursor is either sequestered in a hairpin structure or presented as a single-stranded region accessible for maturation. In the presented work, the regulatory potential of such artificial constructs is tested on different forms of eukaryotic RNase P enzymes-two protein-only RNase P enzymes (PRORP1 and PRORP2) from Arabidopsis thaliana and the ribonucleoprotein of Homo sapiens The PRORP enzymes were analyzed in vitro as well as in vivo in a bacterial RNase P complementation system. We also tested in HEK293T cells whether the riboswitches remain functional with human nuclear RNase P. While the regulatory principle of the synthetic riboswitches applies for all tested RNase P enzymes, the results also show differences in the substrate requirements of the individual enzyme versions. Hence, such designed RNase P riboswitches represent a novel tool to investigate the impact of the structural composition of the 5'-leader on substrate recognition by different types of RNase P enzymes.

Mutations in Drosophila tRNA processing factors cause phenotypes similar to Pontocerebellar Hypoplasia

Mature transfer (t)RNAs are generated by multiple RNA processing events, which can include the excision of intervening sequences. The tRNA splicing endonuclease (TSEN) complex is responsible for cleaving these intron-containing pre-tRNA transcripts. In humans, TSEN copurifies with CLP1, an RNA kinase. Despite extensive work on CLP1, its in vivo connection to tRNA splicing remains unclear. Interestingly, mutations in CLP1 or TSEN genes cause neurological diseases in humans that are collectively termed Pontocerebellar Hypoplasia (PCH). In mice, loss of Clp1 kinase activity results in premature death, microcephaly and progressive loss of motor function. To determine if similar phenotypes are observed in Drosophila, we characterized mutations in crowded-by-cid (cbc), the CLP1 ortholog, as well as in the fly ortholog of human TSEN54. Analyses of organismal viability, larval locomotion and brain size revealed that mutations in both cbc and Tsen54 phenocopy those in mammals in several details. In addition to an overall reduction in brain lobe size, we also found increased cell death in mutant larval brains. Ubiquitous or tissue-specific knockdown of cbc in neurons and muscles reduced viability and locomotor function. These findings indicate that we can successfully model PCH in a genetically-tractable invertebrate.

Mutations in Drosophila tRNA processing factors cause phenotypes similar to Pontocerebellar Hypoplasia

Mature transfer (t)RNAs are generated by multiple RNA processing events, which can include the excision of intervening sequences. The tRNA splicing endonuclease (TSEN) complex is responsible for cleaving these intron-containing pre-tRNA transcripts. In humans, TSEN copurifies with CLP1, an RNA kinase. Despite extensive work on CLP1, its in vivo connection to tRNA splicing remains unclear. Interestingly, mutations in CLP1 or TSEN genes cause neurological diseases in humans that are collectively termed Pontocerebellar Hypoplasia (PCH). In mice, loss of Clp1 kinase activity results in premature death, microcephaly and progressive loss of motor function. To determine if similar phenotypes are observed in Drosophila, we characterized mutations in crowded-by-cid (cbc), the CLP1 ortholog, as well as in the fly ortholog of human TSEN54. Analyses of organismal viability, larval locomotion and brain size revealed that mutations in both cbc and Tsen54 phenocopy those in mammals in several details. In addition to an overall reduction in brain lobe size, we also found increased cell death in mutant larval brains. Ubiquitous or tissue-specific knockdown of cbc in neurons and muscles reduced viability and locomotor function. These findings indicate that we can successfully model PCH in a genetically-tractable invertebrate.

Secondary structure prediction for RNA sequences including N6-methyladenosine

There is increasing interest in the roles of covalently modified nucleotides in RNA. There has been, however, an inability to account for modifications in secondary structure prediction because of a lack of software and thermodynamic parameters. We report the solution for these issues for N⁶-methyladenosine (m⁶A), allowing secondary structure prediction for an alphabet of A, C, G, U, and m⁶A. The RNAstructure software now works with user-defined nucleotide alphabets of any size. We also report a set of nearest neighbor parameters for helices and loops containing m⁶A, using experiments. Interestingly, N⁶-methylation decreases folding stability for adenosines in the middle of a helix, has little effect on folding stability for adenosines at the ends of helices, and increases folding stability for unpaired adenosines stacked on a helix. We demonstrate predictions for an N⁶-methylation-activated protein recognition site from MALAT1 and human transcriptome-wide effects of N⁶-methylation on the probability of adenosine being buried in a helix.

RNA folding free energy nearest neighbor parameters were determined for sequences with the nucleotide m⁶A. The RNAstructure software package can accommodate modified nucleotides, enabling secondary structure prediction of sequences with m⁶A.

RNA modifications stabilize the tertiary structure of tRNAfMet by locally increasing conformational dynamics

A plethora of modified nucleotides extends the chemical and conformational space for natural occurring RNAs. tRNAs constitute the class of RNAs with the highest modification rate. The extensive modification modulates their overall stability, the fidelity and efficiency of translation. However, the impact of nucleotide modifications on the local structural dynamics is not well characterized. Here we show that the incorporation of the modified nucleotides in tRNAfMet from Escherichia coli leads to an increase in the local conformational dynamics, ultimately resulting in the stabilization of the overall tertiary structure. Through analysis of the local dynamics by NMR spectroscopic methods we find that, although the overall thermal stability of the tRNA is higher for the modified molecule, the conformational fluctuations on the local level are increased in comparison to an unmodified tRNA. In consequence, the melting of individual base pairs in the unmodified tRNA is determined by high entropic penalties compared to the modified. Further, we find that the modifications lead to a stabilization of long-range interactions harmonizing the stability of the tRNA's secondary and tertiary structure. Our results demonstrate that the increase in chemical space through introduction of modifications enables the population of otherwise inaccessible conformational substates.

Making Sense of "Nonsense" and More: Challenges and Opportunities in the Genetic Code Expansion, in the World of tRNA Modifications

The evolutional development of the RNA translation process that leads to protein synthesis based on naturally occurring amino acids has its continuation via synthetic biology, the so-called rational bioengineering. Genetic code expansion (GCE) explores beyond the natural translational processes to further enhance the structural properties and augment the functionality of a wide range of proteins. Prokaryotic and eukaryotic ribosomal machinery have been proven to accept engineered tRNAs from orthogonal organisms to efficiently incorporate noncanonical amino acids (ncAAs) with rationally designed side chains. These side chains can be reactive or functional groups, which can be extensively utilized in biochemical, biophysical, and cellular studies. Genetic code extension offers the contingency of introducing more than one ncAA into protein through frameshift suppression, multi-site-specific incorporation of ncAAs, thereby increasing the vast number of possible applications. However, different mediating factors reduce the yield and efficiency of ncAA incorporation into synthetic proteins. In this review, we comment on the recent advancements in genetic code expansion to signify the relevance of systems biology in improving ncAA incorporation efficiency. We discuss the emerging impact of tRNA modifications and metabolism in protein design. We also provide examples of the latest successful accomplishments in synthetic protein therapeutics and show how codon expansion has been employed in various scientific and biotechnological applications.