RNA modification mapping with JACUSA2

RNA editing regulates glutamatergic synapses in the frontal cortex of a molecular subtype of Amyotrophic Lateral Sclerosis

BACKGROUND

Amyotrophic Lateral Sclerosis (ALS) is a highly heterogenous neurodegenerative disorder that primarily affects upper and lower motor neurons, affecting additional cell types and brain regions. Underlying molecular mechanisms are still elusive, in part due to disease heterogeneity. Molecular disease subtyping through integrative analyses including RNA editing profiling is a novel approach for identification of molecular networks involved in pathogenesis.

METHODS

We aimed to highlight the role of RNA editing in ALS, focusing on the frontal cortex and the prevalent molecular disease subtype (ALS-Ox), previously determined by transcriptomic profile stratification. We established global RNA editing (editome) and gene expression (transcriptome) profiles in control and ALS-Ox cases, utilizing publicly available RNA-seq data (GSE153960) and an in-house analysis pipeline. Functional annotation and pathway analyses identified molecular processes affected by RNA editing alterations. Pearson correlation analyses assessed RNA editing effects on expression. Similar analyses on additional ALS-Ox and control samples (GSE124439) were performed for verification. Targeted re-sequencing and qRT-PCR analysis targeting CACNA1C, were performed using frontal cortex tissue from ALS and control samples (n = 3 samples/group).

RESULTS

We identified reduced global RNA editing in the frontal cortex of ALS-Ox cases. Differentially edited transcripts are enriched in synapses, particularly in the glutamatergic synapse pathway. Bioinformatic analyses on additional ALS-Ox and control RNA-seq data verified these findings. We identified increased recoding at the Q621R site in the GRIK2 transcript and determined positive correlations between RNA editing and gene expression alterations in ionotropic receptor subunits GRIA2, GRIA3 and the CACNA1C transcript, which encodes the pore forming subunit of a post-synaptic L-type calcium channel. Experimental data verified RNA editing alterations and editing-expression correlation in CACNA1C, highlighting CACNA1C as a target for further study.

CONCLUSIONS

We provide evidence on the involvement of RNA editing in the frontal cortex of an ALS molecular subtype, highlighting a modulatory role mediated though recoding and gene expression regulation on glutamatergic synapse related transcripts. We report RNA editing effects in disease-related transcripts and validated editing alterations in CACNA1C. Our study provides targets for further functional studies that could shed light in underlying disease mechanisms enabling novel therapeutic approaches.

Current progress in strategies to profile transcriptomic m6A modifications

Various methods have been developed so far for detecting N 6-methyladenosine (m6A). The total m6A level or the m6A status at individual positions on mRNA can be detected and quantified through some sequencing-independent biochemical methods, such as LC/MS, SCARLET, SELECT, and m6A-ELISA. However, the m6A-detection techniques relying on high-throughput sequencing have more effectively advanced the understanding about biological significance of m6A-containing mRNA and m6A pathway at a transcriptomic level over the past decade. Various SGS-based (Second Generation Sequencing-based) methods with different detection principles have been widely employed for this purpose. These principles include m6A-enrichment using antibodies, discrimination of m6A from unmodified A-base by nucleases, a fusion protein strategy relying on RNA-editing enzymes, and marking m6A with chemical/biochemical reactions. Recently, TGS-based (Third Generation Sequencing-based) methods have brought a new trend by direct m6A-detection. This review first gives a brief introduction of current knowledge about m6A biogenesis and function, and then comprehensively describes m6A-profiling strategies including their principles, procedures, and features. This will guide users to pick appropriate methods according to research goals, give insights for developing novel techniques in varying areas, and continue to expand our boundary of knowledge on m6A.

Conquering new grounds: plant organellar C-to-U RNA editing factors can be functional in the plant cytosol

