Post-transcriptional nucleotide modification and alternative folding of RNA

RNA architecture influences plant biology

The majority of the genome is transcribed to RNA in living organisms. RNA transcripts can form astonishing arrays of secondary and tertiary structures via Watson-Crick, Hoogsteen, or wobble base pairing. In vivo, RNA folding is not a simple thermodynamic event of minimizing free energy. Instead, the process is constrained by transcription, RNA-binding proteins, steric factors, and the microenvironment. RNA secondary structure (RSS) plays myriad roles in numerous biological processes, such as RNA processing, stability, transportation, and translation in prokaryotes and eukaryotes. Emerging evidence has also implicated RSS in RNA trafficking, liquid-liquid phase separation, and plant responses to environmental variations such as temperature and salinity. At molecular level, RSS is correlated with splicing, polyadenylation, protein synthesis, and miRNA biogenesis and functions. In this review, we summarize newly reported methods for probing RSS in vivo and functions and mechanisms of RSS in plant physiology.

A novel SHAPE reagent enables the analysis of RNA structure in living cells with unprecedented accuracy

Due to the mounting evidence that RNA structure plays a critical role in regulating almost any physiological as well as pathological process, being able to accurately define the folding of RNA molecules within living cells has become a crucial need. We introduce here 2-aminopyridine-3-carboxylic acid imidazolide (2A3), as a general probe for the interrogation of RNA structures in vivo. 2A3 shows moderate improvements with respect to the state-of-the-art selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) reagent NAI on naked RNA under in vitro conditions, but it significantly outperforms NAI when probing RNA structure in vivo, particularly in bacteria, underlining its increased ability to permeate biological membranes. When used as a restraint to drive RNA structure prediction, data derived by SHAPE-MaP with 2A3 yields more accurate predictions than NAI-derived data. Due to its extreme efficiency and accuracy, we can anticipate that 2A3 will rapidly take over conventional SHAPE reagents for probing RNA structures both in vitro and in vivo.

tRNA-modifying enzyme mutations induce codon-specific mistranslation and protein aggregation in yeast

Protein synthesis rate and accuracy are tightly controlled by the cell and are essential for proteome homoeostasis (proteostasis); however, the full picture of how mRNA translational factors maintain protein synthesis accuracy and co-translational protein folding are far from being fully understood. To address this question, we evaluated the role of 70 yeast tRNA-modifying enzyme genes on protein aggregation and used mass spectrometry to identify the aggregated proteins. We show that modification of uridine at anticodon position 34 (U34) by the tRNA-modifying enzymes Elp1, Elp3, Sml3 and Trm9 is critical for proteostasis, the mitochondrial tRNA-modifying enzyme Slm3 plays a fundamental role in general proteostasis and that stress response proteins whose genes are enriched in codons decoded by tRNAs lacking mcm5U34, mcm5s2U34, ncm5U34, ncm5Um34, modifications are overrepresented in protein aggregates of the ELP1, SLM3 and TRM9 KO strains. Increased rates of amino acid misincorporation were also detected in these strains at protein sites that specifically mapped to the codons sites that are decoded by the hypomodified tRNAs, demonstrating that U34 tRNA modifications safeguard the proteome from translational errors, protein misfolding and proteotoxic stress.

Deciphering Epitranscriptome: Modification of mRNA Bases Provides a New Perspective for Post-transcriptional Regulation of Gene Expression

Gene regulation depends on dynamic and reversibly modifiable biological and chemical information in the epigenome/epitranscriptome. Accumulating evidence suggests that messenger RNAs (mRNAs) are generated in flashing bursts in the cells in a precisely regulated manner. However, the different aspects of the underlying mechanisms are not fully understood. Cellular RNAs are post-transcriptionally modified at the base level, which alters the metabolism of mRNA. The current understanding of epitranscriptome in the animal system is far ahead of that in plants. The accumulating evidence indicates that the epitranscriptomic changes play vital roles in developmental processes and stress responses. Besides being non-genetically encoded, they can be of reversible nature and involved in fine-tuning the expression of gene. However, different aspects of base modifications in mRNAs are far from adequate to assign the molecular basis/functions to the epitranscriptomic changes. Advances in the chemogenetic RNA-labeling and high-throughput next-generation sequencing techniques are enabling functional analysis of the epitranscriptomic modifications to reveal their roles in mRNA biology. Mapping of the common mRNA modifications, including N 6-methyladenosine (m6A), and 5-methylcytidine (m5C), have enabled the identification of other types of modifications, such as N 1-methyladenosine. Methylation of bases in a transcript dynamically regulates the processing, cellular export, translation, and stability of the mRNA; thereby influence the important biological and physiological processes. Here, we summarize the findings in the field of mRNA base modifications with special emphasis on m6A, m5C, and their roles in growth, development, and stress tolerance, which provide a new perspective for the regulation of gene expression through post-transcriptional modification. This review also addresses some of the scientific and technical issues in epitranscriptomic study, put forward the viewpoints to resolve the issues, and discusses the future perspectives of the research in this area.

Conservation of Archaeal C/D Box sRNA-Guided RNA Modifications

Post-transcriptional modifications fulfill many important roles during ribosomal RNA maturation in all three domains of life. Ribose 2'-O-methylations constitute the most abundant chemical rRNA modification and are, for example, involved in RNA folding and stabilization. In archaea, these modification sites are determined by variable sets of C/D box sRNAs that guide the activity of the rRNA 2'-O-methyltransferase fibrillarin. Each C/D box sRNA contains two guide sequences that can act in coordination to bridge rRNA sequences. Here, we will review the landscape of archaeal C/D box sRNA genes and their target sites. One focus is placed on the apparent accelerated evolution of guide sequences and the varied pairing of the two individual guides, which results in different rRNA modification patterns and RNA chaperone activities.

Alteration of mRNA 5-Methylcytosine Modification in Neurons After OGD/R and Potential Roles in Cell Stress Response and Apoptosis

Epigenetic modifications play an important role in central nervous system disorders. As a widespread posttranscriptional RNA modification, the role of the m5C modification in cerebral ischemia-reperfusion injury (IRI) remains poorly defined. Here, we successfully constructed a neuronal oxygen-glucose deprivation/reoxygenation (OGD/R) model and obtained an overview of the transcriptome-wide m5C profiles using RNA-BS-seq. We discovered that the distribution of neuronal m5C modifications was highly conserved, significantly enriched in CG-rich regions and concentrated in the mRNA translation initiation regions. After OGD/R, modification level of m5C increased, whereas the number of methylated mRNA genes decreased. The amount of overlap of m5C sites with the binding sites of most RNA-binding proteins increased significantly, except for that of the RBM3-binding protein. Moreover, hypermethylated genes in neurons were significantly enriched in pathological processes, and the hub hypermethylated genes RPL8 and RPS9 identified by the protein-protein interaction network were significantly related to cerebral injury. Furthermore, the upregulated transcripts with hypermethylated modification were enriched in the processes involved in response to stress and regulation of apoptosis, and these processes were not identified in hypomethylated transcripts. In final, we verified that OGD/R induced neuronal apoptosis in vitro using TUNEL and western blot assays. Our study identified novel m5C mRNAs associated with ischemia-reperfusion in neurons, providing valuable perspectives for future studies on the role of the RNA methylation in cerebral IRI.

