The RNA modification landscape in human disease

RNA modification: mechanisms and therapeutic targets

RNA modifications are dynamic and reversible chemical modifications on substrate RNA that are regulated by specific modifying enzymes. They play important roles in the regulation of many biological processes in various diseases, such as the development of cancer and other diseases. With the help of advanced sequencing technologies, the role of RNA modifications has caught increasing attention in human diseases in scientific research. In this review, we briefly summarized the basic mechanisms of several common RNA modifications, including m6A, m5C, m1A, m7G, Ψ, A-to-I editing and ac4C. Importantly, we discussed their potential functions in human diseases, including cancer, neurological disorders, cardiovascular diseases, metabolic diseases, genetic and developmental diseases, as well as immune disorders. Through the "writing-erasing-reading" mechanisms, RNA modifications regulate the stability, translation, and localization of pivotal disease-related mRNAs to manipulate disease development. Moreover, we also highlighted in this review all currently available RNA-modifier-targeting small molecular inhibitors or activators, most of which are designed against m6A-related enzymes, such as METTL3, FTO and ALKBH5. This review provides clues for potential clinical therapy as well as future study directions in the RNA modification field. More in-depth studies on RNA modifications, their roles in human diseases and further development of their inhibitors or activators are needed for a thorough understanding of epitranscriptomics as well as diagnosis, treatment, and prognosis of human diseases.

Transposable elements as essential elements in the control of gene expression

Interspersed repetitions called transposable elements (TEs), commonly referred to as mobile elements, make up a significant portion of the genomes of higher animals. TEs contribute in controlling the expression of genes locally and even far away at the transcriptional and post-transcriptional levels, which is one of their significant functional effects on gene function and genome evolution. There are different mechanisms through which TEs control the expression of genes. First, TEs offer cis-regulatory regions in the genome with their inherent regulatory features for their own expression, making them potential factors for controlling the expression of the host genes. Promoter and enhancer elements contain cis-regulatory sites generated from TE, which function as binding sites for a variety of trans-acting factors. Second, a significant portion of miRNAs and long non-coding RNAs (lncRNAs) have been shown to have TEs that encode for regulatory RNAs, revealing the TE origin of these RNAs. Furthermore, it was shown that TE sequences are essential for these RNAs' regulatory actions, which include binding to the target mRNA. By being a member of cis-regulatory and regulatory RNA sequences, TEs therefore play essential regulatory roles. Additionally, it has been suggested that TE-derived regulatory RNAs and cis-regulatory regions both contribute to the evolutionary novelty of gene regulation. Additionally, these regulatory systems arising from TE frequently have tissue-specific functions. The objective of this review is to discuss TE-mediated gene regulation, with a particular emphasis on the processes, contributions of various TE types, differential roles of various tissue types, based mostly on recent studies on humans.

N 2-methylguanosine modifications on human tRNAs and snRNA U6 are important for cell proliferation, protein translation and pre-mRNA splicing

Modified nucleotides in non-coding RNAs, such as tRNAs and snRNAs, represent an important layer of gene expression regulation through their ability to fine-tune mRNA maturation and translation. Dysregulation of such modifications and the enzymes installing them have been linked to various human pathologies including neurodevelopmental disorders and cancers. Several methyltransferases (MTases) are regulated allosterically by human TRMT112 (Trm112 in Saccharomyces cerevisiae), but the interactome of this regulator and targets of its interacting MTases remain incompletely characterized. Here, we have investigated the interaction network of human TRMT112 in intact cells and identify three poorly characterized putative MTases (TRMT11, THUMPD3 and THUMPD2) as direct partners. We demonstrate that these three proteins are active N2-methylguanosine (m2G) MTases and that TRMT11 and THUMPD3 methylate positions 10 and 6 of tRNAs, respectively. For THUMPD2, we discovered that it directly associates with the U6 snRNA, a core component of the catalytic spliceosome, and is required for the formation of m2G, the last 'orphan' modification in U6 snRNA. Furthermore, our data reveal the combined importance of TRMT11 and THUMPD3 for optimal protein synthesis and cell proliferation as well as a role for THUMPD2 in fine-tuning pre-mRNA splicing.

End-to-end RT-PCR of long RNA and highly structured RNA

RNA molecules play important roles in numerous normal cellular processes and disease states, from protein coding to gene regulation. RT-PCR, applying the power of polymerase chain reaction (PCR) to RNA by coupling reverse transcription with PCR, is one of the most important techniques to characterize RNA transcripts and monitor gene expression. The ability to analyze full-length RNA transcripts and detect their expression is critical to decipher their biological functions. However, due to the low processivity of retroviral reverse transcriptases (RTs), we can only monitor a small fraction of long RNA transcripts, especially those containing stable secondary and tertiary structures. The full-length sequences can only be deduced by computational analysis, which is often misleading. Group II intron-encoded RTs are a new type of RT enzymes. They have evolved specialized structural elements that unwind template structures and maintain close contact with the RNA template. Therefore, group II intron-encoded RTs are more processive than the retroviral RTs. The discovery, optimization and deployment of processive group II intron RTs provide us the opportunity to analyze RNA transcripts with single molecule resolution. MarathonRT, the most processive group II intron RT, has been extensively optimized for processive reverse transcription. In this chapter, we use MarathonRT to deliver a general protocol for long amplicon generation by RT-PCR, and also provide guidance for troubleshooting and further optimization.

DeepASDPred: a CNN-LSTM-based deep learning method for Autism spectrum disorders risk RNA identification

BACKGROUND

Autism spectrum disorders (ASD) are a group of neurodevelopmental disorders characterized by difficulty communicating with society and others, behavioral difficulties, and a brain that processes information differently than normal. Genetics has a strong impact on ASD associated with early onset and distinctive signs. Currently, all known ASD risk genes are able to encode proteins, and some de novo mutations disrupting protein-coding genes have been demonstrated to cause ASD. Next-generation sequencing technology enables high-throughput identification of ASD risk RNAs. However, these efforts are time-consuming and expensive, so an efficient computational model for ASD risk gene prediction is necessary.

RESULTS

In this study, we propose DeepASDPerd, a predictor for ASD risk RNA based on deep learning. Firstly, we use K-mer to feature encode the RNA transcript sequences, and then fuse them with corresponding gene expression values to construct a feature matrix. After combining chi-square test and logistic regression to select the best feature subset, we input them into a binary classification prediction model constructed by convolutional neural network and long short-term memory for training and classification. The results of the tenfold cross-validation proved our method outperformed the state-of-the-art methods. Dataset and source code are available at https://github.com/Onebear-X/DeepASDPred is freely available.

CONCLUSIONS

Our experimental results show that DeepASDPred has outstanding performance in identifying ASD risk RNA genes.

Biological roles of A-to-I editing: implications in innate immunity, cell death, and cancer immunotherapy

Adenosine-to-inosine (A-to-I) editing, a key RNA modification widely found in eukaryotes, is catalyzed by adenosine deaminases acting on RNA (ADARs). Such RNA editing destabilizes endogenous dsRNAs, which are subsequently recognized by the sensors of innate immune and other proteins as autologous dsRNAs. This prevents the activation of innate immunity and type I interferon-mediated responses, thereby reducing the downstream cell death induced by the activation of the innate immune sensing system. ADARs-mediated editing can also occur in mRNAs and non-coding RNAs (ncRNAs) in different species. In mRNAs, A-to-I editing may lead to missense mutations and the selective splicing of coding regions. Meanwhile, in ncRNAs, A-to-I editing may affect targeting and disrupt ncRNAs maturation, leading to anomalous cell proliferation, invasion, and responses to immunotherapy. This review highlights the biological functions of A-to-I editing, its role in regulating innate immunity and cell death, and its potential molecular significance in tumorigenesis and cancer targeted therapy and immunotherapy.

