Transcriptome-wide identification of adenosine-to-inosine editing using the ICE-seq method

ICLAMP: a novel technique to explore adenosine deamination via inosine chemical labeling and affinity molecular purification

Recent developments in sequencing and bioinformatics have advanced our understanding of adenosine-to-inosine (A-to-I) RNA editing. Surprisingly, recent analyses have revealed the capability of adenosine deaminase acting on RNA (ADAR) to edit DNA:RNA hybrid strands. However, edited inosines in DNA remain largely unexplored. A precise biochemical method could help uncover these potentially rare DNA editing sites. We explore maleimide as a scaffold for inosine labeling. With fluorophore-conjugated maleimide, we were able to label inosine in RNA or DNA. Moreover, with biotin-conjugated maleimide, we purified RNA and DNA containing inosine. Our novel technique of inosine chemical labeling and affinity molecular purification offers substantial advantages and provides a versatile platform for further discovery of A-to-I editing sites in RNA and DNA.

Epitranscriptomic orchestrations: Unveiling the regulatory paradigm of m6A, A-to-I editing, and m5C in breast cancer via long noncoding RNAs and microRNAs

Breast cancer (BC) poses a persistent global health challenge, particularly in countries with elevated human development indices linked to factors such as increased life expectancy, education, and wealth. Despite therapeutic progress, challenges persist, and the role of epitranscriptomic RNA modifications in BC remains inadequately understood. The epitranscriptome, comprising diverse posttranscriptional modifications on RNA molecules, holds the potential to intricately modulate RNA function and regulation, implicating dysregulation in various diseases, including BC. Noncoding RNAs (ncRNAs), acting as posttranscriptional regulators, influence physiological and pathological processes, including cancer. RNA modifications in long noncoding RNAs (lncRNAs) and microRNAs (miRNAs) add an extra layer to gene expression control. This review delves into recent insights into epitranscriptomic RNA modifications, such as N-6-methyladenosine (m6A), adenine-to-inosine (A-to-I) editing, and 5-methylcytosine (m5C), specifically in the context of lncRNA and miRNAs in BC, highlighting their potential implications in BC development and progression. Understanding this intricate regulatory landscape is vital for deciphering the molecular mechanisms underlying BC and identifying potential therapeutic targets.

High-throughput virtual screening to identify potential small molecule inhibitors of the Zα domain of the adenosine deaminases acting on RNA 1(ADAR1)

Changes in RNA editing are closely associated with diseases such as cancer, viral infections, and autoimmune disorders. Adenosine deaminase (ADAR1), which acts on RNA 1, plays a key role in adenosine to inosine editing and is a potential therapeutic target for these various diseases. The p150 subtype of ADAR1 is the only one that contains a Zα domain that binds to both Z-DNA and Z-RNA. The Zα domain modulates immune responses and may be suitable targets for antiviral therapy and cancer immunotherapy. In this study, we attempted to utilize molecular docking to identify potential inhibitors that bind to the ADAR1 Zα domain. The virtual docking method screened the potential activity of more than 100,000 compounds on the Zα domain of ADAR1 and filtered to obtain the highest scoring results.We identified 71 compounds promising to bind to ADAR1 and confirmed that two of them, lithospermic acid and Regaloside B, interacts with the ADAR1 Zα domain by surface plasmonic resonance technique. The molecular dynamics calculation of the complex of lithospermic acid and ADAR1 also showed that the binding effect of lithospermic acid to ADAR1 was stable.This study provides a new perspective for the search of ADAR1 inhibitors, and further studies on the anti-ADAR11 activity of these compounds have broad prospects.

Transfer RNA modifications and cellular thermotolerance

RNA molecules are modified post-transcriptionally to acquire their diverse functions. Transfer RNA (tRNA) has the widest variety and largest numbers of RNA modifications. tRNA modifications are pivotal for decoding the genetic code and stabilizing the tertiary structure of tRNA molecules. Alternation of tRNA modifications directly modulates the structure and function of tRNAs and regulates gene expression. Notably, thermophilic organisms exhibit characteristic tRNA modifications that are dynamically regulated in response to varying growth temperatures, thereby bolstering fitness in extreme environments. Here, we review the history and latest findings regarding the functions and biogenesis of several tRNA modifications that contribute to the cellular thermotolerance of thermophiles.

Quantitative profiling of pseudouridylation landscape in the human transcriptome

Pseudouridine (Ψ) is an abundant post-transcriptional RNA modification in ncRNA and mRNA. However, stoichiometric measurement of individual Ψ sites in human transcriptome remains unaddressed. Here we develop 'PRAISE', via selective chemical labeling of Ψ by bisulfite to induce nucleotide deletion signature during reverse transcription, to realize quantitative assessment of the Ψ landscape in the human transcriptome. Unlike traditional bisulfite treatment, our approach is based on quaternary base mapping and revealed an ~10% median modification level for 2,209 confident Ψ sites in HEK293T cells. By perturbing pseudouridine synthases, we obtained differential mRNA targets of PUS1, PUS7, TRUB1 and DKC1, with TRUB1 targets showing the highest modification stoichiometry. In addition, we quantified known and new Ψ sites in mitochondrial mRNA catalyzed by PUS1. Collectively, we provide a sensitive and convenient method to measure transcriptome-wide Ψ; we envision this quantitative approach would facilitate emerging efforts to elucidate the function and mechanism of mRNA pseudouridylation.

Evaluation and development of deep neural networks for RNA 5-Methyluridine classifications using autoBioSeqpy

Post-transcriptionally RNA modifications, also known as the epitranscriptome, play crucial roles in the regulation of gene expression during development. Recently, deep learning (DL) has been employed for RNA modification site prediction and has shown promising results. However, due to the lack of relevant studies, it is unclear which DL architecture is best suited for some pyrimidine modifications, such as 5-methyluridine (m⁵U). To fill this knowledge gap, we first performed a comparative evaluation of various commonly used DL models for epigenetic studies with the help of autoBioSeqpy. We identified optimal architectural variations for m⁵U site classification, optimizing the layer depth and neuron width. Second, we used this knowledge to develop Deepm5U, an improved convolutional-recurrent neural network that accurately predicts m⁵U sites from RNA sequences. We successfully applied Deepm5U to transcriptomewide m⁵U profiling data across different sequencing technologies and cell types. Third, we showed that the techniques for interpreting deep neural networks, including LayerUMAP and DeepSHAP, can provide important insights into the internal operation and behavior of models. Overall, we offered practical guidance for the development, benchmark, and analysis of deep learning models when designing new algorithms for RNA modifications.

