Detecting the epitranscriptome

Deep-m5U: a deep learning-based approach for RNA 5-methyluridine modification prediction using optimized feature integration

BACKGROUND

RNA 5-methyluridine (m5U) modifications play a crucial role in biological processes, making their accurate identification a key focus in computational biology. This paper introduces Deep-m5U, a robust predictor designed to enhance the prediction of m5U modifications. The proposed method, named Deep-m5U, utilizes a hybrid pseudo-K-tuple nucleotide composition (PseKNC) for sequence formulation, a Shapley Additive exPlanations (SHAP) algorithm for discriminant feature selection, and a deep neural network (DNN) as the classifier.

RESULTS

The model was evaluated using two benchmark datasets, i.e., Full Transcript and Mature mRNA. Deep-m5U achieved overall accuracies of 91.47% and 95.86% for the Full Transcript and Mature mRNA datasets with 10-fold cross-validation, and for independent samples, the model attained 92.94% and 95.17% accuracy.

CONCLUSION

Compared to existing models, Deep-m5U showed approximately 5.23% and 3.73% higher accuracy on the training data and 3.95% and 3.26% higher accuracy on independent samples for the Full Transcript and Mature mRNA datasets, respectively. The reliability and effectiveness of Deep-m5U make it a valuable tool for scientists and a potential asset in pharmaceutical design and research.

Analysis of RNA and Its Modifications

Ribonucleic acids (RNAs) are key biomolecules responsible for the transmission of genetic information, the synthesis of proteins, and modulation of many biochemical processes. They are also often the key components of viruses. Synthetic RNAs or oligoribonucleotides are becoming more widely used as therapeutics. In many cases, RNAs will be chemically modified, either naturally via enzymatic systems within a cell or intentionally during their synthesis. Analytical methods to detect, sequence, identify, and quantify RNA and its modifications have demands that far exceed requirements found in the DNA realm. Two complementary platforms have demonstrated their value and utility for the characterization of RNA and its modifications: mass spectrometry and next-generation sequencing. This review highlights recent advances in both platforms, examines their relative strengths and weaknesses, and explores some alternative approaches that lie at the horizon.

In silico λ-dynamics predicts protein binding specificities to modified RNAs

RNA modifications shape gene expression through a smorgasbord of chemical changes to canonical RNA bases. Although numbering in the hundreds, only a few RNA modifications are well characterized, in part due to the absence of methods to identify modification sites. Antibodies remain a common tool to identify modified RNA and infer modification sites through straightforward applications. However, specificity issues can result in off-target binding and confound conclusions. This work utilizes in silico λ-dynamics to efficiently estimate binding free energy differences of modification-targeting antibodies between a variety of naturally occurring RNA modifications. Crystal structures of inosine and N6-methyladenosine (m6A) targeting antibodies bound to their modified ribonucleosides were determined and served as structural starting points. λ-Dynamics was utilized to predict RNA modifications that permit or inhibit binding to these antibodies. In vitro RNA-antibody binding assays supported the accuracy of these in silico results. High agreement between experimental and computed binding propensities demonstrated that λ-dynamics can serve as a predictive screen for antibody specificity against libraries of RNA modifications. More importantly, this strategy is an innovative way to elucidate how hundreds of known RNA modifications interact with biological molecules without the limitations imposed by in vitro or in vivo methodologies.

Optimized chemical labeling method for isolation of 8-oxoG-modified RNA, ChLoRox-Seq, identifies mRNAs enriched in oxidation and transcriptome-wide distribution biases of oxidation events post environmental stress