Small Nucleolar (Sno)RNA: Therapy Lays in Translation

The ribosome is one of the largest complexes in the cell. Adding to its complexity are more than 200 RNA modification sites present on ribosomal RNAs (rRNAs) in a single human ribosome. These modifications occur in functionally important regions of the rRNA molecule, and they are vital for ribosome function and proper gene expression. Until recent technological advancements, the study of rRNA modifications and their profiles has been extremely laborious, leaving many questions unanswered. Small nucleolar RNAs (snoRNAs) are non-coding RNAs that facilitate and dictate the specificity of rRNA modification deposition, making them an attractive target for ribosome modulation. Here, we propose that through the mapping of rRNA modification profiles, we can identify cell-specific modifications with high therapeutic potential. We also describe the challenges of achieving the targeting specificity needed to implement snoRNAs as therapeutic targets in cancers.

Covalent S-Adenosylhomocysteine-Based DNA Methyltransferase 2 Inhibitors with a New Type of Aryl Warhead

The DNA methyltransferase 2 (DNMT2) is an RNA modifying enzyme associated with pathophysiological processes, such as mental and metabolic disorders or cancer. Although the development of methyltransferase inhibitors remains challenging, DNMT2 is not only a promising target for drug discovery, but also for the development of activity-based probes. Here, we present covalent SAH-based DNMT2 inhibitors decorated with a new type of aryl warhead. Based on a noncovalent DNMT2 inhibitor with N-benzyl substituent, the Topliss scheme was followed for optimization. The results showed that electron-deficient benzyl moieties highly increased affinity. By decorating the structures with strong electron-withdrawing moieties and leaving groups, we adjusted the electrophilicity to create covalent DNMT2 inhibitors. A 4-bromo-3-nitrophenylsulfonamide-decorated SAH derivative (80) turned out to be the most potent (IC50 = 1.2 ± 0.1 μM) and selective inhibitor. Protein mass spectrometry confirmed the covalent reaction with the catalytically active cysteine-79.

circRNA: a promising all-around star in the future

circRNA is a type of RNA molecule with a circular structure. From initially thinking it was useless to now discovering that it has important and complex functions, the importance of circRNA is increasingly being recognized. circRNA was first discovered as an ncRNA with a molecular sponge function. Clinically, due to its special molecular structure, researchers are generally interested in its potential as a biomarker. Recently, circRNA has been proven to have many functions other than encoding proteins. In the clinical setting, circRNA also has strong potential for application in vaccine preparation and targeted therapy. This article discusses the synthesis of circRNA, introduces its functions and discusses its future development prospects.

Vital roles of m5C RNA modification in cancer and immune cell biology

RNA modification plays an important role in epigenetics at the posttranscriptional level, and 5-methylcytosine (m5C) has attracted increasing attention in recent years due to the improvement in RNA m5C site detection methods. By influencing transcription, transportation and translation, m5C modification of mRNA, tRNA, rRNA, lncRNA and other RNAs has been proven to affect gene expression and metabolism and is associated with a wide range of diseases, including malignant cancers. RNA m5C modifications also substantially impact the tumor microenvironment (TME) by targeting different groups of immune cells, including B cells, T cells, macrophages, granulocytes, NK cells, dendritic cells and mast cells. Alterations in immune cell expression, infiltration and activation are highly linked to tumor malignancy and patient prognosis. This review provides a novel and holistic examination of m5C-mediated cancer development by examining the exact mechanisms underlying the oncogenicity of m5C RNA modification and summarizing the biological effects of m5C RNA modification on tumor cells as well as immune cells. Understanding methylation-related tumorigenesis can provide useful insights for the diagnosis as well as the treatment of cancer.

Identification of m7G-Related miRNA Signatures Associated with Prognosis, Oxidative Stress, and Immune Landscape in Lung Adenocarcinoma

The role of N7-methylguanosine(m7G)-related miRNAs in lung adenocarcinoma (LUAD) remains unclear. We used LUAD data from The Cancer Genome Atlas (TCGA) to establish a risk model based on the m7G-related miRNAs, and divided patients into high-risk or low-risk subgroups. A nomogram for predicting overall survival (OS) was then constructed based on the independent risk factors. In addition, we performed a functional enrichment analysis and defined the oxidative stress-related genes, immune landscape as well as a drug response profile in the high-risk and low-risk subgroups. This study incorporated 28 m7G-related miRNAs into the risk model. The data showed a significant difference in the OS between the high-risk and low-risk subgroups. The receiver operating characteristic curve (ROC) predicted that the area under the curve (AUC) of one-year, three-year and five-year OS was 0.781, 0.804 and 0.853, respectively. The C-index of the prognostic nomogram for predicting OS was 0.739. We then analyzed the oxidative stress-related genes and immune landscape in the high-risk and low-risk subgroups. The data demonstrated significant differences in the expression of albumin (ALB), estimated score, immune score, stromal score, immune cell infiltration and functions between the high-risk and low-risk subgroups. In addition, the drug response analysis showed that low-risk subgroups may be more sensitive to tyrosine kinase inhibitor (TKI) and histone deacetylase (HDAC) inhibitors. We successfully developed a novel risk model based on m7G-related miRNAs in this study. The model can predict clinical prognosis and guide therapeutic regimens in patients with LUAD. Our data also provided new insights into the molecular mechanisms of m7G in LUAD.

MRT-ModSeq - Rapid detection of RNA modifications with MarathonRT

Chemical modifications are essential regulatory elements that modulate the behavior and function of cellular RNAs. Despite recent advances in sequencing-based RNA modification mapping, methods combining accuracy and speed are still lacking. Here, we introduce MRT- ModSeq for rapid, simultaneous detection of multiple RNA modifications using MarathonRT. MRT-ModSeq employs distinct divalent cofactors to generate 2-D mutational profiles that are highly dependent on nucleotide identity and modification type. As a proof of concept, we use the MRT fingerprints of well-studied rRNAs to implement a general workflow for detecting RNA modifications. MRT-ModSeq rapidly detects positions of diverse modifications across a RNA transcript, enabling assignment of m1acp3Y, m1A, m3U, m7G and 2'-OMe locations through mutation-rate filtering and machine learning. m1A sites in sparsely modified targets, such as MALAT1 and PRUNE1 could also be detected. MRT-ModSeq can be trained on natural and synthetic transcripts to expedite detection of diverse RNA modification subtypes across targets of interest.

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.

Enzymatic synthesis of RNA standards for mapping and quantifying RNA modifications in sequencing analysis