Plant mitochondrial and chloroplast transcripts are subject to numerous events of specific cytidine-to-uridine (C-to-U) RNA editing to correct genetic information. Key protein factors for this process are specific RNA-binding pentatricopeptide repeat (PPR) proteins, which are encoded in the nucleus and post-translationally imported into the two endosymbiotic organelles. Despite hundreds of C-to-U editing sites in the plant organelles, no comparable editing has been found for nucleo-cytosolic mRNAs raising the question why plant RNA editing is restricted to chloroplasts and mitochondria. Here, we addressed this issue in the model moss Physcomitrium patens, where all PPR-type RNA editing factors comprise specific RNA-binding and cytidine deamination functionalities in single proteins. To explore whether organelle-type RNA editing can principally also take place in the plant cytosol, we expressed PPR56, PPR65 and PPR78, three editing factors recently shown to also function in a bacterial setup, together with cytosolic co-transcribed native targets in Physcomitrium. While we obtained unsatisfying results upon their constitutive expression, we found strong cytosolic RNA editing under hormone-inducible expression. Moreover, RNA-Seq analyses revealed varying numbers of up to more than 900 off-targets in other cytosolic transcripts. We conclude that PPR-mediated C-to-U RNA editing is not per se incompatible with the plant cytosol but that its limited target specificity has restricted its occurrence to the much less complex transcriptomes of mitochondria and chloroplast in the course of evolution.

A-to-I Editing Is Subtype-Specific in Non-Hodgkin Lymphomas

Cancer is a complex and heterogeneous disease, in which a number of genetic and epigenetic changes occur in tumor onset and progression. Recent studies indicate that changes at the RNA level are also involved in tumorigenesis, such as adenosine-to-inosine (A-to-I) RNA editing. Here, we systematically investigate transcriptome-wide A-to-I editing events in a large number of samples from Non-Hodgkin lymphomas (NHLs). Using a computational pipeline that determines significant differences in editing level between NHL and normal samples at known A-to-I editing sites, we identify a number of differentially edited editing sites between NHL subtypes and normal samples. Most of the differentially edited sites are located in non-coding regions, and many such sites show a strong correlation between gene expression level and editing efficiency, indicating that RNA editing might have direct consequences for the cancer cell's aberrant gene regulation status in these cases. Moreover, we establish a strong link between RNA editing and NHL by demonstrating that NHL and normal samples and even NHL subtypes can be distinguished based on genome-wide RNA editing profiles alone. Our study establishes a strong link between RNA editing, cancer and aberrant gene regulation in NHL.

Divergent landscapes of A-to-I editing in postmortem and living human brain

Adenosine-to-inosine (A-to-I) editing is a prevalent post-transcriptional RNA modification within the brain. Yet, most research has relied on postmortem samples, assuming it is an accurate representation of RNA biology in the living brain. We challenge this assumption by comparing A-to-I editing between postmortem and living prefrontal cortical tissues. Major differences were found, with over 70,000 A-to-I sites showing higher editing levels in postmortem tissues. Increased A-to-I editing in postmortem tissues is linked to higher ADAR and ADARB1 expression, is more pronounced in non-neuronal cells, and indicative of postmortem activation of inflammation and hypoxia. Higher A-to-I editing in living tissues marks sites that are evolutionarily preserved, synaptic, developmentally timed, and disrupted in neurological conditions. Common genetic variants were also found to differentially affect A-to-I editing levels in living versus postmortem tissues. Collectively, these discoveries offer more nuanced and accurate insights into the regulatory mechanisms of RNA editing in the human brain.

Divergent landscapes of A-to-I editing in postmortem and living human brain

Adenosine-to-inosine (A-to-I) editing is a prevalent post-transcriptional RNA modification within the brain. Yet, most research has relied on postmortem samples, assuming it is an accurate representation of RNA biology in the living brain. We challenge this assumption by comparing A-to-I editing between postmortem and living prefrontal cortical tissues. Major differences were found, with over 70,000 A-to-I sites showing higher editing levels in postmortem tissues. Increased A-to-I editing in postmortem tissues is linked to higher ADAR1 and ADARB1 expression, is more pronounced in non-neuronal cells, and indicative of postmortem activation of inflammation and hypoxia. Higher A-to-I editing in living tissues marks sites that are evolutionarily preserved, synaptic, developmentally timed, and disrupted in neurological conditions. Common genetic variants were also found to differentially affect A-to-I editing levels in living versus postmortem tissues. Collectively, these discoveries illuminate the nuanced functions and intricate regulatory mechanisms of RNA editing within the human brain.

txtools: an R package facilitating analysis of RNA modifications, structures, and interactions