Bulk increases in nucleobase oxidation, most commonly manifesting as the guanine (G) nucleobase modification 8-oxo-7,8-dihydroguanine (8-oxoG), have been linked to several disease pathologies. Elucidating the effects of RNA oxidation on cellular homoeostasis is limited by a lack of effective tools for detecting specific regions modified with 8-oxoG. Building on a previously published method for studying 8-oxoG in DNA, we developed ChLoRox-Seq, which works by covalently functionalizing 8-oxoG sites in RNA with biotin. Importantly, this method enables antibody-free enrichment of 8-oxoG-containing RNA fragments for Next Generation Sequencing-based detection of modified regions transcriptome-wide. We demonstrate the high specificity of ChLoRox-Seq for functionalizing 8-oxoG over unmodified nucleobases in RNA and benchmark this specificity to a commonly used antibody-based approach. Key advantages of ChLoRox-Seq include: (1) heightened resolution of RNA oxidation regions (e.g. exon-level) and (2) lower experimental costs. By applying ChLoRox-Seq to mRNA extracted from human lung epithelial cells (BEAS-2B) after exposure to environmentally relevant stress, we observe that 8-oxoG modifications tend to cluster in regions that are G-rich and within mRNA transcripts possessing longer 5' UTR and CDS regions. These findings provide new insight into the complex mechanisms that bias the accumulation of RNA oxidation across the transcriptome. Notably, our analysis suggests the possibility that most mRNA oxidation events are probabilistically driven and that mRNAs that possess more favourable intrinsic properties are prone to incur oxidation events at elevated rates. ChLoRox-Seq can be readily applied in future studies to identify regions of elevated RNA oxidation in any cellular model of interest.

Viral oncogenes, viruses, and cancer: a third-generation sequencing perspective on viral integration into the human genome

The link between viruses and cancer has intrigued scientists for decades. Certain viruses have been shown to be vital in the development of various cancers by integrating viral DNA into the host genome and activating viral oncogenes. These viruses include the Human Papillomavirus (HPV), Hepatitis B and C Viruses (HBV and HCV), Epstein-Barr Virus (EBV), and Human T-Cell Leukemia Virus (HTLV-1), which are all linked to the development of a myriad of human cancers. Third-generation sequencing technologies have revolutionized our ability to study viral integration events at unprecedented resolution in recent years. They offer long sequencing capabilities along with the ability to map viral integration sites, assess host gene expression, and track clonal evolution in cancer cells. Recently, researchers have been exploring the application of Oxford Nanopore Technologies (ONT) nanopore sequencing and Pacific BioSciences (PacBio) single-molecule real-time (SMRT) sequencing in cancer research. As viral integration is crucial to the development of cancer via viruses, third-generation sequencing would provide a novel approach to studying the relationship interlinking viral oncogenes, viruses, and cancer. This review article explores the molecular mechanisms underlying viral oncogenesis, the role of viruses in cancer development, and the impact of third-generation sequencing on our understanding of viral integration into the human genome.

Changes in m6A in Steatotic Liver Disease

Fatty liver disease is one of the major causes of morbidity and mortality worldwide. Fatty liver includes non-alcoholic fatty liver disease (NAFLD) and non-alcoholic steatohepatitis (NASH), now replaced by a consensus group as metabolic dysfunction-associated steatotic liver disease (MASLD). While excess nutrition and obesity are major contributors to fatty liver, the underlying mechanisms remain largely unknown and therapeutic interventions are limited. Reversible chemical modifications in RNA are newly recognized critical regulators controlling post-transcriptional gene expression. Among these modifications, N6-methyladenosine (m6A) is the most abundant and regulates transcript abundance in fatty liver disease. Modulation of m6A by readers, writers, and erasers (RWE) impacts mRNA processing, translation, nuclear export, localization, and degradation. While many studies focus on m6A RWE expression in human liver pathologies, limitations of technology and bioinformatic methods to detect m6A present challenges in understanding the epitranscriptomic mechanisms driving fatty liver disease progression. In this review, we summarize the RWE of m6A and current methods of detecting m6A in specific genes associated with fatty liver disease.

Disruption of the mouse liver epitranscriptome by long-term aroclor 1260 exposure

