RMVar: an updated database of functional variants involved in RNA modifications

Spectral libraries from nucleobases and deoxyribonucleosides facilitate the identification of ribonucleosides by nano-flow liquid chromatography-tandem mass spectrometry

RATIONALE

The study addresses the challenge of identifying RNA post-transcriptional modifications when commercial standards are not available to generate reference spectral libraries. It proposes employing homologous nucleobases and deoxyribonucleosides as alternative reference spectral libraries to aid in identifying modified ribonucleosides and distinguishing them from their positional isomers when the standards are unavailable.

METHODS

Complete sets of ribonucleoside, deoxyribonucleoside and nucleobase standards were analyzed using high-performance nano-flow liquid chromatography coupled to an Orbitrap Eclipse Tribrid mass spectrometer. Spectral libraries were constructed from homologous nucleobases and deoxyribonucleosides using targeted MS2 and neutral-loss-triggered MS3 methods, and collision energies were optimized. The feasibility of using these libraries for identifying modified ribonucleosides and their positional isomers was assessed through comparison of spectral fragmentation patterns.

RESULTS

Our analysis reveals that both MS2 and neutral-loss-triggered MS3 methods yielded rich spectra with similar fragmentation patterns across ribonucleosides, deoxyribonucleosides and nucleobases. Moreover, we demonstrate that spectra from nucleobases and deoxyribonucleosides, generated at optimized collision energies, exhibited sufficient similarity to those of modified ribonucleosides to enable their use as reference spectra for accurate identification of positional isomers within ribonucleoside families.

CONCLUSIONS

The study demonstrates the efficacy of utilizing homologous nucleobases and deoxyribonucleosides as interchangeable reference spectral libraries for identifying modified ribonucleosides and their positional isomers. This approach offers a valuable solution for overcoming limitations posed by the unavailability of commercial standards, enhancing the analysis of RNA post-transcriptional modifications via mass spectrometry.

Gene and protein sequence features augment HLA class I ligand predictions

The sensitivity of malignant tissues to T cell-based immunotherapies depends on the presence of targetable human leukocyte antigen (HLA) class I ligands. Peptide-intrinsic factors, such as HLA class I affinity and proteasomal processing, have been established as determinants of HLA ligand presentation. However, the role of gene and protein sequence features as determinants of epitope presentation has not been systematically evaluated. We perform HLA ligandome mass spectrometry to evaluate the contribution of 7,135 gene and protein sequence features to HLA sampling. This analysis reveals that a number of predicted modifiers of mRNA and protein abundance and turnover, including predicted mRNA methylation and protein ubiquitination sites, inform on the presence of HLA ligands. Importantly, integration of such "hard-coded" sequence features into a machine learning approach augments HLA ligand predictions to a comparable degree as experimental measures of gene expression. Our study highlights the value of gene and protein features for HLA ligand predictions.

DUSP5 regulated by YTHDF1-mediated m6A modification promotes epithelial-mesenchymal transition and EGFR-TKI resistance via the TGF-β/Smad signaling pathway in lung adenocarcinoma

BACKGROUND

Lung adenocarcinoma (LUAD) patients have a dismal survival rate because of cancer metastasis and drug resistance. The study aims to identify the genes that concurrently modulate EMT, metastasis and EGFR-TKI resistance, and to investigate the underlying regulatory mechanisms.

METHODS

Cox regression and Kaplan-Meier analyses were applied to identify prognostic oncogenes in LUAD. Gene set enrichment analysis (GSEA) was used to indicate the biological functions of the gene. Wound-healing and Transwell assays were used to detect migratory and invasive ability. EGFR-TKI sensitivity was evaluated by assessing the proliferation, clonogenic survival and metastatic capability of cancer cells with treatment with gefitinib. Methylated RNA immunoprecipitation (MeRIP) and RNA immunoprecipitation (RIP) analyses established the level of m6A modification present on the target gene and the protein's capability to interact with RNA, respectively. Single-sample gene set enrichment (ssGSEA) algorithm used to investigate levels of immune cell infiltration.

RESULTS

Our study identified dual-specificity phosphatase 5 (DUSP5) as a novel and powerful predictor of adverse outcomes for LUAD by using public datasets. Functional enrichment analysis found that DUSP5 was positively enriched in EMT and transforming growth factor-beta (TGF-β) signaling pathway, a prevailing pathway involved in the induction of EMT. As expected, DUSP5 knockdown suppressed EMT via inhibiting the canonical TGF-β/Smad signaling pathway in in vitro experiments. Consistently, knockdown of DUSP5 was first found to inhibit migratory ability and invasiveness of LUAD cells in in vitro and prevent lung metastasis in in vivo. DUSP5 knockdown re-sensitized gefitinib-resistant LUAD cells to gefitinib, accompanying reversion of EMT progress. In LUAD tissue samples, we found 14 cytosine-phosphate-guanine (CpG) sites of DUSP5 that were negatively associated with DUSP5 gene expression. Importantly, 5'Azacytidine (AZA), an FDA-approved DNA methyltransferase inhibitor, restored DUSP5 expression. Moreover, RIP experiments confirmed that YTH N6-methyladenosine RNA binding protein 1 (YTHDF1), a m6A reader protein, could bind DUSP5 mRNA. YTHDF1 promoted DUSP5 expression and the malignant phenotype of LUAD cells. In addition, the DUSP5-derived genomic model revealed the two clusters with distinguishable immune features and tumor mutational burden (TMB).

CONCLUSIONS

Briefly, our study discovered DUSP5 which was regulated by epigenetic modification, might be a potential therapeutic target, especially in LUAD patients with acquired EGFR-TKI resistance.

Unveiling blood pressure-associated genes in aortic cells through integrative analysis of GWAS and RNA modification-associated variants

Background

Genome-wide association studies (GWAS) have identified more than a thousand loci for blood pressure (BP). Functional genes in these loci are cell-type specific. The aim of this study was to elucidate potentially functional genes associated with BP in the aorta through the utilization of RNA modification-associated single-nucleotide polymorphisms (RNAm-SNPs).

Methods

Utilizing large-scale genetic data of 757,601 individuals from the UK Biobank and International Consortium of Blood Pressure consortium, we identified associations between RNAm-SNPs and BP. The association between RNAm-SNPs, gene expression, and BP were examined.

Results

A total of 355 RNAm-SNPs related to m6A, m1A, m5C, m7G, and A-to-I modification were associated with BP. The related genes were enriched in the pancreatic secretion pathway and renin secretion pathway. The BP GWAS signals were significantly enriched with m6A-SNPs, highlighting the potential functional relevance of m6A in physiological processes influencing BP. Notably, m6A-SNPs in CYP11B1, PDE3B, HDAC7, ACE, SLC4A7, PDE1A, FRK, MTHFR, NPPA, CACNA1D, and HDAC9 were identified. Differential methylation and differential expression of the BP genes in FTO-overexpression and METTL14-knockdown vascular smooth muscle cells were detected. RNAm-SNPs were associated with ascending and descending aorta diameter and the genes showed differential methylation between aortic dissection (AD) cases and controls. In scRNA-seq study, we identified ARID5A, HLA-DPB1, HLA-DRA, IRF1, LINC01091, MCL1, MLF1, MLXIPL, NAA16, NADK, RERG, SRM, and USP53 as differential expression genes for AD in aortic cells.

Conclusion

The present study identified RNAm-SNPs in BP loci and elucidated the associations between the RNAm-SNPs, gene expression, and BP. The identified BP-associated genes in aortic cells were associated with AD.

Comprehensive review on lipid metabolism and RNA methylation: Biological mechanisms, perspectives and challenges

Adipose tissue plays a crucial role in maintaining energy balance, regulating hormones, and promoting metabolic health. To address disorders related to obesity and develop effective therapies, it is essential to have a deep understanding of adipose tissue biology. In recent years, RNA methylation has emerged as a significant epigenetic modification involved in various cellular functions and metabolic pathways. Particularly in the realm of adipogenesis and lipid metabolism, extensive research is ongoing to uncover the mechanisms and functional importance of RNA methylation. Increasing evidence suggests that RNA methylation plays a regulatory role in adipocyte development, metabolism, and lipid utilization across different organs. This comprehensive review aims to provide an overview of common RNA methylation modifications, their occurrences, and regulatory mechanisms, focusing specifically on their intricate connections to fat metabolism. Additionally, we discuss the research methodologies used in studying RNA methylation and highlight relevant databases that can aid researchers in this rapidly advancing field.

