Epigenetics: Roles and therapeutic implications of non-coding RNA modifications in human cancers

Epigenetic deregulation in breast cancer microenvironment: Implications for tumor progression and therapeutic strategies

Breast cancer comprises a substantial proportion of cancer diagnoses in women and is a primary cause of cancer-related mortality. While hormone-responsive cases generally have a favorable prognosis, the aggressive nature of triple-negative breast cancer presents challenges, with intrinsic resistance to established treatments being a persistent issue. The complexity intensifies with the emergence of acquired resistance, further complicating the management of breast cancer. Epigenetic changes, encompassing DNA methylation, histone and RNA modifications, and non-coding RNAs, are acknowledged as crucial contributors to the heterogeneity of breast cancer. The unique epigenetic landscape harbored by each cellular component within the tumor microenvironment (TME) adds great diversity to the intricate regulations which influence therapeutic responses. The TME, a sophisticated ecosystem of cellular and non-cellular elements interacting with tumor cells, establishes an immunosuppressive microenvironment and fuels processes such as tumor growth, angiogenesis, and extracellular matrix remodeling. These factors contribute to challenging conditions in cancer treatment by fostering a hypoxic environment, inducing metabolic stress, and creating physical barriers to drug delivery. This article delves into the complex connections between breast cancer treatment response, underlying epigenetic changes, and vital interactions within the TME. To restore sensitivity to treatment, it emphasizes the need for combination therapies considering epigenetic changes specific to individual members of the TME. Recognizing the pivotal role of epigenetics in drug resistance and comprehending the specificities of breast TME is essential for devising more effective therapeutic strategies. The development of reliable biomarkers for patient stratification will facilitate tailored and precise treatment approaches.

Pulling the trigger: Noncoding RNAs in white adipose tissue browning

White adipose tissue (WAT) serves as the primary site for energy storage and endocrine regulation in mammals, while brown adipose tissue (BAT) is specialized for thermogenesis and energy expenditure. The conversion of white adipocytes to brown-like fat cells, known as browning, has emerged as a promising therapeutic strategy for reversing obesity and its associated co-morbidities. Noncoding RNAs (ncRNAs) are a class of transcripts that do not encode proteins but exert regulatory functions on gene expression at various levels. Recent studies have shed light on the involvement of ncRNAs in adipose tissue development, differentiation, and function. In this review, we aim to summarize the current understanding of ncRNAs in adipose biology, with a focus on their role and intricate mechanisms in WAT browning. Also, we discuss the potential applications and challenges of ncRNA-based therapies for overweight and its metabolic disorders, so as to combat the obesity epidemic in the future.

Identification of molecular subtypes and a prognostic signature based on m6A/m5C/m1A-related genes in lung adenocarcinoma

Lung cancer, specifically the histological subtype lung adenocarcinoma (LUAD), has the highest global occurrence and fatality rate. Extensive research has indicated that RNA alterations encompassing m6A, m5C, and m1A contribute actively to tumorigenesis, drug resistance, and immunotherapy responses in LUAD. Nevertheless, the absence of a dependable predictive model based on m6A/m5C/m1A-associated genes hinders accurately predicting the prognosis of patients diagnosed with LUAD. In this study, we collected patient data from The Cancer Genome Atlas (TCGA) and identified genes related to m6A/m5C/m1A modifications using the GeneCards database. The "ConsensusClusterPlus" R package was used to produce molecular subtypes by utilizing genes relevant to m6A/m5C/m1A identified through differential expression and univariate Cox analyses. An independent prognostic factor was identified by constructing a prognostic signature comprising six genes (SNHG12, PABPC1, IGF2BP1, FOXM1, CBFA2T3, and CASC8). Poor overall survival and elevated expression of human leukocyte antigens and immune checkpoints were correlated with higher risk scores. We examined the associations between the sets of genes regulated by m6A/m5C/m1A and the risk model, as well as the immune cell infiltration, using algorithms such as ESTIMATE, CIBERSORT, TIMER, ssGSEA, and exclusion (TIDE). Moreover, we compared tumor stemness indices (TSIs) by considering the molecular subtypes related to m6A/m5C/m1A and risk signatures. Analyses were performed based on the risk signature, including stratification, somatic mutation analysis, nomogram construction, chemotherapeutic response prediction, and small-molecule drug prediction. In summary, we developed a prognostic signature consisting of six genes that have the potential for prognostication in patients with LUAD and the design of personalized treatments that could provide new versions of personalized management for these patients.

An Overview of the Epigenetic Modifications in the Brain under Normal and Pathological Conditions

Epigenetic changes are changes in gene expression that do not involve alterations to the DNA sequence. These changes lead to establishing a so-called epigenetic code that dictates which and when genes are activated, thus orchestrating gene regulation and playing a central role in development, health, and disease. The brain, being mostly formed by cells that do not undergo a renewal process throughout life, is highly prone to the risk of alterations leading to neuronal death and neurodegenerative disorders, mainly at a late age. Here, we review the main epigenetic modifications that have been described in the brain, with particular attention on those related to the onset of developmental anomalies or neurodegenerative conditions and/or occurring in old age. DNA methylation and several types of histone modifications (acetylation, methylation, phosphorylation, ubiquitination, sumoylation, lactylation, and crotonylation) are major players in these processes. They are directly or indirectly involved in the onset of neurodegeneration in Alzheimer's or Parkinson's disease. Therefore, this review briefly describes the roles of these epigenetic changes in the mechanisms of brain development, maturation, and aging and some of the most important factors dynamically regulating or contributing to these changes, such as oxidative stress, inflammation, and mitochondrial dysfunction.

Research progress of N1-methyladenosine RNA modification in cancer

N1-methyladenosine (m1A) is a post-transcriptionally modified RNA molecule that plays a pivotal role in the regulation of various biological functions and activities. Especially in cancer cell invasion, proliferation and cell cycle regulation. Over recent years, there has been a burgeoning interest in investigating the m1A modification of RNA. Most studies have focused on the regulation of m1A in cancer enrichment areas and different regions. This review provides a comprehensive overview of the methodologies employed for the detection of m1A modification. Furthermore, this review delves into the key players in m1A modification, known as the "writers," "erasers," and "readers." m1A modification is modified by the m1A methyltransferases, or writers, such as TRMT6, TRMT61A, TRMT61B, TRMT10C, NML, and, removed by the demethylases, or erasers, including FTO and ALKBH1, ALKBH3. It is recognized by m1A-binding proteins YTHDF1, TYHDF2, TYHDF3, and TYHDC1, also known as "readers". Additionally, we explore the intricate relationship between m1A modification and its regulators and their implications for the development and progression of specific types of cancer, we discuss how m1A modification can potentially facilitate the discovery of novel approaches for cancer diagnosis, treatment, and prognosis. Our summary of m1A methylated adenosine modification detection methods and regulatory mechanisms in various cancers provides useful insights for cancer diagnosis, treatment, and prognosis. Video Abstract.

Identification and validation of N7 -methylguanosine-associated gene NCBP1 as prognostic and immune-associated biomarkers in breast cancer patients

We intend to evaluate the importance of N7 -methylguanosine (m7G) for the prognosis of breast cancer (BC). We gained 29 m7G-related genes from the published literature and among them, 16 m7G-related genes were found to have differential expression. Five differentially expressed genes (CYFIP1, EIF4E, EIF4E3, NCBP1 and WDR4) were linked to overall survival. This suggests that m7G-related genes might be prognostic or therapeutic targets for BC patients. We put the five genes to LASSO regression analysis to create a four-gene signature, including EIF4E, EIF4E3, WDR4 and NCBP1, that divides samples into two risky groups. Survival was drastically worsened in a high-risk group (p < 0.001). The signature's predictive capacity was demonstrated using ROC (10-year AUC 0.689; 10-year AUC 0.615; 3-year AUC 0.602). We found that immune status was significantly different between the two risk groups. In particular, NCBP1 also has a poor prognosis, with higher diagnostic value in ROC. NCBP1 also has different immune states according to its high or low expression. Meanwhile, knockdown of NCBP1 suppresses BC malignancy in vitro. Therefore, m7G RNA regulators are crucial participants in BC and four-gene mRNA levels are important predictors of prognosis. NCBP1 plays a critical target of m7G mechanism in BC.

Epitranscriptomic modifications in mesenchymal stem cell differentiation: advances, mechanistic insights, and beyond