We present txtools, an R package that enables the processing, analysis, and visualization of RNA-seq data at the nucleotide-level resolution, seamlessly integrating alignments to the genome with transcriptomic representation. txtools' main inputs are BAM files and a transcriptome annotation, and the main output is a table, capturing mismatches, deletions, and the number of reads beginning and ending at each nucleotide in the transcriptomic space. txtools further facilitates downstream visualization and analyses. We showcase, using examples from the epitranscriptomic field, how a few calls to txtools functions can yield insightful and ready-to-publish results. txtools is of broad utility also in the context of structural mapping and RNA:protein interaction mapping. By providing a simple and intuitive framework, we believe that txtools will be a useful and convenient tool and pave the path for future discovery. txtools is available for installation from its GitHub repository at https://github.com/AngelCampos/txtools.

No evidence for ac4C within human mRNA upon data reassessment

Cytidine acetylation (ac4C) of RNA is a post-transcriptional modification catalyzed by Nat10. Recently, an approach termed RedaC:T was employed to map ac4C in human mRNA, relying on detection of C>T mutations in WT but not in Nat10-KO cells. RedaC:T suggested widespread ac4C presence. Here, we reanalyze RedaC:T data. We find that mismatch signatures are not reproducible, as C>T mismatches are nearly exclusively present in only one of two biological replicates. Furthermore, all mismatch types-not only C>T-are highly enriched in WT samples, inconsistent with an acetylation signature. We demonstrate that the originally observed enrichment in mutations in one of the WT samples is due to its low complexity, resulting in the technical amplification of all classes of mismatch counts. Removal of duplicate reads abolishes the skewed mismatch patterns. These analyses account for the irreproducible mismatch patterns across samples while failing to find evidence for acetylation of RedaC:T sites.

Detection of ac4C in human mRNA is preserved upon data reassessment

We recently reported the distribution of N4-acetylcytidine (ac4C) in HeLa mRNA at base resolution through chemical reduction and the induction of C:T mismatches in sequencing (RedaC:T-seq). Our results contradicted an earlier report from Schwartz and colleagues utilizing a similar method termed ac4C-seq. Here, we revisit both datasets and reaffirm our findings. Through RedaC:T-seq reanalysis, we establish a low basal error rate at unmodified nucleotides that is not skewed to any specific mismatch type and a prominent increase in C:T substitutions as the dominant mismatch type in both treated wild-type replicates, with a high degree of reproducibility across replicates. In contrast, through ac4C-seq reanalysis, we uncover significant data quality issues including insufficient depth, with one wild-type replicate yielding 2.7 million reads, inconsistencies in reduction efficiencies between replicates, and an overall increase in mismatches involving thymine that could obscure ac4C detection. These analyses bolster the detection of ac4C in HeLa mRNA through RedaC:T-seq.

Recent developments, opportunities, and challenges in the study of mRNA pseudouridylation

Pseudouridine is an abundant mRNA modification found in diverse organisms ranging from bacteria and viruses to multicellular plants and humans. New developments in pseudouridine profiling provide quantitative tools to map mRNA pseudouridylation sites. Sparse biochemical studies establish the potential for mRNA pseudouridylation to affect most stages of the mRNA life cycle from birth to death. This recent progress sets the stage for deeper investigations into the molecular and cellular functions of specific mRNA pseudouridines, including in disease.

Recommendations for detection, validation, and evaluation of RNA editing events in cardiovascular and neurological/neurodegenerative diseases

RNA editing, a common and potentially highly functional form of RNA modification, encompasses two different RNA modifications, namely adenosine to inosine (A-to-I) and cytidine to uridine (C-to-U) editing. As inosines are interpreted as guanosines by the cellular machinery, both A-to-I and C-to-U editing change the nucleotide sequence of the RNA. Editing events in coding sequences have the potential to change the amino acid sequence of proteins, whereas editing events in noncoding RNAs can, for example, affect microRNA target binding. With advancing RNA sequencing technology, more RNA editing events are being discovered, studied, and reported. However, RNA editing events are still often overlooked or discarded as sequence read quality defects. With this position paper, we aim to provide guidelines and recommendations for the detection, validation, and follow-up experiments to study RNA editing, taking examples from the fields of cardiovascular and brain disease. We discuss all steps, from sample collection, storage, and preparation, to different strategies for RNA sequencing and editing-sensitive data analysis strategies, to validation and follow-up experiments, as well as potential pitfalls and gaps in the available technologies. This paper may be used as an experimental guideline for RNA editing studies in any disease context.