Posttranscriptional modifications at the 37th position in the anticodon stem-loop of tRNA: structural insights from MD simulations

Transfer RNA (tRNA) is the most diversely modified RNA. Although the strictly conserved purine position 37 in the anticodon stem-loop undergoes modifications that are phylogenetically distributed, we do not yet fully understand the roles of these modifications. Therefore, molecular dynamics simulations are used to provide molecular-level details for how such modifications impact the structure and function of tRNA. A focus is placed on three hypermodified base families that include the parent i6A, t6A, and yW modifications, as well as derivatives. Our data reveal that the hypermodifications exhibit significant conformational flexibility in tRNA, which can be modulated by additional chemical functionalization. Although the overall structure of the tRNA anticodon stem remains intact regardless of the modification considered, the anticodon loop must rearrange to accommodate the bulky, dynamic hypermodifications, which includes changes in the nucleotide glycosidic and backbone conformations, and enhanced or completely new nucleobase-nucleobase interactions compared to unmodified tRNA or tRNA containing smaller (m1G) modifications at the 37th position. Importantly, the extent of the changes in the anticodon loop is influenced by the addition of small functional groups to parent modifications, implying each substituent can further fine-tune tRNA structure. Although the dominant conformation of the ASL is achieved in different ways for each modification, the molecular features of all modified tRNA drive the ASL domain to adopt the functional open-loop conformation. Importantly, the impact of the hypermodifications is preserved in different sequence contexts. These findings highlight the likely role of regulating mRNA structure and translation.

Structural basis for the substrate specificity and catalytic features of pseudouridine kinase from Arabidopsis thaliana

RNA modifications can regulate the stability of RNAs, mRNA-protein interactions, and translation efficiency. Pseudouridine is a prevalent RNA modification, and its metabolic fate after RNA turnover was recently characterized in eukaryotes, in the plant Arabidopsis thaliana. Here, we present structural and biochemical analyses of PSEUDOURIDINE KINASE from Arabidopsis (AtPUKI), the enzyme catalyzing the first step in pseudouridine degradation. AtPUKI, a member of the PfkB family of carbohydrate kinases, is a homodimeric α/β protein with a protruding small β-strand domain, which serves simultaneously as dimerization interface and dynamic substrate specificity determinant. AtPUKI has a unique nucleoside binding site specifying the binding of pseudourine, in particular at the nucleobase, by multiple hydrophilic interactions, of which one is mediated by a loop from the small β-strand domain of the adjacent monomer. Conformational transition of the dimerized small β-strand domains containing active site residues is required for substrate specificity. These dynamic features explain the higher catalytic efficiency for pseudouridine over uridine. Both substrates bind well (similar Km), but only pseudouridine is turned over efficiently. Our studies provide an example for structural and functional divergence in the PfkB family and highlight how AtPUKI avoids futile uridine phosphorylation which in vivo would disturb pyrimidine homeostasis.

Plasma Metabolomics Reveals Metabolic Profiling For Diabetic Retinopathy and Disease Progression

BACKGROUNDS

Diabetic retinopathy (DR), the main retinal vascular complication of DM, is the leading cause of visual impairment and blindness among working-age people worldwide. The aim of this study was to investigate the difference of plasma metabolic profiles in patients with DR to better understand the mechanism of this disease and disease progression.

METHODS

We used ultrahigh-performance liquid Q-Exactive mass spectrometry and multivariate statistical analyses to conduct a comprehensive analysis of plasma metabolites in a population with DR and proliferative DR (PDR). A risk score based on the level of the selected metabolite was established and evaluated using the least absolute shrinkage and selection operator regularization logistic regression (LASSO-LR) based machine learning model.

RESULTS

22 differentially expressed metabolites which belonged to different metabolic pathway were identified and confirmed to be associated with the occurrence of DR. A risk score based on the level of the selected metabolite pseudouridine was established and evaluated to strongly associated with the occurrence of DR. Four circulating plasma metabolites (pseudouridine, glutamate, leucylleucine and N-acetyltryptophan) were identified to be differentially expressed between patients with PDR and other patients, and a risk score formula based on these plasma metabolites was developed and assessed to be significantly related to PDR.

CONCLUSIONS

Our work highlights the possible use of the risk score assessment based on the plasma metabolites not only reveal in the early diagnosis of DR and PDR but also assist in enhancing current therapeutic strategies in the clinic.

MODOMICS: An Operational Guide to the Use of the RNA Modification Pathways Database

MODOMICS is an established database of RNA modifications that provides comprehensive information concerning chemical structures of modified ribonucleosides, their biosynthetic pathways, the location of modified residues in RNA sequences, and RNA-modifying enzymes. This chapter covers the resources available on MODOMICS web server and the basic steps that can be undertaken by the user to explore them. MODOMICS is available at http://www.genesilico.pl/modomics .

Pseudouridines of tRNA Anticodon Stem-Loop Have Unexpected Role in Mutagenesis in Pseudomonas sp

Pseudouridines are known to be important for optimal translation. In this study we demonstrate an unexpected link between pseudouridylation of tRNA and mutation frequency in Pseudomonas species. We observed that the lack of pseudouridylation activity of pseudouridine synthases TruA or RluA elevates the mutation frequency in Pseudomonas putida 3 to 5-fold. The absence of TruA but not RluA elevates mutation frequency also in Pseudomonas aeruginosa. Based on the results of genetic studies and analysis of proteome data, the mutagenic effect of the pseudouridylation deficiency cannot be ascribed to the involvement of error-prone DNA polymerases or malfunctioning of DNA repair pathways. In addition, although the deficiency in TruA-dependent pseudouridylation made P. putida cells more sensitive to antimicrobial compounds that may cause oxidative stress and DNA damage, cultivation of bacteria in the presence of reactive oxygen species (ROS)-scavenging compounds did not eliminate the mutator phenotype. Thus, the elevated mutation frequency in the absence of tRNA pseudouridylation could be the result of a more specific response or, alternatively, of a cumulative effect of several small effects disturbing distinct cellular functions, which remain undetected when studied independently. This work suggests that pseudouridines link the translation machinery to mutation frequency.

Comparative patterns of modified nucleotides in individual tRNA species from a mesophilic and two thermophilic archaea

To improve and complete our knowledge of archaeal tRNA modification patterns, we have identified and compared the modification pattern (type and location) in tRNAs of three very different archaeal species, Methanococcus maripaludis (a mesophilic methanogen), Pyrococcus furiosus (a hyperthermophile thermococcale), and Sulfolobus acidocaldarius (an acidophilic thermophilic sulfolobale). Most abundant isoacceptor tRNAs (79 in total) for each of the 20 amino acids were isolated by two-dimensional gel electrophoresis followed by in-gel RNase digestions. The resulting oligonucleotide fragments were separated by nanoLC and their nucleotide content analyzed by mass spectrometry (MS/MS). Analysis of total modified nucleosides obtained from complete digestion of bulk tRNAs was also performed. Distinct base- and/or ribose-methylations, cytidine acetylations, and thiolated pyrimidines were identified, some at new positions in tRNAs. Novel, some tentatively identified, modifications were also found. The least diversified modification landscape is observed in the mesophilic Methanococcus maripaludis and the most complex one in Sulfolobus acidocaldarius Notable observations are the frequent occurrence of ac⁴C nucleotides in thermophilic archaeal tRNAs, the presence of m⁷G at positions 1 and 10 in Pyrococcus furiosus tRNAs, and the use of wyosine derivatives at position 37 of tRNAs, especially those decoding U1- and C1-starting codons. These results complete those already obtained by others with sets of archaeal tRNAs from Methanocaldococcus jannaschii and Haloferax volcanii.