Concepts and methods for transcriptome-wide prediction of chemical messenger RNA modifications with machine learning

The expanding field of epitranscriptomics might rival the epigenome in the diversity of biological processes impacted. In recent years, the development of new high-throughput experimental and computational techniques has been a key driving force in discovering the properties of RNA modifications. Machine learning applications, such as for classification, clustering or de novo identification, have been critical in these advances. Nonetheless, various challenges remain before the full potential of machine learning for epitranscriptomics can be leveraged. In this review, we provide a comprehensive survey of machine learning methods to detect RNA modifications using diverse input data sources. We describe strategies to train and test machine learning methods and to encode and interpret features that are relevant for epitranscriptomics. Finally, we identify some of the current challenges and open questions about RNA modification analysis, including the ambiguity in predicting RNA modifications in transcript isoforms or in single nucleotides, or the lack of complete ground truth sets to test RNA modifications. We believe this review will inspire and benefit the rapidly developing field of epitranscriptomics in addressing the current limitations through the effective use of machine learning.

Cryo-EM structure of the RADAR supramolecular anti-phage defense complex

RADAR is a two-protein bacterial defense system that was reported to defend against phage by "editing" messenger RNA. Here, we determine cryo-EM structures of the RADAR defense complex, revealing RdrA as a heptameric, two-layered AAA+ ATPase and RdrB as a dodecameric, hollow complex with twelve surface-exposed deaminase active sites. RdrA and RdrB join to form a giant assembly up to 10 MDa, with RdrA docked as a funnel over the RdrB active site. Surprisingly, our structures reveal an RdrB active site that targets mononucleotides. We show that RdrB catalyzes ATP-to-ITP conversion in vitro and induces the massive accumulation of inosine mononucleotides during phage infection in vivo, limiting phage replication. Our results define ATP mononucleotide deamination as a determinant of RADAR immunity and reveal supramolecular assembly of a nucleotide-modifying machine as a mechanism of anti-phage defense.

Recent Advances of RNA m6A Modifications in Cancer Immunoediting and Immunotherapy

Cancer immunotherapy, which modulates immune responses against tumors using immune-checkpoint inhibitors or adoptive cell transfer, has emerged as a novel and promising therapy for tumors. However, only a minority of patients demonstrate durable responses, while the majority of patients are resistant to immunotherapy. The immune system can paradoxically constrain and promote tumor development and progression. This process is referred to as cancer immunoediting. The mechanisms of resistance to immunotherapy seem to be that cancer cells undergo immunoediting to evade recognition and elimination by the immune system. RNA modifications, specifically N6-methyladenosine (m6A) methylation, have emerged as a key regulator of various post-transcriptional gene regulatory processes, such as RNA export, splicing, stability, and degradation, which play unappreciated roles in various physiological and pathological processes, including immune system development and cancer pathogenesis. Therefore, a deeper understanding of the mechanisms by which RNA modifications impact the cancer immunoediting process can provide insight into the mechanisms of resistance to immunotherapies and the strategies that can be used to overcome such resistance. In this chapter, we briefly introduce the background of cancer immunoediting and immunotherapy. We also review and discuss the roles and mechanisms of RNA m6A modifications in fine-tuning the innate and adaptive immune responses, as well as in regulating tumor escape from immunosurveillance. Finally, we summarize the current strategies targeting m6A regulators for cancer immunotherapy.

Quantitative and site-specific detection of inosine modification in RNA by acrylonitrile labeling-mediated elongation stalling

RNA molecules contain diverse modifications that play crucial roles in a wide variety of biological processes. Inosine is one of the most prevalent modifications in RNA and dysregulation of inosine is correlated with many human diseases. Herein, we established an acrylonitrile labeling-mediated elongation stalling (ALES) method for quantitative and site-specific detection of inosine in RNA from biological samples. In ALES method, inosine is selectively cyanoethylated with acrylonitrile to form N1-cyanoethylinosine (ce1I) through a Michael addition reaction. The N1-cyanoethyl group of ce1I compromises the hydrogen bond between ce1I and other nucleobases, leading to the stalling of reverse transcription at original inosine site. This specific property of stalling at inosine site could be evaluated by subsequent real-time quantitative PCR (qPCR). With the proposed ALES method, we found the significantly increased level of inosine at position Chr1:63117284 of Ino80dos RNA of multiple tissues from sleep-deprived mice compared to the control mice. This is the first report on the investigation of inosine modification in sleep-deprived mice, which may open up new direction for deciphering insomnia from RNA modifications. In addition, we found the decreased level of inosine at GluA2 Q/R site (Chr4:157336723) in glioma tissues, indicating the decreased level of inosine at GluA2 Q/R site may serve as potential indicator for the diagnosis of glioma. Taken together, the proposed ALES method is capable of quantitative and site-specific detection of inosine in RNA, which provides a valuable tool to uncover the functions of inosine in human diseases.

N6-Methyladenosine RNA Modification in Normal and Malignant Hematopoiesis

Over 170 nucleotide variants have been discovered in messenger RNAs (mRNAs) and non-coding RNAs so far. However, only a few of them, including N6-methyladenosine (m6A), 5-methylcytidine (m5C), and N1-methyladenosine (m1A), could be mapped in the transcriptome. These RNA modifications appear to be dynamically regulated, with writer, eraser, and reader proteins being identified for each modification. As a result, there is a growing interest in studying their biological impacts on normal bioprocesses and tumorigenesis over the past few years. As the most abundant internal modification in eukaryotic mRNAs, m6A plays a vital role in the post-transcriptional regulation of mRNA fate via regulating almost all aspects of mRNA metabolism, including RNA splicing, nuclear export, RNA stability, and translation. Studies on mRNA m6A modification serve as a great example for exploring other modifications on mRNA. In this chapter, we will review recent advances in the study of biological functions and regulation of mRNA modifications, specifically m6A, in both normal hematopoiesis and malignant hematopoiesis. We will also discuss the potential of targeting mRNA modifications as a treatment for hematopoietic disorders.