RNA modifications, known as the "epitranscriptome", represent a key layer of regulation that influences a wide array of biological processes in mesenchymal stem cells (MSCs). These modifications, catalyzed by specific enzymes, often termed "writers", "readers", and "erasers", can dynamically alter the MSCs' transcriptomic landscape, thereby modulating cell differentiation, proliferation, and responses to environmental cues. These enzymes include members of the classes METTL, IGF2BP, WTAP, YTHD, FTO, NAT, and others. Many of these RNA-modifying agents are active during MSC lineage differentiation. This review provides a comprehensive overview of the current understanding of different RNA modifications in MSCs, their roles in regulating stem cell behavior, and their implications in MSC-based therapies. It delves into how RNA modifications impact MSC biology, the functional significance of individual modifications, and the complex interplay among these modifications. We further discuss how these intricate regulatory mechanisms contribute to the functional diversity of MSCs, and how they might be harnessed for therapeutic applications. The review also highlights current challenges and potential future directions in the study of RNA modifications in MSCs, emphasizing the need for innovative tools to precisely map these modifications and decipher their context-specific effects. Collectively, this work paves the way for a deeper understanding of the role of the epitranscriptome in MSC biology, potentially advancing therapeutic strategies in regenerative medicine and MSC-based therapies.

An m1A/m6A/m5C-associated long non-coding RNA signature: Prognostic and immunotherapeutic insights into cervical cancer

BACKGROUND

Cervical cancer (CC) remains a significant clinical challenge, even though its fatality rate has been declining in recent years. Particularly in developing countries, the prognosis for CC patients continues to be suboptimal despite numerous therapeutic advances.

METHODS

Using The Cancer Genome Atlas database, we extracted CC-related data. From this, 52 methylation-related genes (MRGs) were identified, leading to the selection of a 10 long non-coding RNA (lncRNA) signature co-expressed with these MRGs. R programming was employed to filter out the methylation-associated lncRNAs. Through univariate, least absolute shrinkage and selection operator (i.e. LASSO) and multivariate Cox regression analysis, an MRG-associated lncRNA model was constructed. The established risk model was further assessed via the Kaplan-Meier method, principal component analysis, functional enrichment annotation and a nomogram. Furthermore, we explored the potential of this model with respect to guiding immune therapeutic interventions and predicting drug sensitivities.

RESULTS

The derived 10-lncRNA signature, linked with MRGs, emerged as an independent prognostic factor. Segmenting patients based on their immunotherapy responses allowed for enhanced differentiation between patient subsets. Lastly, we highlighted potential compounds for distinguishing CC subtypes.

CONCLUSIONS

The risk model, associated with MRG-linked lncRNA, holds promise in forecasting clinical outcomes and gauging the efficacy of immunotherapies for CC patients.

Deciphering the Epigenetic Symphony of Cancer: Insights and Epigenetic Therapies Implications

Epigenetic machinery is a cornerstone in normal cell development, orchestrating tissue-specific gene expression in mammalian cells. Aberrations in this intricate landscape drive substantial changes in gene function, emerging as a linchpin in cancer etiology and progression. While cancer was conventionally perceived as solely a genetic disorder, its contemporary definition encompasses genetic alterations intertwined with disruptive epigenetic abnormalities. This review explores the profound impact of DNA methylation, histone modifications, and noncoding RNAs on fundamental cellular processes. When these pivotal epigenetic mechanisms undergo disruption, they intricately guide the acquisition of the 6 hallmark characteristics of cancer within seemingly normal cells. Leveraging the latest advancements in decoding these epigenetic intricacies holds immense promise, heralding a new era in developing targeted and more efficacious treatment modalities against cancers driven by aberrant epigenetic modifications.

The role of RNA modification in urological cancers: mechanisms and clinical potential

RNA modification is a post-transcriptional level of regulation that is widely distributed in all types of RNAs, including mRNA, tRNA, rRNA, miRNA, and lncRNA, where N6-methyladenine (m6A) is the most abundant mRNA methylation modification. Significant evidence has depicted that m6A modifications are closely related to human diseases, especially cancer, and play pivotal roles in RNA transcription, splicing, stabilization, and translation processes. The most common urological cancers include prostate, bladder, kidney, and testicular cancers, accounting for a certain proportion of human cancers, with an ever-increasing incidence and mortality. The recurrence, systemic metastasis, poor prognosis, and drug resistance of urologic tumors have prompted the identification of new therapeutic targets and mechanisms. Research on m6A modifications may provide new solutions to the current puzzles. In this review, we provide a comprehensive overview of the key roles played by RNA modifications, especially m6A modifications, in urologic cancers, as well as recent research advances in diagnostics and molecularly targeted therapies.

Differentially expressed circRNAs in peripheral blood samples as potential biomarkers and therapeutic targets for acute angle-closure glaucoma

Glaucoma is the leading cause of irreversible blindness globally. Circular RNAs (circRNAs) play vital roles in various biological processes as microRNA (miRNA) sponges and, thus, have been investigated as potential biomarkers and therapeutic targets in numerous human diseases. However, the underlying mechanisms of circRNAs in the pathogenesis of glaucoma remain unclear. Therefore, transcriptome sequencing was performed to identify relevant circRNAs in peripheral blood samples from patients with primary angle-closure glaucoma. Bioinformatics analysis was performed to investigate the potential roles of differentially expressed circRNAs (DEcircRNAs) in the pathogenesis of glaucoma. In total, 481 differentially expressed genes in addition to 345 DEcircRNAs were identified in patients with glaucoma. Based on a public database, targeted gene analysis identified 11 DEcircRNAs that potentially regulate the expression of five genes as miRNA sponges in glaucoma. In addition, quantitative reverse transcription PCR analysis verified that expression of the circRNA hsa-circ-0000745 was positively correlated with the expression of NEAT1 as a potential target gene. These results suggest that DEcircRNAs are involved in a gene expression regulatory network related to immune cell function and progression of glaucoma. Thus, DEcircRNAs in peripheral blood are potential biomarkers and therapeutic targets for glaucoma.

Establishment of a prognostic model toward lung squamous cell carcinoma based on m7G-related genes in the cancer genome atlas

Lung squamous cell carcinoma (LUSC) is a non-small cell lung cancer with a poor prognosis owing to late diagnosis. New molecular markers are urgently needed to improve the diagnosis and prognosis of LUSC. 7-Methylguanosine (m7G) modifications, a tRNA modification, are common in eubacteria, eukaryotes, and a few archaea. These modifications promote the turnover and stability of some mRNAs to prevent mRNA decay, improve translation efficiency, and reduce ribosomal pausing but are associated with poor survival in human cancer cells. However, expression of m7G-related genes in LUSC and their association with prognosis remain unclear. In the present study, we identified nine differentially expressed genes related to prognosis by comparing the expression profiles of tumor tissues (502 LUSC reports) with normal tissues (49 adjacent nontumor lung tissue reports). The genes included six upregulated genes (KLK7, LCE3E, AREG, KLK6, ZBED2, and MAPK4) and three downregulated genes (ADH1C, NTS, and ERLIN2). Based on these nine genes, patients with LUSC were classified into low- and high-risk groups to analyze the trends in prognosis. We found that the nine m7G-related genes play important roles in immune regulation, hormone regulation, and drug sensitivity through pathways including antigen processing and presentation, adherent plaques, extracellular matrix receptor interactions, drug metabolism of cytochrome P-450, and metabolism of cytochrome P-450 to xenobiotics; the functions of these genes are likely accomplished in part by m6A modifications. The effect of m7G-related genes on the diagnosis and prognosis of LUSC was further indicated by population analysis.NEW & NOTEWORTHY Based on the differential expression of 7-methylguanosine (m7G) modification-associated genes between normal and lung squamous cell carcinoma (LUSC) tissues, and considering the performance of our m7G-related gene risk profiles as independent risk factors in predicting overall survival, we conclude that m7G modification is closely linked to the development of LUSC. In addition, this study offers a new genetic marker for predicting the prognosis of patients with LUSC and presents a crucial theoretical foundation for future investigations on the relationship between m7G modification-related genes, immunity, and drug sensitivity in LUSC.

METTL3 Regulates the m6A Modification of NEK7 to Inhibit the Formation of Osteoarthritis

OBJECTIVE

Osteoarthritis (OA) is a common degenerative joint disease. The occurrence of OA slowly destroys the soft tissue structure of the patient's joint. Severe cases could lead to disability. Current studies had shown that inhibition of chondrocytes pyroptosis could slow down the progression of OA. Our work aimed to explore the specific mechanisms and ways of regulating this process.

DESIGN