Interrogation of highly structured RNA with multicomponent deoxyribozyme probes at ambient temperatures

Molecular analysis of RNA through hybridization with sequence-specific probes is challenging due to the intrinsic ability of RNA molecules to form stable secondary and tertiary structures. To overcome the energy barrier toward the probe-RNA complex formation, the probes are made of artificial nucleotides, which are more expensive than their natural counterparts and may still be inefficient. Here, we propose the use of a multicomponent probe based on an RNA-cleaving deoxyribozyme for the analysis of highly structured RNA targets. Efficient interrogation of two native RNA from Saccharomyces cerevisiae-a transfer RNA (tRNA) and 18S ribosomal RNA (rRNA)-was achieved at ambient temperature. We achieved detection limits of tRNA down to ∼0.3 nM, which is two orders of magnitude lower than that previously reported for molecular beacon probes. Importantly, no probe annealing to the target was required, with the hybridization assay performed at 37°C. Excess of nonspecific targets did not compromise the performance of the probe, and high interrogation efficiency was maintained by the probes even in complex matrices, such as cell lysate. A linear dynamic range of 0.3-150 nM tRNA was demonstrated. The probe can be adapted for differentiation of a single mismatch in the tRNA-probe complex. Therefore, this study opens a venue toward highly selective, sensitive, robust, and inexpensive assays for the interrogation of biological RNA.

TET2 chemically modifies tRNAs and regulates tRNA fragment levels

The ten-eleven translocation 2 (TET2) protein, which oxidizes 5-methylcytosine in DNA, can also bind RNA; however, the targets and function of TET2-RNA interactions in vivo are not fully understood. Using stringent affinity tags introduced at the Tet2 locus, we purified and sequenced TET2-crosslinked RNAs from mouse embryonic stem cells (mESCs) and found a high enrichment for tRNAs. RNA immunoprecipitation with an antibody against 5-hydroxymethylcytosine (hm5C) recovered tRNAs that overlapped with those bound to TET2 in cells. Mass spectrometry (MS) analyses revealed that TET2 is necessary and sufficient for the deposition of the hm5C modification on tRNA. Tet2 knockout in mESCs affected the levels of several small noncoding RNAs originating from TET2-bound tRNAs that were enriched by hm5C immunoprecipitation. Thus, our results suggest a new function of TET2 in promoting the conversion of 5-methylcytosine to hm5C on tRNA and regulating the processing or stability of different classes of tRNA fragments.

Cryo-EM structure of the highly atypical cytoplasmic ribosome of Euglena gracilis

Ribosomal RNA is the central component of the ribosome, mediating its functional and architectural properties. Here, we report the cryo-EM structure of a highly divergent cytoplasmic ribosome from the single-celled eukaryotic alga Euglena gracilis. The Euglena large ribosomal subunit is distinct in that it contains 14 discrete rRNA fragments that are assembled non-covalently into the canonical ribosome structure. The rRNA is substantially enriched in post-transcriptional modifications that are spread far beyond the catalytic RNA core, contributing to the stabilization of this highly fragmented ribosome species. A unique cluster of five adenosine base methylations is found in an expansion segment adjacent to the protein exit tunnel, such that it is positioned for interaction with the nascent peptide. As well as featuring distinctive rRNA expansion segments, the Euglena ribosome contains four novel ribosomal proteins, localized to the ribosome surface, three of which do not have orthologs in other eukaryotes.

CIGAR-seq, a CRISPR/Cas-based method for unbiased screening of novel mRNA modification regulators

Cellular RNA is decorated with over 170 types of chemical modifications. Many modifications in mRNA, including m6 A and m5 C, have been associated with critical cellular functions under physiological and/or pathological conditions. To understand the biological functions of these modifications, it is vital to identify the regulators that modulate the modification rate. However, a high-throughput method for unbiased screening of these regulators is so far lacking. Here, we report such a method combining pooled CRISPR screen and reporters with RNA modification readout, termed CRISPR integrated gRNA and reporter sequencing (CIGAR-seq). Using CIGAR-seq, we discovered NSUN6 as a novel mRNA m5 C methyltransferase. Subsequent mRNA bisulfite sequencing in HAP1 cells without or with NSUN6 and/or NSUN2 knockout showed that NSUN6 and NSUN2 worked on non-overlapping subsets of mRNA m5 C sites and together contributed to almost all the m5 C modification in mRNA. Finally, using m1 A as an example, we demonstrated that CIGAR-seq can be easily adapted for identifying regulators of other mRNA modification.

Distinct 5-methylcytosine profiles of circular RNA in human hepatocellular carcinoma

Recent studies have indicated that several circular RNAs (circRNAs) can affect the occurrence and development of hepatocellular carcinoma (HCC). Post-transcriptional methylation modifications, including 5-methylcytosine (m5C) modification, are closely related to the tumorigenesis of cancers. However, the map of m5C modification of circRNA in HCC remains to be investigated. In this study, we performed MeRIP-seq to identify m5C sites on circRNA of human HCC tissues and paired adjacent non-tumor tissues. Further, we analyzed the relationship between m5C and HCC. Moreover, we performed a bioinformatics analysis to predict the function of specific methylated transcripts. We found that there was a significant difference in m5C between HCC tissues and paired non-tumor tissues, suggesting potential critical roles of m5C of circRNA in HCC development. In addition, the Gene Ontology (GO) analysis results indicated that the unique distribution pattern of circRNA m5C in HCC was associated with specific metabolism-associated pathways. In conclusion, our findings suggest a possible association between HCC and m5C of circRNA. Additionally, our results provide new insights into a novel function of m5C RNA methylation of circRNA in HCC progression.

Biphasic Liquid Microjunction Extraction for Profiling Neuronal RNA Modifications by Liquid Chromatography-Tandem Mass Spectrometry

RNA modifications are emerging as critical players in the spatiotemporal regulation of gene expression. Although liquid chromatography-tandem mass spectrometry (LC-MS/MS) enables the simultaneous quantification of numerous enzymatically modified RNAs in a biological sample, conventional RNA extraction and enzymatic digestion protocols that are employed prior to analysis have precluded the application of this technique for small-volume samples. In this study, a biphasic liquid microjunction (LMJ) extraction system using coaxial capillaries that direct and aspirate extraction solvents onto a ∼350 μm diameter sample spot was developed and applied for the extraction of RNA from individual cell clusters in the central nervous system of the marine mollusk Aplysia californica. To maximize RNA recoveries, optimized extraction solvents consisting of 10% methanol and chloroform were evaluated under dynamic and static extraction conditions. An MS-compatible RNA digestion buffer was developed to minimize the number of sample-transfer steps and facilitate the direct enzymatic digestion of extracted RNA within the sample collection tube. Compared to RNA extraction using a conventional phenol-chloroform method, the LMJ-based method provided a 3-fold greater coverage of the neuronal epitranscriptome for similar amounts of tissues and also produced mRNA of sufficient purity for reverse transcription polymerase chain reaction amplification. Using this approach, the expression of RNA-modifying enzymes in a given neuronal cell cluster can be characterized and simultaneously correlated with the LC-MS/MS analysis of RNA modifications within the same subset of neurons.