RNA modifications in aging-associated cardiovascular diseases

Cardiovascular disease (CVD) is a leading cause of morbidity and mortality worldwide that bears an enormous healthcare burden and aging is a major contributing factor to CVDs. Functional gene expression network during aging is regulated by mRNAs transcriptionally and by non-coding RNAs epi-transcriptionally. RNA modifications alter the stability and function of both mRNAs and non-coding RNAs and are involved in differentiation, development, and diseases. Here we review major chemical RNA modifications on mRNAs and non-coding RNAs, including N6-adenosine methylation, N1-adenosine methylation, 5-methylcytidine, pseudouridylation, 2′ -O-ribose-methylation, and N7-methylguanosine, in the aging process with an emphasis on cardiovascular aging. We also summarize the currently available methods to detect RNA modifications and the bioinformatic tools to study RNA modifications. More importantly, we discussed the specific implication of the RNA modifications on mRNAs and non-coding RNAs in the pathogenesis of aging-associated CVDs, including atherosclerosis, hypertension, coronary heart diseases, congestive heart failure, atrial fibrillation, peripheral artery disease, venous insufficiency, and stroke.

Quantification of synthetic errors during chemical synthesis of DNA and its suppression by non-canonical nucleosides

Substitutions, insertions, and deletions derived from synthetic oligonucleotides are the hurdles for the synthesis of long DNA such as genomes. We quantified these synthetic errors by next-generation sequencing and revealed that the quality of the enzymatically amplified final combined product depends on the conditions of the preceding solid phase chemical synthesis, which generates the initial pre-amplified fragments. Among all possible substitutions, the G-to-A substitution was the most prominently observed substitution followed by G-to-T, C-to-T, T-to-C, and A-to-G substitutions. The observed error rate for G-to-A substitution was influenced by capping conditions, suggesting that the capping step played a major role in the generation of G-to-A substitution. Because substitutions observed in long DNA were derived from the generation of non-canonical nucleosides during chemical synthesis, non-canonical nucleosides resistant to side reactions could be used as error-proof nucleosides. As an example of such error-proof nucleosides, we evaluated 7-deaza-2´-deoxyguanosine and 8-aza-7-deaza-2´-deoxyguanosine and showed 50-fold decrease in the error rate of G-to-A substitution when phenoxyacetic anhydride was used as capping reagents. This result is the first example that improves the quality of synthesized sequences by using non-canonical nucleosides as error-proof nucleosides. Our results would contribute to the development of highly accurate template DNA synthesis technologies.

Epitranscriptomic mapping of RNA modifications at single-nucleotide resolution using rhodamine sequencing (Rho-seq)

The recent development of epitranscriptomics revealed a new fundamental layer of gene expression, but the mapping of most RNA modifications remains technically challenging. Here, we describe our protocol for Rho-Seq, which enables the mapping of dihydrouridine RNA modification at single-nucleotide resolution. Rho-Seq relies on specific rhodamine-labeling of a subset of modified nucleotides that hinders reverse transcription. Although Rho-Seq was initially applied to the detection of dihydrouridine, we show here that it is applicable to other modifications including 7-methylguanosine or 4-thiouridine. For complete details on the use and execution of this protocol, please refer to Finet et al. (2022).

Inosine and its methyl derivatives: Occurrence, biogenesis, and function in RNA

Inosine is one of the most common post-transcriptional modifications. Since its discovery, it has been noted for its ability to contribute to non-Watson-Crick interactions within RNA. Rapidly accumulating evidence points to the widespread generation of inosine through hydrolytic deamination of adenosine to inosine by different classes of adenosine deaminases. Three naturally occurring methyl derivatives of inosine, i.e., 1-methylinosine, 2'-O-methylinosine and 1,2'-O-dimethylinosine are currently reported in RNA modification databases. These modifications are expected to lead to changes in the structure, folding, dynamics, stability and functions of RNA. The importance of the modifications is indicated by the strong conservation of the modifying enzymes across organisms. The structure, binding and catalytic mechanism of the adenosine deaminases have been well-studied, but the underlying mechanism of the catalytic reaction is not very clear yet. Here we extensively review the existing data on the occurrence, biogenesis and functions of inosine and its methyl derivatives in RNA. We also included the structural and thermodynamic aspects of these modifications in our review to provide a detailed and integrated discussion on the consequences of A-to-I editing in RNA and the contribution of different structural and thermodynamic studies in understanding its role in RNA. We also highlight the importance of further studies for a better understanding of the mechanisms of the different classes of deamination reactions. Further investigation of the structural and thermodynamic consequences and functions of these modifications in RNA should provide more useful information about their role in different diseases.

Technical challenges in defining RNA modifications

It is well established that DNA base modifications play a key role in gene regulation during development and in response to environmental stress. This type of epigenetic control of development and environmental responses has been intensively studied over the past few decades. Similar to DNA, various RNA species also undergo modifications that play important roles in, for example, RNA splicing, protein translation, and the avoidance of immune surveillance by host. More than 160 different types of RNA modifications have been identified. In addition to base modifications, RNA modification also involves splicing of pre-mRNAs, leading to as many as tens of transcript isoforms from a single pre-RNA, especially in higher organisms. However, the function, prevalence and distribution of RNA modifications are poorly understood. The lack of a suitable method for the reliable identification of RNA modifications constitutes a significant challenge to studying their functions. This review focuses on the technologies that enable de novo identification of RNA base modifications and the alternatively spliced mRNA transcripts.

High-throughput screening to identify potential inhibitors of the Zα domain of the adenosine deaminase 1 (ADAR1)