In this work, the level of N6-methyladenosine (m6A) in clinical tissues was detected by ribonucleic acid (RNA) m6A dot blot. qRT-PCR (quantitative real-time polymerase chain reaction) was used to detect the messenger RNA (mRNA) expression level of m6A modified enzyme in clinical tissues. MTT (3-(4,5)-dimethylthiahiazo(-z-y1)-3,5-di-phenytetrazoliumromid) and flow cytometry were used to detect the effect of sh-METTL3 (methyltransferase like 3) and NIMA-related kinase 7 (NEK7) transfection on chondrocytes pyroptosis in OA. Western blot was used to detect the protein expression levels of pyroptosis-related proteins. ELISA (enzyme-linked immunosorbent assay) was used to measure the protein concentration of inflammatory cytokines. The SRAMP online database was used to predict the m6A site of NEK7. HE staining was used to assess the progression of OA in mice.

RESULTS

The level of m6A in clinical samples of OA patients was higher, and METTL3 was significantly higher expressed in clinical samples of OA patients. We provided evidence that low expression of METTL3 inhibited chondrocytes pyroptosis. In addition, Rescue experiments and in vivo experiments had shown that METTL3 in combination with NEK7 inhibited the progression of OA by promoting chondrocytes pyroptosis.

CONCLUSIONS

METTL3 regulates m6A modification of NEK7 and inhibits OA progression.

Distinct N7-methylguanosine profiles of circular RNAs in drug-resistant acute myeloid leukemia

Post-transcriptional methylation modifications, such as the N7-methylguanosine (m7G) modification, are increasingly acknowledged for their role in the development and resistance to chemotherapy in acute myeloid leukemia (AML). This study employed MeRIP-seq technology to investigate the m7G sites within circular RNAs (circRNAs) derived from human AML cells and drug-resistant AML cells, in order to identify these sites more comprehensively. In addition, a detailed analysis of the relationship between m7G and drug-resistant AML was conducted. The bioinformatics analysis was utilized to predict the functions of specific methylated transcripts. The findings revealed a significant difference in m7G level between AML cells and drug-resistant AML cells, suggesting a potentially critical role of m7G in circRNAs in drug-resistant AML development. The methylation of M7G could affect the circRNA-miRNA-mRNA co-expression during the development of AML resistance, which could further influence the regulation of resistance-associated target genes in AML. Furthermore, gene ontology analysis indicated that the distinct distribution pattern of circRNAs with m7G methylation in drug-resistant AML cells was correlated with metabolism-related pathways. These results suggested a potential association between drug-resistant AML and m7G methylation of circRNAs. Moreover, the results revealed a novel role of m7G RNA methylation in circRNAs in the progression of AML chemoresistance.

Significance of methylation-related genes in diagnosis and subtype classification of renal interstitial fibrosis

BACKGROUND

RNA methylation modifications, such as N1-methyladenosine/N6-methyladenosine /N5-methylcytosine (m¹A/m⁶A/m⁵C), are the most common RNA modifications and are crucial for a number of biological processes. Nonetheless, the role of RNA methylation modifications of m¹A/m⁶A/m⁵C in the pathogenesis of renal interstitial fibrosis (RIF) remains incompletely understood.

METHODS

Firstly, we downloaded 2 expression datasets from the GEO database, namely GSE22459 and GSE76882. In a differential analysis of these datasets between patients with and without RIF, we selected 33 methylation-related genes (MRGs). We then applied a PPI network, LASSO analysis, SVM-RFE algorithm, and RF algorithm to identify key MRGs.

RESULTS

We eventually obtained five candidate MRGs (WTAP, ALKBH5, YTHDF2, RBMX, and ELAVL1) to forecast the risk of RIF. We created a nomogram model derived from five key MRGs, which revealed that the nomogram model may be advantageous to patients. Based on the selected five significant MRGs, patients with RIF were classified into two MRG patterns using consensus clustering, and the correlation between the five MRGs, the two MRG patterns, and the genetic pattern with immune cell infiltration was shown. Moreover, we conducted GO and KEGG analyses on 768 DEGs between MRG clusters A and B to look into their different involvement in RIF. To measure the MRG patterns, a PCA algorithm was developed to determine MRG scores for each sample. The MRG scores of the patients in cluster B were higher than those in cluster A.

CONCLUSIONS

Ultimately, we concluded that cluster A in the two MRG patterns identified on these five key m¹A/m⁶A/m⁵C regulators may be associated with RIF.

Comprehensive multi-omics analysis reveals m7G-related signature for evaluating prognosis and immunotherapy efficacy in osteosarcoma

Background

Osteosarcoma is one of the most prevalent bone malignancies with a poor prognosis. The N7-methylguanosine (m7G) modification facilitates the modification of RNA structure and function tightly associated with cancer. Nonetheless, there is a lack of joint exploration of the relationship between m7G methylation and immune status in osteosarcoma.

Methods

With the support of TARGET and GEO databases, we performed consensus clustering to characterize molecular subtypes based on m7G regulators in all osteosarcoma patients. The least absolute shrinkage and selection operator (LASSO) method, Cox regression, and receiver operating characteristic (ROC) curves were employed to construct and validate m7G-related prognostic features and derived risk scores. In addition, GSVA, ssGSEA, CIBERSORT, ESTIMATE, and gene set enrichment analysis were conducted to characterize biological pathways and immune landscapes. We explored the relationship between risk scores and drug sensitivity, immune checkpoints, and human leukocyte antigens by correlation analysis. Finally, the roles of EIF4E3 in cell function were verified through external experiments.

Results

Two molecular isoforms based on regulator genes were identified, which presented significant discrepancies in terms of survival and activated pathways. Moreover, the six m7G regulators most associated with prognosis in osteosarcoma patients were identified as independent predictors for the construction of prognostic signature. The model was well stabilized and outperformed traditional clinicopathological features to reliably predict 3-year (AUC = 0.787) and 5-year (AUC = 0.790) survival in osteosarcoma cohorts. Patients with increased risk scores had a poorer prognosis, higher tumor purity, lower checkpoint gene expression, and were in an immunosuppressive microenvironment. Furthermore, enhanced expression of EIF4E3 indicated a favorable prognosis and affected the biological behavior of osteosarcoma cells.

Conclusions

We identified six prognostic relevant m7G modulators that may provide valuable indicators for the estimation of overall survival and the corresponding immune landscape in patients with osteosarcoma.

The joint role of methylation and immune-related lncRNAs in ovarian cancer: Defining molecular subtypes and developing prognostic signature

INTRODUCTION

Complex outcome of ovarian cancer (OC) stems from the tumor immune microenvironment (TIME) influenced by genetic and epigenetic factors. This study aimed to comprehensively explored the subclasses of OC through lncRNAs related to both N6-methyladenosine (m6A)/N1-methyladenosine (m1A)/N7-methylguanosine (m7G)/5-methylcytosine (m5C) in terms of epigenetic variability and immune molecules and develop a new set of risk predictive systems.

MATERIAL AND METHODS

The lncRNA data of OC were collected from TCGA. Spearman correlation analysis on lncRNA data of OC with immune-related gene expression and with m6A/m5C/m1A/m7G were respectively conducted. The m6A/m5C/m1A/m7G-related m6A/m5C/m1A/m7G related immune lncRNA subtypes were identified on the basis of the prognostic lncRNAs. Heterogeneity among subtypes was evaluated by tumor mutation analysis, tumor microenvironment (TME) component analysis, response to immune checkpoint blocked (ICB) and chemotherapeutic drugs. A risk predictive system was developed based on the results of Cox regression analysis and random survival forest analysis of the differences between each specific cluster and other clusters.

RESULTS

Three m6A/m5C/m1A/m7G-related immune lncRNA subtypes of OC showing distinct differences in prognosis, mutation pattern, TIME components, immunotherapy and chemotherapy response were identified. A set of risk predictive system consisting of 10 lncRNA for OC was developed, according to which the risk score of samples in each OC dataset was calculated and risk type was defined.

CONCLUSIONS

This study classified three m6A/m5C/m1A/m7G-related immune lncRNA subtypes with distinct heterogeneous mutation patterns, TME components, ICB therapy and immune response, and provided a set of risk predictive system consisted of 10 lncRNA for OC.

Identification and verification of an ALYREF-involved 5-methylcytosine based signature for stratification of prostate cancer patients and prediction of clinical outcome and response to therapies

OBJECTIVES

Due to the heterogeneity of PCa, the clinical indicators used for PCa can't satisfy risk prognostication and personalized treatment. It is imperative to develop novel biomarkers for prognosis prediction and therapy response in PCa. Accumulating evidence shows that non-mutational epigenetic reprogramming, independent from genomic instability and mutation, serves as a newly added hallmark in cancer progression.

METHODS