Chronic environmental exposure to polychlorinated biphenyls (PCBs) is associated with non-alcoholic fatty liver disease (NAFLD) and exacerbated by a high fat diet (HFD). Here, chronic (34 wks.) exposure of low fat diet (LFD)-fed male mice to Aroclor 1260 (Ar1260), a non-dioxin-like (NDL) mixture of PCBs, resulted in steatohepatitis and NAFLD. Twelve hepatic RNA modifications were altered with Ar1260 exposure including reduced abundance of 2'-O-methyladenosine (Am) and N(6)-methyladenosine (m6A), in contrast to increased Am in the livers of HFD-fed, Ar1260-exposed mice reported previously. Differences in 13 RNA modifications between LFD- and HFD- fed mice, suggest that diet regulates the liver epitranscriptome. Integrated network analysis of epitranscriptomic modifications identified a NRF2 (Nfe2l2) pathway in the chronic, LFD, Ar1260-exposed livers and an NFATC4 (Nfatc4) pathway for LFD- vs. HFD-fed mice. Changes in protein abundance were validated. The results demonstrate that diet and Ar1260 exposure alter the liver epitranscriptome in pathways associated with NAFLD.

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.

Recent advances in the plant epitranscriptome

Chemical modifications of RNAs, known as the epitranscriptome, are emerging as widespread regulatory mechanisms underlying gene regulation. The field of epitranscriptomics advances recently due to improved transcriptome-wide sequencing strategies for mapping RNA modifications and intensive characterization of writers, erasers, and readers that deposit, remove, and recognize RNA modifications, respectively. Herein, we review recent advances in characterizing plant epitranscriptome and its regulatory mechanisms in post-transcriptional gene regulation and diverse physiological processes, with main emphasis on N⁶-methyladenosine (m⁶A) and 5-methylcytosine (m⁵C). We also discuss the potential and challenges for utilization of epitranscriptome editing in crop improvement.

Interaction between N6-methyladenosine (m6A) modification and environmental chemical-induced diseases in various organ systems

A wide variety of chemicals are ubiquitous in the environment and thus exposure to these environmental chemicals poses a serious threat to public health. Particularly, environmental factors such as air pollution, heavy metals, and endocrine-disrupting chemicals (EDCs) can lead to diseases in various organ systems. Recent research in environmental epigenetics has demonstrated that N6-methyladenosine (m6A) modification is a key mechanism of environment-related diseases. m6A modification is the most abundant chemical modification in mRNAs, which can specifically regulate gene expression by affecting RNA translation, stability, processing, and nuclear export. In this review, we discussed how environmental chemicals affected m6A modification and mediated environment-related disease occurrence by classifying the diseases of various systems. Here, we conclude that environmental chemicals alter the levels of m6A and its modulators, which then participate in the occurrence of diseases in various systems by regulating gene expression and downstream signaling pathways such as METTL3/m6A ZBTB4/YTHDF2/EZH2, Foxo3a/FTO/m6A ephrin-B2/YTHDF2, and HIF1A/METTL3/m6A BIRC5/IGF2BP3/VEGFA. Considering the significant role of m6A and its modulators in response to environmental chemicals, they are expected to be used as biomarkers of environment-related diseases. Additionally, targeting m6A modulators using small molecule inhibitors and activators is expected to be a new method for the treatment of environment-related diseases. This review systematically and comprehensively clarifies the important role of m6A in diseases caused by environmental chemicals, thus establishing a scientific basis for the treatment of diseases in various organ systems.

Emerging Mutual Regulatory Roles between m6A Modification and microRNAs

N⁶-metyladenosine (m⁶A), one of the most common RNA methylation modifications in mammals, has attracted extensive attentions owing to its regulatory roles in a variety of physiological and pathological processes. As a reversible epigenetic modification on RNAs, m⁶A is dynamically mediated by the functional interplay among the regulatory proteins of methyltransferases, demethylases and methyl-binding proteins. In recent years, it has become increasingly clear that m⁶A modification is associated with the production and function of microRNAs (miRNAs). In this review, we summarize the specific kinds of m⁶A modification methyltransferases, demethylases and methyl-binding proteins. In particular, we focus on describing the roles of m⁶A modification and its regulatory proteins in the production and function of miRNAs in a variety of pathological and physiological processes. More importantly, we further discuss the mediating mechanisms of miRNAs in m⁶A modification and its regulatory proteins during the occurrence and development of various diseases.