Adenosine deaminases acting on RNA 1 (ADAR1) are enzymes involved in editing adenosine to inosine in the dsRNAs of cells associated with cancer development. The p150 isoform of ADAR1 is the only isoform containing the Zα domain that binds to both Z-DNA and Z-RNA. The Zα domain is suggested to modulate the immune response and could be a suitable target for antiviral treatment and cancer immunotherapy. In this study, we aimed to identify potential inhibitors for ADAR1 protein that bind the Zα domain using molecular docking and simulation tools. Virtual docking and molecular dynamics simulation approaches were used to screen the potential activity of 2115 FDA-approved compounds on the Zα domain of ADAR1 and filtered for to obtain the top-scoring hits. The top three compounds with the best XP Gscore—namely alendronate (−7.045), etidronate (−6.923), and zoledronate (−6.77)—were subjected to 50 ns simulations to characterize complex stability and identify the fundamental interactions that contribute to inhibition of the ADAR1 Zα domain. The three compounds were shown to interact with Lys169, Lys170, Asn173, and Tyr177 of the Zα domain-like helical backbone of Z-RNA. The study provides a comprehensive and novel insights of repurposes drugs for the inhibition of ADAR1 function.

Recent technical advances in the study of nucleic acid modifications

Enzyme-mediated chemical modifications of nucleic acids are indispensable regulators of gene expression. Our understanding of the biochemistry and biological significance of these modifications has largely been driven by an ever-evolving landscape of technologies that enable accurate detection, mapping, and manipulation of these marks. Here we provide a summary of recent technical advances in the study of nucleic acid modifications with a focus on techniques that allow accurate detection and mapping of these modifications. For each modification discussed (N⁶-methyladenosine, 5-methylcytidine, inosine, pseudouridine, and N⁴-acetylcytidine), we begin by introducing the "gold standard" technique for its mapping and detection, followed by a discussion of techniques developed to address any shortcomings of the gold standard. By highlighting the commonalities and differences of these techniques, we hope to provide a perspective on the current state of the field and to lay out a guideline for development of future technologies.

General Principles for the Detection of Modified Nucleotides in RNA by Specific Reagents

Epitranscriptomics heavily rely on chemical reagents for the detection, quantification, and localization of modified nucleotides in transcriptomes. Recent years have seen a surge in mapping methods that use innovative and rediscovered organic chemistry in high throughput approaches. While this has brought about a leap of progress in this young field, it has also become clear that the different chemistries feature variegated specificity and selectivity. The associated error rates, e.g., in terms of false positives and false negatives, are in large part inherent to the chemistry employed. This means that even assuming technically perfect execution, the interpretation of mapping results issuing from the application of such chemistries are limited by intrinsic features of chemical reactivity. An important but often ignored fact is that the huge stochiometric excess of unmodified over-modified nucleotides is not inert to any of the reagents employed. Consequently, any reaction aimed at chemical discrimination of modified versus unmodified nucleotides has optimal conditions for selectivity that are ultimately anchored in relative reaction rates, whose ratio imposes intrinsic limits to selectivity. Here chemical reactivities of canonical and modified ribonucleosides are revisited as a basis for an understanding of the limits of selectivity achievable with chemical methods.

The occurrence order and cross-talk of different tRNA modifications

Chemical modifications expand the composition of RNA molecules from four standard nucleosides to over 160 modified nucleosides, which greatly increase the complexity and utility of RNAs. Transfer RNAs (tRNAs) are the most heavily modified cellular RNA molecules and contain the largest variety of modifications. Modification of tRNAs is pivotal for protein synthesis and also precisely regulates the noncanonical functions of tRNAs. Defects in tRNA modifications lead to numerous human diseases. Up to now, more than 100 types of modifications have been found in tRNAs. Intriguingly, some modifications occur widely on all tRNAs, while others only occur on a subgroup of tRNAs or even only a specific tRNA. The modification frequency of each tRNA is approximately 7% to 25%, with 5-20 modification sites present on each tRNA. The occurrence and modulation of tRNA modifications are specifically noticeable as plenty of interplays among different sites and modifications have been discovered. In particular, tRNA modifications are responsive to environmental changes, indicating their dynamic and highly organized nature. In this review, we summarized the known occurrence order, cross-talk, and cooperativity of tRNA modifications.

RNA Modifications in Genomic RNA of Influenza A Virus and the Relationship between RNA Modifications and Viral Infection

Recent studies about the transcriptome-wide presence of RNA modifications have revealed their importance in many cellular functions. Nevertheless, information about RNA modifications in viral RNA is scarce, especially for negative-strand RNA viruses. Here we provide a catalog of RNA modifications including m1A, ac4C, m7G, inosine, and pseudouridine on RNA derived from an influenza A virus infected into A549 cells, as studied by RNA immunoprecipitation followed by deep-sequencing. Possible regions with RNA modifications were found in the negative-strand segments of viral genomic RNA. In addition, our analyses of previously published data revealed that the expression levels of the host factors for RNA modifications were affected by an infection with influenza A virus, and some of the host factors likely have a proviral effect. RNA modification is a novel aspect of host-virus interactions leading to the discovery of previously unrecognized viral pathogenicity mechanisms and has the potential to aid the development of novel antivirals.

RNA Modifications in Pathogenic Bacteria: Impact on Host Adaptation and Virulence

RNA modifications are involved in numerous biological processes and are present in all RNA classes. These modifications can be constitutive or modulated in response to adaptive processes. RNA modifications play multiple functions since they can impact RNA base-pairings, recognition by proteins, decoding, as well as RNA structure and stability. However, their roles in stress, environmental adaptation and during infections caused by pathogenic bacteria have just started to be appreciated. With the development of modern technologies in mass spectrometry and deep sequencing, recent examples of modifications regulating host-pathogen interactions have been demonstrated. They show how RNA modifications can regulate immune responses, antibiotic resistance, expression of virulence genes, and bacterial persistence. Here, we illustrate some of these findings, and highlight the strategies used to characterize RNA modifications, and their potential for new therapeutic applications.

Epitranscriptomics: A New Layer of microRNA Regulation in Cancer