Simultaneous nanopore profiling of mRNA m6A and pseudouridine reveals translation coordination

N6-methyladenosine (m6A) and pseudouridine (Ψ) are the two most abundant modifications in mammalian messenger RNA, but the coordination of their biological functions remains poorly understood. We develop a machine learning-based nanopore direct RNA sequencing method (NanoSPA) that simultaneously analyzes m6A and Ψ in the human transcriptome. Applying NanoSPA to polysome profiling, we reveal opposing transcriptomic co-occurrence of m6A and Ψ and synergistic, hierarchical effects of m6A and Ψ on the polysome.

m6A-Atlas v2.0: updated resources for unraveling the N6-methyladenosine (m6A) epitranscriptome among multiple species

N 6-Methyladenosine (m6A) is one of the most abundant internal chemical modifications on eukaryote mRNA and is involved in numerous essential molecular functions and biological processes. To facilitate the study of this important post-transcriptional modification, we present here m6A-Atlas v2.0, an updated version of m6A-Atlas. It was expanded to include a total of 797 091 reliable m6A sites from 13 high-resolution technologies and two single-cell m6A profiles. Additionally, three methods (exomePeaks2, MACS2 and TRESS) were used to identify >16 million m6A enrichment peaks from 2712 MeRIP-seq experiments covering 651 conditions in 42 species. Quality control results of MeRIP-seq samples were also provided to help users to select reliable peaks. We also estimated the condition-specific quantitative m6A profiles (i.e. differential methylation) under 172 experimental conditions for 19 species. Further, to provide insights into potential functional circuitry, the m6A epitranscriptomics were annotated with various genomic features, interactions with RNA-binding proteins and microRNA, potentially linked splicing events and single nucleotide polymorphisms. The collected m6A sites and their functional annotations can be freely queried and downloaded via a user-friendly graphical interface at: http://rnamd.org/m6a.

Experience with German Research Consortia in the Field of Chemical Biology of Native Nucleic Acid Modifications

The chemical biology of native nucleic acid modifications has seen an intense upswing, first concerning DNA modifications in the field of epigenetics and then concerning RNA modifications in a field that was correspondingly rebaptized epitranscriptomics by analogy. The German Research Foundation (DFG) has funded several consortia with a scientific focus in these fields, strengthening the traditionally well-developed nucleic acid chemistry community and inciting it to team up with colleagues from the life sciences and data science to tackle interdisciplinary challenges. This Perspective focuses on the genesis, scientific outcome, and downstream impact of the DFG priority program SPP1784 and offers insight into how it fecundated further consortia in the field. Pertinent research was funded from mid-2015 to 2022, including an extension related to the coronavirus pandemic. Despite being a detriment to research activity in general, the pandemic has resulted in tremendously boosted interest in the field of RNA and RNA modifications as a consequence of their widespread and successful use in vaccination campaigns against SARS-CoV-2. Funded principal investigators published over 250 pertinent papers with a very substantial impact on the field. The program also helped to redirect numerous laboratories toward this dynamic field. Finally, SPP1784 spawned initiatives for several funded consortia that continue to drive the fields of nucleic acid modification.

Adaptive sampling for nanopore direct RNA-sequencing

Nanopore long-read sequencing enables real-time monitoring and controlling of individual nanopores. This allows us to enrich or deplete specific sequences in DNA sequencing in a process called "adaptive sampling." So far, adaptive sampling (AS) was not applicable to the direct sequencing of RNA. Here, we show that AS is feasible and useful for direct RNA sequencing (DRS), which has its specific technical and biological challenges. Using a well-controlled in vitro transcript-based model system, we identify essential characteristics and parameter settings for AS in DRS, as the superior performance of depletion over enrichment. Here, the efficiency of depletion is close to the theoretical maximum. Additionally, we demonstrate that AS efficiently depletes specific transcripts in transcriptome-wide sequencing applications. Specifically, we applied our AS approach to poly(A)-enriched RNA samples from human-induced pluripotent stem cell-derived cardiomyocytes and mouse whole heart tissue and show efficient 2.5- to 2.8-fold depletion of highly abundant mitochondrial-encoded transcripts. Finally, we characterize depletion and enrichment performance for complex transcriptome subsets, that is, at the level of the entire Chromosome 11, proving the general applicability of direct RNA AS. Our analyses provide evidence that AS is especially useful to enable the detection of lowly expressed transcripts and reduce the sequencing of highly abundant disturbing transcripts.