Functional Characterization of a Putative RNA Demethylase ALKBH6 in Arabidopsis Growth and Abiotic Stress Responses

RNA methylation and demethylation, which is mediated by RNA methyltransferases (referred to as “writers”) and demethylases (referred to as “erasers”), respectively, are emerging as a key regulatory process in plant development and stress responses. Although several studies have shown that AlkB homolog (ALKBH) proteins are potential RNA demethylases, the function of most ALKBHs is yet to be determined. The Arabidopsis thaliana genome contains thirteen genes encoding ALKBH proteins, the functions of which are largely unknown. In this study, we characterized the function of a potential eraser protein, ALKBH6 (At4g20350), during seed germination and seedling growth in Arabidopsis under abiotic stresses. The seeds of T-DNA insertion alkbh6 knockdown mutants germinated faster than the wild-type seeds under cold, salt, or abscisic acid (ABA) treatment conditions but not under dehydration stress conditions. Although no differences in seedling and root growth were observed between the alkbh6 mutant and wild-type under normal conditions, the alkbh6 mutant showed a much lower survival rate than the wild-type under salt, drought, or heat stress. Cotyledon greening of the alkbh6 mutants was much higher than that of the wild-type upon ABA application. Moreover, the transcript levels of ABA signaling-related genes, including ABI3 and ABI4, were down-regulated in the alkbh6 mutant compared to wild-type plants. Importantly, the ALKBH6 protein had an ability to bind to both m⁶A-labeled and m⁵C-labeled RNAs. Collectively, these results indicate that the potential eraser ALKBH6 plays important roles in seed germination, seedling growth, and survival of Arabidopsis under abiotic stresses.

tRNA elbow modifications affect the tRNA pseudouridine synthase TruB and the methyltransferase TrmA

tRNAs constitute the most highly modified class of RNA. Every tRNA contains a unique set of modifications, and Ψ55, m⁵U54, and m⁷G46 are frequently found within the elbow of the tRNA structure. Despite the abundance of tRNA modifications, we are only beginning to understand the orchestration of modification enzymes during tRNA maturation. Here, we investigated whether pre-existing modifications impact the binding affinity or catalysis by tRNA elbow modification enzymes. Specifically, we focused on the Escherichia coli enzymes TruB, TrmA, and TrmB which generate Ψ55, m⁵U54, and m⁷G46, respectively. tRNAs containing a single modification were prepared, and the binding and activity preferences of purified E. coli TrmA, TruB, and TrmB were examined in vitro. TruB preferentially binds and modifies unmodified tRNA. TrmA prefers to modify unmodified tRNA, but binds most tightly to tRNA that already contains Ψ55. In contrast, binding and modification by TrmB is insensitive to the tRNA modification status. Our results suggest that TrmA and TruB are likely to act on mostly unmodified tRNA precursors during the early stages of tRNA maturation whereas TrmB presumably acts on later tRNA intermediates that are already partially modified. In conclusion, we uncover the mechanistic basis for the preferred modification order in the E. coli tRNA elbow region.

Naturally occurring modified ribonucleosides

The chemical identity of RNA molecules beyond the four standard ribonucleosides has fascinated scientists since pseudouridine was characterized as the "fifth" ribonucleotide in 1951. Since then, the ever-increasing number and complexity of modified ribonucleosides have been found in viruses and throughout all three domains of life. Such modifications can be as simple as methylations, hydroxylations, or thiolations, complex as ring closures, glycosylations, acylations, or aminoacylations, or unusual as the incorporation of selenium. While initially found in transfer and ribosomal RNAs, modifications also exist in messenger RNAs and noncoding RNAs. Modifications have profound cellular outcomes at various levels, such as altering RNA structure or being essential for cell survival or organism viability. The aberrant presence or absence of RNA modifications can lead to human disease, ranging from cancer to various metabolic and developmental illnesses such as Hoyeraal-Hreidarsson syndrome, Bowen-Conradi syndrome, or Williams-Beuren syndrome. In this review article, we summarize the characterization of all 143 currently known modified ribonucleosides by describing their taxonomic distributions, the enzymes that generate the modifications, and any implications in cellular processes, RNA structure, and disease. We also highlight areas of active research, such as specific RNAs that contain a particular type of modification as well as methodologies used to identify novel RNA modifications. This article is categorized under: RNA Processing > RNA Editing and Modification.

Transcriptome-Wide 5-Methylcytosine Functional Profiling of Long Non-Coding RNA in Hepatocellular Carcinoma

Background

Growing evidence indicates that methylation status is associated with the pathogenesis of numerous types of cancers. Among these, hepatocellular carcinoma (HCC) is a deadly disease threatening global human health. Although 5-methylcytosine (m5C) has been identified as an important regulatory modification, its distribution in solid tumors, including HCC, remains unclear. The present study aimed to explore the distribution of m5C in HCC.

Materials and Methods

Six pairs of human HCC tissues and adjacent non-tumor tissues were collected to analyze the transcriptome-wide m5C methylation of long non-coding RNA (lncRNA). RNA MeRIP-seq was performed to identify m5C peaks on lncRNA and differences in m5C distribution between HCC and adjacent tissues. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment pathway analyses were explored to predict the possible roles of m5C.

Results

Using m5C peak sequencing, we observed that a sequence motif was necessary for m5C methylation in HCC lncRNA. Unsupervised hierarchical cluster analysis confirmed that lncRNA m5C methylation occurred more frequently in HCC than adjacent non-tumor tissues. RNA sequencing data demonstrated that more genes were up-regulated by methylation in HCC, while methylation down-regulated more genes in adjacent non-tumor tissues. GO and KEGG pathway analyses revealed that genes having a significant correlation with m5C sites in lncRNA were involved in HCC signaling pathways.

Conclusion

Our results revealed the substantially different amounts and distributions of m5C in HCC compared to adjacent non-tumor tissue. We further predicted the cellular functions in HCC that m5C may participate in to provide evidence implicating m5C lncRNA epigenetic regulation in the tumorigenesis and progression in HCC.

Overview of distinct 5-methylcytosine profiles of messenger RNA in human hepatocellular carcinoma and paired adjacent non-tumor tissues

BACKGROUND

Post-transcriptional methylation modifications, including 5-methylcytosine (m5C) modification, are closely related to the tumorigenesis of cancers. However, the mRNA profile of m5C modification in hepatocellular carcinoma (HCC) is unknown.

METHODS

Methylated RNA immunoprecipitation sequencing was performed to identify m5C peaks on mRNA of human HCC tissues and adjacent tissues, and differences in m5C between the two groups were analyzed. In addition, we conducted a bioinformatics analysis to predict the function of specific methylated transcripts.

RESULTS

We found that there was a noticeable difference in m5C between HCC and paired non-tumor tissues, suggesting that m5C could play a role in the pathogenesis of HCC. In addition, analyses of gene ontology and the Kyoto Encyclopedia of Genes and Genomes showed that the unique distribution pattern of mRNA m5C in HCC was associated with a wide range of cellular functions.

CONCLUSIONS

Our results revealed different distribution patterns of m5C in HCC and adjacent tissues and provided new insights into a novel function of m5C RNA methylation of mRNA in HCC progression.

Analysis of mitochondrial m1A/G RNA modification reveals links to nuclear genetic variants and associated disease processes