MicroRNAs are pervasive regulators of gene expression at the post-transcriptional level in metazoan, playing key roles in several physiological and pathological processes. Accordingly, these small non-coding RNAs are also involved in cancer development and progression. Furthermore, miRNAs represent valuable diagnostic and prognostic biomarkers in malignancies. In the last twenty years, the role of RNA modifications in fine-tuning gene expressions at several levels has been unraveled. All RNA species may undergo post-transcriptional modifications, collectively referred to as epitranscriptomic modifications, which, in many instances, affect RNA molecule properties. miRNAs are not an exception, in this respect, and they have been shown to undergo several post-transcriptional modifications. In this review, we will summarize the recent findings concerning miRNA epitranscriptomic modifications, focusing on their potential role in cancer development and progression.

RNA modifications in hematopoietic malignancies: A new research frontier

Both protein-coding and noncoding RNAs can be decorated with a wealth of chemical modifications and such modifications coordinately orchestrate gene expression during normal hematopoietic differentiation and development. However, aberrant expression and/or dysfunction of the relevant RNA modification modulators/regulators ("writers", "erasers", and "readers") drive the initiation and progression of hematopoietic malignancies, and targeting these dysregulated modulators holds potent therapeutic potential for the treatment of hematopoietic malignancies. In this review, we summarize current progress in the understanding of the biological functions and underlying mechanisms of RNA modifications in normal and malignant hematopoiesis, with a focus on the N6-methyladenosine (m6A) modification, and discuss the therapeutic potential of targeting RNA modifications for the treatment of hematopoietic malignancies, especially acute myeloid leukemia (AML).

Long Non-Coding RNA Epigenetics

Long noncoding RNAs exceeding a length of 200 nucleotides play an important role in ensuring cell functions and proper organism development by interacting with cellular compounds such as miRNA, mRNA, DNA and proteins. However, there is an additional level of lncRNA regulation, called lncRNA epigenetics, in gene expression control. In this review, we describe the most common modified nucleosides found in lncRNA, 6-methyladenosine, 5-methylcytidine, pseudouridine and inosine. The biosynthetic pathways of these nucleosides modified by the writer, eraser and reader enzymes are important to understanding these processes. The characteristics of the individual methylases, pseudouridine synthases and adenine-inosine editing enzymes and the methods of lncRNA epigenetics for the detection of modified nucleosides, as well as the advantages and disadvantages of these methods, are discussed in detail. The final sections are devoted to the role of modifications in the most abundant lncRNAs and their functions in pathogenic processes.

Detecting the epitranscriptome

RNA modifications and their corresponding epitranscriptomic writer and eraser enzymes regulate gene expression. Altered RNA modification levels, dysregulated writers, and sequence changes that disrupt epitranscriptomic marks have been linked to mitochondrial and neurological diseases, cancer, and multifactorial disorders. The detection of epitranscriptomics marks is challenging, but different next generation sequencing (NGS)-based and mass spectrometry-based approaches have been used to identify and quantitate the levels of individual and groups of RNA modifications. NGS and mass spectrometry-based approaches have been coupled with chemical, antibody or enzymatic methodologies to identify modifications in most RNA species, mapped sequence contexts and demonstrated the dynamics of specific RNA modifications, as well as the collective epitranscriptome. While epitranscriptomic analysis is currently limited to basic research applications, specific approaches for the detection of individual RNA modifications and the epitranscriptome have potential biomarker applications in detecting human conditions and diseases. This article is categorized under: RNA Structure and Dynamics > Influence of RNA Structure in Biological Systems RNA Processing > tRNA Processing RNA in Disease and Development > RNA in Disease.

Detecting and Characterizing A-To-I microRNA Editing in Cancer

Adenosine to inosine (A-to-I) editing consists of an RNA modification where single adenosines along the RNA sequence are converted into inosines. Such a biochemical transformation is catalyzed by enzymes belonging to the family of adenosine deaminases acting on RNA (ADARs) and occurs either co- or post-transcriptionally. The employment of powerful, high-throughput detection methods has recently revealed that A-to-I editing widely occurs in non-coding RNAs, including microRNAs (miRNAs). MiRNAs are a class of small regulatory non-coding RNAs (ncRNAs) acting as translation inhibitors, known to exert relevant roles in controlling cell cycle, proliferation, and cancer development. Indeed, a growing number of recent researches have evidenced the importance of miRNA editing in cancer biology by exploiting various detection and validation methods. Herein, we briefly overview early and currently available A-to-I miRNA editing detection and validation methods and discuss the significance of A-to-I miRNA editing in human cancer.

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.

Instrumental analysis of RNA modifications

Organisms from all domains of life invest a substantial amount of energy for the introduction of RNA modifications into nearly all transcripts studied to date. Instrumental analysis of RNA can focus on the modified residues and reveal the function of these epitranscriptomic marks. Here, we will review recent advances and breakthroughs achieved by NMR spectroscopy, sequencing, and mass spectrometry of the epitranscriptome.

Analysis of RNA Modifications by Second- and Third-Generation Deep Sequencing: 2020 Update

The precise mapping and quantification of the numerous RNA modifications that are present in tRNAs, rRNAs, ncRNAs/miRNAs, and mRNAs remain a major challenge and a top priority of the epitranscriptomics field. After the keystone discoveries of massive m⁶A methylation in mRNAs, dozens of deep sequencing-based methods and protocols were proposed for the analysis of various RNA modifications, allowing us to considerably extend the list of detectable modified residues. Many of the currently used methods rely on the particular reverse transcription signatures left by RNA modifications in cDNA; these signatures may be naturally present or induced by an appropriate enzymatic or chemical treatment. The newest approaches also include labeling at RNA abasic sites that result from the selective removal of RNA modification or the enhanced cleavage of the RNA ribose-phosphate chain (perhaps also protection from cleavage), followed by specific adapter ligation. Classical affinity/immunoprecipitation-based protocols use either antibodies against modified RNA bases or proteins/enzymes, recognizing RNA modifications. In this survey, we review the most recent achievements in this highly dynamic field, including promising attempts to map RNA modifications by the direct single-molecule sequencing of RNA by nanopores.

Discovering A-to-I RNA Editing Through Chemical Methodology "ICE-seq"