Pathological and Therapeutic Advances in Parkinson's Disease: Mitochondria in the Interplay

Parkinson's disease (PD) is the second most common neurodegenerative illness majorly affecting the population between the ages of 55 to 65 years. Progressive dopaminergic neuronal loss and the collective assemblage of misfolded alpha-synuclein in the substantia nigra, remain notable neuro-pathological hallmarks of the disease. Multitudes of mechanistic pathways have been proposed in attempts to unravel the pathogenesis of PD but still, it remains elusive. The convergence of PD pathology is found in organelle dysfunction where mitochondria remain a major contributor. Mitochondrial processes like bioenergetics, mitochondrial dynamics, and mitophagy are under strict regulation by the mitochondrial genome and nuclear genome. These processes aggravate neurodegenerative activities upon alteration through neuroinflammation, oxidative damage, apoptosis, and proteostatic stress. Therefore, the mitochondria have grabbed a central position in the patho-mechanistic exploration of neurodegenerative diseases like PD. The management of PD remains a challenge to physicians to date, due to the variable therapeutic response of patients and the limitation of conventional chemical agents which only offer symptomatic relief with minimal to no disease-modifying effect. This review describes the patho-mechanistic pathways involved in PD not only limited to protein dyshomeostasis and oxidative stress, but explicit attention has been drawn to exploring mechanisms like organelle dysfunction, primarily mitochondria and mitochondrial genome influence, while delineating the newer exploratory targets such as GBA1, GLP, LRRK2, and miRNAs and therapeutic agents targeting them.

The host tRNA epitranscriptome: A new player in RNA virus infections

Viruses completely depend on the host translation machineries to express the viral proteins. Recent data reveal an unprecedented interaction of positive strand RNA ((+)RNA) viruses with the host tRNA epitranscriptome to favor viral protein expression via a specific reprogramming of codon optimality that ultimately favors decoding of the viral codons. We propose that this feature is shared by multiple RNA viruses and that the involved tRNA modifying enzymes represent promising novel targets for the development of broad-spectrum antivirals.

The Transition from Cancer "omics" to "epi-omics" through Next- and Third-Generation Sequencing

Deciphering cancer etiopathogenesis has proven to be an especially challenging task since the mechanisms that drive tumor development and progression are far from simple. An astonishing amount of research has revealed a wide spectrum of defects, including genomic abnormalities, epigenomic alterations, disturbance of gene transcription, as well as post-translational protein modifications, which cooperatively promote carcinogenesis. These findings suggest that the adoption of a multidimensional approach can provide a much more precise and comprehensive picture of the tumor landscape, hence serving as a powerful tool in cancer research and precision oncology. The introduction of next- and third-generation sequencing technologies paved the way for the decoding of genetic information and the elucidation of cancer-related cellular compounds and mechanisms. In the present review, we discuss the current and emerging applications of both generations of sequencing technologies, also referred to as massive parallel sequencing (MPS), in the fields of cancer genomics, transcriptomics and proteomics, as well as in the progressing realms of epi-omics. Finally, we provide a brief insight into the expanding scope of sequencing applications in personalized cancer medicine and pharmacogenomics.

Geographic encoding of transcripts enabled high-accuracy and isoform-aware deep learning of RNA methylation

As the most pervasive epigenetic mark present on mRNA and lncRNA, N⁶-methyladenosine (m⁶A) RNA methylation regulates all stages of RNA life in various biological processes and disease mechanisms. Computational methods for deciphering RNA modification have achieved great success in recent years; nevertheless, their potential remains underexploited. One reason for this is that existing models usually consider only the sequence of transcripts, ignoring the various regions (or geography) of transcripts such as 3′UTR and intron, where the epigenetic mark forms and functions. Here, we developed three simple yet powerful encoding schemes for transcripts to capture the submolecular geographic information of RNA, which is largely independent from sequences. We show that m⁶A prediction models based on geographic information alone can achieve comparable performances to classic sequence-based methods. Importantly, geographic information substantially enhances the accuracy of sequence-based models, enables isoform- and tissue-specific prediction of m⁶A sites, and improves m⁶A signal detection from direct RNA sequencing data. The geographic encoding schemes we developed have exhibited strong interpretability, and are applicable to not only m⁶A but also N¹-methyladenosine (m¹A), and can serve as a general and effective complement to the widely used sequence encoding schemes in deep learning applications concerning RNA transcripts.