In this study, we integrated multi-center cohorts (N > 1300) to develop a RNA 5-methylcytosine regulator-based signature, the m5C score. We performed unsupervised clustering and LASSO regression to identify novel m5C-related subtypes and calculate the m5C score. Then we assessed the role of m5C cluster and m5C score in several clinical aspects such as prognosis in various molecular subtypes, responses to chemotherapy, androgen receptor signaling inhibitor (ARSI) therapy and immunotherapy in PCa. Finally, we validated the cancer-promoting performance of ALYREF through clinical data analysis and experiments in vivo and in vitro.

RESULTS

The investigation revealed that the m5C score could accurately predict the biochemical recurrence (BCR) in different subtypes (the PAM50 subtypes and immunophenotypes) and the responses to chemotherapy, ARSI therapy, and immunotherapy (PD1/PD-L1). A high m5C score indicated a poor BCR prognosis in every subtype of PCa, unfavorable responses in ARSI therapy and immunotherapy (PD1/PD-L1). Moreover, the m5C reader gene termed ALYREF, yielding the highest weighed coefficient, promoted PCa progression through in silico analysis and experimental validations (in vivo and in vitro).

CONCLUSIONS

The m5C signature can function in many aspects of PCa, such as the development and prognosis of the disease, and multiple therapy responses. Further, the m5C reader, ALYREF, was identified as a prognostic biomarker and a potential therapeutic target for PCa. The m5C signature could act as a brand-new tool for predicting the prognosis of patients in different molecular subtypes and patients' therapy responses and promoting customized treatments.

Expression pattern, tumor immune landscape, and prognostic value of N7‑methylguanosine regulators in bladder urothelial carcinoma

N7-Methylguanosine (m7G) modification is important in post-transcriptional regulation. dysregulation of m7G RNA modification has been reported to be markedly associated with cancer. However, its importance in bladder urothelial carcinoma (BLCA) remains poorly characterized. The present study systematically analyzed mRNA gene expression data and clinical information from The Cancer Genome Atlas and further constructed robust risk signatures for the four regulators of m7G RNA modification (nudix hydrolase 11, gem nuclear organelle-associated protein 5, eukaryotic translation initiation factor 3 subunit D and cytoplasmic FMR1 interacting protein 1). The differential expression and cell function of m7G-related genes in bladder cancer cells were verified by reverse transcription-quantitative PCR, Cell Counting Kit-8 and colony formation assays. The four-gene-based model could accurately predict the prognosis of BLCA. Nomogram-based clinical decisions had a higher net benefit compared with that of individual predictors. Through immune infiltration analysis, it was found that immune cell infiltration affected the prognosis of patients with BLCA. Finally, the present study identified potential therapeutics that differ between high and low-risk groups based on four genes. In summary, the current findings revealed an essential role for m7G RNA modification regulators in BLCA, and developed risk signatures as promising prognostic markers in patients with BLCA.

N6-methyladenosine related gene expression signatures for predicting the overall survival and immune responses of patients with colorectal cancer

N6-methyladenosine (m6A) modification has been demonstrated to exhibit a crucial prognostic effect on colorectal cancer (CRC). Nonetheless, potential mechanism of m6A in survival rate and immunotherapeutic response remains unknown. Here we investigated the genes associated with m6A regulators and developed a risk score for predicting the overall survival (OS) of CRC patients. RNA-seq transcriptomic profiling data of COAD/READ samples were obtained from The Cancer Genome Atlas (TCGA) database. Absolute Shrinkage and Selection Operator (LASSO)- Cox regression analysis was conducted to identify the m6A-related gene expression signatures and the selected genes were inputted into stepwise regression to develop a prognostic risk score in TCGA, and its predictive performance of CRC survival was further validated in Gene Expression Omnibus (GEO) datasets. According to our results, the risk score comprising 18 m6A-related mRNAs was significantly associated with CRC survival in both TCGA and GEO datasets. And the stratified analysis also confirmed that high-risk score acted as a poor factor in different age, sex, T stage, and tumour, node, metastasis (TNM) stages. The m6A-related prognostic score in combination with clinical characteristics yielded time-dependent area under the receiver operating characteristic curve (AUCs) of 0.85 (95%CI: 0.79-0.91), 0.84 (95%CI: 0.79-0.90) and 0.80 (95%CI: 0.71-0.88) for the prediction of the 1-, 3-, 5-year OS of CRC in TCGA cohort. Furthermore, mutation of oncogenes occurred more frequently in the high-risk group and the composition of immune cells in tumour microenvironment (TME) was significantly distinct between the low- and high-risk groups. The low-risk group had a lower microsatellite instability (MSI) score, T-cell exclusion score and dysfunction score, implying that low-risk patients may have a better immunotherapy response than high-risk patients. In summary, a prognostic risk score derived from m6A-related gene expression signatures could serve as a potential prognostic predictor for CRC survival and indicator for predicting immunotherapy response in CRC patients.

m7G-related gene NUDT4 as a novel biomarker promoting cancer cell proliferation in lung adenocarcinoma

Background

Lung cancer is the leading cause of mortality in cancer patients. N7-methylguanosine (m7G) modification as a translational regulation pattern has been reported to participate in multiple types of cancer progression, but little is known in lung cancer. This study attempts to explore the role of m7G-related proteins in genetic and epigenetic variations in lung adenocarcinoma, and its relationship with clinical prognosis, immune infiltration, and immunotherapy.

Methods

Sequencing data were obtained from the Genomic Data Commons (GDC) Data Portal and Gene Expression Omnibus (GEO) databases. Consensus clustering was utilized to distinguish m7G clusters, and responses to immunotherapy were also evaluated. Moreover, univariate and multivariate Cox and Least absolute shrinkage and selection operator LASSO Cox regression analyses were used to screen independent prognostic factors and generated risk scores for constructing a survival prediction model. Multiple cell types such as epithelial cells and immune cells were identified to verify the bulk RNA results. Short hairpin RNA (shRNA) Tet-on plasmids, Clustered Regularly Interspaced Short Palindromic Repeats CRISPR/Cas9 for knockout plasmids, and nucleoside diphosphate linked to moiety X-type motif 4 (NUDT4) overexpression plasmids were constructed to inhibit or promote tumor cell NUDT4 expression, then RT-qPCR, Cell Counting Kit-8 CCK8 proliferation assay, and Transwell assay were used to observe tumor cell biological functions.

Results

Fifteen m7G-related genes were highly expressed in tumor samples, and 12 genes were associated with poor prognosis. m7G cluster-B had lower immune infiltration level, worse survival, and samples that predicted poor responses to immunotherapy. The multivariate Cox model showed that NUDT4 and WDR4 (WD repeat domain 4) were independent risk factors. Single-cell m7G gene set variation analysis (GSVA) scores also had a negative correlation tendency with immune infiltration level and T-cell Programmed Death-1 PD-1 expression, but the statistics were not significant. Knocking down and knocking out the NUDT4 expression significantly inhibited cell proliferation capability in A549 and H1299 cells. In contrast, overexpressing NUDT4 promoted tumor cell proliferation. However, there was no difference in migration capability in the knockdown, knockout, or overexpression groups.

Conclusions

Our study revealed that m7G modification-related proteins are closely related to the tumor microenvironment, immune cell infiltration, responses to immunotherapy, and patients' prognosis in lung adenocarcinoma and could be useful biomarkers for the identification of patients who could benefit from immunotherapy. The m7G modification protein NUDT4 may be a novel biomarker in promoting the progression of lung cancer.

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.

Molecular Characterization and Clinical Characteristics of m5C-Based RNA Methylation in Spinal Cord Injury: Validated by qPCR