RNA editing of adenosines to inosines contributes to a wide range of biological processes by regulating gene expression post-transcriptionally. To understand the effect, accurate mapping of inosines is necessary. The most conventional method to identify an editing site is to compare the cDNA sequence with its corresponding genomic sequence. However, this method has a high false discovery rate because guanosine signals, due to experimental errors or noise in the obtained sequences, contaminate genuine inosine signals detected as guanosine. To ensure high accuracy, we developed the Inosine Chemical Erasing (ICE) method to accurately and biochemically identify inosines in RNA strands utilizing inosine cyanoethylation and reverse transcription-PCR. Furthermore, we applied this technique to next-generation sequencing technology, called ICE-seq, to conduct an unbiased genome-wide screening of A-to-I editing sites in the transcriptome.

iMRM: a platform for simultaneously identifying multiple kinds of RNA modifications

MOTIVATION

RNA modifications play critical roles in a series of cellular and developmental processes. Knowledge about the distributions of RNA modifications in the transcriptomes will provide clues to revealing their functions. Since experimental methods are time consuming and laborious for detecting RNA modifications, computational methods have been proposed for this aim in the past five years. However, there are some drawbacks for both experimental and computational methods in simultaneously identifying modifications occurred on different nucleotides.

RESULTS

To address such a challenge, in this article, we developed a new predictor called iMRM, which is able to simultaneously identify m6A, m5C, m1A, ψ and A-to-I modifications in Homo sapiens, Mus musculus and Saccharomyces cerevisiae. In iMRM, the feature selection technique was used to pick out the optimal features. The results from both 10-fold cross-validation and jackknife test demonstrated that the performance of iMRM is superior to existing methods for identifying RNA modifications.

AVAILABILITY AND IMPLEMENTATION

A user-friendly web server for iMRM was established at http://www.bioml.cn/XG_iRNA/home. The off-line command-line version is available at https://github.com/liukeweiaway/iMRM.

CONTACT

greatchen@ncst.edu.cn.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

circMeta: a unified computational framework for genomic feature annotation and differential expression analysis of circular RNAs

MOTIVATION

Circular RNAs (circRNAs), a class of non-coding RNAs generated from non-canonical back-splicing events, have emerged to play key roles in many biological processes. Though numerous tools have been developed to detect circRNAs from rRNA-depleted RNA-seq data based on back-splicing junction-spanning reads, computational tools to identify critical genomic features regulating circRNA biogenesis are still lacking. In addition, rigorous statistical methods to perform differential expression (DE) analysis of circRNAs remain under-developed.

RESULTS

We present circMeta, a unified computational framework for circRNA analyses. circMeta has three primary functional modules: (i) a pipeline for comprehensive genomic feature annotation related to circRNA biogenesis, including length of introns flanking circularized exons, repetitive elements such as Alu elements and SINEs, competition score for forming circulation and RNA editing in back-splicing flanking introns; (ii) a two-stage DE approach of circRNAs based on circular junction reads to quantitatively compare circRNA levels and (iii) a Bayesian hierarchical model for DE analysis of circRNAs based on the ratio of circular reads to linear reads in back-splicing sites to study spatial and temporal regulation of circRNA production. Both proposed DE methods without and with considering host genes outperform existing methods by obtaining better control of false discovery rate and comparable statistical power. Moreover, the identified DE circRNAs by the proposed two-stage DE approach display potential biological functions in Gene Ontology and circRNA-miRNA-mRNA networks that are not able to be detected using existing mRNA DE methods. Furthermore, top DE circRNAs have been further validated by RT-qPCR using divergent primers spanning back-splicing junctions.

AVAILABILITY AND IMPLEMENTATION

The software circMeta is freely available at https://github.com/lichen-lab/circMeta.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

CRISPR Tools for Physiology and Cell State Changes: Potential of Transcriptional Engineering and Epigenome Editing

Given the large amount of genome-wide data that have been collected during the last decades, a good understanding of how and why cells change during development, homeostasis, and disease might be expected. Unfortunately, the opposite is true; triggers that cause cellular state changes remain elusive, and the underlying molecular mechanisms are poorly understood. Although genes with the potential to influence cell states are known, the historic dependency on methods that manipulate gene expression outside the endogenous chromatin context has prevented us from understanding how cells organize, interpret, and protect cellular programs. Fortunately, recent methodological innovations are now providing options to answer these outstanding questions, by allowing to target and manipulate individual genomic and epigenomic loci. In particular, three experimental approaches are now feasible due to DNA targeting tools, namely, activation and/or repression of master transcription factors in their endogenous chromatin context; targeting transcription factors to endogenous, alternative, or inaccessible sites; and finally, functional manipulation of the chromatin context. In this article, we discuss the molecular basis of DNA targeting tools and review the potential of these new technologies before we summarize how these have already been used for the manipulation of cellular states and hypothesize about future applications.

Role of RNA modifications in cancer

Specific chemical modifications of biological molecules are an efficient way of regulating molecular function, and a plethora of downstream signalling pathways are influenced by the modification of DNA and proteins. Many of the enzymes responsible for regulating protein and DNA modifications are targets of current cancer therapies. RNA epitranscriptomics, the study of RNA modifications, is the new frontier of this arena. Despite being known since the 1970s, eukaryotic RNA modifications were mostly identified on transfer RNA and ribosomal RNA until the last decade, when they have been identified and characterized on mRNA and various non-coding RNAs. Increasing evidence suggests that RNA modification pathways are also misregulated in human cancers and may be ideal targets of cancer therapy. In this Review we highlight the RNA epitranscriptomic pathways implicated in cancer, describing their biological functions and their connections to the disease.

RNA Modifications in Cancer: Functions, Mechanisms, and Therapeutic Implications

Over 170 chemical modifications have been identified in protein-coding and noncoding RNAs and shown to exhibit broad impacts on gene expression. Dysregulation of RNA modifications caused by aberrant expression of or mutations in RNA modifiers aberrantly reprograms the epitranscriptome and skews global gene expression, which in turn leads to tumorigenesis and drug resistance. Here we review current knowledge of the functions and underlying mechanisms of aberrant RNA modifications in human cancers, particularly several common RNA modifications, including N6-methyladenosine (m6A), A-to-I editing, pseudouridine (ψ), 5-methylcytosine (m5C), 5-hydroxymethylcytosine (hm5C), N1-methyladenosine (m1A), and N4-acetylcytidine (ac4C), providing insights into therapeutic implications of targeting RNA modifications and the associated machineries for cancer therapy.