Synthesis of RNA standards that contain an internal site-specific modification is important for mapping and quantification of the modified nucleotide in sequencing analysis. While RNA containing a site-specific modification can be readily synthesized by solid-state coupling for less than 100-mer nucleotides, longer RNA must be synthesized by enzymatic ligation in the presence of a DNA splint. However, long RNAs have structural heterogeneity, and those generated by in vitro transcription have 3'-end sequence heterogeneity, which together substantially reduce the yield of ligation. Here we describe a method of 3-part splint ligation that joins an in vitro transcribed left-arm RNA, an in vitro transcribed right-arm RNA, and a chemically synthesized modification-containing middle RNA, with an efficiency higher than previously reported. We report that the improved efficiency is largely attributed to the inclusion of a pair of DNA disruptors proximal to the ligation sites, and to a lesser extent to the homogeneous processing of the 3'-end of the left-arm RNA. The yields of the ligated long RNA are sufficiently high to afford purification to homogeneity for practical RNA research. We also verify the sequence accuracy at each ligation junction by nanopore sequencing.

Methylated guanosine and uridine modifications in S. cerevisiae mRNAs modulate translation elongation

Chemical modifications to protein encoding messenger RNAs (mRNAs) influence their localization, translation, and stability within cells. Over 15 different types of mRNA modifications have been observed by sequencing and liquid chromatography coupled to tandem mass spectrometry (LC-MS/MS) approaches. While LC-MS/MS is arguably the most essential tool available for studying analogous protein post-translational modifications, the high-throughput discovery and quantitative characterization of mRNA modifications by LC-MS/MS has been hampered by the difficulty of obtaining sufficient quantities of pure mRNA and limited sensitivities for modified nucleosides. We have overcome these challenges by improving the mRNA purification and LC-MS/MS pipelines. The methodologies we developed result in no detectable non-coding RNA modifications signals in our purified mRNA samples, quantify 50 ribonucleosides in a single analysis, and provide the lowest limit of detection reported for ribonucleoside modification LC-MS/MS analyses. These advancements enabled the detection and quantification of 13 S. cerevisiae mRNA ribonucleoside modifications and reveal the presence of four new S. cerevisiae mRNA modifications at low to moderate levels (1-methyguanosine, N2-methylguanosine, N2,N2-dimethylguanosine, and 5-methyluridine). We identified four enzymes that incorporate these modifications into S. cerevisiae mRNAs (Trm10, Trm11, Trm1, and Trm2, respectively), though our results suggest that guanosine and uridine nucleobases are also non-enzymatically methylated at low levels. Regardless of whether they are incorporated in a programmed manner or as the result of RNA damage, we reasoned that the ribosome will encounter the modifications that we detect in cells. To evaluate this possibility, we used a reconstituted translation system to investigate the consequences of modifications on translation elongation. Our findings demonstrate that the introduction of 1-methyguanosine, N2-methylguanosine and 5-methyluridine into mRNA codons impedes amino acid addition in a position dependent manner. This work expands the repertoire of nucleoside modifications that the ribosome must decode in S. cerevisiae. Additionally, it highlights the challenge of predicting the effect of discrete modified mRNA sites on translation de novo because individual modifications influence translation differently depending on mRNA sequence context.

WDR4 gene polymorphisms and Wilms tumor susceptibility in Chinese children: A five-center case-control study

Wilms tumor is the most common embryonal renal malignancy in children. WDR4 is an indispensable noncatalytic subunit of the RNA N7-methylguanosine (m7G) methyltransferase complex and plays an essential role in tumorigenesis. However, the relationship between polymorphisms in the WDR4 gene and susceptibility to Wilms tumor remains to be fully investigated. We performed a large case-control study involving 414 patients and 1199 cancer-free controls to investigate whether single nucleotide polymorphisms (SNPs) in the WDR4 gene are associated with Wilms tumor susceptibility. WDR4 gene polymorphisms (rs2156315 C > T, rs2156316 C > G, rs6586250 C > T, rs15736 G > A, and rs2248490 C > G) were genotyped using the TaqMan assay. In addition, unconditioned logistic regression analysis was performed, odds ratios (ORs) and 95% confidence intervals (CIs) were used to assess the association between WDR4 gene SNPs and Wilms tumor susceptibility as well as the strength of the associations. We found that only the rs6586250 C>T polymorphism was significantly associated with an increased risk of Wilms tumor (adjusted OR=2.99, 95% CI = 1.28-6.97, P = 0.011 for the rs6586250 TT genotype; adjusted OR=3.08, 95% CI = 1.33-7.17, P = 0.009 for the rs6586250 CC/CT genotype). Furthermore, the stratification analysis revealed that patients with the rs6586250 TT genotype and carriers with 1-5 risk genotypes exhibited statistically significant associations with increased Wilms tumor risk in specific subgroups. However, the rs2156315 CT/TT genotype was identified as having a protective effect against Wilms tumor in the age >18 months subgroup compared with the rs2156315 CC genotype. In brief, our study demonstrated that the rs6586250 C > T polymorphism of the WDR4 gene was significantly associated with Wilms tumor. This finding may contribute to the understanding of the genetic mechanism of Wilms tumor.

m6Aminer: Predicting the m6Am Sites on mRNA by Fusing Multiple Sequence-Derived Features into a CatBoost-Based Classifier

As one of the most important post-transcriptional modifications, m6Am plays a fairly important role in conferring mRNA stability and in the progression of cancers. The accurate identification of the m6Am sites is critical for explaining its biological significance and developing its application in the medical field. However, conventional experimental approaches are time-consuming and expensive, making them unsuitable for the large-scale identification of the m6Am sites. To address this challenge, we exploit a CatBoost-based method, m6Aminer, to identify the m6Am sites on mRNA. For feature extraction, nine different feature-encoding schemes (pseudo electron-ion interaction potential, hash decimal conversion method, dinucleotide binary encoding, nucleotide chemical properties, pseudo k-tuple composition, dinucleotide numerical mapping, K monomeric units, series correlation pseudo trinucleotide composition, and K-spaced nucleotide pair frequency) were utilized to form the initial feature space. To obtain the optimized feature subset, the ExtraTreesClassifier algorithm was adopted to perform feature importance ranking, and the top 300 features were selected as the optimal feature subset. With different performance assessment methods, 10-fold cross-validation and independent test, m6Aminer achieved average AUC of 0.913 and 0.754, demonstrating a competitive performance with the state-of-the-art models m6AmPred (0.905 and 0.735) and DLm6Am (0.897 and 0.730). The prediction model developed in this study can be used to identify the m6Am sites in the whole transcriptome, laying a foundation for the functional research of m6Am.

m5U-SVM: identification of RNA 5-methyluridine modification sites based on multi-view features of physicochemical features and distributed representation

BACKGROUND

RNA 5-methyluridine (m5U) modifications are obtained by methylation at the C₅ position of uridine catalyzed by pyrimidine methylation transferase, which is related to the development of human diseases. Accurate identification of m5U modification sites from RNA sequences can contribute to the understanding of their biological functions and the pathogenesis of related diseases. Compared to traditional experimental methods, computational methods developed based on machine learning with ease of use can identify modification sites from RNA sequences in an efficient and time-saving manner. Despite the good performance of these computational methods, there are some drawbacks and limitations.

RESULTS

In this study, we have developed a novel predictor, m5U-SVM, based on multi-view features and machine learning algorithms to construct predictive models for identifying m5U modification sites from RNA sequences. In this method, we used four traditional physicochemical features and distributed representation features. The optimized multi-view features were obtained from the four fused traditional physicochemical features by using the two-step LightGBM and IFS methods, and then the distributed representation features were fused with the optimized physicochemical features to obtain the new multi-view features. The best performing classifier, support vector machine, was identified by screening different machine learning algorithms. Compared with the results, the performance of the proposed model is better than that of the existing state-of-the-art tool.

CONCLUSIONS