RNA modifications affect the stability and function of RNA species, regulating important downstream processes. Modification levels are often dynamic, varying between tissues and individuals, although it is not always clear what modulates this or what impact it has on biological systems. Here, we quantify variation in m1A/G RNA modification levels at functionally important positions in the human mitochondrial genome across 11,552 samples from 39 tissue/cell types and find that modification levels are associated with mitochondrial transcript processing. We identify links between mitochondrial RNA modification levels and genetic variants in the nuclear genome, including a missense mutation in LONP1, and find that genetic variants within MRPP3 and TRMT61B are associated with RNA modification levels across a large number of tissues. Genetic variants linked to RNA modification levels are associated with multiple disease/disease-related phenotypes, including blood pressure, breast cancer and psoriasis, suggesting a role for mitochondrial RNA modification in complex disease.

Ali et al. analyze publicly available RNA-seq data from different tissues to quantify variation in m1A/G methylation levels in mitochondrial RNAs. They show a link between mitochondrial m1A/G modification levels and nuclear genetic variants, many of which are associated with disease.

A Mini-review of the Computational Methods Used in Identifying RNA 5-Methylcytosine Sites

RNA 5-methylcytosine (m⁵C) is one of the pillars of post-transcriptional modification (PTCM). A growing body of evidence suggests that m⁵C plays a vital role in RNA metabolism. Accurate localization of RNA m⁵C sites in tissue cells is the premise and basis for the in-depth understanding of the functions of m⁵C. However, the main experimental methods of detecting m⁵C sites are limited to varying degrees. Establishing a computational model to predict modification sites is an excellent complement to wet experiments for identifying m⁵C sites. In this review, we summarized some available m⁵C predictors and discussed the characteristics of these methods.

The guide sRNA sequence determines the activity level of box C/D RNPs

2'-O-rRNA methylation, which is essential in eukaryotes and archaea, is catalysed by the Box C/D RNP complex in an RNA-guided manner. Despite the conservation of the methylation sites, the abundance of site-specific modifications shows variability across species and tissues, suggesting that rRNA methylation may provide a means of controlling gene expression. As all Box C/D RNPs are thought to adopt a similar structure, it remains unclear how the methylation efficiency is regulated. Here, we provide the first structural evidence that, in the context of the Box C/D RNP, the affinity of the catalytic module fibrillarin for the substrate-guide helix is dependent on the RNA sequence outside the methylation site, thus providing a mechanism by which both the substrate and guide RNA sequences determine the degree of methylation. To reach this result, we develop an iterative structure-calculation protocol that exploits the power of integrative structural biology to characterize conformational ensembles.

Identification of the 3-amino-3-carboxypropyl (acp) transferase enzyme responsible for acp3U formation at position 47 in Escherichia coli tRNAs

tRNAs from all domains of life contain modified nucleotides. However, even for the experimentally most thoroughly characterized model organism Escherichia coli not all tRNA modification enzymes are known. In particular, no enzyme has been found yet for introducing the acp³U modification at position 47 in the variable loop of eight E. coli tRNAs. Here we identify the so far functionally uncharacterized YfiP protein as the SAM-dependent 3-amino-3-carboxypropyl transferase catalyzing this modification and thereby extend the list of known tRNA modification enzymes in E. coli. Similar to the Tsr3 enzymes that introduce acp modifications at U or m¹Ψ nucleotides in rRNAs this protein contains a DTW domain suggesting that acp transfer reactions to RNA nucleotides are a general function of DTW domain containing proteins. The introduction of the acp³U-47 modification in E. coli tRNAs is promoted by the presence of the m⁷G-46 modification as well as by growth in rich medium. However, a deletion of the enzymes responsible for the modifications at position 46 and 47 in the variable loop of E. coli tRNAs did not lead to a clearly discernible phenotype suggesting that these two modifications play only a minor role in ensuring the proper function of tRNAs in E. coli.

Ribonucleic Acid Sequence Characterization by Negative Electron Transfer Dissociation Mass Spectrometry

Modified oligonucleotides represent a promising avenue for drug development, with small interfering RNAs (siRNA) and microRNAs gaining traction in the therapeutic market. Mass spectrometry (MS)-based analysis offers many benefits for characterizing modified nucleic acids. Negative electron transfer dissociation (NETD) has proven valuable in sequencing oligonucleotide anions, particularly because it can retain modifications while generating sequence-informative fragments. We show that NETD can be successfully implemented on a widely available quadrupole-Orbitrap-linear ion trap mass spectrometer that uses a front-end glow discharge source to generate radical fluoranthene reagent cations. We characterize both unmodified and modified ribonucleic acids and present the first application of activated-ion negative electron transfer dissociation (AI-NETD) to nucleic acids. AI-NETD achieved 100% sequence coverage for both a 6-mer (5'-rGmUrArCmUrG-3') with 2'-O-methyl modifications and a 21-mer (5'-rCrArUrCrCrUrCrUrArGrArGrGrArUrArGrArArUrG-3'), the luciferase antisense siRNA. Both NETD and AI-NETD afforded complete sequence coverage of these molecules while maintaining a relatively low degree of undesired base-loss products and internal products relative to collision-based methods.

Ribosomal RNA Methyltransferase RsmC Moonlights as an RNA Chaperone

Ribosomes are ribonucleoprotein particles that are essential for protein biosynthesis in all forms of life. During ribosome biogenesis, transcription, folding, modification, and processing of rRNA are coupled to the assembly of proteins. Various assembly factors are required to synchronize all different processes that occur during ribosome biogenesis. Herein, the RNA chaperone and RNA strand annealing activity of rRNA modification enzyme ribosome small subunit methyltransferase C (RsmC), which modifies guanine to 2-methylguanosine (m² G) at position 1207 of 16S rRNA (Escherichia coli nucleotide numbering) located at helix 34 (h34), are reported. A 25-fold increase in the h34 RNA strand annealing rates is observed in the presence of RsmC. Single-molecule FRET experiments confirmed the ability of protein RsmC to denature a non-native structure formed by one of the two h34 strands and to form a native-like duplex. This observed RNA chaperone activity of protein RsmC might play a vital role in the rapid generation of functional ribosomes.

mRNA modification orchestrates cancer stem cell fate decisions

Despite their small numbers, cancer stem cells play a central role in driving cancer cell growth, chemotherapeutic resistance, and distal metastasis. Previous studies mainly focused on how DNA or histone modification determines cell fate in cancer. However, it is still largely unknown how RNA modifications orchestrate cancer cell fate decisions. More than 170 distinct RNA modifications have been identified in the RNA world, while only a few RNA base modifications have been found in mRNA. Growing evidence indicates that three mRNA modifications, inosine, 5-methylcytosine, and N⁶-methyladenosine, are essential for the regulation of spatiotemporal gene expression during cancer stem cell fate transition. Furthermore, transcriptome-wide mapping has found that the aberrant deposition of mRNA modification, which can disrupt the gene regulatory network and lead to uncontrollable cancer cell growth, is widespread across different cancers. In this review, we try to summarize the recent advances of these three mRNA modifications in maintaining the stemness of cancer stem cells and discuss the underlying molecular mechanisms, which will shed light on the development of novel therapeutic approaches for eradicating cancer stem cells.

RF-PseU: A Random Forest Predictor for RNA Pseudouridine Sites