NSUN6 mediates 5-methylcytosine modification of METTL3 and promotes colon adenocarcinoma progression

Colon adenocarcinoma (COAD) is a common and fatal malignant tumor of digestive system with complex etiology. 5-Methylcytosine (m5C) modification of RNA by the NSUN gene family (NSUN1-NSUN7) and DNMT2 reshape cell biology and regulate tumor development. However, the expression profile, prognostic significance and function of these m5C modifiers in COAD remain largely unclear. By mining multiple integrated tumor databases, we found that NSUN1, NSUN2, NSUN5, and NSUN6 were overexpressed in COAD tumor samples relative to normal samples. Clinically, high expression of NSUN6 was significantly associated with shorter survival (including both disease-free survival and overall survival) in COAD patients. NSUN6 was further confirmed to be upregulated at both tissue and cellular levels of COAD, suggesting that NSUN6 plays a critical role in disease progression. Through comprehensive gene enrichment analysis and cell-based functional validation, it was revealed that NSUN6 promoted the cell cycle progression and cell proliferation of COAD. Mechanistically, NSUN6 upregulates the expression of oncogenic METTL3 and catalyzes its m5C modification in COAD cells. Overexpression of METTL3 significantly relieved the cell cycle inhibition of COAD caused by NSUN6 deficiency. Furthermore, NSUN6 was negatively associated with the abundance of infiltrating immune cells in COAD tumors, such as activated B cells, natural killer cells, effector memory CD8 T cells, and regulatory T cells. Importantly, pan-cancer analysis further uncovered that NSUN6 was dysregulated and heterogeneous in various tumors. Thus our findings extend the role of m5C transferase in COAD and suggest that NSUN6 is a potential biomarker and target for this malignancy.

N6-methyladenosine-modified SENP1, identified by IGF2BP3, is a novel molecular marker in acute myeloid leukemia and aggravates progression by activating AKT signal via de-SUMOylating HDAC2

BACKGROUND

Elevated evidence suggests that the SENPs family plays an important role in tumor progression. However, the role of SENPs in AML remains unclear.

METHODS

We evaluated the expression pattern of SENP1 based on RNA sequencing data obtained from OHSU, TCGA, TARGET, and MILE datasets. Clinical samples were used to verify the expression of SENP1 in the AML cells. Lentiviral vectors shRNA and sgRNA were used to intervene in SENP1 expression in AML cells, and the effects of SENP1 on AML proliferation and anti-apoptosis were detected using in vitro and in vivo models. Chip-qPCR, MERIP-qPCR, CO-IP, RNA pulldown, and dual-luciferase reporter gene assays were used to explore the regulatory mechanisms of SNEP1 in AML.

RESULTS

SENP1 was significantly upregulated in high-risk AML patients and closely related to poor prognosis. The AKT/mTOR signaling pathway is a key downstream pathway that mediates SENP1's regulation of AML proliferation and anti-apoptosis. Mechanistically, the CO-IP assay revealed binding between SENP1 and HDAC2. SUMO and Chip-qPCR assays suggested that SENP1 can desumoylate HDAC2, which enhances EGFR transcription and activates the AKT pathway. In addition, we found that IGF2BP3 expression was upregulated in high-risk AML patients and was positively correlated with SENP1 expression. MERIP-qPCR and RIP-qPCR showed that IGF2BP3 binds SENP1 3-UTR in an m6A manner, enhances SENP1 expression, and promotes AKT pathway conduction.

CONCLUSIONS

Our findings reveal a distinct mechanism of SENP1-mediated HDAC2-AKT activation and establish the critical role of the IGF2BP3/SENP1signaling axis in AML development.

Homozygous EPRS1 missense variant causing hypomyelinating leukodystrophy-15 alters variant-distal mRNA m6A site accessibility

Hypomyelinating leukodystrophy (HLD) is an autosomal recessive disorder characterized by defective central nervous system myelination. Exome sequencing of two siblings with severe cognitive and motor impairment and progressive hypomyelination characteristic of HLD revealed homozygosity for a missense single-nucleotide variant (SNV) in EPRS1 (c.4444 C > A; p.Pro1482Thr), encoding glutamyl-prolyl-tRNA synthetase, consistent with HLD15. Patient lymphoblastoid cell lines express markedly reduced EPRS1 protein due to dual defects in nuclear export and cytoplasmic translation of variant EPRS1 mRNA. Variant mRNA exhibits reduced METTL3 methyltransferase-mediated writing of N6-methyladenosine (m6A) and reduced reading by YTHDC1 and YTHDF1/3 required for efficient mRNA nuclear export and translation, respectively. In contrast to current models, the variant does not alter the sequence of m6A target sites, but instead reduces their accessibility for modification. The defect was rescued by antisense morpholinos predicted to expose m6A sites on target EPRS1 mRNA, or by m6A modification of the mRNA by METTL3-dCas13b, a targeted RNA methylation editor. Our bioinformatic analysis predicts widespread occurrence of SNVs associated with human health and disease that similarly alter accessibility of distal mRNA m6A sites. These results reveal a new RNA-dependent etiologic mechanism by which SNVs can influence gene expression and disease, consequently generating opportunities for personalized, RNA-based therapeutics targeting these disorders.

RNA modifications in pulmonary diseases

Threatening public health, pulmonary disease (PD) encompasses diverse lung injuries like chronic obstructive PD, pulmonary fibrosis, asthma, pulmonary infections due to pathogen invasion, and fatal lung cancer. The crucial involvement of RNA epigenetic modifications in PD pathogenesis is underscored by robust evidence. These modifications not only shape cell fates but also finely modulate the expression of genes linked to disease progression, suggesting their utility as biomarkers and targets for therapeutic strategies. The critical RNA modifications implicated in PDs are summarized in this review, including N6-methylation of adenosine, N1-methylation of adenosine, 5-methylcytosine, pseudouridine (5-ribosyl uracil), 7-methylguanosine, and adenosine to inosine editing, along with relevant regulatory mechanisms. By shedding light on the pathology of PDs, these summaries could spur the identification of new biomarkers and therapeutic strategies, ultimately paving the way for early PD diagnosis and treatment innovation.

Cigarette smoking combined with genetic variation regulates the m6A methylation of CRNKL1 and is associated with bladder cancer risk