Multi-omics analysis reveals the molecular response to heat stress in a "red tide" dinoflagellate

BACKGROUND

"Red tides" are harmful algal blooms caused by dinoflagellate microalgae that accumulate toxins lethal to other organisms, including humans via consumption of contaminated seafood. These algal blooms are driven by a combination of environmental factors including nutrient enrichment, particularly in warm waters, and are increasingly frequent. The molecular, regulatory, and evolutionary mechanisms that underlie the heat stress response in these harmful bloom-forming algal species remain little understood, due in part to the limited genomic resources from dinoflagellates, complicated by the large sizes of genomes, exhibiting features atypical of eukaryotes.

RESULTS

We present the de novo assembled genome (~ 4.75 Gbp with 85,849 protein-coding genes), transcriptome, proteome, and metabolome from Prorocentrum cordatum, a globally abundant, bloom-forming dinoflagellate. Using axenic algal cultures, we study the molecular mechanisms that underpin the algal response to heat stress, which is relevant to current ocean warming trends. We present the first evidence of a complementary interplay between RNA editing and exon usage that regulates the expression and functional diversity of biomolecules, reflected by reduction in photosynthesis, central metabolism, and protein synthesis. These results reveal genomic signatures and post-transcriptional regulation for the first time in a pelagic dinoflagellate.

CONCLUSIONS

Our multi-omics analyses uncover the molecular response to heat stress in an important bloom-forming algal species, which is driven by complex gene structures in a large, high-G+C genome, combined with multi-level transcriptional regulation. The dynamics and interplay of molecular regulatory mechanisms may explain in part how dinoflagellates diversified to become some of the most ecologically successful organisms on Earth.

Detection of queuosine and queuosine precursors in tRNAs by direct RNA sequencing

Queuosine (Q) is a complex tRNA modification found in bacteria and eukaryotes at position 34 of four tRNAs with a GUN anticodon, and it regulates the translational efficiency and fidelity of the respective codons that differ at the Wobble position. In bacteria, the biosynthesis of Q involves two precursors, preQ0 and preQ1, whereas eukaryotes directly obtain Q from bacterial sources. The study of queuosine has been challenging due to the limited availability of high-throughput methods for its detection and analysis. Here, we have employed direct RNA sequencing using nanopore technology to detect the modification of tRNAs with Q and Q precursors. These modifications were detected with high accuracy on synthetic tRNAs as well as on tRNAs extracted from Schizosaccharomyces pombe and Escherichia coli by comparing unmodified to modified tRNAs using the tool JACUSA2. Furthermore, we present an improved protocol for the alignment of raw sequence reads that gives high specificity and recall for tRNAs ex cellulo that, by nature, carry multiple modifications. Altogether, our results show that 7-deazaguanine-derivatives such as queuosine are readily detectable using direct RNA sequencing. This advancement opens up new possibilities for investigating these modifications in native tRNAs, furthering our understanding of their biological function.

MePMe-seq: antibody-free simultaneous m6A and m5C mapping in mRNA by metabolic propargyl labeling and sequencing

Internal modifications of mRNA have emerged as widespread and versatile regulatory mechanism to control gene expression at the post-transcriptional level. Most of these modifications are methyl groups, making S-adenosyl-L-methionine (SAM) a central metabolic hub. Here we show that metabolic labeling with a clickable metabolic precursor of SAM, propargyl-selenohomocysteine (PSH), enables detection and identification of various methylation sites. Propargylated A, C, and G nucleosides form at detectable amounts via intracellular generation of the corresponding SAM analogue. Integration into next generation sequencing enables mapping of N6-methyladenosine (m6A) and 5-methylcytidine (m5C) sites in mRNA with single nucleotide precision (MePMe-seq). Analysis of the termination profiles can be used to distinguish m6A from 2'-O-methyladenosine (Am) and N1-methyladenosine (m1A) sites. MePMe-seq overcomes the problems of antibodies for enrichment and sequence-motifs for evaluation, which was limiting previous methodologies. Metabolic labeling via clickable SAM facilitates the joint evaluation of methylation sites in RNA and potentially DNA and proteins.