One of the ubiquitous chemical modifications in RNA, pseudouridine modification is crucial for various cellular biological and physiological processes. To gain more insight into the functional mechanisms involved, it is of fundamental importance to precisely identify pseudouridine sites in RNA. Several useful machine learning approaches have become available recently, with the increasing progress of next-generation sequencing technology; however, existing methods cannot predict sites with high accuracy. Thus, a more accurate predictor is required. In this study, a random forest-based predictor named RF-PseU is proposed for prediction of pseudouridylation sites. To optimize feature representation and obtain a better model, the light gradient boosting machine algorithm and incremental feature selection strategy were used to select the optimum feature space vector for training the random forest model RF-PseU. Compared with previous state-of-the-art predictors, the results on the same benchmark data sets of three species demonstrate that RF-PseU performs better overall. The integrated average leave-one-out cross-validation and independent testing accuracy scores were 71.4% and 74.7%, respectively, representing increments of 3.63% and 4.77% versus the best existing predictor. Moreover, the final RF-PseU model for prediction was built on leave-one-out cross-validation and provides a reliable and robust tool for identifying pseudouridine sites. A web server with a user-friendly interface is accessible at http://148.70.81.170:10228/rfpseu.

Transcriptome Analysis of ppdnmt2 and Identification of Superoxide Dismutase as a Novel Interactor of DNMT2 in the Moss Physcomitrella patens

DNMT2 is a DNA/tRNA cytosine methyltransferase that is highly conserved in structure and function in eukaryotes. In plants however, limited information is available on the function of this methyltransferase. We have previously reported that in the moss Physcomitrella patens, DNMT2 plays a crucial role in stress recovery and tRNA^Asp transcription/stability under salt stress. To further investigate the role of PpDNMT2 at genome level, in this study we have performed RNA sequencing of ppdnmt2. Transcriptome analysis reveals a number of genes and pathways to function differentially and suggests a close link between PpDNMT2 function and osmotic and ionic stress tolerance. We propose PpDNMT2 to play a pivotal role in regulating salt tolerance by affecting molecular networks involved in stress perception and signal transduction that underlie maintenance of ion homeostasis in cells. We also examined interactome of PpDNMT2 using affinity purification (AP) coupled to mass spectrometry (AP-MS). Quantitative proteomic analysis reveals several chloroplast proteins involved in light reactions and carbon assimilation and proteins involved in stress response and some not implicated in stress to co-immunoprecipitate with PpDNMT2. Comparison between transcriptome and interactome datasets has revealed novel association between PpDNMT2 activity and the antioxidant enzyme Superoxide dismutase (SOD), protein turnover mediated by the Ubiquitin-proteasome system and epigenetic gene regulation. PpDNMT2 possibly exists in complex with CuZn-SODs in vivo and the two proteins also directly interact in the yeast nucleus as observed by yeast two-hybrid assay. Taken together, the work presented in this study sheds light on diverse roles of PpDNMT2 in maintaining molecular and physiological homeostasis in P. patens. This is a first report describing transcriptome and interactome of DNMT2 in any land plant.

Targeting the Ribosome Biogenesis Key Molecule Fibrillarin to Avoid Chemoresistance

Background:

Inherent or acquired chemo resistance in cancer patients has been a perpetual limitation in cancer treatment. Expanding knowledge on essential cellular processes opens a new window for therapeutic targeting. Ribosome biogenesis is a process that shows potential due to its fundamental role in cell development and contribution to tumorigenesis as a result of its upregulation. Inhibiting components of ribosome biogenesis has been explored and has shown interesting results. Yet, an important key component, methyltransferase Fibrillarin (FBL), which influences both the abundance and composition of ribosomes, has not been exploited thus far.

Methods:

In this literature review, we describe relevant aspects of ribosome biogenesis in cancer to emphasize the potential of FBL as a therapeutic target, in order to lower the genotoxic effects of anti-cancer treatment.

Results:

Remarkably, the amplification of the 19q13 cytogenetic band, including the gene coding for FBL, correlated to cell viability and resistance in pancreatic cells as well as to a trend toward a shorter survival in pancreatic cancer patients.

:

Targeting ribosome biogenesis, more specifically compared to the secondary effects of chemotherapeutics such as 5-fluorouracil or oxaliplatin, has been achieved by compound CX-5461. The cell dependent activity of this Pol I inhibitor has been reported in ovarian cancer, melanoma and leukemia models with active or mutated p53 status, presenting a promising mechanism to evade p53 resistance.

Conclusion:

Targeting critical ribosome biogenesis components in order to decrease the genotoxic activity in cancer cell looks promising. Hence, we believe that targeting key protein rRNA methyltransferase FBL shows great potential, due to its pivotal role in ribosome biogenesis, its correlation to an improved survival rate at low expression in breast cancer patients and its association with p53.

RNA modifications and the link to human disease

Ribonucleic acid (RNA) is involved in translation and transcription, which are the mechanisms in which cells express genes (Alberts et al., 2002). The three classes of RNA discussed are transfer RNA (tRNA), messenger RNA (mRNA), and ribosomal RNA (rRNA). mRNA is the transcript encoded from DNA, rRNA is associated with ribosomes, and tRNA is associated with amino acids and is used to read mRNA transcripts to make proteins (Lodish, Berk, Zipursky, et al., 2000). Interestingly, the function of tRNA, rRNA, and mRNA can be significantly altered by chemical modifications at the co-transcriptional and post-transcriptional levels, and there are over 171 of these modifications identified thus far (Boccaletto et al., 2018; Modomics-Modified bases, 2017). Several of these modifications are linked to diseases such as cancer, diabetes, and neurological disorders. In this review, we will introduce a few RNA modifications with biological functions and how dysregulation of these RNA modifications is linked to human disease.

Deducing putative ancestral forms of GNRA/receptor interactions from the ribosome

Stable RNAs rely on a vast repertoire of long-range interactions to assist in the folding of complex cellular machineries such as the ribosome. The universally conserved L39/H89 interaction is a long-range GNRA-like/receptor interaction localized in proximity to the peptidyl transferase center of the large subunit of the ribosome. Because of its central location, L39/H89 likely originated at an early evolutionary stage of the ribosome and played a significant role in its early function. However, L39/H89 self-assembly is impaired outside the ribosomal context. Herein, we demonstrate that structural modularity principles can be used to re-engineer L39/H89 to self-assemble in vitro. The new versions of L39/H89 improve affinity and loop selectivity by several orders of magnitude and retain the structural and functional features of their natural counterparts. These versions of L39/H89 are proposed to be ancestral forms of L39/H89 that were capable of assembling and folding independently from proteins and post-transcriptional modifications. This work demonstrates that novel RNA modules can be rationally designed by taking advantage of the modular syntax of RNA. It offers the prospect of creating new biochemical models of the ancestral ribosome and increases the tool kit for RNA nanotechnology and synthetic biology.

Positive-sense RNA viruses reveal the complexity and dynamics of the cellular and viral epitranscriptomes during infection