m5U-SVM provides an effective tool that successfully captures sequence-related attributes of modifications and can accurately predict m5U modification sites from RNA sequences. The identification of m5U modification sites helps to understand and delve into the related biological processes and functions.

Biosensors integrated 3D organoid/organ-on-a-chip system: A real-time biomechanical, biophysical, and biochemical monitoring and characterization

As a full-fidelity simulation of human cells, tissues, organs, and even systems at the microscopic scale, Organ-on-a-Chip (OOC) has significant ethical advantages and development potential compared to animal experiments. The need for the design of new drug high-throughput screening platforms and the mechanistic study of human tissues/organs under pathological conditions, the evolving advances in 3D cell biology and engineering, etc., have promoted the updating of technologies in this field, such as the iteration of chip materials and 3D printing, which in turn facilitate the connection of complex multi-organs-on-chips for simulation and the further development of technology-composite new drug high-throughput screening platforms. As the most critical part of organ-on-a-chip design and practical application, verifying the success of organ model modeling, i.e., evaluating various biochemical and physical parameters in OOC devices, is crucial. Therefore, this paper provides a logical and comprehensive review and discussion of the advances in organ-on-a-chip detection and evaluation technologies from a broad perspective, covering the directions of tissue engineering scaffolds, microenvironment, single/multi-organ function, and stimulus-based evaluation, and provides a more comprehensive review of the progress in the significant organ-on-a-chip research areas in the physiological state.

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

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

Construction of a Diagnostic m7G Regulator-Mediated Scoring Model for Identifying the Characteristics and Immune Landscapes of Osteoarthritis

With the increasingly serious burden of osteoarthritis (OA) on modern society, it is urgent to propose novel diagnostic biomarkers and differentiation models for OA. 7-methylguanosine (m7G), as one of the most common base modification forms in post transcriptional regulation, through which the seventh position N of guanine (G) of messenger RNA is modified by methyl under the action of methyltransferase; it has been found that it plays a crucial role in different diseases. Therefore, we explored the relationship between OA and m7G. Based on the expression level of 18 m7G-related regulators, we identified nine significant regulators. Then, via a series of methods of machine learning, such as support vector machine recursive feature elimination, random forest and lasso-cox regression analysis, a total of four significant regulators were further identified (DCP2, EIF4E2, LARP1 and SNUPN). Additionally, according to the expression level of the above four regulators, two different m7G-related clusters were divided via consensus cluster analysis. Furthermore, via immune infiltration, differential expression analysis and enrichment analysis, we explored the characteristic of the above two different clusters. An m7G-related scoring model was constructed via the PCA algorithm. Meanwhile, there was a different immune status and correlation for immune checkpoint inhibitors between the above two clusters. The expression difference of the above four regulators was verified via real-time quantitative polymerase chain reaction. Overall, a total of four biomarkers were identified and two different m7G-related subsets of OA with different immune microenvironment were obtained. Meanwhile, the construction of m7G-related Scoring model may provide some new strategies and insights for the therapy and diagnosis of OA patients.

Expression patterns of eight RNA-modified regulators correlating with immune infiltrates during the progression of osteoarthritis

Background

RNA modifications in eukaryotic cells have emerged as an exciting but under-explored area in recent years and are considered to be associated with many human diseases. While several studies have been published relating to m6A in osteoarthritis (OA), we only have limited knowledge of other kinds of RNA modifications. Our study investigated eight RNA modifiers' specific roles in OA including A-to-I, APA, m5C, m6A, m7G, mcm5s2U, Nm and Ψ together with their relationship with immune infiltration.

Methods

RNA modification patterns in OA samples were identified based on eight-type RNA modifiers and their correlation with the degree of immune infiltration was also methodically investigated. Receiver operating characteristic curves (ROC) and qRT-PCR was performed to confirm the abnormal expression of hub genes. The RNA modification score (Rmscore) was generated by the applications of principal component analysis (PCA) algorithm in order to quantify RNA modification modes in individual OA patients.

Results

We identified 21 differentially-expressed RNA modification related genes between OA and healthy samples. For example, CFI, CBLL1 and ALKBH8 were expressed at high levels in OA (P<0.001), while RPUSD4, PUS1, NUDT21, FBL and WDR4 were expressed at low levels (P<0.001). Two candidate RNA modification regulators (WDR4 and CFI) were screened out utilizing a random forest machine learning model. We then identified two distinctive RNA modification modes in OA which were found to display distinctive biological features. High Rmscore, characterized by increased immune cell infiltration, indicated an inflamed phenotype.

Conclusions

Our study was the first to systematically reveal the crosstalk and dysregulations eight-type of RNA modifications in OA. Assessing individuals' RNA modification patterns will be conductive to enhance our understanding of the properties of immune infiltration, provide novel diagnostic and prognostic biomarkers, and guide more effective immunotherapy strategies in the future.

Contemporary Progress and Opportunities in RNA-Targeted Drug Discovery

The surprising discovery that RNAs are the predominant gene products to emerge from the human genome catalyzed a renaissance in RNA biology. It is now well-understood that RNAs act as more than just a messenger and comprise a large and diverse family of ribonucleic acids of differing sizes, structures, and functions. RNAs play expansive roles in the cell, contributing to the regulation and fine-tuning of nearly all aspects of gene expression and genome architecture. In line with the significance of these functions, we have witnessed an explosion in discoveries connecting RNAs with a variety of human diseases. Consequently, the targeting of RNAs, and more broadly RNA biology, has emerged as an untapped area of drug discovery, making the search for RNA-targeted therapeutics of great interest. In this Microperspective, I highlight contemporary learnings in the field and present my views on how to catapult us toward the systematic discovery of RNA-targeted medicines.

The lncRNA epigenetics: The significance of m6A and m5C lncRNA modifications in cancer

Most of our transcribed RNAs are represented by non-coding sequences. Long non-coding RNAs (lncRNAs) are transcripts with no or very limited protein coding ability and a length &gt;200nt. They can be epigenetically modified. N6-methyladenosine (m6A), N1-methyladenosine (m1A), 5-methylcytosine (m5C), 7-methylguanosine (m7G) and 2’-O-methylation (Nm) are some of the lncRNAs epigenetic modifications. The epigenetic modifications of RNA are controlled by three classes of enzymes, each playing a role in a specific phase of the modification. These enzymes are defined as “writers”, “readers” and “erasers”. m6A and m5C are the most studied epigenetic modifications in RNA. These modifications alter the structure and properties, thus modulating the functions and interactions of lncRNAs. The aberrant expression of several lncRNAs is linked to the development of a variety of cancers and the epigenetic signatures of m6A- or m5C-related lncRNAs are increasingly recognized as potential biomarkers of prognosis, predictors of disease stage and overall survival. In the present manuscript, the most up to date literature is reviewed with the focus on m6A and m5C modifications of lncRNAs and their significance in cancer.

Current and potential roles of RNA modification-mediated autophagy dysregulation in cancer