Comprehensive mapping of exon junction complex binding sites reveals universal EJC deposition in Drosophila

BACKGROUND

The exon junction complex (EJC) is involved in most steps of the mRNA life cycle, ranging from splicing to nonsense-mediated mRNA decay (NMD). It is assembled by the splicing machinery onto mRNA in a sequence-independent manner. A fundamental open question is whether the EJC is deposited onto all exon‒exon junctions or only on a subset of them. Several previous studies have made observations supportive of the latter, yet these have been limited by methodological constraints.

RESULTS

In this study, we sought to overcome these limitations via the integration of two different approaches for transcriptome-wide mapping of EJCs. Our results revealed that nearly all, if not all, internal exons consistently harbor an EJC in Drosophila, demonstrating that EJC presence is an inherent consequence of the splicing reaction. Furthermore, our study underscores the limitations of eCLIP methods in fully elucidating the landscape of RBP binding sites. Our findings highlight how highly specific (low false positive) methodologies can lead to erroneous interpretations due to partial sensitivity (high false negatives).

CONCLUSIONS

This study contributes to our understanding of EJC deposition and its association with pre-mRNA splicing. The universal presence of EJC on internal exons underscores its significance in ensuring proper mRNA processing. Additionally, our observations highlight the need to consider both specificity and sensitivity in RBP mapping methodologies.

Nanopore-based direct sequencing of RNA transcripts with 10 different modified nucleotides reveals gaps in existing technology

RNA undergoes complex posttranscriptional processing including chemical modifications of the nucleotides. The resultant-modified nucleotides are an integral part of RNA sequences that must be considered in studying the biology of RNA and in the design of RNA therapeutics. However, the current "RNA-sequencing" methods primarily sequence complementary DNA rather than RNA itself, which means that the modifications present in RNA are not captured in the sequencing results. Emerging direct RNA-sequencing technologies, such as those offered by Oxford Nanopore, aim to address this limitation. In this study, we synthesized and used Nanopore technology to sequence RNA transcripts consisting of canonical nucleotides and 10 different modifications in various concentrations. The results show that direct RNA sequencing still has a baseline error rate of >10%, and although some modifications can be detected, many remain unidentified. Thus, there is a need to develop sequencing technologies and analysis methods that can comprehensively capture the total complexity of RNA. The RNA sequences obtained through this project are made available for benchmarking analysis methods.

FLARE: a fast and flexible workflow for identifying RNA editing foci

BACKGROUND

Fusion of RNA-binding proteins (RBPs) to RNA base-editing enzymes (such as APOBEC1 or ADAR) has emerged as a powerful tool for the discovery of RBP binding sites. However, current methods that analyze sequencing data from RNA-base editing experiments are vulnerable to false positives due to off-target editing, genetic variation and sequencing errors.

RESULTS

We present FLagging Areas of RNA-editing Enrichment (FLARE), a Snakemake-based pipeline that builds on the outputs of the SAILOR edit site discovery tool to identify regions statistically enriched for RNA editing. FLARE can be configured to analyze any type of RNA editing, including C to U and A to I. We applied FLARE to C-to-U editing data from a RBFOX2-APOBEC1 STAMP experiment, to show that our approach attains high specificity for detecting RBFOX2 binding sites. We also applied FLARE to detect regions of exogenously introduced as well as endogenous A-to-I editing.

CONCLUSIONS

FLARE is a fast and flexible workflow that identifies significantly edited regions from RNA-seq data. The FLARE codebase is available at https://github.com/YeoLab/FLARE .

Benchmarking RNA Editing Detection Tools