m6A mRNA methylation: A pleiotropic regulator of cancer

In recent days, RNA modifications are gaining the interest of biologist worldwide. Till date, a total of 171 RNA modifications has been reported, and the number may increase with advancing technologies. The mRNA undergoes modifications like m⁵M, hm⁵C, m¹A, m⁶A and pseudouridine, collectively called as epitranscriptomic alterations, each of them has their functional significance. m⁶A modification is the most common one which occurs at the motif of RRm⁶AACH in mRNA. The altered profiles of these epitranscriptomic changes are reported in multiple cancers. The present review discusses the dynamic nature of functional enzymes called methyltransferase (writer), demethylase (erasers) and m⁶A binding proteins (readers) and importance of the balance between these proteins for the homeostasis of our body functions like metabolism, circadian rhythm, immune response, viral replications, embryogenesis and cancer development. Nevertheless, the main focus has been on cancer development and progression. The understanding of such differential modifications are at infancy and may provide bring about a paradigm shift in our understanding of cancer for management and treatment.

Single-nucleotide variants in human RNA: RNA editing and beyond

Through analysis of paired high-throughput DNA-Seq and RNA-Seq data, researchers quickly recognized that RNA-Seq can be used for more than just gene expression quantification. The alternative applications of RNA-Seq data are abundant, and we are particularly interested in its usefulness for detecting single-nucleotide variants, which arise from RNA editing, genomic variants and other RNA modifications. A stunning discovery made from RNA-Seq analyses is the unexpectedly high prevalence of RNA-editing events, many of which cannot be explained by known RNA-editing mechanisms. Over the past 6-7 years, substantial efforts have been made to maximize the potential of RNA-Seq data. In this review we describe the controversial history of mining RNA-editing events from RNA-Seq data and the corresponding development of methodologies to identify, predict, assess the quality of and catalog RNA-editing events as well as genomic variants.

Identification of Adenosine-to-Inosine RNA Editing with Acrylonitrile Reagents

New chemical probes have been designed to facilitate the identification of adenosine-to-inosine (A-to-I) edited RNAs. These reagents combine a conjugate acceptor for selective inosine covalent modification with functional groups for bioorthogonal biotinylation. The resulting biotinylated RNA was enriched and verified with RT-qPCR. This powerful chemical approach provides new opportunities to identify and quantify A-to-I editing sites.

Formation of tRNA Wobble Inosine in Humans Is Disrupted by a Millennia-Old Mutation Causing Intellectual Disability

The formation of inosine at the wobble position of eukaryotic tRNAs is an essential modification catalyzed by the ADAT2/ADAT3 complex. In humans, a valine-to-methionine mutation (V144M) in ADAT3 that originated ∼1,600 years ago is the most common cause of autosomal recessive intellectual disability (ID) in Arabia. While the mutation is predicted to affect protein structure, the molecular and cellular effects of the V144M mutation are unknown.

The formation of inosine at the wobble position of eukaryotic tRNAs is an essential modification catalyzed by the ADAT2/ADAT3 complex. In humans, a valine-to-methionine mutation (V144M) in ADAT3 that originated ∼1,600 years ago is the most common cause of autosomal recessive intellectual disability (ID) in Arabia. While the mutation is predicted to affect protein structure, the molecular and cellular effects of the V144M mutation are unknown. Here, we show that cell lines derived from ID-affected individuals expressing only ADAT3-V144M exhibit decreased wobble inosine in certain tRNAs. Moreover, extracts from the same cell lines of ID-affected individuals display a severe reduction in tRNA deaminase activity. While ADAT3-V144M maintains interactions with ADAT2, the purified ADAT2/3-V144M complexes exhibit defects in activity. Notably, ADAT3-V144M exhibits an increased propensity to form aggregates associated with cytoplasmic chaperonins that can be suppressed by ADAT2 overexpression. These results identify a key role for ADAT2-dependent folding of ADAT3 in wobble inosine modification and indicate that proper formation of an active ADAT2/3 complex is crucial for proper neurodevelopment.

Patterns of Ancestral Animal Codon Usage Bias Revealed through Holozoan Protists

Choanoflagellates and filastereans are the closest known single celled relatives of Metazoa within Holozoa and provide insight into how animals evolved from their unicellular ancestors. Codon usage bias has been extensively studied in metazoans, with both natural selection and mutation pressure playing important roles in different species. The disparate nature of metazoan codon usage patterns prevents the reconstruction of ancestral traits. However, traits conserved across holozoan protists highlight characteristics in the unicellular ancestors of Metazoa. Presented here are the patterns of codon usage in the choanoflagellates Monosiga brevicollis and Salpingoeca rosetta, as well as the filasterean Capsaspora owczarzaki. Codon usage is shown to be remarkably conserved. Highly biased genes preferentially use GC-ending codons, however there is limited evidence this is driven by local mutation pressure. The analyses presented provide strong evidence that natural selection, for both translational accuracy and efficiency, dominates codon usage bias in holozoan protists. In particular, the signature of selection for translational accuracy can be detected even in the most weakly biased genes. Biased codon usage is shown to have coevolved with the tRNA species, with optimal codons showing complementary binding to the highest copy number tRNA genes. Furthermore, tRNA modification is shown to be a common feature for amino acids with higher levels of degeneracy and highly biased genes show a strong preference for using modified tRNAs in translation. The translationally optimal codons defined here will be of benefit to future transgenics work in holozoan protists, as their use should maximise protein yields from edited transgenes.

The juvenility-associated long noncoding RNA Gm14230 maintains cellular juvenescence