Aberrant patterns of 5-methylcytosine (m5C)-based ribonucleic acid (RNA) methylation have critical roles in various human diseases, but their importance in spinal cord injury (SCI) is largely unknown. We explore the expression patterns and potential roles of m5C-based regulators of RNA modification after SCI. We analyzed 16 m5C-based regulators of RNA modification in tissues with SCI and normal rats from the Gene Expression Omnibus database. We constructed a "gene signature" of m5C-based regulators of RNA modification to predict the prognosis of SCI using least absolute shrinkage and selection operator regression and random-forest strategy. We found that the m5C-related genes, deoxyribonucleic acid (DNA) methyltransferase1 (Dnmt1), methyl-CpG binding domain protein 2 (Mbd2), ubiquitin-like with PHD and ring finger domains 1 (Uhrf1), uracil-N-glycosylase (Ung), and zinc finger and BTB(brica-brac, tramtrack, and broad) domain containing 38 (Zbtb38) had high expression, and zinc finger and BTB domain containing 4 (Zbtb4) had low expression in SCI. Analysis of the correlation between the gene sets of m5C-based regulators of RNA modification and immune-cell infiltration and immune response revealed Dnmt1, DNA methyltransferases 3A (Dnmt3a), Mbd2, and Ung to be positive regulators of the immune microenvironment, and Zbtb4 may negatively regulate the immune environment. Then, two molecular subtypes were identified based on 16 m5C-regulated genes. Functional-enrichment analysis of differentially expressed genes between different patterns of m5C-based modification was undertaken. Through the creation of a protein-protein interaction network, we screened 11 hub genes. We demonstrated their importance between SCI group and sham group using real-time reverse transcription-quantitative polymerase chain reaction in rat model. Expression of hub genes did not correlate with mitophagy but was positively correlated with endoplasmic reticulum stress (ERS), which suggested that there may be differences in ERS between different patterns of m5C-based modification. This present study explored and discovered the close link between m5C regulators-related genes and SCI. We also hope our findings may contribute to further mechanistic and therapeutic research on the role of key m5C regulators after SCI.

Recent advances in cancer fusion transcript detection

Extensive investigation of gene fusions in cancer has led to the discovery of novel biomarkers and therapeutic targets. To date, most studies have neglected chromosomal rearrangement-independent fusion transcripts and complex fusion structures such as double or triple-hop fusions, and fusion-circRNAs. In this review, we untangle fusion-related terminology and propose a classification system involving both gene and transcript fusions. We highlight the importance of RNA-level fusions and how long-read sequencing approaches can improve detection and characterization. Moreover, we discuss novel bioinformatic tools to identify fusions in long-read sequencing data and strategies to experimentally validate and functionally characterize fusion transcripts.

M7G-related LncRNAs: A comprehensive analysis of the prognosis and immunity in glioma

Today, numerous international researchers have demonstrated that N7-methylguanosine (m7G) related long non-coding RNAs (m7G-related lncRNAs) are closely linked to the happenings and developments of various human beings’ cancers. However, the connection between m7G-related lncRNAs and glioma prognosis has not been investigated. We did this study to look for new potential biomarkers and construct an m7G-related lncRNA prognostic signature for glioma. We identified those lncRNAs associated with DEGs from glioma tissue sequences as m7G-related lncRNAs. First, we used Pearson’s correlation analysis to identify 28 DEGs by glioma and normal brain tissue gene sequences and predicated 657 m7G-related lncRNAs. Then, eight lncRNAs associated with prognosis were obtained and used to construct the m7G risk score model by lasso and Cox regression analysis methods. Furthermore, we used Kaplan-Meier analysis, time-dependent ROC, principal component analysis, clinical variables, independent prognostic analysis, nomograms, calibration curves, and expression levels of lncRNAs to determine the model’s accuracy. Importantly, we validated the model with external and internal validation methods and found it has strong predictive power. Finally, we performed functional enrichment analysis (GSEA, aaGSEA enrichment analyses) and analyzed immune checkpoints, associated pathways, and drug sensitivity based on predictors. In conclusion, we successfully constructed the formula of m7G-related lncRNAs with powerful predictive functions. Our study provides instructional value for analyzing glioma pathogenesis and offers potential research targets for glioma treatment and scientific research.

Scoring model based on the signature of non-m6A-related neoantigen-coding lncRNAs assists in immune microenvironment analysis and TCR-neoantigen pair selection in gliomas

Background

Small peptides encoded by long non-coding RNAs (lncRNAs) have attracted attention for their various functions. Recent studies indicate that these small peptides participate in immune responses and antigen presentation. However, the significance of RNA modifications remains unclear.

Methods

Thirteen non-m6A-related neoantigen-coding lncRNAs were selected for analysis from the TransLnc database. Next, a neoantigen activation score (NAS) model was established based on the characteristics of the lncRNAs. Machine learning was employed to expand the model to two additional RNA-seq and two single-cell sequencing datasets for further validation. The DLpTCR algorithm was used to predict T cell receptor (TCR)-peptide binding probability.

Results

The non-m6A-related NAS model predicted patients’ overall survival outcomes more precisely than the m6A-related NAS model. Furthermore, the non-m6A-related NAS was positively correlated with tumor cells’ evolutionary level, immune infiltration, and antigen presentation. However, high NAS gliomas also showed more PD-L1 expression and high mutation frequencies of T-cell positive regulators. Interestingly, results of intercellular communication analysis suggest that T cell-high neoplastic cell interaction is weaker in both of the NAS groups which might arise from decreased IFNGR1 expression. Moreover, we identified unique TCR-peptide pairs present in all glioma samples based on peptides encoded by the 13 selected lncRNAs. And increased levels of neoantigen-active TCR patterns were found in high NAS gliomas.

Conclusions

Our work suggests that non-m6A-related neoantigen-coding lncRNAs play an essential role in glioma progression and that screened TCR clonotypes might provide potential avenues for chimeric antigen receptor T cell (CAR-T) therapy for gliomas.

A novel m7G methylation–related signature associated with chromosome homeostasis in patients with lung adenocarcinoma

Lung adenocarcinoma (LUAD) is a malignant tumor of the respiratory system with poor prognosis. Recent studies have revealed that N7-methylguanosine (m7G) methylation is a widespread modification occurring in RNA. But the expression of m7G methylation–related genes in LUAD and their correlations with prognosis are still unclear. In this study, we found 12 m7G methylation–related regulators with differential expression between LUAD and normal lung tissues. According to differentially expressed genes (DEGs), all LUAD cases were separated into two subtypes. The prognostic value of each m7G methylation–related gene for survival was evaluated to construct a multigene signature using The Cancer Genome Atlas (TCGA) cohort. Finally, an m7G methylation–related prognostic signature based on three genes was built to classify LUAD patients into two risk groups. Patients in the high-risk group showed significantly reduced overall survival (OS) when compared with patients in the low-risk group (p &lt; 0.05). The receiver operating characteristic (ROC) curve analysis confirmed the predictive capacity of the signature. The Gene Ontology (GO) functional annotation analysis disclosed that chromosome homeostasis plays an important role in this process. The gene set enrichment analysis (ssGSEA) implied that the immune status was decreased in the high-risk group. To sum up, m7G methylation–related genes play a vital role in tumor immunity and the related signature is a reliable predictor for LUAD prognosis.

The role of RNA modification in hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is a highly mortal type of primary liver cancer. Abnormal epigenetic modifications are present in HCC, and RNA modification is dynamic and reversible and is a key post-transcriptional regulator. With the in-depth study of post-transcriptional modifications, RNA modifications are aberrantly expressed in human cancers. Moreover, the regulators of RNA modifications can be used as potential targets for cancer therapy. In RNA modifications, N6-methyladenosine (m6A), N7-methylguanosine (m7G), and 5-methylcytosine (m5C) and their regulators have important regulatory roles in HCC progression and represent potential novel biomarkers for the confirmation of diagnosis and treatment of HCC. This review focuses on RNA modifications in HCC and the roles and mechanisms of m6A, m7G, m5C, N1-methyladenosine (m1A), N3-methylcytosine (m3C), and pseudouridine (ψ) on its development and maintenance. The potential therapeutic strategies of RNA modifications are elaborated for HCC.

The pattern of expression and prognostic value of key regulators for m7G RNA methylation in hepatocellular carcinoma

N7-methylguanosine (m7G) modification on internal RNA positions plays a vital role in several biological processes. Recent research shows m7G modification is associated with multiple cancers. However, in hepatocellular carcinoma (HCC), its implications remain to be determined. In this place, we need to interrogate the mRNA patterns for 29 key regulators of m7G RNA modification and assess their prognostic value in HCC. Initial, the details from The Cancer Genome Atlas (TCGA) database concerning transcribed gene data and clinical information of HCC patients were inspected systematically. Second, according to the mRNA profiles of 29 m7G RNA methylation regulators, two clusters (named 1 and 2, respectively) were identified by consensus clustering. Furthermore, robust risk signature for seven m7G RNA modification regulators was constructed. Last, we used the Gene Expression Omnibus (GEO) dataset to validate the prognostic associations of the seven-gene risk signature. We figured out that 24/29 key regulators of m7G RNA modification varied remarkably in their grades of expression between the HCC and the adjacent tumor control tissues. Cluster one compared with cluster two had a substandard prognosis and was also positively correlated with T classification (T), pathological stage, and vital status (fustat) significantly. Consensus clustering results suggested the expression pattern of m7G RNA modification regulators was correlated with the malignancy of HCC strongly. In addition, cluster one was extensively enriched in metabolic-related pathways. Seven optimal genes (METTL1, WDR4, NSUN2, EIF4E, EIF4E2, NCBP1, and NCBP2) were selected to establish the risk model for HCC. Indicating by further analyses and validation, the prognostic model has fine anticipating command and this probability signature might be a self supporting presage factor for HCC. Finally, a new prognostic nomogram based on age, gender, pathological stage, histological grade, and prospects were established to forecast the prognosis of HCC patients accurately. In essence, we detected association of HCC severity and expression levels of m7G RNA modification regulators, and developed a risk score model for predicting prognosis of HCC patients’ progression.