RNA, like DNA and proteins, can undergo modifications. To date, over 170 RNA modifications have been identified, leading to the emergence of a new research area known as epitranscriptomics. RNA editing is the most frequent RNA modification in mammalian transcriptomes, and two types have been identified: (1) the most frequent, adenosine to inosine (A-to-I); and (2) the less frequent, cysteine to uracil (C-to-U) RNA editing. Unlike other epitranscriptomic marks, RNA editing can be readily detected from RNA sequencing (RNA-seq) data without any chemical conversions of RNA before sequencing library preparation. Furthermore, analyzing RNA editing patterns from transcriptomic data provides an additional layer of information about the epitranscriptome. As the significance of epitranscriptomics, particularly RNA editing, gains recognition in various fields of biology and medicine, there is a growing interest in detecting RNA editing sites (RES) by analyzing RNA-seq data. To cope with this increased interest, several bioinformatic tools are available. However, each tool has its advantages and disadvantages, which makes the choice of the most appropriate tool for bench scientists and clinicians difficult. Here, we have benchmarked bioinformatic tools to detect RES from RNA-seq data. We provide a comprehensive view of each tool and its performance using previously published RNA-seq data to suggest recommendations on the most appropriate for utilization in future studies.

PIE-seq: identifying RNA-binding protein targets by dual RNA-deaminase editing and sequencing

RNA-binding proteins (RBPs) are essential for gene regulation, but it remains a challenge to identify their RNA targets across cell types. Here we present PIE-Seq to investigate Protein-RNA Interaction with dual-deaminase Editing and Sequencing by conjugating C-to-U and A-to-I base editors to RBPs. We benchmark PIE-Seq and demonstrate its sensitivity in single cells, its application in the developing brain, and its scalability with 25 human RBPs. Bulk PIE-Seq identifies canonical binding features for RBPs such as PUM2 and NOVA1, and nominates additional target genes for most tested RBPs such as SRSF1 and TDP-43/TARDBP. Homologous RBPs frequently edit similar sequences and gene sets in PIE-Seq while different RBP families show distinct targets. Single-cell PIE-PUM2 uncovers comparable targets to bulk samples and applying PIE-PUM2 to the developing mouse neocortex identifies neural-progenitor- and neuron-specific target genes such as App. In summary, PIE-Seq provides an orthogonal approach and resource to uncover RBP targets in mice and human cells.

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.

Comprehensive identification of diverse ribosomal RNA modifications by targeted nanopore direct RNA sequencing and JACUSA2

Ribosomal RNAs are decorated by numerous post-transcriptional modifications whose exact roles in ribosome biogenesis, function, and human pathophysiology remain largely unknown. Here, we report a targeted direct rRNA sequencing approach involving a substrate selection step and demonstrate its suitability to identify differential modification sites in combination with the JACUSA2 software. We compared JACUSA2 to other tools designed for RNA modification detection and show that JACUSA2 outperforms other software with regard to detection of base modifications such as methylation, acetylation and aminocarboxypropylation. To illustrate its widespread usability, we applied our method to a collection of CRISPR-Cas9 engineered colon carcinoma cells lacking specific enzymatic activities responsible for particular rRNA modifications and systematically compared them to isogenic wild-type RNAs. Besides the numerous 2'-O methylated riboses and pseudouridylated residues, our approach was suitable to reliably identify differential base methylation and acetylation events. Importantly, our method does not require any prior knowledge of modification sites or the need to train complex models. We further report for the first time detection of human rRNA modifications by direct RNA-sequencing on Flongle flow cells, the smallest-scale nanopore flow cell available to date. The use of these smaller flow cells reduces RNA input requirements, making our workflow suitable for the analysis of samples with limited availability and clinical work.

DirectRMDB: a database of post-transcriptional RNA modifications unveiled from direct RNA sequencing technology

With advanced technologies to map RNA modifications, our understanding of them has been revolutionized, and they are seen to be far more widespread and important than previously thought. Current next-generation sequencing (NGS)-based modification profiling methods are blind to RNA modifications and thus require selective chemical treatment or antibody immunoprecipitation methods for particular modification types. They also face the problem of short read length, isoform ambiguities, biases and artifacts. Direct RNA sequencing (DRS) technologies, commercialized by Oxford Nanopore Technologies (ONT), enable the direct interrogation of any given modification present in individual transcripts and promise to address the limitations of previous NGS-based methods. Here, we present the first ONT-based database of quantitative RNA modification profiles, DirectRMDB, which includes 16 types of modification and a total of 904,712 modification sites in 25 species identified from 39 independent studies. In addition to standard functions adopted by existing databases, such as gene annotations and post-transcriptional association analysis, we provide a fresh view of RNA modifications, which enables exploration of the epitranscriptome in an isoform-specific manner. The DirectRMDB database is freely available at: http://www.rnamd.org/directRMDB/.