More than 140 post-transcriptional modifications (PTMs) are known to decorate cellular RNAs, but their incidence, identity and significance in viral RNA are still largely unknown. We have developed an agnostic analytical approach to comprehensively survey PTMs on viral and cellular RNAs. Specifically, we used mass spectrometry to analyze PTMs on total RNA isolated from cells infected with Zika virus, Dengue virus, hepatitis C virus (HCV), poliovirus and human immunodeficiency virus type 1. All five RNA viruses significantly altered global PTM landscapes. Examination of PTM profiles of individual viral genomes isolated by affinity capture revealed a plethora of PTMs on viral RNAs, which far exceeds the handful of well-characterized modifications. Direct comparison of viral epitranscriptomes identified common and virus-specific PTMs. In particular, specific dimethylcytosine modifications were only present in total RNA from virus-infected cells, and in intracellular HCV RNA, and viral RNA from Zika and HCV virions. Moreover, dimethylcytosine abundance during viral infection was modulated by the cellular DEAD-box RNA helicase DDX6. By opening the Pandora's box on viral PTMs, this report presents numerous questions and hypotheses on PTM function and strongly supports PTMs as a new tier of regulation by which RNA viruses subvert the host and evade cellular surveillance systems.

Double methylation of tRNA-U54 to 2'-O-methylthymidine (Tm) synergistically decreases immune response by Toll-like receptor 7

Sensing of nucleic acids for molecular discrimination between self and non-self is a challenging task for the innate immune system. RNA acts as a potent stimulus for pattern recognition receptors including in particular human Toll-like receptor 7 (TLR7). Certain RNA modifications limit potentially harmful self-recognition of endogenous RNA. Previous studies had identified the 2′-O-methylation of guanosine 18 (Gm18) within tRNAs as an antagonist of TLR7 leading to an impaired immune response. However, human tRNA^Lys₃ was non-stimulatory despite lacking Gm18. To identify the underlying molecular principle, interferon responses of human peripheral blood mononuclear cells to differentially modified tRNA^Lys₃ were determined. The investigation of synthetic modivariants allowed attributing a significant part of the immunosilencing effect to the 2′-O-methylthymidine (m⁵Um) modification at position 54. The effect was contingent upon the synergistic presence of both methyl groups at positions C5 and 2’O, as shown by the fact that neither Um54 nor m⁵U54 produced any effect alone. Testing permutations of the nucleobase at ribose-methylated position 54 suggested that the extent of silencing and antagonism of the TLR7 response was governed by hydrogen patterns and lipophilic interactions of the nucleobase. The results identify a new immune-modulatory endogenous RNA modification that limits TLR7 activation by RNA.

Structure of HIV-1 RT/dsRNA initiation complex prior to nucleotide incorporation

RT is a key enzyme in the life cycle of HIV, and is targeted by multiple antiviral drugs. Although for most of its function RT binds a dsDNA or RNA–DNA template–primer substrate, initiation of reverse transcription involves binding of dsRNA. The current study presents a structure of an RT/dsRNA complex that has the basic components of a reverse transcription initiation complex (RTIC). The unique structural features help understand the significantly slower rate of nucleotide incorporation by an RTIC compared with a catalytically efficient reverse transcription elongation complex. This complex may help in designing new experiments for understanding the intricate and slow process of reverse transcription initiation.

The initiation phase of HIV reverse transcription has features that are distinct from its elongation phase. The first structure of a reverse transcription initiation complex (RTIC) that trapped the complex after incorporation of one ddCMP nucleotide was published recently [Larsen KP, et al. (2018) Nature 557:118–122]. Here we report a crystal structure of a catalytically active HIV-1 RT/dsRNA complex that mimics the state of the RTIC before the first nucleotide incorporation. The structure reveals that the dsRNA-bound conformation of RT is closer to that of RT bound to a nonnucleoside RT inhibitor (NNRTI) and dsDNA; a hyperextended thumb conformation helps to accommodate the relatively wide dsRNA duplex. The RNA primer 3′ end is positioned 5 Å away from the polymerase site; however, unlike in an NNRTI-bound state in which structural elements of RT restrict the movement of the primer, the primer terminus of dsRNA is not blocked from reaching the active site of RT. The observed structural changes and energetic cost of bringing the primer 3′ end to the priming site are hypothesized to explain the slower nucleotide incorporation rate of the RTIC. An unusual crystal lattice interaction of dsRNA with its symmetry mate is reminiscent of the RNA architecture within the extended vRNA–tRNA^Lys3 in the RTIC. This RT/dsRNA complex captures the key structural characteristics and components of the RTIC, including the RT conformational changes and interactions with the dsRNA primer-binding site region, and these features have implications for better understanding of RT initiation.

Importance of m N6-methyladenosine (m6A) RNA modification in cancer

RNA methylation, which was identified back in 1970s, has gained remarkable interest in recent years as it was shown to be a reversible modification involved in many cellular processes like mRNA and miRNA processing, mRNA localisation, translation suppression, or activation. These, in turn, affect important bioprocesses such as tissue development, sex determination, and DNA damage response. Important group of proteins are responsible for adding, recognizing, and removing the methyl group to and from the RNA molecules, which are referred as writers, readers, and erasers, respectively. If any of the processes is not strictly controlled, this can cause abnormalities in gene expression, which result in diseases including cancers such as lung, pancreas, glioblastoma, and breast cancer. Mechanisms of RNA methylation and its role in various cancer types and diagnostic methods for RNA methylation are discussed in this article.

The Elongator subunit Elp3 is a non-canonical tRNA acetyltransferase

The Elongator complex catalyzes posttranscriptional tRNA modifications by attaching carboxy-methyl (cm⁵) moieties to uridine bases located in the wobble position. The catalytic subunit Elp3 is highly conserved and harbors two individual subdomains, a radical S-adenosyl methionine (rSAM) and a lysine acetyltransferase (KAT) domain. The details of its modification reaction cycle and particularly the substrate specificity of its KAT domain remain elusive. Here, we present the co-crystal structure of bacterial Elp3 (DmcElp3) bound to an acetyl-CoA analog and compare it to the structure of a monomeric archaeal Elp3 from Methanocaldococcus infernus (MinElp3). Furthermore, we identify crucial active site residues, confirm the importance of the extended N-terminus for substrate recognition and uncover the specific induction of acetyl-CoA hydrolysis by different tRNA species. In summary, our results establish the clinically relevant Elongator subunit as a non-canonical acetyltransferase and genuine tRNA modification enzyme.

Elp3 is the catalytic subunit of the eukaryotic Elongator complex that catalyzes posttranscriptional tRNA modifications. Here the authors present the crystal structures of an acetyl-CoA analog bound bacterial Elp3 and a monomeric archaeal Elp3 and show that Elp3 functions as a tRNA modification enzyme in all domains of life.

Evidence for Hidden Involvement of N3-Protonated Guanine in RNA Structure and Function

Charged nucleobases have been found to occur in several known RNA molecules and are considered essential for their structure and function. The mechanism of their involvement is however not yet fully understood. Revelation of the role of N7-protonated guanine, in modulating the geometry and stability of noncanonical base pairs formed through its unprotonated edges [Watson-Crick (WC) and sugar], has triggered the need to evaluate the feasibility of similar roles of other protonated nucleobases [Halder et al., Phys Chem Chem Phys, 2015, 17, 26249]. In this context, N3 protonation of guanine makes an interesting case as its influence on the charge distribution of the WC edge is similar to that of N7 protonation, though its thermodynamic cost of protonation is significantly higher. In this work, we have carried out structural bioinformatics analyses and quantum mechanics-based calculations to show that N3 protonation of guanine may take place in a cellular environment, at least in the G:C W:W Trans and G:G W:H Cis base pairs. Our results provide a reasonable starting point for future investigations in order to address the larger mechanistic question.