Regulatory roles of non-coding RNAs and m6A modification in trophoblast functions and the occurrence of its related adverse pregnancy outcomes

Adverse pregnancy outcomes, such as preeclampsia, gestational diabetes mellitus, fetal growth restriction, and recurrent miscarriage, occur frequently in pregnant women and might further induce morbidity and mortality for both mother and fetus. Increasing studies have shown that dysfunctions of human trophoblast are related to these adverse pregnancy outcomes. Recent studies also showed that environmental toxicants could induce trophoblast dysfunctions. Moreover, non-coding RNAs (ncRNAs) have been reported to play important regulatory roles in various cellular processes. However, the roles of ncRNAs in the regulation of trophoblast dysfunctions and the occurrence of adverse pregnancy outcomes still need to be further investigated, especially with exposure to environmental toxicants. In this review, we analyzed the regulatory mechanisms of ncRNAs and m6A methylation modification in the dysfunctions of trophoblast cells and the occurrence of adverse pregnancy outcomes and also summarized the harmful effects of environmental toxicants. In addition to DNA replication, mRNA transcription, and protein translation, ncRNAs and m6A modification might be considered as the fourth and fifth elements that regulate the genetic central dogma, respectively. Environmental toxicants might also affect these processes. In this review, we expect to provide a deeper scientific understanding of the occurrence of adverse pregnancy outcomes and to discover potential biomarkers for the diagnosis and treatment of these outcomes.

Molecular characterization, clinical relevance and immune feature of m7G regulator genes across 33 cancer types

Over 170 RNA modifications have been identified after transcriptions, involving in regulation of RNA splicing, processing, translation and decay. Growing evidence has unmasked the crucial role of N6-methyladenosine (m6A) in cancer development and progression, while, as a relative newly found RNA modification, N7-methylguanosine (m7G) is also certified to participate in tumorigenesis via different catalytic machinery from that of m6A. However, system analysis on m7G RNA modification-related regulator genes is lack. In this study, we first investigated the genetic alteration of m7G related regulator genes in 33 cancers, and found mRNA expression levels of most regulator genes were positively correlated with copy number variation (CNV) and negatively correlated with methylation in most cancers. We built a m7G RNA modification model based on the enrichment of the regulator gene scores to evaluate the m7G modification levels in 33 cancers, and investigated the connections of m7G scores to clinical outcomes. Furthermore, we paid close attention to the role of m7G in immunology due to the widely used immune checkpoint blockade therapy. Our results showed the higher m7G scores related to immunosuppression of tumor cells. Further confirmation with phase 3 clinical data with application of anti-PDL1/PDL indicated the impact of m7G modification level on immunotherapy effect. Relevance of m7G regulator genes and drug sensitivity was also evaluated to provide a better treatment choice when treating cancers. In summary, our study uncovered the profile of m7G RNA modification through various cancers, and figured out the connection of m7G modification levels with therapeutical outcomes, providing potential better options of cancer treatment.

A novel prognostic 7-methylguanosine signature reflects immune microenvironment and alternative splicing in glioma based on multi-omics analysis

Glioma is the most common type of central nervous system tumor with increasing incidence. 7-methylguanosine (m7G) is one of the diverse RNA modifications that is known to regulate RNA metabolism and its dysregulation was associated with various cancers. However, the expression pattern of m7G regulators and their roles in regulating tumor immune microenvironments (TIMEs) as well as alternative splicing events (ASEs) in glioma has not been reported. In this study, we showed that m7G regulators displayed a close correlation with each other and most of them were differentially expressed between normal and glioma tissues. Two m7G signatures were then constructed to predict the overall survival of both GBM and LGG patients with moderate predictive performance. The risk score calculated from the regression coefficient and expression level of signature genes was proved to be an independent prognostic factor for patients with LGG, thus, a nomogram was established on the risk score and other independent clinical parameters to predict the survival probability of LGG patients. We also investigated the correlation of m7G signatures with TIMEs in terms of immune scores, expression levels of HLA and immune checkpoint genes, immune cell composition, and immune-related functions. While exploring the correlation between signature genes and the ASEs in glioma, we found that EIF4E1B was a key regulator and might play dual roles depending on glioma grade. By incorporating spatial transcriptomic data, we found a cluster of cells featured by high expression of PTN exhibited the highest m7G score and may communicate with adjacent cancer cells via SPP1 and PTN signaling pathways. In conclusion, our work brought novel insights into the roles of m7G modification in TIMEs and ASEs in glioma, suggesting that evaluation of m7G in glioma could predict prognosis. Moreover, our data suggested that blocking SPP1 and PTN pathways might be a strategy for combating glioma.

Liver-specific deletion of miR-181ab1 reduces liver tumour progression via upregulation of CBX7

MiR-181 expression levels increased in hepatocellular carcinoma (HCC) compared to non-cancerous tissues. MiR-181 has been widely reported as a possible driver of tumourigenesis but also acts as a tumour suppressor. In addition, the miR-181 family regulates the development and function of immune and vascular cells, which play vital roles in the progression of tumours. More complicatedly, many genes have been identified as miR-181 targets to mediate the effects of miR-181. However, the role of miR-181 in the development of primary tumours remains largely unexplored. We aimed to examine the function of miR-181 and its vital mediators in the progression of diethylnitrosamine-induced primary liver cancers in mice. The size of liver tumours was significantly reduced by 90% in global (GKO) or liver-specific (LKO) 181ab1 knockout mice but not in hematopoietic and endothelial lineage-specific knockout mice, compared to WT mice. In addition, the number of tumours was significantly reduced by 50% in GKO mice. Whole-genome RNA-seq analysis and immunohistochemistry showed that epithelial-mesenchymal transition was partially reversed in GKO tumours compared to WT tumours. The expression of CBX7, a confirmed miR-181 target, was up-regulated in GKO compared to WT tumours. Stable CBX7 expression was achieved with an AAV/Transposase Hybrid-Vector System and up-regulated CBX7 expression inhibited liver tumour progression in WT mice. Hepatic CBX7 deletion restored the progression of LKO liver tumours. MiR-181a expression was the lowest and CBX7 expression the highest in iClust2 and 3 subclasses of human HCC compared to iClust1. Gene expression profiles of GKO tumours overlapped with low-proliferative peri-portal-type HCCs. Liver-specific loss of miR-181ab1 inhibited primary liver tumour progression via up-regulating CBX7 expression, but tumour induction requires both hepatic and non-hepatic miR-181. Also, miR-181ab1-deficient liver tumours may resemble low-proliferative periportal-type human HCC.

The role of N6-methyladenosine-modified non-coding RNAs in the pathological process of human cancer

Non-coding RNAs (ncRNAs) account for the majority of the widespread transcripts of mammalian genomes. They rarely encode proteins and peptides, but their regulatory role is crucial in numerous physiological and pathological processes. The m6A (N6-methyladenosine) modification is one of the most common internal RNA modifications in eukaryotes and is associated with all aspects of RNA metabolism. Accumulating researches have indicated a close association between m6A modification and ncRNAs, and suggested m6A-modified ncRNAs played a crucial role in tumor progression. The correlation between m6A modification and ncRNAs offers a novel perspective for investigating the potential mechanisms of cancer pathological processes, which suggests that both m6A modification and ncRNAs are critical prognostic markers and therapeutic targets in numerous malignancies. In the present report, we summarized the interaction between m6A modification and ncRNA, emphasizing how their interaction regulates pathological processes in cancer.

Methyladenosine Modification in RNAs: From Regulatory Roles to Therapeutic Implications in Cancer

Simple Summary

Cancer remains a burden to the public health all over the world. An increasing number of studies have concentrated on the role of methyladenosine modifications on cancers. Methyladenosine modifications mainly include N6-methyladenosine (m⁶A), N1-methyladenosine (m¹A), and 2’-O-methyladenosine (m⁶A_m), of which dynamic changes could modulate the metabolism of RNAs in eukaryotic cells. Mounting evidence has confirmed the crucial role of methyladenosine modification in cancer, offering possibilities for cancer therapy. In this review, we discussed the regulatory role of methyladenosine modification on cancer, as well as their potential for treatment.