Spatiotemporal and genetic regulation of A-to-I editing throughout human brain development

Posttranscriptional RNA modifications by adenosine-to-inosine (A-to-I) editing are abundant in the brain, yet elucidating functional sites remains challenging. To bridge this gap, we investigate spatiotemporal and genetically regulated A-to-I editing sites across prenatal and postnatal stages of human brain development. More than 10,000 spatiotemporally regulated A-to-I sites were identified that occur predominately in 3' UTRs and introns, as well as 37 sites that recode amino acids in protein coding regions with precise changes in editing levels across development. Hyper-edited transcripts are also enriched in the aging brain and stabilize RNA secondary structures. These features are conserved in murine and non-human primate models of neurodevelopment. Finally, thousands of cis-editing quantitative trait loci (edQTLs) were identified with unique regulatory effects during prenatal and postnatal development. Collectively, this work offers a resolved atlas linking spatiotemporal variation in editing levels to genetic regulatory effects throughout distinct stages of brain maturation.

Plant mitochondrial RNA editing factors can perform targeted C-to-U editing of nuclear transcripts in human cells

RNA editing processes are strikingly different in animals and plants. Up to thousands of specific cytidines are converted into uridines in plant chloroplasts and mitochondria whereas up to millions of adenosines are converted into inosines in animal nucleo-cytosolic RNAs. It is unknown whether these two different RNA editing machineries are mutually incompatible. RNA-binding pentatricopeptide repeat (PPR) proteins are the key factors of plant organelle cytidine-to-uridine RNA editing. The complete absence of PPR mediated editing of cytosolic RNAs might be due to a yet unknown barrier that prevents its activity in the cytosol. Here, we transferred two plant mitochondrial PPR-type editing factors into human cell lines to explore whether they could operate in the nucleo-cytosolic environment. PPR56 and PPR65 not only faithfully edited their native, co-transcribed targets but also different sets of off-targets in the human background transcriptome. More than 900 of such off-targets with editing efficiencies up to 91%, largely explained by known PPR-RNA binding properties, were identified for PPR56. Engineering two crucial amino acid positions in its PPR array led to predictable shifts in target recognition. We conclude that plant PPR editing factors can operate in the entirely different genetic environment of the human nucleo-cytosol and can be intentionally re-engineered towards new targets.

ADAR1 averts fatal type I interferon induction by ZBP1

Mutations of the ADAR1 gene encoding an RNA deaminase cause severe diseases associated with chronic activation of type I interferon (IFN) responses, including Aicardi–Goutières syndrome and bilateral striatal necrosis^1–3. The IFN-inducible p150 isoform of ADAR1 contains a Zα domain that recognizes RNA with an alternative left-handed double-helix structure, termed Z-RNA^4,5. Hemizygous ADAR1 mutations in the Zα domain cause type I IFN-mediated pathologies in humans^2,3 and mice^6–8; however, it remains unclear how the interaction of ADAR1 with Z-RNA prevents IFN activation. Here we show that Z-DNA-binding protein 1 (ZBP1), the only other protein in mammals known to harbour Zα domains⁹, promotes type I IFN activation and fatal pathology in mice with impaired ADAR1 function. ZBP1 deficiency or mutation of its Zα domains reduced the expression of IFN-stimulated genes and largely prevented early postnatal lethality in mice with hemizygous expression of ADAR1 with mutated Zα domain (Adar1^mZα/– mice). Adar1^mZα/– mice showed upregulation and impaired editing of endogenous retroelement-derived complementary RNA reads, which represent a likely source of Z-RNAs activating ZBP1. Notably, ZBP1 promoted IFN activation and severe pathology in Adar1^mZα/– mice in a manner independent of RIPK1, RIPK3, MLKL-mediated necroptosis and caspase-8-dependent apoptosis, suggesting a novel mechanism of action. Thus, ADAR1 prevents endogenous Z-RNA-dependent activation of pathogenic type I IFN responses by ZBP1, suggesting that ZBP1 could contribute to type I interferonopathies caused by ADAR1 mutations.

ADAR1 prevents Z-RNA-dependent activation of pathogenic type I interferon responses by ZBP1, whose activity may contribute to pathology in type I interferonopathies with ADAR1 mutations.