More than a duologue: In-depth insights into epitranscriptomics and ferroptosis

Beyond transcription, RNA molecules are enzymatically modified to influence the biological functions of living organisms. The term “epitranscriptomics” describes the changes in RNA strands aside from altering the innate sequences. Modifications on adenosine (A) are the most widely characterized epitranscriptomic modification, including N⁶-methyladenosine (m⁶A), N¹-methyladenosine (m¹A), polyadenylation, and adenosine-to-inosine (A-to-I) RNA editing, and modifications on other nucleotides seem to be fewer, such as N⁷-methylguanosine (m⁷G), 5-methylcytosine (m⁵C), and pseudouridine (Ψ). These changes on the RNA strand surface, exclusively by their RNA-modifying proteins (RMPs), are reported in various biological phenomena, including programmed cell death (PCD). One necro-biological phenomenon that has been observed for long but has started to gain heed in recent years is “ferroptosis.” The phospholipid peroxidation by polyunsaturated-fatty-acid-containing-phospholipid hydroperoxyl (PLOOH) radicals destroys membrane integrity due to a series of mechanisms. The Fenton reaction, constituting the final Haber–Weiss reaction that is less recognized, collaboratively leading to the conversion of polyunsaturated fatty acid (PUFA) to PLOOH, is the etymological origin of ferroptosis. However, it is with increasing evidence that ferroptotic signaling is also intervened by epitranscriptomic modifications, although the truth is still ambiguous. We attempted to delineate some up-to-date discoveries on both epitranscriptomics and ferroptosis, bringing up the fundamentals to address any potential connection between the two. Next, we discussed whether a duologal relationship, or more, exists between the two, taking the ROS level and iron status into consideration. Lastly, we surveyed future perspectives that would favor the understanding of these topics.

Programmable RNA targeting by bacterial Argonaute nucleases with unconventional guide binding and cleavage specificity

Argonaute proteins are programmable nucleases that have defense and regulatory functions in both eukaryotes and prokaryotes. All known prokaryotic Argonautes (pAgos) characterized so far act on DNA targets. Here, we describe a new class of pAgos that uniquely use DNA guides to process RNA targets. The biochemical and structural analysis of Pseudooceanicola lipolyticus pAgo (PliAgo) reveals an unusual organization of the guide binding pocket that does not rely on divalent cations and the canonical set of contacts for 5'-end interactions. Unconventional interactions of PliAgo with the 5'-phosphate of guide DNA define its new position within pAgo and shift the site of target RNA cleavage in comparison with known Argonautes. The specificity for RNA over DNA is defined by ribonucleotide residues at the cleavage site. The analysed pAgos sense mismatches and modifications in the RNA target. The results broaden our understanding of prokaryotic defense systems and extend the spectrum of programmable nucleases with potential use in RNA technology.

Atomistic Simulation of Molecules Interacting with Biological Nanopores: From Current Understanding to Future Directions

Biological nanopores have been at the focus of numerous studies due to their role in many biological processes as well as their (prospective) technological applications. Among many other topics, recent studies on nanopores have addressed two key areas: antibiotic permeation through bacterial channels and sensing of analytes. Although the two areas are quite far apart in terms of their objectives, in both cases atomistic simulations attempt to understand the solute dynamics and the solute-protein interactions within the channel lumen. While decades of studies on various channels have culminated in an improved understanding of the key molecular factors and led to practical applications in some cases, successful utilization is limited. In this Perspective we summarize recent progress in understanding key issues in molecular simulations of antibiotic translocation and in the development of nanopore sensors. Moreover, we comment on possible advancements in computational algorithms that can potentially resolve some of the issues.