Abstract

Methyladenosine modifications are the most abundant RNA modifications, including N6-methyladenosine (m⁶A), N1-methyladenosine (m¹A), and 2’-O-methyladenosine (m⁶A_m). As reversible epigenetic modifications, methyladenosine modifications in eukaryotic RNAs are not invariable. Drastic alterations of m⁶A are found in a variety of diseases, including cancers. Dynamic changes of m⁶A modification induced by abnormal methyltransferase, demethylases, and readers can regulate cancer progression via interfering with the splicing, localization, translation, and stability of mRNAs. Meanwhile, m⁶A, m¹A, and m⁶A_m modifications also exert regulatory effects on noncoding RNAs in cancer progression. In this paper, we reviewed recent findings concerning the underlying biomechanism of methyladenosine modifications in oncogenesis and metastasis and discussed the therapeutic potential of methyladenosine modifications in cancer treatments.

The Value of m5C-Related lncRNAs in the Prognostic Assessment and Immunotherapy of Stomach Adenocarcinoma

Long noncoding RNAs (lncRNAs) are closely associated with a variety of tumors, including stomach adenocarcinoma (STAD). However, the role of 5-methylcytosine- (m5C-) related lncRNAs in STAD is still uncertain. This study investigated the value of m5C-related lncRNAs in prognostic evaluation and immunotherapy of STAD. STAD transcriptome sequencing data were downloaded from The Cancer Genome Atlas (TCGA) database, and m5C-related lncRNAs were screened by Pearson correlation analysis. A prognostic m5C-related lncRNA signature (m5CRLSig) associated with STAD was established using univariate and multivariate Cox regression analysis. We constructed a prognostic risk model for STAD with six m5C-related lncRNAs. The receiver operating characteristic (ROC) curve was used to examine the predictive efficacy. Univariate and multifactorial Cox regression analysis and principal component analysis (PCA) validated m5CRLSig as an independent factor of STAD prognosis. The clinicopathological characteristics of the model showed higher risk scores for stages II-IV, grade 3, N1-3, and death status. The calibration curve of a nomogram revealed that the nomogram had an excellent predictive function for survival risk. Furthermore, the expression of six m5C-related lncRNAs in STAD and paracancerous tissues was detected by quantitative real-time PCR (qRT-PCR), which verified the feasibility of the m5CRLSig as a prognostic marker for STAD. m5C-related lncRNAs were linked to multiple immune-associated pathways, according to gene set enrichment analysis (GSEA). CIBERSORT analysis indicated that m5CRLSig was involved in immune cell infiltration. Risk score was associated with immune checkpoint gene expression, immune function scores, and chemotherapeutic drug sensitivity. Therefore, m5CRLSig can efficiently assess the prognosis of STAD patients and can be used as a biological marker for immunotherapy.

YTHDF3 Facilitates eIF2AK2 and eIF3A Recruitment on mRNAs to Regulate Translational Processes in Oxaliplatin-Resistant Colorectal Cancer

Oxaliplatin, as a first-line drug, frequently causes chemo-resistance in colorectal cancer (CRC). The role of N6-methyladenosine (m6A) modification in multiple biological functions has been well studied. However, the molecular mechanisms underlying m6A methylation in modulating anti-cancer drug resistance in CRC remain obscure. In the present study, we found that YTH m6A RNA-binding protein 3 (YTHDF3) was highly expressed in oxaliplatin-resistant (OXAR) CRC tissues and cells. Moreover, we observed that YTHDF3 could recognize the 5' untranslated region of significantly m6A-methylated RNAs, which were associated with tumor resistance and recruit eukaryotic translation initiation factor 3 subunit A (eIF3A) to facilitate the translation of these target genes. Furthermore, we determined that eukaryotic translation initiation factor 2 alpha kinase 2 (eIF2AK2) bridged YTHDF3 and eIF3A, enhancing the stability of the YTHDF3/eIF3A complex in OXAR CRC cells. Taken together, our data identified YTHDF3 as a novel hallmark and revealed the molecular mechanism of YTHDF3 on gene translation via coordination with eIF2AK2 in OXAR CRC cells.

The Key Role of RNA Modification in Breast Cancer

The modulation of the function and expression of epigenetic regulators of RNA modification has gradually become the hotspot of cancer research. Studies have shown that alteration of epigenetic modifications can promote the development and metastasis of breast cancer. This review highlights the progress in characterization of the link between RNA modification and the prognosis, carcinogenesis and treatment of breast cancer, which may provide a new theoretical basis for development of effective strategies for monitoring of breast cancer based on epigenetics.

Impact of non-coding RNAs on cancer directed immune therapies: Now then and forever

Accumulating evidence demonstrates that the host genome's epigenetic modifications are essential for living organisms to adapt to extreme conditions. DNA methylation, covalent modifications of histone, and inter-association of non-coding RNAs facilitate the cellular manifestation of epigenetic changes in the genome. Out of various factors involved in the epigenetic programming of the host, non-coding RNAs (ncRNAs) such as microRNA (miRNA), long non-coding RNA (lncRNA), circular RNA, snoRNA and piRNA are new generation non-coding molecules that influence a variety of cellular processes like immunity, cellular differentiation, and tumor development. During tumor development, temporal changes in miRNA/LncRNA rheostat influence sterile inflammatory responses accompanied by the changes in the carcinogenic signaling in the host. At the cellular level, this is manifested by the up-regulation of Inflammasome and inflammatory pathways, which promotes cancer-related inflammation. Given this, we discuss the potential of lncRNAs, miRNAs, circular RNA, snoRNA and piRNA in regulating inflammation and tumor development in the host. This article is protected by copyright. All rights reserved.

Role of main RNA modifications in cancer: N6-methyladenosine, 5-methylcytosine, and pseudouridine

Cancer is one of the major diseases threatening human life and health worldwide. Epigenetic modification refers to heritable changes in the genetic material without any changes in the nucleic acid sequence and results in heritable phenotypic changes. Epigenetic modifications regulate many biological processes, such as growth, aging, and various diseases, including cancer. With the advancement of next-generation sequencing technology, the role of RNA modifications in cancer progression has become increasingly prominent and is a hot spot in scientific research. This review studied several common RNA modifications, such as N⁶-methyladenosine, 5-methylcytosine, and pseudouridine. The deposition and roles of these modifications in coding and noncoding RNAs are summarized in detail. Based on the RNA modification background, this review summarized the expression, function, and underlying molecular mechanism of these modifications and their regulators in cancer and further discussed the role of some existing small-molecule inhibitors. More in-depth studies on RNA modification and cancer are needed to broaden the understanding of epigenetics and cancer diagnosis, treatment, and prognosis.

Role of epigenetic transgenerational inheritance in generational toxicology

Many environmental toxicants have been shown to be associated with the transgenerational inheritance of increased disease susceptibility. This review describes the generational toxicity of some of these chemicals and their role in the induction of epigenetic transgenerational inheritance of disease. Epigenetic factors include DNA methylation, histone modifications, retention of histones in sperm, changes to chromatin structure, and expression of non-coding RNAs. For toxicant-induced epigenetic transgenerational inheritance to occur, exposure to a toxicant must result in epigenetic changes to germ cells (sperm or eggs) since it is the germ cells that carry molecular information to subsequent generations. In addition, the epigenetic changes induced in transgenerational generation animals must cause alterations in gene expression in these animals' somatic cells. In some cases of generational toxicology, negligible changes are seen in the directly exposed generations, but increased disease rates are seen in transgenerational descendants. Governmental policies regulating toxicant exposure should take generational effects into account. A new approach that takes into consideration generational toxicity will be needed to protect our future populations.

Single-molecule long-read sequencing reveals the potential impact of posttranscriptional regulation on gene dosage effects on the avian Z chromosome

Background

Mammalian sex chromosomes provide dosage compensation, but avian lack a global mechanism of dose compensation. Herein, we employed nanopore sequencing to investigate the genetic basis of gene expression and gene dosage effects in avian Z chromosomes at the posttranscriptional level.

Results

In this study, the gonad and head skin of female and male duck samples (n = 4) were collected at 16 weeks of age for Oxford nanopore sequencing. Our results revealed a dosage effect and local regulation of duck Z chromosome gene expression. Additionally, AS and APA achieve tissue-specific gene expression, and male-biased lncRNA regulates its Z-linked target genes, with a positive regulatory role for gene dosage effects on the duck Z chromosome. In addition, GO enrichment and KEGG pathway analysis showed that the dosage effects of Z-linked genes were mainly associated with the cellular response to hormone stimulus, melanin biosynthetic, metabolic pathways, and melanogenesis, resulting in sex differences.