Autophagy, a cellular lysosomal degradation and survival pathway, supports nutrient recycling and adaptation to metabolic stress and participates in various stages of tumor development, including tumorigenesis, metastasis, and malignant state maintenance. Among the various factors contributing to the dysregulation of autophagy in cancer, RNA modification can regulate autophagy by directly affecting the expression of core autophagy proteins. We propose that autophagy disorder mediated by RNA modification is an important mechanism for cancer development. Therefore, this review mainly discusses the role of RNA modification-mediated autophagy regulation in tumorigenesis. We summarize the molecular basis of autophagy and the core proteins and complexes at different stages of autophagy, especially those involved in cancer development. Moreover, we describe the crosstalk of RNA modification and autophagy and review the recent advances and potential role of the RNA modification/autophagy axis in the development of multiple cancers. Furthermore, the dual role of the RNA modification/autophagy axis in cancer drug resistance is discussed. A comprehensive understanding and extensive exploration of the molecular crosstalk of RNA modifications with autophagy will provide important insights into tumor pathophysiology and provide more options for cancer therapeutic strategies.

RNA out of the mist

RNA has long been regarded primarily as the intermediate between genes and proteins. It was a surprise then to discover that eukaryotic genes are mosaics of mRNA sequences interrupted by large tracts of transcribed but untranslated sequences, and that multicellular organisms also express many long 'intergenic' and antisense noncoding RNAs (lncRNAs). The identification of small RNAs that regulate mRNA translation and half-life did not disturb the prevailing view that animals and plant genomes are full of evolutionary debris and that their development is mainly supervised by transcription factors. Gathering evidence to the contrary involved addressing the low conservation, expression, and genetic visibility of lncRNAs, demonstrating their cell-specific roles in cell and developmental biology, and their association with chromatin-modifying complexes and phase-separated domains. The emerging picture is that most lncRNAs are the products of genetic loci termed 'enhancers', which marshal generic effector proteins to their sites of action to control cell fate decisions during development.

The emergent role of mitochondrial RNA modifications in metabolic alterations

Mitochondrial epitranscriptomics refers to the modifications occurring in all the different RNA types of mitochondria. Although the number of mitochondrial RNA modifications is less than those in cytoplasm, substantial evidence indicates that they play a critical role in accurate protein synthesis. Recent evidence supported those modifications in mitochondrial RNAs also have crucial implications in mitochondrial-related diseases. In the light of current knowledge about the involvement, the association between mitochondrial RNA modifications and diseases arises from studies focusing on mutations in both mitochondrial and nuclear DNA genes encoding enzymes involved in such modifications. Here, we review the current evidence available for mitochondrial RNA modifications and their role in metabolic disorders, and we also explore the possibility of using them as promising targets for prevention and early detection. Finally, we discuss future directions of mitochondrial epitranscriptomics in these metabolic alterations, and how these RNA modifications may offer a new diagnostic and theragnostic avenue for preventive purposes. This article is categorized under: RNA Processing > RNA Editing and Modification.

The chemical landscape of the human ribosome at 1.67 Å resolution

The ability of ribosomes to translate the genetic code into protein requires a finely tuned ion and solvent ecosystem. However, the lack of high-resolution structures has precluded accurate positioning of all the functional elements of the ribosome and limited our understanding of the specific role of ribosomal RNA chemical modifications in modulating ribosome function in health and disease. Here, using a new sample preparation methodology based on functionalised pristine graphene-coated grids, we solve the cryo-EM structure of the human large ribosomal subunit to a resolution of 1.67 Å. The accurate assignment of water molecules, magnesium and potassium ions in our model highlights the fundamental biological role of ribosomal RNA methylation in harnessing unconventional carbon-oxygen hydrogen bonds to establish chemical interactions with the environment and fine-tune the functional interplay with tRNA. In addition, the structures of three translational inhibitors bound to the human large ribosomal subunit at better than 2 Å resolution provide mechanistic insights into how three key druggable pockets of the ribosome are targeted and illustrate the potential of this methodology to accelerate high-throughput structure-based design of anti-cancer therapeutics.

Expression pattern and prognostic value of key regulators for N7-methylguanosine RNA modification in prostate cancer

Alterations in the regulators of RNA methylation modifications, such as N7-methylguanosine (m7G), have been implicated in a variety of diseases. Therefore, the analysis and identification of disease-related m7G modification regulators will accelerate advances in understanding disease pathogenesis. However, the implications of alterations in the regulators of m7G modifications remain poorly understood in prostate adenocarcinoma. In the present study, we analyze the expression patterns of 29 m7G RNA modification regulators in prostate adenocarcinoma using The Cancer Genome Atlas (TCGA) and perform consistent clustering analysis of differentially expressed genes (DEGs). We find that 18 m7G-related genes are differentially expressed in tumor and normal tissues. In different cluster subgroups, DEGs are mainly enriched in tumorigenesis and tumor development. Furthermore, immune analyses demonstrate that patients in cluster 1 have significantly higher scores for stromal and immune cells, such as B cells, T cells, and macrophages. Then, a TCGA-related risk model is developed and successfully validated using a Gene Expression Omnibus external dataset. Two genes ( EIF4A1 and NCBP2) are determined to be prognostically significant. Most importantly, we construct tissue microarrays from 26 tumor specimens and 20 normal specimens, and further confirm that EIF4A1 and NCBP2 are associated with tumor progression and Gleason score. Therefore, we conclude that the m7G RNA methylation regulators may be involved in the poor prognosis of patients with prostate adenocarcinoma. The results of this study may provide support for exploring the underlying molecular mechanisms of m7G regulators, especially EIF4A1 and NCBP2.

Hierarchical DNA branch assembly-encoded fluorescent nanoladders for single-cell transcripts imaging

Spatial visualization of single-cell transcripts is limited by signal specificity and multiplexing. Here, we report hierarchical DNA branch assembly-encoded fluorescent nanoladders, which achieve denoised and highly multiplexed signal amplification for single-molecule transcript imaging. This method first offers independent RNA-primed rolling circle amplification without nonspecific amplification based on circular DNAzyme. It then executes programmable DNA branch assembly on these amplicons to encode virtual signals for visualizing numbers of targets by FISH. In theory, more virtual signals can be encoded via the increase of detection spectral channels and repeats of the same sequences on barcode. Our method almost eliminates nonspecific amplification in fixed cells (reducing nonspecific spots of single cells from 16 to nearly zero), and achieves simultaneous quantitation of nine transcripts by using only two detection spectral channels. We demonstrate accurate RNA profiling in different cancer cells, and reveal diverse localization patterns for spatial regulation of transcripts.

m5C-dependent cross-regulation between nuclear reader ALYREF and writer NSUN2 promotes urothelial bladder cancer malignancy through facilitating RABL6/TK1 mRNAs splicing and stabilization

The significance of 5-methylcytosine (m⁵C) methylation in human malignancies has become an increasing focus of investigation. Here, we show that m⁵C regulators including writers, readers and erasers, are predominantly upregulated in urothelial carcinoma of the bladder (UCB) derived from Sun Yat-sen University Cancer Center and The Cancer Genome Atlas cohort. In addition, NOP2/Sun RNA methyltransferase family member 2 (NSUN2) as a methyltransferase and Aly/REF export factor (ALYREF) as a nuclear m⁵C reader, are frequently coexpressed in UCB. By applying patient-derived organoids model and orthotopic xenograft mice model, we demonstrate that ALYREF enhances proliferation and invasion of UCB cells in an m⁵C-dependent manner. Integration of tanscriptome-wide RNA bisulphite sequencing (BisSeq), RNA-sequencing (RNA-seq) and RNA Immunoprecipitation (RIP)-seq analysis revealed that ALYREF specifically binds to hypermethylated m⁵C site in RAB, member RAS oncogene family like 6 (RABL6) and thymidine kinase 1 (TK1) mRNA via its K171 domain. ALYREF controls UCB malignancies through promoting hypermethylated RABL6 and TK1 mRNA for splicing and stabilization. Moreover, ALYREF recognizes hypermethylated m⁵C site of NSUN2, resulting in NSUN2 upregulation in UCB. Clinically, the patients with high coexpression of ALYREF/RABL6/TK1 axis had the poorest overall survival. Our study unveils an m⁵C dependent cross-regulation between nuclear reader ALYREF and m⁵C writer NSUN2 in activation of hypermethylated m⁵C oncogenic RNA through promoting splicing and maintaining stabilization, consequently leading to tumor progression, which provides profound insights into therapeutic strategy for UCB.