Epitranscriptomics in fibroblasts and fibrosis

Fibroblasts play an important role in the pathogenic mechanisms of several socially significant diseases, including pulmonary and cardiovascular fibrosis, liver cirrhosis, systemic sclerosis, progressive kidney disease. The alterations of the epitranscriptome, including more than 170 distinct post-transcriptional RNA modifications or editing events, justified their investigation as an important modulator of fibrosis. Recent development of high-throughput methods allows the identification of RNA modification sites and their mechanistic aspect in the fibrosis development. The most common RNA modification is methylation of N⁶-adenosine deposited by the m⁶A methyltransferase complex (METTL3/14/16, WTAP, KIAA1429, and RBM15/15B), erased by demethylases (FTO and ALKBH5), and recognized by binding proteins (e.g., YTHDF1/2/3, YTHDC1/2, IGF2BP1/2/3, etc.). Adenosine to inosine (A-to-I) RNA editing is another abundant editing event converting adenosine to inosine in double-stranded RNA regions through the action of the adenosine deaminase (ADAR) proteins. Last, but not least, 5-methylcytosine (m⁵C) regulates the stability and translation of mRNAs. All those RNA modifications have been observed in mRNA as well as the non-coding regions of pre-mRNA and ncRNAs, and demonstrate to be involved in fibrosis in different cellular and animal models. This Mini-Review focuses on the latest research on epitranscriptomic marks related to fibroblast biology and fibrosis as well as elucidates the future research directions in this context.

Epitranscriptomics of SARS-CoV-2 Infection

Recent studies on the epitranscriptomic code of SARS-CoV-2 infection have discovered various RNA modifications, such as N6-methyladenosine (m6A), pseudouridine (Ψ), and 2′-O-methylation (Nm). The effects of RNA methylation on SARS-CoV-2 replication and the enzymes involved in this mechanism are emerging. In this review, we summarize the advances in this emerging field and discuss the role of various players such as readers, writers, and erasers in m6A RNA methylation, the role of pseudouridine synthase one and seven in epitranscriptomic modification Ψ, an isomer of uridine, and role of nsp16/nsp10 heterodimer in 2′-O-methylation of the ribose sugar of the first nucleotide of SARS-CoV-2 mRNA. We also discuss RNA expression levels of various enzymes involved in RNA modifications in blood cells of SARS-CoV-2 infected individuals and their impact on host mRNA modification. In conclusion, these observations will facilitate the development of novel strategies and therapeutics for targeting RNA modification of SARS-CoV-2 RNA to control SARS-CoV-2 infection.

The Chloroplast Epitranscriptome: Factors, Sites, Regulation, and Detection Methods

Modifications in nucleic acids are present in all three domains of life. More than 170 distinct chemical modifications have been reported in cellular RNAs to date. Collectively termed as epitranscriptome, these RNA modifications are often dynamic and involve distinct regulatory proteins that install, remove, and interpret these marks in a site-specific manner. Covalent nucleotide modifications-such as methylations at diverse positions in the bases, polyuridylation, and pseudouridylation and many others impact various events in the lifecycle of an RNA such as folding, localization, processing, stability, ribosome assembly, and translational processes and are thus crucial regulators of the RNA metabolism. In plants, the nuclear/cytoplasmic epitranscriptome plays important roles in a wide range of biological processes, such as organ development, viral infection, and physiological means. Notably, recent transcriptome-wide analyses have also revealed novel dynamic modifications not only in plant nuclear/cytoplasmic RNAs related to photosynthesis but especially in chloroplast mRNAs, suggesting important and hitherto undefined regulatory steps in plastid functions and gene expression. Here we report on the latest findings of known plastid RNA modifications and highlight their relevance for the post-transcriptional regulation of chloroplast gene expression and their role in controlling plant development, stress reactions, and acclimation processes.