Cigarette smoking was known to accelerate the occurrence and development of bladder cancer by regulating RNA modification. However, the association between the combination of cigarette smoking and RNA modification-related single nucleotide polymorphisms (RNAm-SNPs) and bladder cancer risk remains unclear. In this study, 1681 participants, including 580 cases and 1101 controls, were recruited for genetic association analysis. In total, 1 287 990 RNAm-SNPs involving nine RNA modifications (m6A, m1A, m6Am, 2'-O-Me, m5C, m7G, A-to-I, m5U, and pseudouridine modification) were obtained from the RMVar database. The interactive effect of cigarette smoking and RNAm-SNPs on bladder cancer risk was assessed through joint analysis. The susceptibility analysis revealed that 89 RNAm-SNPs involving m6A, m1A, and A-to-I modifications were associated with bladder cancer risk. Among them, m6A-related rs2273058 in CRNKL1 was associated with bladder cancer risk (odds ratios (OR) = 1.35, padj = 1.78 × 10-4), and CRNKL1 expression was increased in bladder cancer patients (p = 0.035). Cigarette smoking combined with the A allele of rs2273058 increased bladder cancer risk compared with nonsmokers with the G allele of rs2273058 (OR = 2.40, padj = 3.11 × 10-9). Mechanistically, the A allele of rs2273058 endowed CRNKL1 with an additional m6A motif, facilitating recognition by m6A reader IGF2BP1, thereby promoting CRNKL1 expression under cigarette smoking (r = 0.142, p = 0.017). Moreover, elevated CRNKL1 expression may accelerate cell cycle and proliferation, thereby increasing bladder cancer risk. In summary, our study demonstrated that cigarette smoking combined with RNAm-SNPs contributes to bladder cancer risk, which provides a potential target for bladder cancer prevention.

A LATS2 and ALKBH5 positive feedback loop supports their oncogenic roles

N(6)-methyladenosine (m6A) critically regulates RNA dynamics in various biological processes. The m6A demethylase ALKBH5 promotes tumorigenesis of glioblastoma, while the intricate web that orchestrates its regulation remains enigmatic. Here, we discover that cell density affects ALKBH5 subcellular localization and m6A dynamics. Mechanistically, ALKBH5 is phosphorylated by the large tumor suppressor kinase 2 (LATS2), preventing its nuclear export and enhancing protein stability. Furthermore, phosphorylated ALKBH5 reciprocally erases m6A from LATS2 mRNA, thereby stabilizing this transcript. Unexpectedly, LATS2 depletion suppresses glioblastoma stem cell self-renewal independent of yes-associated protein activation. Additionally, deficiency in either LATS2 or ALKBH5 phosphorylation impedes tumor progression in mouse xenograft models. Moreover, high levels of LATS2 expression and ALKBH5 phosphorylation are associated with tumor malignancy in patients with gliomas. Collectively, our study unveils an oncogenic positive feedback loop between LATS2 and ALKBH5, revealing a non-canonical branch of the Hippo pathway for RNA processing and suggesting potential anti-cancer interventions.

Serum micro-RNAs with mutation-targeted RNA modification: a potent cancer detection tool constructed using an optimized machine learning workflow

RNA modifications affect fundamental biological processes and diseases and are a research hotspot. Several micro-RNAs (miRNAs) exhibit genetic variant-targeted RNA modifications that can greatly alter their biofunctions and influence their effect on cancer. Therefore, the potential role of these miRNAs in cancer can be implicated in new prevention and treatment strategies. In this study, we determined whether RMvar-related miRNAs were closely associated with tumorigenesis and identified cancer-specific signatures based on these miRNAs with variants targeting RNA modifications using an optimized machine learning workflow. An effective machine learning workflow, combining least absolute shrinkage and selection operator analyses, recursive feature elimination, and nine types of machine learning algorithms, was used to screen candidate miRNAs from 504 serum RMvar-related miRNAs and construct a diagnostic signature for cancer detection based on 43,047 clinical samples (with an area under the curve value of 0.998, specificity of 93.1%, and sensitivity of 99.3% in the validation cohort). This signature demonstrated a satisfactory diagnostic performance for certain cancers and different conditions, including distinguishing early-stage tumors. Our study revealed the close relationship between RMvar-related miRNAs and tumors and proposed an effective cancer screening tool.

N6-methyladenosine-associated genetic variants in NECTIN2 and HPCAL1 are risk factors for abdominal aortic aneurysm

Although N6-methyladenosine (m6A) modification has been implicated in the pathogenesis of abdominal aortic aneurysm (AAA), the relationship between m6A-associated single nucleotide polymorphisms (m6A-SNPs) and AAA remains unknown. This study used integrative multi-omics analysis and clinical validation approaches to systematically identify potential m6A-SNPs connected with AAA risk. We found that rs6859 and rs10198139 could modulate the expression of local genes, NECTIN2 and HPCAL1, respectively, which exhibited upregulation in AAA tissues, and their risk variants were significantly correlated with an increased susceptibility to AAA. Incorporating rs6859 and rs10198139 improved the efficiency of AAA risk prediction compared to the model considering only conventional risk factors. Additionally, these two SNPs were predicted to be located within the regulatory sequences, and rs6859 showed a substantial impact on m6A modification levels. Our findings suggest that m6A-SNPs rs6859 and rs10198139 confer an elevated risk of AAA, possibly by promoting local gene expression through an m6A-mediated manner.

Regulatory role of m6A epitranscriptomic modifications in normal development and congenital malformations during embryogenesis

The discovery of N6-methyladenosine (m6A) methylation and its role in translation has led to the emergence of a new field of research. Despite accumulating evidence suggesting that m6A methylation is essential for the pathogenesis of cancers and aging diseases by influencing RNA stability, localization, transformation, and translation efficiency, its role in normal and abnormal embryonic development remains unclear. An increasing number of studies are addressing the development of the nervous and gonadal systems during embryonic development, but only few are assessing that of the immune, hematopoietic, urinary, and respiratory systems. Additionally, these studies are limited by the requirement for reliable embryonic animal models and the difficulty in collecting tissue samples of fetuses during development. Multiple studies on the function of m6A methylation have used suitable cell lines to mimic the complex biological processes of fetal development or the early postnatal phase; hence, the research is still in the primary stage. Herein, we discuss current advances in the extensive biological functions of m6A methylation in the development and maldevelopment of embryos/fetuses and conclude that m6A modification occurs extensively during fetal development. Aberrant expression of m6A regulators is probably correlated with single or multiple defects in organogenesis during the intrauterine life. This comprehensive review will enhance our understanding of the pivotal role of m6A modifications involved in fetal development and examine future research directions in embryogenesis.

m6A-TCPred: a web server to predict tissue-conserved human m6A sites using machine learning approach

BACKGROUND

N6-methyladenosine (m6A) is the most prevalent post-transcriptional modification in eukaryotic cells that plays a crucial role in regulating various biological processes, and dysregulation of m6A status is involved in multiple human diseases including cancer contexts. A number of prediction frameworks have been proposed for high-accuracy identification of putative m6A sites, however, none have targeted for direct prediction of tissue-conserved m6A modified residues from non-conserved ones at base-resolution level.

RESULTS

We report here m6A-TCPred, a computational tool for predicting tissue-conserved m6A residues using m6A profiling data from 23 human tissues. By taking advantage of the traditional sequence-based characteristics and additional genome-derived information, m6A-TCPred successfully captured distinct patterns between potentially tissue-conserved m6A modifications and non-conserved ones, with an average AUROC of 0.871 and 0.879 tested on cross-validation and independent datasets, respectively.

CONCLUSION

Our results have been integrated into an online platform: a database holding 268,115 high confidence m6A sites with their conserved information across 23 human tissues; and a web server to predict the conserved status of user-provided m6A collections. The web interface of m6A-TCPred is freely accessible at: www.rnamd.org/m6ATCPred .

N6-methyladenosine-induced METTL1 promotes tumor proliferation via CDK4

N6-methyladenosine (m6A) and N7-methylguanosine (m7G) modification of RNA represent two major intracellular post-transcriptional regulation modes of gene expression. However, the crosstalk of these two epigenetic modifications in tumorigenesis remain poorly understood. Here, we show that m6A methyltransferase METTL3-mediated METTL1 promotes cell proliferation of head and neck squamous cell carcinoma (HNSC) through m7G modification of the cell-cycle regulator CDK4. By mining the database GEPIA, METTL1 was shown to be up-regulated in a broad spectrum of human cancers and correlated with patient clinical outcomes, particularly in HNSC. Mechanistically, METTL3 methylates METTL1 mRNA and mediates its elevation in HNSC via m6A. Functionally, over-expression of METTL1 enhances HNSC cell growth and facilitates cell-cycle progress, while METTL1 knockdown represses these biological behaviors. Moreover, METTL1 physically binds to CDK4 transcript and regulates its m7G modification level to stabilize CDK4. Importantly, the inhibitory effects of METTL1 knockdown on the proliferation of HNSC, esophageal cancer (ESCA), stomach adenocarcinoma (STAD), and colon adenocarcinoma (COAD) were significantly mitigated by over-expression of CDK4. Taken together, this study expands the understanding of epigenetic mechanisms involved in tumorigenesis and identifies the METTL1/CDK4 axis as a potential therapeutic target for digestive system tumors.

ALKHB5-demethylated lncRNA SNHG15 promotes myeloma tumorigenicity by increasing chromatin accessibility and recruiting H3K36me3 modifier SETD2

Chromatin instability plays a crucial role in multiple myeloma (MM) relapse and progression, but its mechanism remains obscure. Here, we uncovered that m6A-demethylase ALKBH5 upregulated and stabilized long noncoding RNA (lncRNA) small nucleolar RNA host gene 15 (SNHG15), which was elevated in MM and positively correlated with unfavorable clinical prognosis factors. ALKBH5-SNHG15 axis participated in viability and migration/invasion of myeloma cell lines and MM-xenografted SCID/NOD mice. Mechanically, ALKBH5 promoted the expression of trimethylated histone H3 at lysine 36 (H3K36me3) methyltransferase SETD2 through lncRNA SNHG15-mediated protein stability. ALKBH5-SNHG15 axis increased chromatin accessibility and altered the H3K36me3 enrichment at the gene body, which is responsible for transcription elongation. Our study suggested a novel epigenetically interaction of N6-methyladenosine (m6A) methylation, lncRNA SNHG15, and histone SETD2/H3K36me3 modifications in myeloma progression, indicating that ALKBH5 and lncRNA SNHG15 could serve as potential novel therapeutic targets for MM treatment.NEW & NOTEWORTHY To our knowledge, this study first demonstrated the prognostic significance and biological function of long noncoding RNA (lncRNA) small nucleolar RNA host gene 15 (SNHG15) in multiple myeloma (MM), and indicated a novel revelation on the effect of N6-methyladenosine (m6A)-regulated lncRNA on MM tumorigenicity. Moreover, the novel chromatin-regulatory mechanism of lncRNA by interacting with epigenetic modifiers including m6A demethylase ALKBH5 and H3K36me3 methyltransferase SETD2 in myeloma progression elucidated intricate mechanism of tumor pathogenesis.

AStruct: detection of allele-specific RNA secondary structure in structuromic probing data

BACKGROUND

Uncovering functional genetic variants from an allele-specific perspective is of paramount importance in advancing our understanding of gene regulation and genetic diseases. Recently, various allele-specific events, such as allele-specific gene expression, allele-specific methylation, and allele-specific binding, have been explored on a genome-wide scale due to the development of high-throughput sequencing methods. RNA secondary structure, which plays a crucial role in multiple RNA-associated processes like RNA modification, translation and splicing, has emerged as an essential focus of relevant research. However, tools to identify genetic variants associated with allele-specific RNA secondary structures are still lacking.

RESULTS

Here, we develop a computational tool called 'AStruct' that enables us to detect allele-specific RNA secondary structure (ASRS) from RT-stop based structuromic probing data. AStruct shows robust performance in both simulated datasets and public icSHAPE datasets. We reveal that single nucleotide polymorphisms (SNPs) with higher AStruct scores are enriched in coding regions and tend to be functional. These SNPs are highly conservative, have the potential to disrupt sites involved in m6A modification or protein binding, and are frequently associated with disease.

CONCLUSIONS

AStruct is a tool dedicated to invoke allele-specific RNA secondary structure events at heterozygous SNPs in RT-stop based structuromic probing data. It utilizes allelic variants, base pairing and RT-stop information under different cell conditions to detect dynamic and functional ASRS. Compared to sequence-based tools, AStruct considers dynamic cell conditions and outperforms in detecting functional variants. AStruct is implemented in JAVA and is freely accessible at: https://github.com/canceromics/AStruct .

A novel serum m7G-harboring microRNA signature for cancer detection

Background: Emerging evidence points to the exceptional importance and value of m7G alteration in the diagnosis and prognosis of cancers. Nonetheless, a biomarker for precise screening of various cancer types has not yet been developed based on serum m7G-harboring miRNAs. Methods: A total of 20,702 serum samples, covering 12 cancer types and consisting of 7,768 cancer samples and 12,934 cancer-free samples were used in this study. A m7G target miRNA diagnostic signature (m7G-miRDS) was established through the least absolute shrinkage and selection operator (LASSO) analyses in a training dataset (n = 10,351), and validated in a validation dataset (n = 10,351). Results: The m7G-miRDS model, a 12 m7G-target-miRNAs signature, demonstrated high accuracy and was qualified for cancer detection. In the training and validation cohort, the area under the curve (AUC) reached 0.974 (95% CI 0.971-0.977) and 0.972 (95% CI 0.969-0.975), respectively. The m7G-miRDS showed superior sensitivity in each cancer type and had a satisfactory AUC in identifying bladder cancer, lung cancer and esophageal cancer. Additionally, the diagnostic performance of m7G-miRDS was not interfered by the gender, age and benign disease. Conclusion: Our results greatly extended the value of serum circulating miRNAs and m7G in cancer detection, and provided a new direction and strategy for the development of novel biomarkers with high accuracy, low cost and less invasiveness for mass cancer screening, such as ncRNA modification.

Exploring N6-methyladenosine (m6A) modification in tree species: opportunities and challenges

N 6-methyladenosine (m6A) in eukaryotes is the most common and widespread internal modification in mRNA. The modification regulates mRNA stability, translation efficiency, and splicing, thereby fine-tuning gene regulation. In plants, m6A is dynamic and critical for various growth stages, embryonic development, morphogenesis, flowering, stress response, crop yield, and biomass. Although recent high-throughput sequencing approaches have enabled the rapid identification of m6A modification sites, the site-specific mechanism of this modification remains unclear in trees. In this review, we discuss the functional significance of m6A in trees under different stress conditions and discuss recent advancements in the quantification of m6A. Quantitative and functional insights into the dynamic aspect of m6A modification could assist researchers in engineering tree crops for better productivity and resistance to various stress conditions.

PRMD: an integrated database for plant RNA modifications

The scope and function of RNA modifications in model plant systems have been extensively studied, resulting in the identification of an increasing number of novel RNA modifications in recent years. Researchers have gradually revealed that RNA modifications, especially N6-methyladenosine (m6A), which is one of the most abundant and commonly studied RNA modifications in plants, have important roles in physiological and pathological processes. These modifications alter the structure of RNA, which affects its molecular complementarity and binding to specific proteins, thereby resulting in various of physiological effects. The increasing interest in plant RNA modifications has necessitated research into RNA modifications and associated datasets. However, there is a lack of a convenient and integrated database with comprehensive annotations and intuitive visualization of plant RNA modifications. Here, we developed the Plant RNA Modification Database (PRMD; http://bioinformatics.sc.cn/PRMD and http://rnainformatics.org.cn/PRMD) to facilitate RNA modification research. This database contains information regarding 20 plant species and provides an intuitive interface for displaying information. Moreover, PRMD offers multiple tools, including RMlevelDiff, RMplantVar, RNAmodNet and Blast (for functional analyses), and mRNAbrowse, RNAlollipop, JBrowse and Integrative Genomics Viewer (for displaying data). Furthermore, PRMD is freely available, making it useful for the rapid development and promotion of research on plant RNA modifications.

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

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

m7GHub V2.0: an updated database for decoding the N7-methylguanosine (m7G) epitranscriptome

With recent progress in mapping N7-methylguanosine (m7G) RNA methylation sites, tens of thousands of experimentally validated m7G sites have been discovered in various species, shedding light on the significant role of m7G modification in regulating numerous biological processes including disease pathogenesis. An integrated resource that enables the sharing, annotation and customized analysis of m7G data will greatly facilitate m7G studies under various physiological contexts. We previously developed the m7GHub database to host mRNA m7G sites identified in the human transcriptome. Here, we present m7GHub v.2.0, an updated resource for a comprehensive collection of m7G modifications in various types of RNA across multiple species: an m7GDB database containing 430 898 putative m7G sites identified in 23 species, collected from both widely applied next-generation sequencing (NGS) and the emerging Oxford Nanopore direct RNA sequencing (ONT) techniques; an m7GDiseaseDB hosting 156 206 m7G-associated variants (involving addition or removal of an m7G site), including 3238 disease-relevant m7G-SNPs that may function through epitranscriptome disturbance; and two enhanced analysis modules to perform interactive analyses on the collections of m7G sites (m7GFinder) and functional variants (m7GSNPer). We expect that m7Ghub v.2.0 should serve as a valuable centralized resource for studying m7G modification. It is freely accessible at: www.rnamd.org/m7GHub2.

A review on the role of RNA methylation in aging-related diseases

Senescence is the underlying mechanism of organism aging and is robustly regulated at the post-transcriptional level. This regulation involves the chemical modifications, of which the RNA methylation is the most common. Recently, a rapidly growing number of studies have demonstrated that methylation is relevant to aging and aging-associated diseases. Owing to the rapid development of detection methods, the understanding on RNA methylation has gone deeper. In this review, we summarize the current understanding on the influence of RNA modification on cellular senescence, with a focus on mRNA methylation in aging-related diseases, and discuss the emerging potential of RNA modification in diagnosis and therapy.

Global RNA modifications to the MALAT1 triple helix differentially affect thermostability and weaken binding to METTL16

Therapeutic mRNAs are generated using modified nucleotides, namely N1-methylpseudouridine (m1Ψ) triphosphate, so that the mRNA evades detection by the immune system. RNA modifications, even at a single-nucleotide position, perturb RNA structure, although it is not well understood how structure and function is impacted by globally modified RNAs. Therefore, we examined the metastasis-associated lung adenocarcinoma transcript 1 triple helix, a highly structured stability element that includes single-, double-, and triple-stranded RNA, globally modified with N6-methyladenosine (m6A), pseudouridine (Ψ), or m1Ψ. UV thermal denaturation assays showed that m6A destabilizes both the Hoogsteen and Watson-Crick faces of the RNA by ∼20 °C, Ψ stabilizes the Hoogsteen and Watson-Crick faces of the RNA by ∼12 °C, and m1Ψ has minimal effect on the stability of the Hoogsteen face of the RNA but increases the stability of the Watson-Crick face by ∼9 °C. Native gel-shift assays revealed that binding of the methyltransferase-like protein 16 to the metastasis-associated lung adenocarcinoma transcript 1 triple helix was weakened by at least 8-, 99-, and 23-fold, respectively, when RNA is globally modified with m6A, Ψ, or m1Ψ. These results demonstrate that a more thermostable RNA structure does not lead to tighter RNA-protein interactions, thereby highlighting the regulatory power of RNA modifications by multiple means.

Identifying prognostic characteristics of m6A-related glycolysis gene and predicting the immune infiltration landscape in bladder cancer

BACKGROUNDS

Glucose metabolism is associated with the development of cancers, and m6A RNA methylation regulator-related genes play vital roles in bladder urothelial carcinoma (BLCA). However, the role of m6A-related glucose metabolism genes in BLCA occurrence and development has not yet been reported. Our study aims to integrate m6A- and glycolysis-related genes and find potential gene targets for clinical diagnosis and prognosis of BLCA patients.

METHODS

Sequencing data and clinical information on BLCA were extracted from common databases. Univariate Cox analysis was used to screen prognosis-related m6A glucose metabolism genes; BLCA subtypes were distinguished using consensus clustering analysis. Subsequently, genes associated with BLCA occurrence and development were identified using the "limma" R package. The risk score was then calculated, and a nomogram was constructed to predict survival rate of BLCA patients. Functional and immune microenvironment analyses were performed to explore potential functions and mechanisms of the different risk groups.

RESULTS

Based on 70 prognosis-related m6A glucose metabolism genes, BLCA was classified into two subtypes, and 34 genes associated with its occurrence and development were identified. Enrichment analysis revealed an association of genes in high-risk groups with tricarboxylic acid cycle function and glycolysis. Moreover, significantly higher levels of seven immune checkpoints, 14 immune checkpoint inhibitors, and 32 immune factors were found in high-risk score groups.

CONCLUSIONS

This study identified two biomarkers associated with BLCA prognosis; these findings may deepen our understanding of the role of m6A-related glucose metabolism genes in BLCA development. We constructed a m6A-related glucose metabolism- and immune-related gene risk model, which could effectively predict patient prognosis and immunotherapy response and guide individualized immunotherapy.

m5U-GEPred: prediction of RNA 5-methyluridine sites based on sequence-derived and graph embedding features

5-Methyluridine (m5U) is one of the most common post-transcriptional RNA modifications, which is involved in a variety of important biological processes and disease development. The precise identification of the m5U sites allows for a better understanding of the biological processes of RNA and contributes to the discovery of new RNA functional and therapeutic targets. Here, we present m5U-GEPred, a prediction framework, to combine sequence characteristics and graph embedding-based information for m5U identification. The graph embedding approach was introduced to extract the global information of training data that complemented the local information represented by conventional sequence features, thereby enhancing the prediction performance of m5U identification. m5U-GEPred outperformed the state-of-the-art m5U predictors built on two independent species, with an average AUROC of 0.984 and 0.985 tested on human and yeast transcriptomes, respectively. To further validate the performance of our newly proposed framework, the experimentally validated m5U sites identified from Oxford Nanopore Technology (ONT) were collected as independent testing data, and in this project, m5U-GEPred achieved reasonable prediction performance with ACC of 91.84%. We hope that m5U-GEPred should make a useful computational alternative for m5U identification.

The m6A reader YTHDF2 alleviates the inflammatory response by inhibiting IL-6R/JAK2/STAT1 pathway-mediated high-mobility group box-1 release

Background

Sepsis is a common severe complication in major burn victims and is characterized by a dysregulated systemic response to inflammation. YTH domain family 2 (YTHDF2), a well-studied N6-methyladenosine (m6A) reader that specifically recognizes and binds to m6A-modified transcripts to mediate their degradation, is connected to pathogenic and physiological processes in eukaryotes, but its effect on sepsis is still unknown. We aimed to discover the effects and mechanisms of YTHDF2 in sepsis.

Methods

Quantitative reverse transcription-polymerase chain reaction (qRT-PCR) and western blot analyses were used to measure the expression of YTHDF2, the interleukin 6 receptor (IL-6R), high-mobility group box-1 (HMGB1), Janus kinase 2 (JAK2) and signal transducer and activator of transcription 1 (STAT1) under different in vitro conditions. Enzyme-linked immunosorbent assays were utilized to evaluate the expression of HMGB1, IL-6, IL-1β and tumor necrosis factor-α. To confirm that YTHDF2 specifically targets IL-6R mRNA, RNA immunoprecipitation and dual-luciferase reporter assays were performed. Finally, we utilized a mouse model of lipopolysaccharide (LPS)-induced sepsis to verify the effects of YTHDF2 in vivo.

Results

According to our findings, YTHDF2 was expressed at a low level in peripheral blood mononuclear cells from septic mice and patients as well as in LPS-induced RAW264.7 cells. Overexpression of YTHDF2 alleviated the inflammatory response by inhibiting HMGB1 release and JAK2/STAT1 signalling in LPS-stimulated cells. Mechanistically, YTHDF2 suppressed JAK2/STAT1 signalling by directly recognizing the m6A-modified site in IL-6R and decreasing the stability of IL-6R mRNA, thereby inhibiting HMGB1 release. In vivo experiments showed that YTHDF2 played a protective role in septic mice by suppressing the IL-6R/JAK2/STAT1/HMGB1 axis.

Conclusions

In summary, these findings demonstrate that YTHDF2 plays an essential role as an inhibitor of inflammation to reduce the release of HMGB1 by inhibiting the IL-6R/JAK2/STAT1 pathway, indicating that YTHDF2 is a novel target for therapeutic interventions in sepsis.

Shining a spotlight on m6A and the vital role of RNA modification in endometrial cancer: a review

RNA modifications are mostly dynamically reversible post-transcriptional modifications, of which m6A is the most prevalent in eukaryotic mRNAs. A growing number of studies indicate that RNA modification can finely tune gene expression and modulate RNA metabolic homeostasis, which in turn affects the self-renewal, proliferation, apoptosis, migration, and invasion of tumor cells. Endometrial carcinoma (EC) is the most common gynecologic tumor in developed countries. Although it can be diagnosed early in the onset and have a preferable prognosis, some cases might develop and become metastatic or recurrent, with a worse prognosis. Fortunately, immunotherapy and targeted therapy are promising methods of treating endometrial cancer patients. Gene modifications may also contribute to these treatments, as is especially the case with recent developments of new targeted therapeutic genes and diagnostic biomarkers for EC, even though current findings on the relationship between RNA modification and EC are still very limited, especially m6A. For example, what is the elaborate mechanism by which RNA modification affects EC progression? Taking m6A modification as an example, what is the conversion mode of methylation and demethylation for RNAs, and how to achieve selective recognition of specific RNA? Understanding how they cope with various stimuli as part of in vivo and in vitro biological development, disease or tumor occurrence and development, and other processes is valuable and RNA modifications provide a distinctive insight into genetic information. The roles of these processes in coping with various stimuli, biological development, disease, or tumor development in vivo and in vitro are self-evident and may become a new direction for cancer in the future. In this review, we summarize the category, characteristics, and therapeutic precis of RNA modification, m6A in particular, with the purpose of seeking the systematic regulation axis related to RNA modification to provide a better solution for the treatment of EC.

Combined signature of N7-methylguanosine regulators with their related genes and the tumor microenvironment: a prognostic and therapeutic biomarker for breast cancer

Background

Identifying predictive markers for breast cancer (BC) prognosis and immunotherapeutic responses remains challenging. Recent findings indicate that N7-methylguanosine (m7G) modification and the tumor microenvironment (TME) are critical for BC tumorigenesis and metastasis, suggesting that integrating m7G modifications and TME cell characteristics could improve the predictive accuracy for prognosis and immunotherapeutic responses.

Methods

We utilized bulk RNA-sequencing data from The Cancer Genome Atlas Breast Cancer Cohort and the GSE42568 and GSE146558 datasets to identify BC-specific m7G-modification regulators and associated genes. We used multiple m7G databases and RNA interference to validate the relationships between BC-specific m7G-modification regulators (METTL1 and WDR4) and related genes. Single-cell RNA-sequencing data from GSE176078 confirmed the association between m7G modifications and TME cells. We constructed an m7G-TME classifier, validated the results using an independent BC cohort (GSE20685; n = 327), investigated the clinical significance of BC-specific m7G-modifying regulators by reverse transcription-quantitative polymerase chain reaction (RT-qPCR) analysis, and performed tissue-microarray assays on 192 BC samples.

Results

Immunohistochemistry and RT-qPCR results indicated that METTL1 and WDR4 overexpression in BC correlated with poor patient prognosis. Moreover, single-cell analysis revealed relationships between m7G modification and TME cells, indicating their potential as indicators of BC prognosis and treatment responses. The m7G-TME classifier enabled patient subgrouping and revealed significantly better survival and treatment responses in the m7Glow+TMEhigh group. Significant differences in tumor biological functions and immunophenotypes occurred among the different subgroups.

Conclusions

The m7G-TME classifier offers a promising tool for predicting prognosis and immunotherapeutic responses in BC, which could support personalized therapeutic strategies.

A New Understanding of Long Non-Coding RNA in Hepatocellular Carcinoma-From m6A Modification to Blood Biomarkers

With recent advancements in biological research, long non-coding RNAs (lncRNAs) with lengths exceeding 200 nucleotides have emerged as pivotal regulators of gene expression and cellular phenotypic modulation. Despite initial skepticism due to their low sequence conservation and expression levels, their significance in various biological processes has become increasingly apparent. We provided an overview of lncRNAs and discussed their defining features and modes of operation. We then explored their crucial function in the hepatocarcinogenesis process, elucidating their complex involvement in hepatocellular carcinoma (HCC). The influential role of lncRNAs within the HCC tumor microenvironment is emphasized, illustrating their potential as key modulators of disease dynamics. We also investigated the significant influence of N6-methyladenosine (m6A) modification on lncRNA function in HCC, enhancing our understanding of both their roles and their upstream regulators. Additionally, the potential of lncRNAs as promising biomarkers was discussed in liver cancer diagnosis, suggesting a novel avenue for future research and clinical application. Finally, our work underscored the dual potential of lncRNAs as both contributors to HCC pathogenesis and innovative tools for its diagnosis. Existing challenges and prospective trajectories in lncRNA research are also discussed, emphasizing their potential in advancing liver cancer research.

RNA modifications in cancer

Currently, more than 170 modifications have been identified on RNA. Among these RNA modifications, various methylations account for two-thirds of total cases and exist on almost all RNAs. Roles of RNA modifications in cancer are garnering increasing interest. The research on m6A RNA methylation in cancer is in full swing at present. However, there are still many other popular RNA modifications involved in the regulation of gene expression post-transcriptionally besides m6A RNA methylation. In this review, we focus on several important RNA modifications including m1A, m5C, m7G, 2'-O-Me, Ψ and A-to-I editing in cancer, which will provide a new perspective on tumourigenesis by peeking into the complex regulatory network of epigenetic RNA modifications, transcript processing, and protein translation.

m6A-TSHub: Unveiling the Context-specific m6A Methylation and m6A-affecting Mutations in 23 Human Tissues

As the most pervasive epigenetic marker present on mRNAs and long non-coding RNAs (lncRNAs), N6-methyladenosine (m6A) RNA methylation has been shown to participate in essential biological processes. Recent studies have revealed the distinct patterns of m6A methylome across human tissues, and a major challenge remains in elucidating the tissue-specific presence and circuitry of m6A methylation. We present here a comprehensive online platform, m6A-TSHub, for unveiling the context-specific m6A methylation and genetic mutations that potentially regulate m6A epigenetic mark. m6A-TSHub consists of four core components, including (1) m6A-TSDB, a comprehensive database of 184,554 functionally annotated m6A sites derived from 23 human tissues and 499,369 m6A sites from 25 tumor conditions, respectively; (2) m6A-TSFinder, a web server for high-accuracy prediction of m6A methylation sites within a specific tissue from RNA sequences, which was constructed using multi-instance deep neural networks with gated attention; (3) m6A-TSVar, a web server for assessing the impact of genetic variants on tissue-specific m6A RNA modifications; and (4) m6A-CAVar, a database of 587,983 The Cancer Genome Atlas (TCGA) cancer mutations (derived from 27 cancer types) that were predicted to affect m6A modifications in the primary tissue of cancers. The database should make a useful resource for studying the m6A methylome and the genetic factors of epitranscriptome disturbance in a specific tissue (or cancer type). m6A-TSHub is accessible at www.xjtlu.edu.cn/biologicalsciences/m6ats.

Epigenetic drugs as new emerging therapeutics: What is the scale's orientation of application and challenges?

Epigenetics is the study of heritable changes in gene expression or function without altering the DNA sequence. Important factors are part of epigenetic events, such as methylation, DNA histone rearrangements, nucleosome transposition, and non-coding RNAs. Dysregulated epigenetic mechanics are associated with various cancers' initiation, development, and metastasis. It is known that the occurrence and development of cancer can be controlled by regulating unexpected epigenetic events. Epi-drugs are used singly or in combination with chemotherapy and enhance antitumor activity, reduce drug resistance, and stimulate the host immune response. Despite these benefits, epigenetic therapy as a single therapy or in combination with other drugs leads to adverse effects. This review article introduces and compares the advantages, disadvantages, and side effects of using these drugs for the first time since their introduction. Also, this article describes the mechanism of action of various epigenetic drugs. Recommendations for future use of epigenetic drugs as cancer therapeutics are suggested as an overall conclusion.

Crosstalk between m6A and coding/non-coding RNA in cancer and detection methods of m6A modification residues

N6-methyladenosine (m6A) is one of the most common and well-known internal RNA modifications that occur on mRNAs or ncRNAs. It affects various aspects of RNA metabolism, including splicing, stability, translocation, and translation. An abundance of evidence demonstrates that m6A plays a crucial role in various pathological and biological processes, especially in tumorigenesis and tumor progression. In this article, we introduce the potential functions of m6A regulators, including "writers" that install m6A marks, "erasers" that demethylate m6A, and "readers" that determine the fate of m6A-modified targets. We have conducted a review on the molecular functions of m6A, focusing on both coding and noncoding RNAs. Additionally, we have compiled an overview of the effects noncoding RNAs have on m6A regulators and explored the dual roles of m6A in the development and advancement of cancer. Our review also includes a detailed summary of the most advanced databases for m6A, state-of-the-art experimental and sequencing detection methods, and machine learning-based computational predictors for identifying m6A sites.

Novel insight into RNA modifications in tumor immunity: Promising targets to prevent tumor immune escape

An immunosuppressive state is a typical feature of the tumor microenvironment. Despite the dramatic success of immune checkpoint inhibitor (ICI) therapy in preventing tumor cell escape from immune surveillance, primary and acquired resistance have limited its clinical use. Notably, recent clinical trials have shown that epigenetic drugs can significantly improve the outcome of ICI therapy in various cancers, indicating the importance of epigenetic modifications in immune regulation of tumors. Recently, RNA modifications (N⁶-methyladenosine [m⁶A], N¹-methyladenosine [m¹A], 5-methylcytosine [m⁵C], etc.), novel hotspot areas of epigenetic research, have been shown to play crucial roles in protumor and antitumor immunity. In this review, we provide a comprehensive understanding of how m⁶A, m¹A, and m⁵C function in tumor immunity by directly regulating different immune cells as well as indirectly regulating tumor cells through different mechanisms, including modulating the expression of immune checkpoints, inducing metabolic reprogramming, and affecting the secretion of immune-related factors. Finally, we discuss the current status of strategies targeting RNA modifications to prevent tumor immune escape, highlighting their potential.

Rescuing biologically relevant consensus regions across replicated samples

BACKGROUND

Protein-DNA binding sites of ChIP-seq experiments are identified where the binding affinity is significant based on a given threshold. The choice of the threshold is a trade-off between conservative region identification and discarding weak, but true binding sites.

RESULTS

We rescue weak binding sites using MSPC, which efficiently exploits replicates to lower the threshold required to identify a site while keeping a low false-positive rate, and we compare it to IDR, a widely used post-processing method for identifying highly reproducible peaks across replicates. We observe several master transcription regulators (e.g., SP1 and GATA3) and HDAC2-GATA1 regulatory networks on rescued regions in K562 cell line.

CONCLUSIONS

We argue the biological relevance of weak binding sites and the information they add when rescued by MSPC. An implementation of the proposed extended MSPC methodology and the scripts to reproduce the performed analysis are freely available at https://genometric.github.io/MSPC/ ; MSPC is distributed as a command-line application and an R package available from Bioconductor ( https://doi.org/doi:10.18129/B9.bioc.rmspc ).

A positive feedback circuit between RN7SK snRNA and m6A readers is essential for tumorigenesis

N6-Methyladenosine (m6A) RNA modification, methylation at the N6 position of adenosine, plays critical roles in tumorigenesis. m6A readers recognize m6A modifications and thus act as key executors for the biological consequences of RNA methylation. However, knowledge about the regulatory mechanism(s) of m6A readers is extremely limited. In this study, RN7SK was identified as a small nuclear RNA that interacts with m6A readers. m6A readers recognized and facilitated secondary structure formation of m6A-modified RN7SK, which in turn prevented m6A reader mRNA degradation from exonucleases. Thus, a positive feedback circuit between RN7SK and m6A readers is established in tumor cells. From findings on the interaction with RN7SK, new m6A readers, such as EWS RNA binding protein 1 (EWSR1) and KH RNA binding domain containing, signal transduction-associated 1 (KHDRBS1), were identified and shown to boost Wnt/β-catenin signaling and tumorigenesis by suppressing translation of Cullin1 (CUL1). Moreover, several Food and Drug Administration-approved small molecules were demonstrated to reduce RN7SK expression and inhibit tumorigenesis. Together, these findings reveal a common regulatory mechanism of m6A readers and indicate that targeting RN7SK has strong potential for tumor treatment.

Integrating machine learning and bioinformatics analysis to m6A regulator-mediated methylation modification models for predicting glioblastoma patients' prognosis and immunotherapy response

BACKGROUND

Epigenetic regulations of immune responses are essential for cancer development and growth. As a critical step, comprehensive and rigorous explorations of m6A methylation are important to determine its prognostic significance, tumor microenvironment (TME) infiltration characteristics and underlying relationship with glioblastoma (GBM).

METHODS

To evaluate m6A modification patterns in GBM, we conducted unsupervised clustering to determine the expression levels of GBM-related m6A regulatory factors and performed differential analysis to obtain m6A-related genes. Consistent clustering was used to generate m6A regulators cluster A and B. Machine learning algorithms were implemented for identifying TME features and predicting the response of GBM patients receiving immunotherapy.

RESULTS

It is found that the m6A regulatory factor significantly regulates the mutation of GBM and TME. Based on Europe, America, and China data, we established m6Ascore through the m6A model. The model accurately predicted the results of 1206 GBM patients from the discovery cohort. Additionally, a high m6A score was associated with poor prognoses. Significant TME features were found among the different m6A score groups, which demonstrated positive correlations with biological functions (i.e., EMT2) and immune checkpoints.

CONCLUSIONS

m6A modification was important to characterize the tumorigenesis and TME infiltration in GBM. The m6Ascore provided GBM patients with valuable and accurate prognosis and prediction of clinical response to various treatment modalities, which could be useful to guide patient treatments.

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.

Ythdf2 promotes pulmonary hypertension by suppressing Hmox1-dependent anti-inflammatory and antioxidant function in alveolar macrophages

Pulmonary hypertension (PH) is a devastating disease characterized by irreversible pulmonary vascular remodeling (PVR) that causes right ventricular failure and death. The early alternative activation of macrophages is a critical event in the development of PVR and PH, but the underlying mechanisms remain elusive. Previously we have shown that N⁶-methyladenosine (m⁶A) modifications of RNA contribute to phenotypic switching of pulmonary artery smooth muscle cells and PH. In the current study, we identify Ythdf2, an m⁶A reader, as an important regulator of pulmonary inflammation and redox regulation in PH. In a mouse model of PH, the protein expression of Ythdf2 was increased in alveolar macrophages (AMs) during the early stages of hypoxia. Mice with a myeloid specific knockout of Ythdf2 (Ythdf2^{Lyz2 Cre}) were protected from PH with attenuated right ventricular hypertrophy and PVR compared to control mice and this was accompanied by decreased macrophage polarization and oxidative stress. In the absence of Ythdf2, heme oxygenase 1 (Hmox1) mRNA and protein expression were significantly elevated in hypoxic AMs. Mechanistically, Ythdf2 promoted the degradation of Hmox1 mRNA in a m⁶A dependent manner. Furthermore, an inhibitor of Hmox1 promoted macrophage alternative activation, and reversed the protection from PH seen in Ythdf2^{Lyz2 Cre} mice under hypoxic exposure. Together, our data reveal a novel mechanism linking m⁶A RNA modification with changes in macrophage phenotype, inflammation and oxidative stress in PH, and identify Hmox1 as a downstream target of Ythdf2, suggesting that Ythdf2 may be a therapeutic target in PH.

Genome-wide identification of RNA modification-related single nucleotide polymorphisms associated with rheumatoid arthritis

BACKGROUND

RNA modification plays important roles in many biological processes, such as gene expression control. The aim of this study was to identify single nucleotide polymorphisms related to RNA modification (RNAm-SNPs) for rheumatoid arthritis (RA) as putative functional variants.

METHODS

We examined the association of RNAm-SNPs with RA in summary data from a genome-wide association study of 19,234 RA cases and 61,565 controls. We performed eQTL and pQTL analyses for the RNAm-SNPs to find associated gene expression and protein levels. Furthermore, we examined the associations of gene expression and circulating protein levels with RA using two-sample Mendelian randomization analysis methods.

RESULTS

A total of 160 RNAm-SNPs related to m⁶A, m¹A, A-to-I, m⁷G, m⁵C, m⁵U and m⁶Am modifications were identified to be significantly associated with RA. These RNAm-SNPs were located in 62 protein-coding genes, which were significantly enriched in immune-related pathways. RNAm-SNPs in important RA susceptibility genes, such as PADI2, SPRED2, PLCL2, HLA-A, HLA-B, HLA-DRB1, HLA-DPB1, TRAF1 and TXNDC11, were identified. Most of these RNAm-SNPs showed eQTL effects, and the expression levels of 26 of the modifiable genes (e.g., PADI2, TRAF1, HLA-A, HLA-DRB1, HLA-DPB1 and HLA-B) in blood cells were associated with RA. Circulating protein levels, such as CFB, GZMA, HLA-DQA2, IL21, LRPAP1 and TFF3, were affected by RNAm-SNPs and were associated with RA.

CONCLUSION

The present study identified RNAm-SNPs in the reported RA susceptibility genes and suggested that RNAm-SNPs may affect RA risk by affecting the expression levels of corresponding genes and proteins.

The m6A methyltransferase METTL16 negatively regulates MCP1 expression in mesenchymal stem cells during monocyte recruitment

Mesenchymal stem cells (MSCs) possess strong immunoregulatory functions, one aspect of which is recruiting monocytes from peripheral vessels to local tissue by secreting monocyte chemoattractant protein 1 (MCP1). However, the regulatory mechanisms of MCP1 secretion in MSCs are still unclear. Recently, the N6-methyladenosine (m6A) modification was reported to be involved in the functional regulation of MSCs. In this study, we demonstrated that methyltransferase-like 16 (METTL16) negatively regulated MCP1 expression in MSCs through the m6A modification. Specifically, the expression of METTL16 in MSCs decreased gradually and was negatively correlated with the expression of MCP1 after coculture with monocytes. Knocking down METTL16 markedly enhanced MCP1 expression and the ability to recruit monocytes. Mechanistically, knocking down METTL16 decreased MCP1 mRNA degradation, which was mediated by the m6A reader YTH N6-methyladenosine RNA-binding protein 2 (YTHDF2). We further revealed that YTHDF2 specifically recognized m6A sites on MCP1 mRNA in the CDS region and thus negatively regulated MCP1 expression. Moreover, an in vivo assay showed that MSCs transfected with METTL16 siRNA showed greater ability to recruit monocytes. These findings reveal a potential mechanism by which the m6A methylase METTL16 regulates MCP1 expression through YTHDF2-mediated mRNA degradation and suggest a potential strategy to manipulate MCP1 expression in MSCs.

METTL3 inhibits autoreactive Th17 cell responses in experimental autoimmune uveitis via stabilizing ASH1L mRNA

Methyltransferase like 3 (METTL3), a primary N6-methyladenosine (m6A) methyltransferase, has been implicated in various biological and pathological processes including immune responses. However, the functions and mechanisms of METTL3 in pathogenic T helper (Th)17 cells are poorly understood. Here we found significantly decreased METTL3 expression along with reduced m6A levels in eyeballs and T cells of experimental autoimmune uveitis (EAU). Overexpression of METTL3 ameliorated the development of EAU and suppressed pathogenic Th17 cell responses in vivo and in vitro. Mechanistically, METTL3 promoted the expression of absent, small, or homeotic-like 1 (ASH1L) via enhancing its stability in a YT521-B homology domain containing 2 (YTHDC2)-dependent manner, which further decreased the expression of IL-17 and IL-23 receptor (IL-23R), resulting in reduced pathogenic Th17 responses. Together, our data reveal a pivotal role of METTL3 in regulating pathogenic Th17 responses, which may contribute to human uveitis therapy.

RBM3 suppresses stemness remodeling of prostate cancer in bone microenvironment by modulating N6-methyladenosine on CTNNB1 mRNA

Bone metastasis is the most happened metastatic event in prostate cancer (PCa) and needs a large effort in treatment. When PCa metastasizes to the bone, the new microenvironment can induce the epigenome reprogramming and stemness remodeling of cancer cells, thereby increasing the adaptability of cancer cells to the bone microenvironment, and this even leads to the occurrence of secondary tumor metastasis. Our group has previously found that RNA binding motif 3 (RBM3) affects the stem cell-like properties of PCa by interfering with alternative splicing of CD44. However, whether RBM3, as a stress-response protein, can resist microenvironmental remodeling of PCa particularly in bone metastasis remains unknown. By co-culturing PCa cells with osteoblasts to mimic PCa bone metastases, we found that RBM3 upregulates the N6-methyladenosine (m⁶A) methylation on the mRNA of catenin beta 1 (CTNNB1) in a manner dependent on methyltransferase 3 (METTL3), an N6-adenosine-methyltransferase complex catalytic subunit. Consequently, this modification results in a decreased stability of CTNNB1 mRNA and a followed inactivation of Wnt signaling, which ultimately inhibits the stemness remodeling of PCa cells by osteoblasts. Thus, the present study may extend our understanding of the inhibitory role of RBM3 on particularly bone metastasis of PCa.

Occurrence of RNA post-transcriptional modifications in plant viruses and viroids and their correlation with structural and functional features

Post-transcriptional modifications of RNA bases are widespread across all the tree of life and have been linked to RNA maturation, stability, and molecular interactions. RNA modifications have been extensively described in endogenous eukaryotic mRNAs, however, little is known about the presence of RNA modifications in plant viral and subviral RNAs. Here, we used a computational approach to infer RNA modifications in plant-pathogenic viruses and viroids using high-throughput annotation of modified ribonucleotides (HAMR), a software that predicts modified ribonucleotides using high-throughput RNA sequencing data. We analyzed datasets from representative members of different plant viruses and viroids and compared them to plant-endogenous mRNAs. Our approach was able to predict potential RNA chemical modifications (RCMs) in all analyzed pathogens. We found that both DNA and RNA viruses presented a wide range of RCM proportions while viroids had lowest values. Furthermore, we found that for viruses with segmented genomes, some genomic RNAs had a higher proportion of RCM. Interestingly, nuclear-replicating viroids showed most of the predicted modifications located in the pathogenesis region, pointing towards a possible functional role of RCMs in their infectious cycle. Thus, our results strongly suggest that plant viral and subviral RNAs might contain a variety of previously unreported RNA modifications, thus opening a new perspective in the multifaceted process of plant-pathogen interactions.

Inhibition of IGF2BP1 attenuates renal injury and inflammation by alleviating m6A modifications and E2F1/MIF pathway

Septic acute kidney injury (AKI) is characterized by inflammation. Pyroptosis often occurs during AKI and is associated with the development of septic AKI. This study found that induction of insulin-like growth factor 2 mRNA binding protein 1 (IGF2BP1) to a higher level can induce pyroptosis in renal tubular cells. Meanwhile, macrophage migration inhibitory factor (MIF), a subunit of NLRP3 inflammasomes, was essential for IGF2BP1-induced pyroptosis. A putative m6A recognition site was identified at the 3'-UTR region of E2F transcription factor 1 (E2F1) mRNA via bioinformatics analyses and validated using mutation and luciferase experiments. Further actinomycin D (Act D) chase experiments showed that IGF2BP1 stabilized E2F1 mRNA dependent on m6A. Electrophoretic mobility shift assay (EMSA) and chromatin immunoprecipitation (ChIP) indicated that E2F1 acted as a transcription factor to promote MIF expression. Thus, IGF2BP1 upregulated MIF through directly upregulating E2F1 expression via m6A modification. Experiments on mice with cecum ligation puncture (CLP) surgery verified the relationships between IGF2BP1, E2F1, and MIF and demonstrated the significance of IGF2BP1 in MIF-associated pyroptosis in vivo. In conclusion, IGF2BP1 was a potent pyroptosis inducer in septic AKI through targeting the MIF component of NLRP3 inflammasomes. Inhibiting IGF2BP1 could be an alternate pyroptosis-based treatment for septic AKI.

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.

RNA Epigenetics in Chronic Lung Diseases

Chronic lung diseases are highly prevalent worldwide and cause significant mortality. Lung cancer is the end stage of many chronic lung diseases. RNA epigenetics can dynamically modulate gene expression and decide cell fate. Recently, studies have confirmed that RNA epigenetics plays a crucial role in the developing of chronic lung diseases. Further exploration of the underlying mechanisms of RNA epigenetics in chronic lung diseases, including lung cancer, may lead to a better understanding of the diseases and promote the development of new biomarkers and therapeutic strategies. This article reviews basic information on RNA modifications, including N⁶ methylation of adenosine (m⁶A), N¹ methylation of adenosine (m¹A), N⁷-methylguanosine (m⁷G), 5-methylcytosine (m⁵C), 2'O-methylation (2'-O-Me or Nm), pseudouridine (5-ribosyl uracil or Ψ), and adenosine to inosine RNA editing (A-to-I editing). We then show how they relate to different types of lung disease. This paper hopes to summarize the mechanisms of RNA modification in chronic lung disease and finds a new way to develop early diagnosis and treatment of chronic lung disease.

Predicting N6-Methyladenosine Sites in Multiple Tissues of Mammals through Ensemble Deep Learning

N6-methyladenosine (m⁶A) is the most abundant within eukaryotic messenger RNA modification, which plays an essential regulatory role in the control of cellular functions and gene expression. However, it remains an outstanding challenge to detect mRNA m⁶A transcriptome-wide at base resolution via experimental approaches, which are generally time-consuming and expensive. Developing computational methods is a good strategy for accurate in silico detection of m⁶A modification sites from the large amount of RNA sequence data. Unfortunately, the existing computational models are usually only for m⁶A site prediction in a single species, without considering the tissue level of species, while most of them are constructed based on low-confidence level data generated by an m⁶A antibody immunoprecipitation (IP)-based sequencing method, thereby restricting reliability and generalizability of proposed models. Here, we review recent advances in computational prediction of m⁶A sites and construct a new computational approach named im6APred using ensemble deep learning to accurately identify m⁶A sites based on high-confidence level data in multiple tissues of mammals. Our model im6APred builds upon a comprehensive evaluation of multiple classification methods, including four traditional classification algorithms and three deep learning methods and their ensembles. The optimal base-classifier combinations are then chosen by five-fold cross-validation test to achieve an effective stacked model. Our model im6APred can produce the area under the receiver operating characteristic curve (AUROC) in the range of 0.82-0.91 on independent tests, indicating that our model has the ability to learn general methylation rules on RNA bases and generalize to m⁶A transcriptome-wide identification. Moreover, AUROCs in the range of 0.77-0.96 were achieved using cross-species/tissues validation on the benchmark dataset, demonstrating differences in predictive performance at the tissue level and the need for constructing tissue-specific models for m⁶A site prediction.