Juvenile animals possess distinct properties that are missing in adults. These properties include capabilities for higher growth, faster wound healing, plasticity and regeneration. However, the molecular mechanisms underlying these juvenile physiological properties are not fully understood. To obtain insight into the distinctiveness of juveniles from adults at the molecular level, we assessed long noncoding RNAs (lncRNAs) that are highly expressed selectively in juvenile cells. The noncoding elements of the transcriptome were investigated in hepatocytes and cardiomyocytes isolated from juvenile and adult mice. Here, we identified 62 juvenility-associated lncRNAs (JAlncs), which are selectively expressed in both hepatocytes and cardiomyocytes from juvenile mice. Among these common (shared) JAlncs, Gm14230 is evolutionarily conserved and is essential for cellular juvenescence. Loss of Gm14230 impairs cell growth and causes cellular senescence. Gm14230 safeguards cellular juvenescence through recruiting the histone methyltransferase Ezh2 to Tgif2, thereby repressing the functional role of Tgif2 in cellular senescence. Thus, we identify Gm14230 as a juvenility-selective lncRNA required to maintain cellular juvenescence.

Mapping the dsRNA World

Long double-stranded RNAs (dsRNAs) are abundantly expressed in animals, in which they frequently occur in introns and 3' untranslated regions of mRNAs. Functions of long, cellular dsRNAs are poorly understood, although deficiencies in adenosine deaminases that act on RNA, or ADARs, promote their recognition as viral dsRNA and an aberrant immune response. Diverse dsRNA-binding proteins bind cellular dsRNAs, hinting at additional roles. Understanding these roles is facilitated by mapping the genomic locations that express dsRNA in various tissues and organisms. ADAR editing provides a signature of dsRNA structure in cellular transcripts. In this review, we detail approaches to map ADAR editing sites and dsRNAs genome-wide, with particular focus on high-throughput sequencing methods and considerations for their successful application to the detection of editing sites and dsRNAs.

Epitranscriptomics: RNA Modifications in Bacteria and Archaea

The increasingly complex functionality of RNA is contrasted by its simple chemical composition. RNA is generally built from only four different nucleotides (adenine, guanine, cytosine, and uracil). To date, >160 chemical modifications are known to decorate RNA molecules and thereby alter their function or stability. Many RNA modifications are conserved throughout bacteria, archaea, and eukaryotes, while some are unique to each branch of life. Most known modifications occur at internal positions, while there is limited diversity at the termini. The dynamic nature of RNA modifications and newly discovered regulatory functions of some of these RNA modifications gave birth to a new field, now often referred to as "epitranscriptomics." This review highlights the major developments in this field and summarizes detection principles for internal as well as 5'-terminal mRNA modifications in prokaryotes and archaea to investigate their biological significance.

Methods for RNA Modification Mapping Using Deep Sequencing: Established and New Emerging Technologies

New analytics of post-transcriptional RNA modifications have paved the way for a tremendous upswing of the biological and biomedical research in this field. This especially applies to methods that included RNA-Seq techniques, and which typically result in what is termed global scale modification mapping. In this process, positions inside a cell’s transcriptome are receiving a status of potential modification sites (so called modification calling), typically based on a score of some kind that issues from the particular method applied. The resulting data are thought to represent information that goes beyond what is contained in typical transcriptome data, and hence the field has taken to use the term “epitranscriptome”. Due to the high rate of newly published mapping techniques, a significant number of chemically distinct RNA modifications have become amenable to mapping, albeit with variegated accuracy and precision, depending on the nature of the technique. This review gives a brief overview of known techniques, and how they were applied to modification calling.

Existence of Internal N7-Methylguanosine Modification in mRNA Determined by Differential Enzyme Treatment Coupled with Mass Spectrometry Analysis

The recent discovery of reversible chemical modifications on mRNA has opened a new era of post-transcriptional gene regulation in eukaryotes. Among the 15 types of modifications identified in mRNA of eukaryotes, N7-methylguanosine (m⁷G) is unique owing to its presence in the 5' cap structure. It remains unknown whether m⁷G is also present internally in mRNA, and this is largely attributed to the lack of an appropriate analytical method to differentiate internal m⁷G in mRNA from that in the 5' cap. To address this analytical challenge, we developed a novel strategy of combining differential enzymatic digestion with liquid chromatography-tandem mass spectrometry analysis to quantify the levels of these two types of m⁷G modifications in mRNA. In particular, we found that S1 nuclease and phosphodiesterase I exhibit differential activities toward internal and 5'-terminal m⁷G. By using this method, we found that internal m⁷G was present in mRNA of cultured human cells as well as plants and rat tissue. In addition, our results showed that plants contain higher levels of internal m⁷G in mRNA than mammals. We also observed that exposure of rice to cadmium (Cd) stimulated marked diminution in the levels of m⁷G at both the 5' cap and internal positions of mRNA, which was correlated with the Cd-induced elevated expression of m⁷G-decapping enzymes. Taken together, we reported here a strategy to distinguish internal and 5'-terminal m⁷G in mRNA, and by using this method, we demonstrated the prevalence of internal m⁷G modification in mRNA, which we believe will stimulate future functional studies of m⁷G on post-transcriptional gene regulation in eukaryotes.

Transcriptome-wide identification of A-to-I RNA editing sites using ICE-seq

In A-to-I RNA editing, adenosine is converted to inosine in double-stranded regions of RNAs. Inosine, an abundant epitranscriptomic mark, contributes to a wide range of biological processes by regulating gene expression post-transcriptionally. To understand the effect of A-to-I RNA editing on regulation of the epitranscriptome, accurate mapping of inosines is necessary. To this end, we established a biochemical method called inosine chemical erasing sequencing (ICE-seq) that enables unbiased and reliable identification of A-to-I RNA editing sites throughout the transcriptome. Here, we describe our updated protocol for ICE-seq in the human transcriptome.

Nucleoside analogs in the study of the epitranscriptome

Over 150 unique RNA modifications are now known including several nonstandard nucleotides present in the body of messenger RNAs. These modifications can alter a transcript's function and are collectively referred to as the epitrancriptome. Chemically modified nucleoside analogs are poised to play an important role in the study of these epitranscriptomic marks. Introduced chemical features on nucleic acid strands provide unique structures or reactivity that can be used for downstream detection or quantification. Three methods are used in the field to synthesize RNA containing chemically modified nucleoside analogs. Nucleoside analogs can be introduced by metabolic labeling, via polymerases with modified nucleotide triphosphates or via phosphoramidite-based chemical synthesis. In this review, these methods for incorporation of nucleoside analogs will be discussed with specific recently published examples pertaining to the study of the epitranscriptome.