m6A Modification-Association with Oxidative Stress and Implications on Eye Diseases

Oxidative stress (OS) refers to a state of imbalance between oxidation and antioxidation. OS is considered to be an important factor leading to aging and a range of diseases. The eyes are highly oxygen-consuming organs. Due to its continuous exposure to ultraviolet light, the eye is particularly vulnerable to the impact of OS, leading to eye diseases such as corneal disease, cataracts, glaucoma, etc. The N6-methyladenosine (m6A) modification is the most investigated RNA post-transcriptional modification and participates in a variety of cellular biological processes. In this study, we review the role of m6A modification in oxidative stress-induced eye diseases and some therapeutic methods to provide a relatively overall understanding of m6A modification in oxidative stress-related eye diseases.

RNA Modification Detection Using Nanopore Direct RNA Sequencing and nanoDoc2

RNA modifications regulate multiple aspects of cellular function including RNA splicing, translation, export, decay, stability, and phase separation. One of the comprehensive ways to detect such modifications is by the recent advancement of direct RNA sequencing from Oxford Nanopore Technologies (ONT). However, this method obtains a large amount of data with high complexity in the form of raw current signal that poses a new informatics challenge to accurately detect those modifications. Here, we provide nanoDoc2, a software to detect multiple types of RNA modification from nanopore direct RNA sequencing data. The nanoDoc2 includes a novel signal segmentation algorithm based on the trace value-a base probability feature that is added by the Guppy basecalling program from ONT during processing of the raw signal. The core of nanoDoc2 includes a machine learning algorithm in which a 6-mer segmented raw current signal is analyzed by deep one-class classification using a WaveNet-based neural network. As an output, an RNA modification is detected by a statistical score in each candidate position. Herein, we describe the detailed instructions on how to use nanoDoc2 for signal segmentation, train/test the neural network, and finally predict RNA modifications present in nanopore direct RNA sequencing data.

Comprehensive analysis of transcriptomic profiling of 5-methylcytosin modification in placentas from preeclampsia and normotensive pregnancies

Increasing evidence suggests that RNA m5C modification and its regulators have been confirmed to be associated with the pathogenesis of many diseases. However, the distribution and biological functions of m5C in mRNAs of placental tissues remain unknown. we collected placentae from normotensive pregnancies (CTR) and preeclampsia patients (PE) to analyze the transcriptomic profiling of m5C RNA methylation through m5C RNA immunoprecipitation (UMI-MeRIP-Seq). we discovered that overall m5C methylation peaks were decreased in placental tissues from PE patients. And, 2844 aberrant m5C peaks were identified, of which respectively 1304 m5C peaks were upregulated and 1540 peaks were downregulated. The distribution of m5C peaks were mainly located in CDS (coding sequences) regions in placental tissues of both groups, but compared with the CTR group, the m5C peak in PE group before the stop code of CDS was significantly increased and even higher than the peak value after start code in CDS. Differentially methylated genes were mainly enriched in MAPK/cAMP signaling pathway. Moreover, the up-regulated genes with hypermethylated modification were enriched in the processes of hypoxia, inflammation/immune response. Finally, through analyzing the mRNA expression levels of m5C RNA methylation regulators, we found only DNMT3B and TET3 were significantly upregulated in PE samples than in control group. And they are not only negatively correlated with each other, but also closely related to those differentially expressed genes modified by differential methylation.Our findings provide new insights regarding alterations of m5C RNA modification into the pathogenic mechanisms of PE.

Self-attention enabled deep learning of dihydrouridine (D) modification on mRNAs unveiled a distinct sequence signature from tRNAs

Dihydrouridine (D) is a modified pyrimidine nucleotide universally found in viral, prokaryotic, and eukaryotic species. It serves as a metabolic modulator for various pathological conditions, and its elevated levels in tumors are associated with a series of cancers. Precise identification of D sites on RNA is vital for understanding its biological function. A number of computational approaches have been developed for predicting D sites on tRNAs; however, none have considered mRNAs. We present here DPred, the first computational tool for predicting D on mRNAs in yeast from the primary RNA sequences. Built on a local self-attention layer and a convolutional neural network (CNN) layer, the proposed deep learning model outperformed classic machine learning approaches (random forest, support vector machines, etc.) and achieved reasonable accuracy and reliability with areas under the curve of 0.9166 and 0.9027 in jackknife cross-validation and on an independent testing dataset, respectively. Importantly, we showed that distinct sequence signatures are associated with the D sites on mRNAs and tRNAs, implying potentially different formation mechanisms and putative divergent functionality of this modification on the two types of RNA. DPred is available as a user-friendly Web server.

The Repertoire of RNA Modifications Orchestrates a Plethora of Cellular Responses

Although a plethora of DNA modifications have been extensively investigated in the last decade, recent breakthroughs in molecular biology, including high throughput sequencing techniques, have enabled the identification of post-transcriptional marks that decorate RNAs; hence, epitranscriptomics has arisen. This recent scientific field aims to decode the regulatory layer of the transcriptome and set the ground for the detection of modifications in ribose nucleotides. Until now, more than 170 RNA modifications have been reported in diverse types of RNA that contribute to various biological processes, such as RNA biogenesis, stability, and transcriptional and translational accuracy. However, dysfunctions in the RNA-modifying enzymes that regulate their dynamic level can lead to human diseases and cancer. The present review aims to highlight the epitranscriptomic landscape in human RNAs and match the catalytic proteins with the deposition or deletion of a specific mark. In the current review, the most abundant RNA modifications, such as N6-methyladenosine (m⁶A), N5-methylcytosine (m⁵C), pseudouridine (Ψ) and inosine (I), are thoroughly described, their functional and regulatory roles are discussed and their contributions to cellular homeostasis are stated. Ultimately, the involvement of the RNA modifications and their writers, erasers, and readers in human diseases and cancer is also discussed.

Selective Chemical Labeling Strategy for Oligonucleotides Determination: A First Application to Full-Range Profiling of Transfer RNA Modifications