Prediction of intracellular targets of a small compound by analyzing peptides presented on MHC class I

In chemical biology, the elucidation of chemical target is crucial for successful drug development. Because MHC class I molecules present peptides from intracellular damaged proteins, it might be possible to identify targets of a chemical by analyzing peptide sequences on MHC class I. Therefore, we treated cells with the autophagy-inducing chemical TMD-457 and identified the peptides presented on MHC class I. Many of the peptides were derived from molecules involved in ER trafficking and ER stress, which were confirmed by morphological and biochemical analyses. Therefore, our results demonstrate that analyzing MHC class I peptides is useful for the detection of chemical targets.

ncRNA Editing: Functional Characterization and Computational Resources

Noncoding RNAs (ncRNAs) have received much attention due to their central role in gene expression and translational regulation as well as due to their involvement in several biological processes and disease development. Small noncoding RNAs (sncRNAs), such as microRNAs and piwiRNAs, have been thoroughly investigated and functionally characterized. Long noncoding RNAs (lncRNAs), known to play an important role in chromatin-interacting transcription regulation, posttranscriptional regulation, cell-to-cell signaling, and protein regulation, are also being investigated to further elucidate their functional roles.Next-generation sequencing (NGS) technologies have greatly aided in characterizing the ncRNAome. Moreover, the coupling of NGS technology together with bioinformatics tools has been essential to the genome-wide detection of RNA modifications in ncRNAs. RNA editing, a common human co-transcriptional and posttranscriptional modification, is a dynamic biological phenomenon able to alter the sequence and the structure of primary transcripts (both coding and noncoding RNAs) during the maturation process, consequently influencing the biogenesis, as well as the function, of ncRNAs. In particular, the dysregulation of the RNA editing machineries have been associated with the onset of human diseases.In this chapter we discuss the potential functions of ncRNA editing and describe the knowledge base and bioinformatics resources available to investigate such phenomenon.

Nucleotide Modification Alters MicroRNA-Dependent Silencing of MicroRNA Switches

mRNA therapeutics hold great promise for the treatment of human diseases. While incorporating naturally occurring modified nucleotides during synthesis has greatly increased their potency and safety, challenges in selective expression have hindered clinical applications. MicroRNA (miRNA)-regulated in vitro-transcribed mRNAs, called miRNA switches, have been used to control the expression of exogenous mRNA in a cell-selective manner. However, the effect of nucleotide modifications on miRNA-dependent silencing has not been examined. Here we show that the incorporation of pseudouridine, N1-methylpseudourdine, or pseudouridine and 5-methylcytidine, which increases translation, tends to decrease the regulation of miRNA switches. Moreover, pseudouridine and 5-methylcytidine modification enables one miRNA target site at the 3′ UTR to be as effective as four target sites. We also demonstrate that the effects of pseudouridine, pseudouridine and 5-methylcytidine, and N1-methylpseudourdine modification are miRNA switch specific and do not correlate with the proportion of modified nucleotides in the miRNA target site. Furthermore, modified miRNA switches containing seed-complementary target sites are poorly regulated by miRNA. We also show that placing the miRNA target site in the 5′ UTR of the miRNA switch abolishes the effect of nucleotide modification on miRNA-dependent silencing. This work provides insights into the influence of nucleotide modifications on miRNA-dependent silencing and informs the design of optimal miRNA switches.

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

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

Identification of a 2'-O-Methyluridine Nucleoside Hydrolase Using the Metagenomic Libraries

Ribose methylation is among the most ubiquitous modifications found in RNA. 2'-O-methyluridine is found in rRNA, snRNA, snoRNA and tRNA of Archaea, Bacteria, and Eukaryota. Moreover, 2'-O-methylribonucleosides are promising starting materials for the production of nucleic acid-based drugs. Despite the countless possibilities of practical use for the metabolic enzymes associated with methylated nucleosides, there are very few reports regarding the metabolic fate and enzymes involved in the metabolism of 2'-O-alkyl nucleosides. The presented work focuses on the cellular degradation of 2'-O-methyluridine. A novel enzyme was found using a screening strategy that employs Escherichia coli uracil auxotroph and the metagenomic libraries. A 2'-O-methyluridine hydrolase (RK9NH) has been identified together with an aldolase (RK9DPA)-forming a part of a probable gene cluster that is involved in the degradation of 2'-O-methylated nucleosides. The RK9NH is functional in E. coli uracil auxotroph and in vitro. The RK9NH nucleoside hydrolase could be engineered to enzymatically produce 2'-O-methylated nucleosides that are of great demand as raw materials for production of nucleic acid-based drugs. Moreover, RK9NH nucleoside hydrolase converts 5-fluorouridine, 5-fluoro-2'-deoxyuridine and 5-fluoro-2'-O-methyluridine into 5-fluorouracil, which suggests it could be employed in cancer therapy.

RNA Modifications: Reversal Mechanisms and Cancer

An emerging molecular understanding of RNA alkylation and its removal is transforming our knowledge of RNA biology and its interplay with cancer chemotherapy responses. DNA modifications are known to perform critical functions depending on the genome template, including gene expression, DNA replication timing, and DNA damage protection, yet current results suggest that the chemical diversity of DNA modifications pales in comparison to those on RNA. More than 150 RNA modifications have been identified to date, and their complete functional implications are still being unveiled. These include intrinsic roles such as proper processing and RNA maturation; emerging evidence has furthermore uncovered RNA modification "readers", seemingly analogous to those identified for histone modifications. These modification recognition factors may regulate mRNA stability, localization, and interaction with translation machinery, affecting gene expression. Not surprisingly, tumors differentially modulate factors involved in expressing these marks, contributing to both tumorigenesis and responses to alkylating chemotherapy. Here we describe the current understanding of RNA modifications and their removal, with a focus primarily on methylation and alkylation as functionally relevant changes to the transcriptome. Intriguingly, some of the same RNA modifications elicited by physiological processes are also produced by alkylating agents, thus blurring the lines between what is a physiological mark and a damage-induced modification. Furthermore, we find that a high level of gene expression of enzymes with RNA dealkylation activity is a sensitive readout for poor survival in four different cancer types, underscoring the likely importance of examining RNA dealkylation mechanisms to cancer biology and for cancer treatment and prognosis.

RNA modifications in structure prediction - Status quo and future challenges

Chemical modifications of RNA nucleotides change their identity and characteristics and thus alter genetic and structural information encoded in the genomic DNA. tRNA and rRNA are probably the most heavily modified genes, and often depend on derivatization or isomerization of their nucleobases in order to correctly fold into their functional structures. Recent RNomics studies, however, report transcriptome wide RNA modification and suggest a more general regulation of structuredness of RNAs by this so called epitranscriptome. Modification seems to require specific substrate structures, which in turn are stabilized or destabilized and thus promote or inhibit refolding events of regulatory RNA structures. In this review, we revisit RNA modifications and the related structures from a computational point of view. We discuss known substrate structures, their properties such as sub-motifs as well as consequences of modifications on base pairing patterns and possible refolding events. Given that efficient RNA structure prediction methods for canonical base pairs have been established several decades ago, we review to what extend these methods allow the inclusion of modified nucleotides to model and study epitranscriptomic effects on RNA structures.