Conclusions

Our data suggested that post transcriptional regulation (AS, APA and lncRNA) has a potential impact on the gene expression effects of avian Z chromosomes. Our study provides a new view of gene regulation underlying the dose effects in avian Z chromosomes at the RNA post transcriptional level.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08360-8.

LncRNA LINC02253 activates KRT18/MAPK/ERK pathway by mediating N6-methyladenosine modification of KRT18 mRNA in gastric cancer

Long non-coding RNAs (lncRNAs) play a crucial role in gastric cancer (GC) progression. And understanding the role of N6-methyladenosine (m6A) in tumorigenesis is an emerging field in cancer research. Here, we identified a novel oncogene, lncRNA LINC02253, in GC. LINC02253 expression was found to be significantly increased in GC. And LINC02253 expression was closely correlated with tumor size, lymph node metastasis and TNM stage of GC. Besides, GC patients with higher LINC02253 expression had worse 5-year overall survival. Additionally, LINC02253 promoted GC cell growth, migration and invasion both in vitro and in vivo. Mechanistically, we determined that LINC02253 increased KRT18 expression through enhancing the stability of KRT18 mRNA. Furthermore, LINC02253 increased m6A modification of KRT18 mRNA to stabilize KRT18 mRNA by recruiting m6A writer METTL3. And, rescue experiments revealed that KRT18 mediated the effects of LINC02253 on growth, migration and invasion of GC cells through activating MAPK/ERK signaling pathway. In conclusion, we demonstrates that oncogenic lncRNA LINC02253 positively regulates GC growth and metastasis via increasing METTL3-mediated mRNA stability of KRT18, extending the understanding of GC pathogenesis regulated by lncRNAs.

Epi-miRNAs: Regulators of the Histone Modification Machinery in Human Cancer

Cancer is a leading cause of death and disability worldwide. Epigenetic deregulation is one of the most critical mechanisms in carcinogenesis and can be classified into effects on DNA methylation and histone modification. MicroRNAs are small noncoding RNAs involved in fine-tuning their target genes after transcription. Various microRNAs control the expression of histone modifiers and are involved in a variety of cancers. Therefore, overexpression or downregulation of microRNAs can alter cell fate and cause malignancies. In this review, we discuss the role of microRNAs in regulating the histone modification machinery in various cancers, with a focus on the histone-modifying enzymes such as acetylases, deacetylases, methyltransferases, demethylases, kinases, phosphatases, desumoylases, ubiquitinases, and deubiquitinases. Understanding of microRNA-related aberrations underlying histone modifiers in pathogenesis of different cancers can help identify novel therapeutic targets or early detection approaches that allow better management of patients or monitoring of treatment response.

Identification and Validation of HOTAIRM1 as a Novel Biomarker for Oral Squamous Cell Carcinoma

ORAL squamous cell carcinoma (OSCC) is a malignant tumor with the highest incidence among tumors involving the oral cavity maxillofacial region, and is notorious for its high recurrence and metastasis potential. Long non-coding RNAs (lncRNAs), which regulate the genesis and evolution of cancers, are potential prognostic biomarkers. This study identified HOTAIRM1 as a novel significantly upregulated lncRNA in OSCC, which is strongly associated with unfavorable prognosis of OSCC. Systematic bioinformatics analyses demonstrated that HOTAIRM1 was closely related to tumor stage, overall survival, genome instability, the tumor cell stemness, the tumor microenvironment, and immunocyte infiltration. Using biological function prediction methods, including Weighted gene co-expression network analysis (WGCNA), Gene set enrichment analysis (GSEA), and Gene set variation analysis (GSVA), HOTAIRM1 plays a pivotal role in OSCC cell proliferation, and is mainly involved in the regulation of the cell cycle. In vitro, cell loss-functional experiments confirmed that HOTAIRM1 knockdown significantly inhibited the proliferation of OSCC cells, and arrested the cell cycle in G1 phase. At the molecular level, PCNA and CyclinD1 were obviously reduced after HOTAIRM1 knockdown. The expression of p53 and p21 was upregulated while CDK4 and CDK6 expression was decreased by HOTAIRM1 knockdown. In vivo, knocking down HOTAIRM1 significantly inhibited tumor growth, including the tumor size, weight, volume, angiogenesis, and hardness, monitored by ultrasonic imaging and magnetic resonance imaging In summary, our study reports that HOTAIRM1 is closely associated with tumorigenesis of OSCC and promotes cell proliferation by regulating cell cycle. HOTAIRM1 could be a potential prognostic biomarker and a therapeutic target for OSCC.

Atranorin driven by nano materials SPION lead to ferroptosis of gastric cancer stem cells by weakening the mRNA 5-hydroxymethylcytidine modification of the Xc-/GPX4 axis and its expression

Gastric cancer is a highly malignant tumor. Gastric cancer stem cells (GCSCs) are the main causes of drug resistance, metastasis, recurrence, and poor prognosis. As a secondary metabolite of lichen, Atranorin has a variety of biological effects, such as antibacterial, anti-inflammatory, analgesic, and wound healing; however, its killing effect on GCSCs has not been reported. In this study, we constructed Atranorin complexes comprising superparamagnetic iron oxide nanoparticles (SPION) (Atranorin@SPION). In vitro and in vivo experiments confirmed that Atranorin@SPION could significantly inhibit the proliferation, invasion, angiogenesis, and tumorigenicity of CD44+/ CD24+ GCSCs, and induce oxidative stress injury, Fe₂+ accumulation, and ferroptosis. Quantitative real-time reverse transcription PCR and western blotting results showed that Atranorin@SPION not only reduced the expression levels of GCSC stem cell markers and cell proliferation and division markers, but also significantly inhibited the expression levels of key molecules in the cystine/glutamate transporter (Xc-)/glutathione peroxidase 4 (GPX4) and Tet methylcytosine dioxygenase (TET) family proteins. The results of high performance liquid chromatography-mass spectrometry and Dot blotting showed that Atranorin@SPION significantly inhibited the mRNA 5‑hydroxymethylcytidine modification of GCSCs. Meanwhile, the results of RNA immunoprecipitation-PCR also indicated that Atranorin@SPIONs significantly reduced the 5-hydroxymethylcytidine modification level of GPX4 and SLC7A11 mRNA 3' untranslated region in GCSCs, resulting in a decrease in their stability, shortening their half-lives and reducing translation activity. Therefore, this study revealed that Atranorin@SPIONs induced ferroptosis of GCSCs by weakening the expression of the Xc-/GPX4 axis and the 5-hydroxymethylcytidine modification of mRNAs in the pathway, thereby achieving their therapeutic effect on gastric cancer.

5-Methylcytosine RNA Methyltransferases-Related Long Non-coding RNA to Develop and Validate Biochemical Recurrence Signature in Prostate Cancer

The effects of 5-methylcytosine in RNA (m5C) in various human cancers have been increasingly studied recently; however, the m5C regulator signature in prostate cancer (PCa) has not been well established yet. In this study, we identified and characterized a series of m5C-related long non-coding RNAs (lncRNAs) in PCa. Univariate Cox regression analysis and least absolute shrinkage and selector operation (LASSO) regression analysis were implemented to construct a m5C-related lncRNA prognostic signature. Consequently, a prognostic m5C-lnc model was established, including 17 lncRNAs: MAFG-AS1, AC012510.1, AC012065.3, AL117332.1, AC132192.2, AP001160.2, AC129510.1, AC084018.2, UBXN10-AS1, AC138956.2, ZNF32-AS2, AC017100.1, AC004943.2, SP2-AS1, Z93930.2, AP001486.2, and LINC01135. The high m5C-lnc score calculated by the model significantly relates to poor biochemical recurrence (BCR)-free survival (p < 0.0001). Receiver operating characteristic (ROC) curves and a decision curve analysis (DCA) further validated the accuracy of the prognostic model. Subsequently, a predictive nomogram combining the prognostic model with clinical features was created, and it exhibited promising predictive efficacy for BCR risk stratification. Next, the competing endogenous RNA (ceRNA) network and lncRNA–protein interaction network were established to explore the potential functions of these 17 lncRNAs mechanically. In addition, functional enrichment analysis revealed that these lncRNAs are involved in many cellular metabolic pathways. Lastly, MAFG-AS1 was selected for experimental validation; it was upregulated in PCa and probably promoted PCa proliferation and invasion in vitro. These results offer some insights into the m5C's effects on PCa and reveal a predictive model with the potential clinical value to improve the prognosis of patients with PCa.