To date, the extremely high polarity and poor signal intensity of macromolecular nucleic acids are greatly impeding the progress of mass spectrometry technology in the quality control of nucleic acid drugs and the characterization of DNA oxidation and RNA modifications. We recently described a general N-(tert-butyldimethylsilyl)-N-methyl-trifluoroacetamide (MTBSTFA) labeling method for oligonucleotide determination and applied it to the full-range profiling of tRNA in vitro and in vivo studies for the first time. The primary advantages of this method include strong retention, no observable byproducts, predictable and easily interpreted MS² data, and the circumvention of instrument harmful reagents that were necessary in previous methods. Selective labeling of N-(tert-butyldimethylsilyl)-N-methyl-trifluoroacetamide to the terminal phosphate groups of oligonucleotides endows it broadly applicable for DNA/RNA profiling. Moreover, the improvement of sequence coverage was achieved in yeast tRNA^phe(GAA) analysis owing to this method's good detection capability of 1-12 nucleotides in length. We also extended this strategy to determine the abundance of modified bases and discover new modifications via digesting RNA into single-nucleotide products, promoting the comprehensive mapping of RNA. The easy availability of derivatization reagent and the simple, rapid one-step reaction render it easy to operate for researchers. When applied in characterizing tRNAs in HepG2 cells and rats with nonalcoholic fatty liver disease, a fragment of U[m¹G][m²G], specific for tRNA^Asn(QUU) in cells, was significantly upregulated, indicating a possible clue to nonalcoholic fatty liver disease pathogenesis.

Nucleos'ID: A New Search Engine Enabling the Untargeted Identification of RNA Post-transcriptional Modifications from Tandem Mass Spectrometry Analyses of Nucleosides

As RNA post-transcriptional modifications are of growing interest, several methods were developed for their characterization. One of them established for their identification, at the nucleosidic level, is the hyphenation of separation methods, such as liquid chromatography or capillary electrophoresis, to tandem mass spectrometry. However, to our knowledge, no software is yet available for the untargeted identification of RNA post-transcriptional modifications from MS/MS data-dependent acquisitions. Thus, very long and tedious manual data interpretations are required. To meet the need of easier and faster data interpretation, a new user-friendly search engine, called Nucleos'ID, was developed for CE-MS/MS and LC-MS/MS users. Performances of this new software were evaluated on CE-MS/MS data from nucleoside analyses of already well-described Saccharomyces cerevisiae transfer RNA and Bos taurus total tRNA extract. All samples showed great true positive, true negative, and false discovery rates considering the database size containing all modified and unmodified nucleosides referenced in the literature. The true positive and true negative rates obtained were above 0.94, while the false discovery rates were between 0.09 and 0.17. To increase the level of sample complexity, untargeted identification of several RNA modifications from Pseudomonas aeruginosa 70S ribosome was achieved by the Nucleos'ID search following CE-MS/MS analysis.

A five-pseudouridylation-associated-LncRNA classifier for primary prostate cancer prognosis prediction

Background: Prostate cancer (PCa) is one of the most common cancers in males around the globe, and about one-third of patients with localized PCa will experience biochemical recurrence (BCR) after radical prostatectomy or radiation therapy. Reportedly, a proportion of patients with BCR had a poor prognosis. Cumulative studies have shown that RNA modifications participate in the cancer-related transcriptome, but the role of pseudouridylation occurring in lncRNAs in PCa remains opaque. Methods: Spearman correlation analysis and univariate Cox regression were utilized to determine pseudouridylation-related lncRNAs with prognostic value in PCa. Prognostic pseudouridylation-related lncRNAs were included in the LASSO (least absolute shrinkage and selection operator) regression algorithm to develop a predictive model. KM (Kaplan-Meier) survival analysis and ROC (receiver operating characteristic) curves were applied to validate the constructed model. A battery of biological cell assays was conducted to confirm the cancer-promoting effects of RP11-468E2.5 in the model. Results: A classifier containing five pseudouridine-related lncRNAs was developed to stratify PCa patients on BCR and named the "ψ-lnc score." KM survival analysis showed patients in the high ψ-lnc score group experienced BCR more than those in the low ψ-lnc score group. ROC curves demonstrated that ψ-lnc score outperformed other clinical indicators in BCR prediction. An external dataset, GSE54460, was utilized to validate the predictive model's efficacy and authenticity. A ceRNA (competitive endogenous RNA) network was constructed to explore the model's potential molecular functions and was annotated through GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway analyses. RP11-468E2.5 was picked for further investigation, including pan-cancer analysis and experimental validation. Preliminarily, RP11-468E2.5 was confirmed as a tumor promoter. Conclusion: We provide some evidence that pseudouridylation in lncRNA played a role in the development of PCa and propose a novel prognostic classifier for clinical practice.

A novel RNA modification prognostic signature for predicting the characteristics of the tumor microenvironment in gastric cancer

Gastric cancer (GC) is one of the most common neoplastic malignancies, which permutes a fourth of cancer-related mortality globally. RNA modification plays a significant role in tumorigenesis, the underlying molecular mechanism of how different RNA modifications directly affect the tumor microenvironment (TME) in GC is unclear. Here, we profiled the genetic and transcriptional alterations of RNA modification genes (RMGs) in GC samples from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) cohorts. Through the unsupervised clustering algorithm, we identified three distinct RNA modification clusters and found that they participate in different biological pathways and starkly correlate with the clinicopathological characteristics, immune cell infiltration, and prognosis of GC patients. Subsequently, univariate Cox regression analysis unveiled 298 of 684 subtype-related differentially expressed genes (DEGs) are tightly interwoven to prognosis. In addition, we conducted the principal component analysis to develop the RM_Score system, which was used to quantify and predict the prognostic value of RNA modification in GC. Our analysis indicated that patients with high RM_Score were characterized by higher tumor mutational burden, mutation frequency, and microsatellite instability which were more susceptible to immunotherapy and had a favorable prognosis. Altogether, our study uncovered RNA modification signatures that may have a potential role in the TME and prediction of clinicopathological characteristics. Identification of these RNA modifications may provide a new understanding of immunotherapy strategies for gastric cancer.

Structures and mechanisms of tRNA methylation by METTL1-WDR4

Specific, regulated modification of RNAs is important for proper gene expression^1,2. tRNAs are rich with various chemical modifications that affect their stability and function^3,4. 7-Methylguanosine (m⁷G) at tRNA position 46 is a conserved modification that modulates steady-state tRNA levels to affect cell growth^5,6. The METTL1-WDR4 complex generates m⁷G46 in humans, and dysregulation of METTL1-WDR4 has been linked to brain malformation and multiple cancers^7-22. Here we show how METTL1 and WDR4 cooperate to recognize RNA substrates and catalyse methylation. A crystal structure of METTL1-WDR4 and cryo-electron microscopy structures of METTL1-WDR4-tRNA show that the composite protein surface recognizes the tRNA elbow through shape complementarity. The cryo-electron microscopy structures of METTL1-WDR4-tRNA with S-adenosylmethionine or S-adenosylhomocysteine along with METTL1 crystal structures provide additional insights into the catalytic mechanism by revealing the active site in multiple states. The METTL1 N terminus couples cofactor binding with conformational changes in the tRNA, the catalytic loop and the WDR4 C terminus, acting as the switch to activate m⁷G methylation. Thus, our structural models explain how post-translational modifications of the METTL1 N terminus can regulate methylation. Together, our work elucidates the core and regulatory mechanisms underlying m⁷G modification by METTL1, providing the framework to understand its contribution to biology and disease.

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.

ALKBH5 ALLEVIATES HYPOXIA POSTCONDITIONING INJURY IN d -GALACTOSE-INDUCED SENESCENT CARDIOMYOCYTES BY REGULATING STAT3

ABSTRACT

Ischemic postconditioning (I/Post) reduces I/R injury by activating endogenous cardioprotection mechanisms, such as the JAK/signal transducer and activator of transcription 3 (STAT3) and PI3K/Akt pathways, which offer a traditional approach to myocardial protection. According to a previous study, cardioprotection by I/Post may be lost in aged mice, and in our previous research, hypoxic postconditioning (H/Post) lacked a protective effect in senescent cardiomyocytes, which was associated with low expression of long noncoding RNA H19. The N6-methyladenosine (m 6 A) modification is a dynamic and reversible process that has been confirmed to play a role in cardiovascular diseases. However, the mechanisms of m 6 A modification in myocardial I/Post remain to be explored. Neonatal cardiomyocytes were isolated from 2-day-old Sprague-Dawley rats, and senescence was induced by d -galactose, followed by stimulation of hypoxia-reoxygenation and H/Post. Hypoxic injury was evaluated by cell viability and the Bcl-2/Bax protein ratio. Total m 6 A levels were measured using a colorimetric m 6 A RNA Methylation Quantification Kit, and the m 6 A modified and differentially expressed mRNA was determined by MeRIP (methylated RNA immunoprecipitation). We found that H/Post increased m 6 A methylation and decreased RNA mA demethylase alkB homolog 5 (ALKBH5) expression in aged cardiomyocytes. Furthermore, ALKBH5 knockdown exacerbated injury following H/Post in senescent cardiomyocytes. In addition, ALKBH5 regulated STAT3 expression by mediating its m 6 A modification and long noncoding RNA H19/miR-124-3p. ALKBH5 also alleviated the H/Post injury induced by the low expression of STAT3 in senescent cardiomyocytes.

Nanopore Direct RNA Sequencing Data Processing and Analysis Using MasterOfPores

This chapter describes MasterOfPores v.2 (MoP2), an open-source suite of pipelines for processing and analyzing direct RNA Oxford Nanopore sequencing data. The MoP2 relies on the Nextflow DSL2 framework and Linux containers, thus enabling reproducible data analysis in transcriptomic and epitranscriptomic studies. We introduce the key concepts of MoP2 and provide a step-by-step fully reproducible and complete example of how to use the workflow for the analysis of S. cerevisiae total RNA samples sequenced using MinION flowcells. The workflow starts with the pre-processing of raw FAST5 files, which includes basecalling, read quality control, demultiplexing, filtering, mapping, estimation of per-gene/transcript abundances, and transcriptome assembly, with support of the GPU computing for the basecalling and read demultiplexing steps. The secondary analyses of the workflow focus on the estimation of RNA poly(A) tail lengths and the identification of RNA modifications. The MoP2 code is available at https://github.com/biocorecrg/MOP2 and is distributed under the MIT license.

Parent tRNA Modification Status Determines the Induction of Functional tRNA-Derived RNA by Respiratory Syncytial Virus Infection

tRNA-derived RNA fragments (tRFs) are a recently discovered family of small noncoding RNAs (sncRNAs). We previously reported that respiratory syncytial virus (RSV) infection induces functional tRFs, which are derived from a limited subset of parent tRNAs, in airway epithelial cells. Such induction is also observed in nasopharyngeal wash samples from RSV patients and correlates to RSV genome copies, suggesting a clinical significance of tRFs in RSV infection. This work also investigates whether the modification of parent tRNAs is changed by RSV to induce tRFs, using one of the most inducible tRFs as a model. We discovered that RSV infection changed the methylation modification of adenine at position 57 in tRNA glutamic acid, with a codon of CTC (tRNA-GluCTC), and the change is essential for its cleavage. AlkB homolog 1, a previously reported tRNA demethylase, appears to remove methyladenine from tRNA-GluCTC, prompting the subsequent production of tRFs from the 5'-end of tRNA-GluCTC, a regulator of RSV replication. This study demonstrates for the first time the importance of post-transcriptional modification of tRNAs in tRF biogenesis following RSV infection, providing critical insights for antiviral strategy development.

Dynamic regulation and key roles of ribonucleic acid methylation

Ribonucleic acid (RNA) methylation is the most abundant modification in biological systems, accounting for 60% of all RNA modifications, and affects multiple aspects of RNA (including mRNAs, tRNAs, rRNAs, microRNAs, and long non-coding RNAs). Dysregulation of RNA methylation causes many developmental diseases through various mechanisms mediated by N ⁶-methyladenosine (m⁶A), 5-methylcytosine (m⁵C), N ¹-methyladenosine (m¹A), 5-hydroxymethylcytosine (hm⁵C), and pseudouridine (Ψ). The emerging tools of RNA methylation can be used as diagnostic, preventive, and therapeutic markers. Here, we review the accumulated discoveries to date regarding the biological function and dynamic regulation of RNA methylation/modification, as well as the most popularly used techniques applied for profiling RNA epitranscriptome, to provide new ideas for growth and development.

Interplay of four types of RNA modification writers revealed distinct tumor microenvironment and biological characteristics in pancreatic cancer

Background

Pancreatic cancer (PC) is one of the most lethal malignancies and carries a dismal mortality and morbidity. Four types of RNA modification (namely m6A, m1A, APA and A-to-I) could be catalyzed by distinct enzymatic compounds ("writers"), mediating numerous epigenetic events in carcinogenesis and immunomodulation. We aim to investigate the interplay mechanism of these writers in immunogenomic features and molecular biological characteristics in PC.

Methods

We first accessed the specific expression pattern and transcriptional variation of 26 RNA modification writers in The Cancer Genome Atlas (TCGA) dataset. Unsupervised consensus clustering was performed to divide patients into two RNA modification clusters. Then, based on the differentially expressed genes (DEGs) among two clusters, RNA modification score (WM_Score) model was established to determine RNA modification-based subtypes and was validated in International Cancer Genome Consortium (ICGC) dataset. What's more, we manifested the unique status of WM_Score in transcriptional and post-transcriptional regulation, molecular biological characteristics, targeted therapies and immunogenomic patterns.

Results

We documented the tight-knit correlations between transcriptional expression and variation of RNA modification writers. We classified patients into two distinct RNA modification patterns (WM_Score_high and _low), The WM_Score_high subgroup was correlated with worse prognosis, Th2/Th17 cell polarization and oncogenic pathways (e.g. EMT, TGF-β, and mTORC1 signaling pathways), whereas the WM_Score_low subgroup associated with favorable survival rate and Th1 cell trend. WM_Score model also proved robust predictive power in interpreting transcriptional and post-transcriptional events. Additionally, the potential targeted compounds with related pathways for the WM_Score model were further identified.

Conclusions

Our research unfolds a novel horizon on the interplay network of four RNA modifications in PC. This WM_Score model demonstrated powerful predictive capacity in epigenetic, immunological and biological landscape, providing a theoretical basis for future clinical judgments of PC.

The emerging role of epigenetic regulation in the progression of silicosis

Silicosis is one of the most severe occupational diseases worldwide and is characterized by silicon nodules and diffuse pulmonary fibrosis. However, specific treatments for silicosis are still lacking at present. Therefore, elucidating the pathogenesis of silicosis plays a significant guiding role for its treatment and prevention. The occurrence and development of silicosis are accompanied by many regulatory mechanisms, including epigenetic regulation. The main epigenetic regulatory mechanisms of silicosis include DNA methylation, non-coding RNA (ncRNA), and histone modifications. In recent years, the expression and regulation of genes related to silicosis have been explored at epigenetic level to reveal its pathogenesis further, and the identification of aberrant epigenetic markers provides new biomarkers for prediction and diagnosis of silicosis. Here, we summarize the studies on the role of epigenetic changes in the pathogenesis of silicosis to give some clues for finding specific therapeutic targets for silicosis.