RNA m6A Modification in Cancers: Molecular Mechanisms and Potential Clinical Applications

Multichrome encoding-based multiplexed, spatially resolved imaging reveals single-cell RNA epigenetic modifications heterogeneity

Understanding the heterogeneity of epigenetic modifications within single cells is pivotal for unraveling the nature of the complexity of gene expression and cellular function. In this study, we have developed a strategy based on multichrome encoding and "AND" Boolean logic recognition for multiplexed, spatially resolved imaging of single-cell RNA epigenetic modifications, termed as PRoximity Exchange-assisted Encoding of Multichrome (PREEM). Through the implementation of this strategy, we can now map the expression and nuclear distribution of multiple site-specific RNA N6-methyladenosine (m6A) modifications at the single-molecule resolution level in single-cells, and reveal the previously unknown heterogeneity. Notably, we demonstrate how these patterns change after treatment with various drugs. Moreover, cyclic imaging with tailed DNA self-assembly further suggest the scalability and adaptability of PREEM's design. As an innovative epigenetic modification imaging tool, PREEM not only broadens the horizons of single-cell epigenetics research, enabling joint analysis of multiple targets beyond the limitations of imaging channels, but also reveals cell-to-cell variability, thereby enhancing our capacity to explore cellular functions.

RNA N6-Methyladenosine-Binding Protein YTHDFs Redundantly Attenuate Cancer Immunity by Downregulating IFN-γ Signaling in Gastric Cancer

Immunotherapy holds potential as a treatment for gastric cancer (GC), though immune checkpoint inhibitor (ICI) resistance remains an obstacle. One resistance mechanism involves defects in interferon-γ (IFN-γ) signaling, in which IFN-γ is linked to improved responsiveness to ICIs. Herein, the roles of RNA N6-methyladenosine (m6A) modifications in regulation of IFN-γ signaling and the responsiveness to ICIs are unveiled. The m6A-binding protein YTH N6-methyladenosine RNA-binding protein F1 (YTHDF1) is overexpressed in GC tissues, correlating with the suppression of cancer immunity and poorer survival rates. YTHDF1 overexpression impaired the responsiveness to IFN-γ in GC cells, and knockdown studies indicated the redundant effects of YTHDF2 and YTHDF3 with YTHDF1 in IFN-γ responsiveness. RNA immunoprecipitation sequencing revealed YTHDFs directly target interferon regulatory factor 1 (IRF1) mRNA, a master regulator of IFN-γ signaling, leading to reduced RNA stability and consequent downregulation of IFN-γ signaling. Furthermore, in mouse syngeneic tumor models, Ythdf1 depletion in cancer cells resulted in reduced tumor growth and increased tumor-infiltrating lymphocytes, which are attributed to the augmentation of IFN-γ signaling. Collectively, these findings highlight how YTHDFs modulate cancer immunity by influencing IFN-γ signaling through IRF1 regulation, suggesting their viability as therapeutic targets in cancer immunotherapy.

The FvABF3-FvALKBH10B-FvSEP3 cascade regulates fruit ripening in strawberry

Fruit ripening is a highly-orchestrated process that requires the fine-tuning and precise control of gene expression, which is mainly governed by phytohormones, epigenetic modifiers, and transcription factors. How these intrinsic regulators coordinately modulate the ripening remains elusive. Here we report the identification and characterization of FvALKBH10B as an N6-methyladenosine (m6A) RNA demethylase necessary for the normal ripening of strawberry (Fragaria vesca) fruit. FvALKBH10B is induced by phytohormone abscisic acid (ABA), and ABA-Responsive Element Binding Factor 3 (FvABF3), a master regulator in ABA signaling, is responsible for this activation. FvALKBH10B mutation leads to a delay in fruit ripening and causes global m6A hypermethylation of 1859 genes. Further analyses show that FvALKBH10B positively modulates the mRNA stability of SEPALLATA3 (FvSEP3) encoding a transcription factor via m6A demethylation. In turn, FvSEP3 targets numerous ripening-related genes including those associated with biosynthesis of ABA and anthocyanin and regulates their expression. Our findings uncover an FvABF3-FvALKBH10B-FvSEP3 cascade in controlling fruit ripening in strawberry and provide insights into the complex regulatory networks involved in this process.

Identification of diagnostic biomarkers and molecular subtype analysis associated with m6A in Tuberculosis immunopathology using machine learning

Tuberculosis (TB), ranking just below COVID-19 in global mortality, is a highly complex infectious disease involving intricate immunological molecules, diverse signaling pathways, and multifaceted immune processes. N6-methyladenosine (m6A), a critical epigenetic modification, regulates various immune-metabolic and pathological pathways, though its precise role in TB pathogenesis remains largely unexplored. This study aims to identify m6A-associated genes implicated in TB, elucidate their mechanistic contributions, and evaluate their potential as diagnostic biomarkers and tools for molecular subtyping. Using TB-related datasets from the GEO database, this study identified differentially expressed genes associated with m6A modification. We applied four machine learning algorithms-Random Forest, Support Vector Machine, Extreme Gradient Boosting, and Generalized Linear Model-to construct diagnostic models focusing on m6A regulatory genes. The Random Forest algorithm was selected as the optimal model based on performance metrics (area under the curve [AUC] = 1.0, p < 0.01), and a clinical predictive model was developed based on these critical genes. Patients were stratified into distinct subtypes according to m6A gene expression profiles, followed by immune infiltration analysis across subtypes. Additionally, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses elucidated the biological functions and pathways associated with the identified genes. Quantitative real-time PCR (RT-qPCR) was used to validate the expression of key m6A regulatory genes. Analysis of the GSE83456 dataset revealed four differentially expressed m6A-related genes-YTHDF1, HNRNPC, LRPPRC, and ELAVL1-identified as critical m6A regulators in TB through the Random Forest model. The diagnostic significance of these genes was further supported by a nomogram, achieving a high predictive accuracy (95% confidence interval [CI]: 0.87-0.94). Consensus clustering classified patients into two m6A subtypes with distinct immune profiles, as principal component analysis (PCA) showed significantly higher m6A scores in Group A than in Group B (p < 0.05). Immune infiltration analysis highlighted significant correlations between key m6A genes and specific immune cell infiltration patterns across subtypes. This study highlights the potential of key m6A regulatory genes as diagnostic biomarkers and immunotherapy targets for TB, supporting their role in TB pathogenesis. Future research should aim to further validate these findings across diverse cohorts to enhance their clinical applicability.

Development of a novel disulfidptosis-correlated m6A/m1A/m5C/m7G gene signature to predict prognosis and therapeutic response for lung adenocarcinoma patients by integrated machine-learning

BACKGROUND

Lung adenocarcinoma (LUAD) represents a significant global health burden, necessitating advanced prognostic tools for improved patient management. RNA modifications (m6A, m1A, m5C, m7G), and disulfidptosis, a novel cell death mechanism, have emerged as promising biomarkers and therapeutic targets in cancer.

METHODS

We systematically compiled disulfidptosis-correlated genes and RNA modification-related genes from existing literature. A novel disulfidptosis-correlated m6A/m1A/m5C/m7G riskscore was computed using integrated machine-learning algorithms. Transcriptomic data from TCGA and GEO databases were downloaded analyzed. Single-cell RNA-sequencing data from the TISCH database was processed using the Seurat package. Genes' protein-protein interaction network was constructed using the String database. Functional phenotype analysis was performed using GSVA, ClusterProfiler, and IOBR packages. Consensus clustering divided patients into two distinct groups. Drug sensitivity predictions were obtained from the GDSC1 database and predicted using the Oncopredict package.

RESULTS

The disulfidptosis-correlated m6A/m1A/m5C/m7G risk score effectively stratified LUAD patients into prognostically distinct groups, demonstrating superior predictive accuracy compared to conventional clinical parameters. Patients in different risk groups exhibited significant molecular and clinical differences. Subsequent analyses identified two molecular subtypes associated with RNA modification and disulfidptosis, revealing differences in immune infiltration and prognosis. Functional enrichment analyses highlighted pathways involving RNA modification and disulfidptosis, underscoring their roles in LUAD pathogenesis. Single-cell analysis revealed distinct features between high- and low-risk status cells.

CONCLUSION

This study introduces a novel disulfidptosis-correlated m6A/m1A/m5C/m7G risk score as a robust prognostic tool for LUAD, integrating insights from RNA modifications and cell death mechanisms. The risk score enhances prognostic stratification and identifies potential targets for personalized therapeutic strategies in LUAD. This comprehensive approach emphasizes the critical roles of RNA modifications and disulfidptosis in LUAD biology, paving the way for future research and clinical applications aimed at improving patient outcomes.

RNA modification in normal hematopoiesis and hematologic malignancies

N6-methyladenosine (m6A) is the most abundant RNA modification in eukaryotic cells. Previous studies have shown that m6A plays a critical role under both normal physiological and pathological conditions. Hematopoiesis and differentiation are highly regulated processes, and recent studies on m6A mRNA methylation have revealed how this modification controls cell fate in both normal and malignant hematopoietic states. However, despite these insights, a comprehensive understanding of its complex roles between normal hematopoietic development and malignant hematopoietic diseases remains elusive. This review first provides an overview of the components and biological functions of m6A modification regulators. Additionally, it highlights the origin, differentiation process, biological characteristics, and regulatory mechanisms of hematopoietic stem cells, as well as the features, immune properties, and self-renewal pathways of leukemia stem cells. Last, the article systematically reviews the latest research advancements on the roles and mechanisms of m6A regulatory factors in normal hematopoiesis and related malignant diseases. More importantly, this review explores how targeting m6A regulators and various signaling pathways could effectively intervene in the development of leukemia, providing new insights and potential therapeutic targets. Targeting m6A modification may hold promise for achieving more precise and effective leukemia treatments.

Identification of potential novel N6-methyladenosine effector-related lncRNA biomarkers for serous ovarian carcinoma: a machine learning-based exploration in the framework of 3P medicine

Background

Serous ovarian carcinoma (SOC) is considered the most lethal gynecological malignancy. The current lack of reliable prognostic biomarkers for SOC reduces the efficacy of predictive, preventive, and personalized medicine (PPPM/3PM) in patients with SOC, leading to unsatisfactory therapeutic outcomes. N6-methyladenosine (m6A) modification-associated long noncoding RNAs (lncRNAs) are effective predictors of SOC. In this study, an effective risk prediction model for SOC was constructed based on m6A modification-associated lncRNAs.

Methods

Transcriptomic data and clinical information of patients with SOC were downloaded from The Cancer Genome Atlas. Candidate lncRNAs were identified using univariate and multivariate and least absolute shrinkage and selection operator-penalized Cox regression analyses. The molecular mechanisms of m6A effector-related lncRNAs were explored via Gene Ontology, pathway analysis, gene set enrichment analysis, and gene set variation analysis (GSVA). The extent of immune cell infiltration was assessed using various algorithms, including CIBERSORT, Microenvironment Cell Populations counter, xCell, European Prospective Investigation into Cancer and Nutrition, and GSVA. The calcPhenotype algorithm was used to predict responses to the drugs commonly used in ovarian carcinoma therapy. In vitro experiments, such as migration and invasion Transwell assays, wound healing assays, and dot blot assays, were conducted to elucidate the functional roles of candidate lncRNAs.

Results

Six m6A effector-related lncRNAs that were markedly associated with prognosis were used to establish an m6A effector-related lncRNA risk model (m6A-LRM) for SOC. Immune microenvironment analysis suggested that the high-risk group exhibited a proinflammatory state and displayed increased sensitivity to immunotherapy. A nomogram was constructed with the m6A effector-related lncRNAs to assess the prognostic value of the model. Sixteen drugs potentially targeting m6A effector-related lncRNAs were identified. Furthermore, we developed an online web application for clinicians and researchers (https://leley.shinyapps.io/OC_m6A_lnc/). Overexpression of the lncRNA RP11-508M8.1 promoted SOC cell migration and invasion. METTL3 is an upstream regulator of RP11-508M8.1. The preliminary regulatory axis METTL3/m6A/RP11-508M8.1/hsa-miR-1270/ARSD underlying SOC was identified via a combination of in vitro and bioinformatic analyses.

Conclusion

In this study, we propose an innovative prognostic risk model and provide novel insights into the mechanism underlying the role of m6A-related lncRNAs in SOC. Incorporating the m6A-LRM into PPPM may help identify high-risk patients and personalize treatment as early as possible.

RNA modification-related genes illuminate prognostic signature and mechanism in esophageal squamous cell carcinoma

Emerging studies have demonstrated the link between RNA modifications and various cancers, while the predictive value and functional mechanisms of RNA modification-related genes (RMGs) in esophageal squamous cell carcinoma (ESCC) remain unclear. Here we established a prognostic signature for ESCC based on five RMGs. The analysis of ESCC clinical samples further verified the prognostic power of the prognostic signature. Moreover, we found that the knockdown of NSUN6 promotes ESCC progression in vitro and in vivo, whereas the overexpression of NSUN6 inhibits the malignant phenotype of ESCC cells. Mechanically, NSUN6 mediated tRNA m5C modifications selectively enhance the translation efficiency of CDH1 mRNA in a codon dependent manner. Rescue assays revealed that E-cadherin is an essential downstream target that mediates NSUN6's function in the regulation of ESCC progression. These findings offer additional insights into the link between ESCC and RMGs, as well as provide potential strategies for ESCC management and therapy.

Non-coding RNA-Mediated N6-Methyladenosine (m6A) deposition: A pivotal regulator of cancer, impacting key signaling pathways in carcinogenesis and therapy response

The emergence of RNA modifications has recently been considered as critical post-transcriptional regulations which governed gene expression. N6-methyladenosine (m6A) modification is the most abundant type of RNA modification which is mediated by three distinct classes of proteins called m6A writers, readers, and erasers. Accumulating evidence has been made in understanding the role of m6A modification of non-coding RNAs (ncRNAs) in cancer. Importantly, aberrant expression of ncRNAs and m6A regulators has been elucidated in various cancers. As the key role of ncRNAs in regulation of cancer hallmarks is well accepted now, it could be accepted that m6A modification of ncRNAs could affect cancer progression. The present review intended to discuss the latest knowledge and importance of m6A epigenetic regulation of ncRNAs including mircoRNAs, long non-coding RNAs, and circular RNAs, and their interaction in the context of cancer. Moreover, the current insight into the underlying mechanisms of therapy resistance and also immune response and escape mediated by m6A regulators and ncRNAs are discussed.

Comprehensive analysis of m6A RNA methylation regulators in esophageal carcinoma

Background

N6-methyladenosine (m6A) is the most pervasive modification of RNA methylation in eukaryotic cells. m6A modification plays a pivotal role in tumorigenesis and progression in many types of cancers. Until now, the role of m6A modification in esophageal carcinoma (ESCA) has remained obscure. The aim of the study was to construct and validate prognostic signatures based on m6A regulators for ESCA.

Methods

Transcriptomic data, somatic mutations and clinical information were obtained from The Cancer Genome Atlas (TCGA). Copy number variations were obtained from the UCSC (University of California, Santa Cruz) Xena database. We curated 21 m6A regulators and performed consensus clustering analysis to quantify the m6A modification pattern.

Results

Of the 184 patients, 23 (12.5%) were genetically altered in m6A regulators, with the highest frequency of mutations in ZC3H13 and LRPPRC. We constructed a m6A score system to investigate the prognosis of ESCA. The m6A score was closely related to immune cell infiltration in the tumor immune microenvironment. Patients with a high m6A score had an unfavorable prognosis. The combination of tumor mutation burden and m6A score would improve the prognostic value.

Conclusions

Our study established and validated a strong prognostic signature based on m6A regulators. This can be used to accurately predict the prognosis of ESCA.

IGF2BP2 promotes cell invasion and epithelial-mesenchymal transition through Src-mediated upregulation of EREG in oral cancer

Insulin-like growth factor 2 mRNA binding protein 2 (IGF2BP2), with high affinity to a myriad of RNA transcripts, has been shown to elicit promotive effects on tumorigenesis and metastasis. Yet, the functional involvement of IGF2BP2 in the progression of oral squamous cell carcinoma (OSCC) remains poorly understood. In this study, we showed that IGF2BP2 was upregulated in head and neck cancer, and high levels of IGF2BP2 were associated with poor survival. In in vitro experiments, IGF2BP2 promoted migration and invasion responses of OSCC cells. Moreover, we identified an IGF2BP2-regulated gene, EREG, which functioned as a modulator of OSCC invasion downstream of IGF2BP2. In addition, EREG expression triggered the epithelia-mesenchymal transition (EMT) in OSCC, as evidenced by the observation that knockdown of EREG weakened the induction of EMT mediated by IFG2BP2, and replenishment of EREG favored the EMT in IGF2BP2-depleted cells. Such IGF2BP2-regulated EREG expression, EMT, and cell invasion were dependent on the activation of FAK/Src signaling pathway. Collectively, these findings suggest that EREG, serving as a functional mediator of IGF2BP2-regulated EMT and cell invasion in oral cancer, may be implicated as a potential target for antimetastatic therapies.

Identification and validation of N6-methyladenosine (m6A)-related lncRNAs signature for predicting the prognosis of laryngeal carcinoma, especially for smoking patients

Laryngeal cancer (LC), a highly fatal tumor in the head and neck region, has been the focus of research in recent years. The study of LC has primarily focused on the role of long non-coding RNAs (lncRNAs) in regulating gene expression, as they have emerged as pivotal factors in this biological process. Additionally, a reversible RNA modification called N6-methyladenosine (m6A) has been observed to have a significant impact on gene expression as well. The purpose of this research is to investigate the impact of m6A-related lncRNAs on the prognosis of laryngeal squamous cell carcinoma (LSCC). Specifically, this investigation analyzed the m6A-related regulators' patterns of expression and mutation, encompassing a total of 15 regulators. Drawing upon the expression levels of prognostic m6A-regulated lncRNAs, two distinct lncRNA clusters were identified. Further analysis revealed differentially expressed lncRNAs between these clusters. In addition to studying the expression of lncRNAs, the researchers also examined the distribution of clinical characteristics and the tumor microenvironment (TME) in relation to the identified lncRNA clusters. This provided valuable insights into potential associations between lncRNA expression patterns and the clinical features of LSCC. Through the establishment of a risk model associated with lncRNAs, we were able to further investigate their clinical features, prognosis, and immune status. Additionally, we conducted a separate analysis of LINC00528, a lncRNA associated with smoking, examining its expression, overall survival time, correlated mRNAs, and conducting enrichment of Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG), as well as determining the sensitivity of related drugs. RT-qPCR results also indicated an increase in LINC00528 expression among smoking LSCC patients. The findings suggest that a high expression level of LINC00528 in LSCC patients may lead to a more favorable prognosis, providing new insights for the management and treatment of LSCC patients, particularly those with high expression of LINC00528. Overall, this research sheds light on the prognostic impact of m6A-regulated lncRNAs in LSCC. The implications of these findings for the advancement of innovative therapeutic approaches for LSCC patients are noteworthy.

Comprehensive analysis of m6A regulators associated with immune infiltration in Hepatitis B virus-related hepatocellular carcinoma

BACKGROUND

N6A methylation (m6A) is a significant epigenetic modification that critically impacts post-transcriptional regulation and tumor occurrence and development. While previous studies have identified a role for epigenetic regulation in hepatocellular carcinoma (HCC), the potential function of the m6A cluster in Hepatitis B virus (HBV)-related HCC remains unclear.

METHODS

The related information was downloaded from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO). Based on the expression of 20 m6A regulators, we comprehensively evaluated the m6A clusters and systematically explored the correlation between these clusters and immune cell infiltration characteristics of the tumor microenvironment (TME). The patients were divided into low- and high-m6A score groups. Then, the immune cell infiltration, chemokines, and cytokines levels, and drug sensitivity were further explored between the two groups.

RESULTS

The m6A cluster predicted a better prognosis that was accompanied by increased immune cell infiltration. Using these results, an m6A score was established that could predict overall survival, immune checkpoints, and clinical treatments for patients with HBV-related HCC. This study demonstrated that m6A modifications affected tumorigenesis, TME, and the prognosis of patients with HBV-related HCC.

CONCLUSION

A comprehensive assessment of m6A patterns could improve the current understanding of immune cell infiltration patterns and inform the development of individualized cancer treatments.

CRISPR-based m6A modification and its potential applications in telomerase regulation

Telomerase determines cell lifespan by controlling chromosome stability and cell viability, m⁶A epigenetic modification plays an important role in the regulation of telomerase activity. Using CRISPR epigenome editing to analyze specific m⁶A modification sites in telomerase will provide an important tool for analyzing the molecular mechanism of m⁶A modification regulating telomerase activity. In this review, we clarified the relevant applications of CRISPR system, paid special attention to the regulation of m⁶A modification in stem cells and cancer cells based on CRISPR system, emphasized the regulation of m⁶A modification on telomerase activity, pointed out that m⁶A modification sites regulate telomerase activity, and discussed strategies based on telomerase activity and disease treatment, which are helpful to promote the research of anti-aging and tumor related diseases.

Underexplored reciprocity between genome-wide methylation status and long non-coding RNA expression reflected in breast cancer research: potential impacts for the disease management in the framework of 3P medicine

Breast cancer (BC) is the most common female malignancy reaching a pandemic scale worldwide. A comprehensive interplay between genetic alterations and shifted epigenetic regions synergistically leads to disease development and progression into metastatic BC. DNA and histones methylations, as the most studied epigenetic modifications, represent frequent and early events in the process of carcinogenesis. To this end, long non-coding RNAs (lncRNAs) are recognized as potent epigenetic modulators in pathomechanisms of BC by contributing to the regulation of DNA, RNA, and histones' methylation. In turn, the methylation status of DNA, RNA, and histones can affect the level of lncRNAs expression demonstrating the reciprocity of mechanisms involved. Furthermore, lncRNAs might undergo methylation in response to actual medical conditions such as tumor development and treated malignancies. The reciprocity between genome-wide methylation status and long non-coding RNA expression levels in BC remains largely unexplored. Since the bio/medical research in the area is, per evidence, strongly fragmented, the relevance of this reciprocity for BC development and progression has not yet been systematically analyzed. Contextually, the article aims at:consolidating the accumulated knowledge on both-the genome-wide methylation status and corresponding lncRNA expression patterns in BC andhighlighting the potential benefits of this consolidated multi-professional approach for advanced BC management. Based on a big data analysis and machine learning for individualized data interpretation, the proposed approach demonstrates a great potential to promote predictive diagnostics and targeted prevention in the cost-effective primary healthcare (sub-optimal health conditions and protection against the health-to-disease transition) as well as advanced treatment algorithms tailored to the individualized patient profiles in secondary BC care (effective protection against metastatic disease). Clinically relevant examples are provided, including mitochondrial health control and epigenetic regulatory mechanisms involved.

RNA N6-methyladenosine modification in regulating cancer stem cells and tumor immune microenvironment and its implication for cancer therapy

Therapy resistance is a well-known phenomenon in cancer treatment. It can be intrinsic or acquired, accountable for frequent tumor relapse and death worldwide. The interplay between cancer cells and their neighboring environment can activate complex signaling mechanisms influencing epigenetic changes and maintain cancer cell survival leading to the malignant phenotype. Cancer stem cells (CSCs) are tumor-initiating cells (TICs) and constitute the primary source of drug resistance and tumor recurrence. Studies have shown that cancer cells exhibit dysregulated RNA N6-methyladenosine (m6A) "writers," "erasers," and "readers" levels after acquiring drug resistance. The present review provides novel insight into the role of m6A modifiers involved in CSC generation, cancer cell proliferation, and therapy resistance. m6A RNA modifications in the cross-talk between CSC and the tumor immune microenvironment (TIME) have also been highlighted. Further, we have discussed the therapeutic potential of targeting m6A machinery for cancer diagnosis and the development of new therapies for cancer treatment.

METTL3 protects METTL14 from STUB1-mediated degradation to maintain m6 A homeostasis

N6 -Methyladenosine (m6 A) is an important RNA modification catalyzed by methyltransferase-like 3 (METTL3) and METTL14. m6 A homeostasis mediated by the methyltransferase (MTase) complex plays key roles in various biological processes. However, the mechanism underlying METTL14 protein stability and its role in m6 A homeostasis remain elusive. Here, we show that METTL14 stability is regulated by the competitive interaction of METTL3 with the E3 ligase STUB1. STUB1 directly interacts with METTL14 to mediate its ubiquitination at lysine residues K148, K156, and K162 for subsequent degradation, resulting in a significant decrease in total m6 A levels. The amino acid regions 450-454 and 464-480 of METTL3 are essential to promote METTL14 stabilization. Changes in STUB1 expression affect METTL14 protein levels, m6 A modification and tumorigenesis. Collectively, our findings uncover an ubiquitination mechanism controlling METTL14 protein levels to fine-tune m6 A homeostasis. Finally, we present evidence that modulating STUB1 expression to degrade METTL14 could represent a promising therapeutic strategy against cancer.

A pan-cancer landscape of IGF2BPs and their association with prognosis, stemness and tumor immune microenvironment

Background

The human insulin-like growth factor 2 mRNA binding proteins 1-3 (IGF2BP1-3, also called IMP1-3) play essential roles in mRNA regulation, including its splicing, translocation, stability, and translation. However, knowledge regarding the involvement of IGF2BPs in tumor immunity and stemness across cancer types is still lacking.

Methods

In this study, we comprehensively analyzed pan-cancer multi-omic data to determine the correlation of IGF2BPs mRNA and protein expression with various cancer parameters such as mutation frequency, prognostic value, the tumor microenvironment (TME), checkpoint blockade, tumor immune infiltration, stemness and drug sensitivity. Validation of the expression of IGF2BPs in cancer samples and glioma cells were performed by quantitative real-time (qRT)-PCR, and immunofluorescence staining. Investigation of the functional role of IGF2BP3 in glioma stem cells(GSCs) were performed by sphere formation, cytotoxicity, transwell, and wound healing assays.

Results

We found that IGF2BP1 and 3 are either absent or expressed at very low levels in most normal tissues. However, IGF2BP1-3 can be re-expressed in a broad range of cancer types and diverse cancer cell lines, where their expression often correlates with poor prognosis. Immunofluorescence staining and qRT-PCR analyses also showed that the expression of IGF2BP2 and IGF2BP3 were higher in cancer tissues than that in adjacent normal tissues. Moreover, IGF2BPs are associated with TME and stemness in human pan-cancer. Remarkably, IGF2BP3 participated in the maintenance and self-renewal of glioma stem cell (GSCs). Knockdown of IGF2BP3 attenuated GSC and glioma cell proliferation, invasion, and migration.

Conclusions

Our systematic pan-cancer study confirmed the identification of IGF2BPs as therapeutic targets and highlighted the need to study their association with stemness, and the TME, which contribute to the cancer drug-discovery research. Especially, preliminary studies demonstrate the IGF2BP3 as a potential negative regulator of glioma tumorigenesis by modulating stemness.

Role of Circulating Tumor DNA in Colorectal Cancer

Colorectal cancer (CRC) is a very common gastrointestinal tumor, ranking second in the global cause of cancer death. Because of the invasive nature of biopsy and cannot reflect the heterogeneity of tumor or monitor the dynamic progress of tumor, it is necessary to induce a novel noninvasive method to improve the current treatment strategies of colorectal cancer. Among all the components of liquid biopsy, circulating tumor DNA (ctDNA) may have the best future. CtDNA maintains the same genomic characteristics as those in matched tumor tissues, so it allows quantitative evaluation and analysis of mutation load in body fluid. Furthermore, because the half-life of ctDNA is from 16 min to several hours in circulation, the circulating ctDNA can be measured repeatedly within a certain period to monitor the response of CRC to treatment, the occurrence of drug resistance, and the diagnosis of recurrence.

Epigenetic alterations in glioblastomas: Diagnostic, prognostic and therapeutic relevance

Glioblastoma, the most common and heterogeneous tumor affecting brain parenchyma, is dismally characterized by a very poor prognosis. Thus, the search of new, more effective treatments is a vital need. Here, we will review the druggable epigenetic features of glioblastomas that are, indeed, currently explored in preclinical studies and in clinical trials for the development of more effective, personalized treatments. In detail, we will review the studies that have led to the identification of epigenetic signatures, IDH mutations, MGMT gene methylation, histone modification alterations, H3K27 mutations and epitranscriptome landscapes of glioblastomas, in each case discussing the corresponding targeted therapies and their potential efficacy. Finally, we will emphasize how recent technological improvements permit to routinely investigate many glioblastoma epigenetic biomarkers in clinical practice, further enforcing the hope that personalized drugs, targeting specific epigenetic features, could be in future a therapeutic option for selected patients.

Epigenetic Regulation in Chromium-, Nickel- and Cadmium-Induced Carcinogenesis

Environmental and occupational exposure to heavy metals, such as hexavalent chromium, nickel, and cadmium, are major health concerns worldwide. Some heavy metals are well-documented human carcinogens. Multiple mechanisms, including DNA damage, dysregulated gene expression, and aberrant cancer-related signaling, have been shown to contribute to metal-induced carcinogenesis. However, the molecular mechanisms accounting for heavy metal-induced carcinogenesis and angiogenesis are still not fully understood. In recent years, an increasing number of studies have indicated that in addition to genotoxicity and genetic mutations, epigenetic mechanisms play critical roles in metal-induced cancers. Epigenetics refers to the reversible modification of genomes without changing DNA sequences; epigenetic modifications generally involve DNA methylation, histone modification, chromatin remodeling, and non-coding RNAs. Epigenetic regulation is essential for maintaining normal gene expression patterns; the disruption of epigenetic modifications may lead to altered cellular function and even malignant transformation. Therefore, aberrant epigenetic modifications are widely involved in metal-induced cancer formation, development, and angiogenesis. Notably, the role of epigenetic mechanisms in heavy metal-induced carcinogenesis and angiogenesis remains largely unknown, and further studies are urgently required. In this review, we highlight the current advances in understanding the roles of epigenetic mechanisms in heavy metal-induced carcinogenesis, cancer progression, and angiogenesis.

Characterization of the Immune Microenvironmental Landscape of Lung Squamous Cell Carcinoma with Immune Cell Infiltration

BACKGROUND

Increasing evidence supports that immune cell infiltration (ICI) patterns play a key role in the tumor progression of lung squamous cell carcinoma (LUSC). However, to date, the immune infiltration picture of LUSC has not been elucidated.

METHOD

TCGA was used to download multiomics data from LUSC samples. At the same time, we included two datasets on lung squamous cell carcinoma, GSE17710 and GSE157010. To reveal the landscape of tumor immune microenvironment (TIME), the ESTIMATE algorithm, ssGSEA approach, and CIBERSORT analysis are used. To quantify the ICI pattern in a single tumor, consistent clustering is used to determine the LUSC subtype based on the ICI pattern, and principal component analysis (PCA) is used to obtain the ICI score. The prognostic value of the Kaplan-Meier curves is confirmed. GSEA (Gene Set Enrichment Analysis) was used to perform functional annotation. To investigate the immunotherapeutic effects of the ICI score, the immunophenotyping score (IPS) is used. Finally, analyze the mutation data with the "maftools" R package.

RESULTS

We identified four different immune infiltration patterns with different prognosis and biological characteristics in 792 LUSC samples. The identification of ICI patterns in individual tumors developed under ICI-related characteristic genes based on the ICI score helps to analyze the biological process, clinical results, immune cell infiltration, immunotherapy effects, and genetic variation. Immune failure is indicated by a high ICI score subtype marked by immunosuppression. Patients with low ICI scores have an abundance of efficient immune cells, which corresponds to the immunological activation phenotype and may have therapeutic benefits. The immunophenotypic score was used as a surrogate indicator of immunotherapy results, and samples with low ICI scores obtained significantly higher immunophenotypic scores. Finally, the relationship between the ICI score and tumor mutation burden (TMB) was proven.

CONCLUSION

This study fully clarified the indispensable role of the ICI model in the complexity and diversity of TIME. The quantitative identification of ICI patterns in a single tumor will help draw the picture of TIME and further optimize precision immunotherapy.

Application of nucleoside or nucleotide analogues in RNA dynamics and RNA-binding protein analysis

Cellular RNAs undergo dynamic changes during RNA biological processes, which are tightly orchestrated by RNA-binding proteins (RBPs). Yet, the investigation of RNA dynamics is hurdled by highly abundant steady-state RNAs, which make the signals of dynamic RNAs less detectable. Notably, the exert of nucleoside or nucleotide analogue-based RNA technologies has provided a remarkable platform for RNA dynamics research, revealing diverse unnoticed features in RNA metabolism. In this review, we focus on the application of two types of analogue-based RNA sequencing, antigen-/antibody- and click chemistry-based methodologies, and summarize the RNA dynamics features revealed. Moreover, we discuss emerging single-cell newly transcribed RNA sequencing methodologies based on nucleoside analogue labeling, which provides novel insights into RNA dynamics regulation at single-cell resolution. On the other hand, we also emphasize the identification of RBPs that interact with polyA, non-polyA RNAs, or newly transcribed RNAs and also their associated RNA-binding domains at genomewide level through ultraviolet crosslinking and mass spectrometry in different contexts. We anticipated that further modification and development of these analogue-based RNA and RBP capture technologies will aid in obtaining an unprecedented understanding of RNA biology. This article is categorized under: RNA Interactions with Proteins and Other Molecules &gt; Protein-RNA Recognition RNA Structure and Dynamics &gt; RNA Structure, Dynamics and Chemistry RNA Methods &gt; RNA Analyses in Cells.

Epigenetic Dysregulations in Arsenic-Induced Carcinogenesis

Arsenic is a crucial environmental metalloid whose high toxicity levels negatively impact human health. It poses significant health concerns to millions of people in developed and developing countries such as the USA, Canada, Bangladesh, India, China, and Mexico by enhancing sensitivity to various types of diseases, including cancers. However, how arsenic causes changes in gene expression that results in heinous conditions remains elusive. One of the proposed essential mechanisms that still has seen limited research with regard to causing disease upon arsenic exposure is the dysregulation of epigenetic components. In this review, we have extensively summarized current discoveries in arsenic-induced epigenetic modifications in carcinogenesis and angiogenesis. Importantly, we highlight the possible mechanisms underlying epigenetic reprogramming through arsenic exposure that cause changes in cell signaling and dysfunctions of different epigenetic elements.

SETD2 regulates gene transcription patterns and is associated with radiosensitivity in lung adenocarcinoma

Lung adenocarcinoma (LUAD) has high morbidity and mortality worldwide, and its prognosis remains unsatisfactory. Identification of epigenetic biomarkers associated with radiosensitivity is beneficial for precision medicine in LUAD patients. SETD2 is important in repairing DNA double-strand breaks and maintaining chromatin integrity. Our studies established a comprehensive analysis pipeline, which identified SETD2 as a radiosensitivity signature. Multi-omics analysis revealed enhanced chromatin accessibility and gene transcription by SETD2. In both LUAD bulk RNA sequencing (RNA-seq) and single-cell RNA sequencing (scRNA-seq), we found that SETD2-associated positive transcription patterns were associated with DNA damage responses. SETD2 knockdown significantly upregulated tumor cell apoptosis, attenuated proliferation and migration of LUAD tumor cells, and enhanced radiosensitivity in vitro. Moreover, SETD2 was a favorably prognostic factor whose effects were antagonized by the m6A-related genes RBM15 and YTHDF3 in LUAD. In brief, SETD2 was a promising epigenetic biomarker in LUAD patients.

N(6)-methyladenosine methylation-regulated polo-like kinase 1 cell cycle homeostasis as a potential target of radiotherapy in pancreatic adenocarcinoma

In pancreatic cancer, methyltransferase-like 3 (METTL3), a N(6)-methyladenosine (m6A) methyltransferase, has a favorable effect on tumors and is a risk factor for patients' prognosis. However, the details of what genes are regulated by METTL3 remain unknown. Several RNAs are methylated, and what genes are favored in pancreatic cancer remains unclear. By epitranscriptomic analysis, we report that polo-like kinase 1 (PLK1) is an important hub gene defining patient prognosis in pancreatic cancer and that RNA methylation is involved in regulating its cell cycle-specific expression. We found that insulin like growth factor 2 mRNA binding protein 2 (IGF2BP2) binds to m6A of PLK1 3' untranslated region and is involved in upregulating PLK1 expression and that demethylation of this site activates the ataxia telangiectasia and Rad3-related protein pathway by replicating stress and increasing mitotic catastrophe, resulting in increased radiosensitivity. This suggests that PLK1 methylation is essential for cell cycle maintenance in pancreatic cancer and is a new therapeutic target.

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.

Characterization of the Immune Cell Infiltration Landscape Uncovers Prognostic and Immunogenic Characteristics in Lung Adenocarcinoma

The immune cell infiltration in TME has been reported to be associated with prognosis and immunotherapy efficiency of lung cancers. However, to date, the immune infiltrative landscape of lung adenocarcinoma (LUAD) has not been elucidated yet. Therefore, this study aimed to identify a new transcriptomic-based TME classification and develop a risk scoring system to predict the clinical outcomes of patients with LUAD. We applied “CIBERSORT” algorithm to analyze the transcriptomic data of LUAD samples and classified LUAD into four discrete subtypes according to the distinct immune cell infiltration patterns. Furthermore, we established a novel predictive tool (TMEscore) to quantify the immune infiltration patterns for each LUAD patient by principal component analysis. The TMEscore displayed as a reliable and independent prognostic biomarker for LUAD, with worse survival in TMEscrore-high patients and better survival in TMEscrore-low patients in both TCGA and other five GEO cohorts. In addition, enriched pathways and genomic alterations were also analyzed and compared in different TMEscore subgroups, and we observed that high TMEscore was significantly correlated with more aggressive molecular changes, while the low TMEscore subgroup enriched in immune active-related pathways. The TMEscore-low subtype showed overexpression of PD-1, CTLA4, and associations of other markers of sensitivity to immunotherapy, including TMB, immunophenoscore (IPS) analysis, and tumor immune dysfunction and exclusion (TIDE) algorithm. Conclusively, TMEscore is a promising and reliable biomarker to distinguish the prognosis, the molecular and immune characteristics, and the benefit from ICIs treatments in LUAD.

N6-methyladenosine methyltransferase KIAA1429 elevates colorectal cancer aerobic glycolysis via HK2-dependent manner

Emerging evidence has emphasized the critical roles played by N⁶-methyladenosine RNA (m⁶A) modification in colorectal carcinoma (CRC) initiation and progression. However, the roles and mechanism of m⁶A and KIAA1429 in CRC progression require further clarification. Here, our research aimed to investigate the functions of KIAA1429 in CRC tumorigenesis. Results indicated that KIAA1429 up-regulation closely correlated to the poor prognosis of CRC patients. Bio-functional assays demonstrated that KIAA1429 promoted the aerobic glycolysis, including glucose uptake, lactate production, ATP generation and extracellular acidification rate (ECAR). Mechanistically, KIAA1429 positively up-regulated HK2 level via increasing its mRNA stability by binding the m⁶A site of HK2 mRNA via m⁶A-independent manner. Collectively, our work indicates that KIAA1429 has the potential to promote CRC carcinogenesis by targeting HK2 via m⁶A-independent manner, providing insight into the critical roles of m⁶A in CRC.

Crosstalk of Histone and RNA Modifications Identified a Stromal-Activated Subtype with Poor Survival and Resistance to Immunotherapy in Gastric Cancer

Emerging evidence has revealed the pivotal role of epigenetic modifications in shaping the tumor microenvironment (TME). However, crosstalk between different modification types and their clinical relevance in cancers remain largely unexplored. In this study, using ChIP/MeRIP-seq data of seven human gastric cell lines, we systematically characterized the crosstalk of four epigenetic modification types including H3K4me1, H3K4me3, H3K27ac, and N6-methyladenosine (m6A) and identified a recurrent subtype with high FTO expression and low HDAC1 expression across three independent gastric cancer (GC) cohorts, which we named the epigenetic-modification-dysregulated (EMD) subtype. Patients of the EMD subtype were featured with poor survival, stromal activation, and immune suppression. Extensive relevance to clinical characteristics was observed in the EMD subtype, including the Lauren classification, MSI status, histological grade, TNM stage, the Asian Cancer Research Group classification, and the immune/fibrotic classification. An EMD score was then constructed using WGCNA and ssGSEA algorithms, to precisely recognize the EMD subtype and indicate prognosis and response to immunotherapy in multiple independent GC cohorts. Correlations of the EMD score with tumor mutation burden, tumor purity, aneuploidy score, tumorigenic pathways, TME characteristics, and FTO/HDAC1 ratio were measured. In vitro experiments were performed to demonstrate the correlation between FTO and the epithelial–mesenchymal transition pathway, which suggested FTO as a targetable vulnerability for GC patients with a high EMD score. Altogether, by comprehensively analyzing the epigenetic modification patterns of 1518 GC patients, we identified a novel stromal-activated subtype with poor survival and resistance to immunotherapy, which might benefit from the combined immune checkpoint inhibition therapy with FTO inhibition.

Improved Methods for Deamination-Based m6A Detection

N⁶-methyladenosine (m⁶A) is a critical regulator of gene expression and cellular function. Much of our knowledge of m⁶A has been enabled by the identification of m⁶A sites transcriptome-wide. However, global m⁶A profiling methods require high amounts of input RNA to accurately identify methylated RNAs, making m⁶A profiling from rare cell types or scarce tissue samples infeasible. To overcome this issue, we previously developed DART-seq, which relies on the expression of a fusion protein consisting of the APOBEC1 cytidine deaminase tethered to the m⁶A-binding YTH domain. APOBEC1-YTH directs C-to-U mutations adjacent to m⁶A sites, therefore enabling single nucleotide-resolution m⁶A mapping. Here, we present an improved version of DART-seq which utilizes a variant of the YTH domain engineered to achieve enhanced m⁶A recognition. In addition, we develop in vitro DART-seq and show that it performs similarly to cellular DART-seq and can map m⁶A in any sample of interest using nanogram amounts of total RNA. Altogether, these improvements to the DART-seq approach will enable better m⁶A detection and will facilitate the mapping of m⁶A in samples not previously amenable to global m⁶A profiling.

Comprehensive Analysis Revealed the Potential Implications of m6A Regulators in Lung Adenocarcinoma

Background: The biological significance of RNA N6-methyladenosine (m6A) decoration in tumorigenicity and progression has been highlighted in recent studies, but whether m6A modification plays a potential role in tumor microenvironment (TME) formation and immune regulation in lung adenocarcinoma (LUAD) remains unclear.

Methods: m6A modification features were evaluated by analyzing the multi-omics features of 17 m6A regulators in over 1900 LUAD samples, and at the same time, the correlation between these modification patterns and TME characteristics was analyzed. An m6A score signature–based principal component analysis (PCA) algorithm was constructed to assess the prognosis and responses of individual patients to immunotherapeutic and targeted therapies.

Results: Three different m6A modification patterns were determined in 1901 LUAD samples, which were found to be related to diverse clinical outcomes via different biological pathways. Based on the m6A score extracted from the m6A-associated signature genes, LUAD patients were separated into high- and low-m6A score groups. It was discovered that patients with high m6A scores had longer survival, lower tumor mutation loads, and low PD-L1/PDCD1/CTLA4/TAG3 expression level. In addition, LUAD patients with high m6A scores displayed lower IC₅₀ to some targeted drugs, including nilotinib, erlotinib, imatinib, and lapatinib.

Conclusion: m6A modification was significantly associated with the TME and clinical outcomes. These findings may help gain more insights into the role of m6A decoration in the molecular mechanism of LUAD, thus facilitating the development of more effective personalized treatment strategies.

N6-Methyladenosine Modification Profile in Bovine Mammary Epithelial Cells Treated with Heat-Inactivated Staphylococcus aureus

The symptoms of mastitis caused by Staphylococcus aureus (S. aureus) in dairy cows are not obvious and difficult to identify, resulting in major economic losses. N6-Methyladenosine (m6A) modification has been reported to be closely associated with the occurrence of many diseases. However, only a few reports have described the role of m6A modification in S. aureus-induced mastitis. In this study, after 24 h of treatment with inactivated S. aureus, MAC-T cells (an immortalized bovine mammary epithelial cell line) showed increased expression levels of the inflammatory factors IL-1β, IL-6, TNF-α, and reactive oxygen species. We found that the mRNA levels of METLL3, METLL14, WTAP, and ALKBH5 were also upregulated. Methylated RNA immunoprecipitation sequencing analysis revealed that 133 genes were m6A hypermethylated, and 711 genes were m6A hypomethylated. Biological functional analysis revealed that the differential m6A methylated genes were mainly related to oxidative stress, lipid metabolism, inflammatory response, and so on. In the present study, we also identified 62 genes with significant changes in m6A modification and mRNA expression levels. These findings elucidated the m6A modification spectrum induced by S. aureus in MAC-T cells and provide the basis for subsequent m6A research on mastitis.

Decreased Expression of Programmed Death Ligand-L1 by Seven in Absentia Homolog 2 in Cholangiocarcinoma Enhances T-Cell-Mediated Antitumor Activity

N6-methyladenosine (m6A) has been reported as an important mechanism of post-transcriptional regulation. Programmed death ligand 1 (PD-L1) is a primary immune inhibitory molecule expressed on tumor cells that promotes immune evasion. In addition, seven in absentia homolog 2 (Siah2), a RING E3 ubiquitin ligase, has been involved in tumorigenesis and cancer progression. However, the role of m6A-METTL14-Siah2-PD-L1 axis in immunotherapy remains to be elucidated. In this study, we showed that METTL14, a component of the m⁶A methyltransferase complex, induced Siah2 expression in cholangiocarcinoma (CCA). METTL14 was shown to enrich m⁶A modifications in the 3'UTR region of the Siah2 mRNA, thereby promoting its degradation in an YTHDF2-dependent manner. Furthermore, co-immunoprecipitation experiments demonstrated that Siah2 interacted with PD-L1 by promoting its K63-linked ubiquitination. We also observed that in vitro and in vivo Siah2 knockdown inhibited T cells expansion and cytotoxicity by sustaining tumor cell PD-L1 expression. The METTL14-Siah2-PD-L1-regulating axis was further confirmed in human CCA specimens. Analysis of specimens from patients receiving anti-PD1 immunotherapy suggested that tumors with low Siah2 levels were more sensitive to anti-PD1 immunotherapy. Taken together, our results evidenced a new regulatory mechanism of Siah2 by METTL14-induced mRNA epigenetic modification and the potential role of Siah2 in cancer immunotherapy.

Crosstalk Between Four Types of RNA Modification Writers Characterizes the Tumor Immune Microenvironment Infiltration Patterns in Skin Cutaneous Melanoma

The “writers” of four types of adenosine (A)-related RNA modifications (N6-methyladenosine, N1-methyladenosine, alternative polyadenylation, as well as A-to-inosine RNA editing) are closely related to the tumorigenesis and progression of many cancer types, including skin cutaneous melanoma (SKCM). However, the potential roles of the crosstalk between these RNA modification “writers” in the tumor microenvironment (TME) remain unclear. The RNA modification patterns were identified using an unsupervised clustering method. Subsequently, based on differentially expressed genes responsible for the aforementioned RNA modification patterns, an RNA modification “writer” scoring model (W_Score) was constructed to quantify the RNA modification-associated subtypes in individual patients. Moreover, a correlation analysis for W_Score and the TME characteristics, clinical features, molecular subtypes, drug sensitivities, immune responses, and prognosis was performed. We identified three RNA modification patterns, corresponding to distinct tumor immune microenvironment characteristics and survival outcomes. Based on the W_Score score, which was extracted from the RNA modification-related signature genes, patients with SKCM were divided into high- and low-W_Score groups. The low-W_Score group was characterized by better survival outcomes and strengthened immunocyte infiltration. Further analysis showed that the low-W_Score group was positively associated with higher tumor mutation burden and PD-L1 expression. Of note, two immunotherapy cohorts demonstrated that patients with low W_Score exhibited long-term clinical benefits and an enhanced immune response. This study is the first to systematically analyze four types of A-related RNA modifications in SKCM, revealing that these “writers” essentially contribute to TME complexity and diversity. We quantitatively evaluated the RNA modification patterns in individual tumors, which could aid in developing personalized immunotherapy strategies for patients.

N6-methyladenosine demethylase FTO promotes growth and metastasis of gastric cancer via m6A modification of caveolin-1 and metabolic regulation of mitochondrial dynamics

Gastric cancer (GC) is the fifth most common tumor and the third most deadly cancer worldwide. N6-methyladenosine (m6A) modification has been reported to play a regulatory role in human cancers. However, the exact role of m6A in GC remains largely unknown, and the dysregulation of m6A on mitochondrial metabolism has never been studied. In the present study, we demonstrated that FTO, a key demethylase for RNA m6A modification, was up-regulated in GC tissues, especially in tissues with liver metastasis. Functionally, FTO acted as a promoter for the proliferation and metastasis in GC. Moreover, FTO enhanced the degradation of caveolin-1 mRNA via its demethylation, which regulated the mitochondrial fission/fusion and metabolism. Collectively, our current findings provided some valuable insights into FTO-mediated m6A demethylation modification and could be used as a new strategy for more careful surveillance and aggressive therapeutic intervention.

Predicting RNA 5-Methylcytosine Sites by Using Essential Sequence Features and Distributions

Methylation is one of the most common and considerable modifications in biological systems mediated by multiple enzymes. Recent studies have shown that methylation has been widely identified in different RNA molecules. RNA methylation modifications have various kinds, such as 5-methylcytosine (m5C). However, for individual methylation sites, their functions still remain to be elucidated. Testing of all methylation sites relies heavily on high-throughput sequencing technology, which is expensive and labor consuming. Thus, computational prediction approaches could serve as a substitute. In this study, multiple machine learning models were used to predict possible RNA m5C sites on the basis of mRNA sequences in human and mouse. Each site was represented by several features derived from $k$ -mers of an RNA subsequence containing such site as center. The powerful max-relevance and min-redundancy (mRMR) feature selection method was employed to analyse these features. The outcome feature list was fed into incremental feature selection method, incorporating four classification algorithms, to build efficient models. Furthermore, the sites related to features used in the models were also investigated.

The Prognostic Role of m6A-Related Genes in Paediatric Neuroblastoma Patients

Methods

The gene expression data were extracted from the “Therapeutically Applicable Research to Generate Effective Treatments” (TARGET) database. The differentially expressed genes (DEGs) were identified, and the relationships between DEGs and m6A genes were explored. Then, the correlations among the m6A genes in neuroblastoma were investigated. Finally, the prognostic role of the m6A genes was studied, and the risk model was constructed.

Results

81 NB patients were extracted from the TARGET database. After comparing the gene expression between unfavorable and favorable outcome groups, 73 DEGs were identified, including 54 downregulated genes and 19 upregulated genes. In NB patients, we found that IGF2BP3, METTL14, and METTL16 are prognostic factors for disease-free survival (DFS) while IGF2BP3, METTL14, and METTL16 are prognostic factors for overall survival (OS). Besides, after the risk model construction, the OS between the two risk groups was drawn (log-rank p = 1.64e − 08, HR = 3.438, 95% CI 2.24-5.278). The 1-, 3-, and 5-year time-dependent receiving operating characteristic (ROC) curves were also illustrated, and the areas under the receiver operating characteristic curves (AUCs) attained 0.75, 0.798, and 0.768, respectively.

Conclusions

IGF2BP3, METTL14, and METTL16 were identified as the significant factors for DFS and OS in NB patients.

Identification of a N6-Methyladenosine (m6A)-Related lncRNA Signature for Predicting the Prognosis and Immune Landscape of Lung Squamous Cell Carcinoma

Background

m6A-related lncRNAs emerged as potential targets for tumor diagnosis and treatment. This study aimed to identify m6A-regulated lncRNAs in lung squamous cell carcinoma (LUSC) patients.

Materials and Methods

RNA sequencing and the clinical data of LUSC patients were downloaded from The Cancer Genome Atlas (TCGA) database. The m6A-related lncRNAs were identified by using Pearson correlation assay. Univariate and multivariate Cox regression analyses were utilized to construct a risk model. The performance of the risk model was validated using Kaplan–Meier survival analysis and receiver operating characteristics (ROC). Immune estimation of LUSC was downloaded from TIMER, and the correlations between the risk score and various immune cells infiltration were analyzed using various methods. Differences in immune functions and expression of immune checkpoint inhibitors and m6A regulators between high-risk and low-risk groups were further explored. Finally, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were utilized to explore the biological functions of AL122125.1.

Results

A total of 351 m6A-related lncRNAs were obtained from TCGA. Seven lncRNAs demonstrated prognostic values. A further multivariate Cox regression assay constructed a risk model consisting of two lncRNAs (AL122125.1 and HORMAD2-AS1). The Kaplan–Meier analysis and area under the curve indicated that this risk model could be used to predict the prognosis of LUSC patients. The m6A-related lncRNAs were immune-associated. There were significant correlations between risk score and immune cell infiltration, immune functions, and expression of immune checkpoint inhibitors. Meanwhile, there were significant differences in the expression of m6A regulators between the high- and low-risk groups. Moreover, GO and KEGG analyses revealed that the upregulated expression of AL122125.1 was tumor-related.

Conclusion

In this study, we constructed an m6A-related lncRNA risk model to predict the survival of LUSC patients. This study could provide a novel insight to the prognosis and treatment of LUSC patients.

Crystal structure of human METTL6, the m3C methyltransferase

In tRNA, the epigenetic m3C modification at position 32 in the anticodon loop is highly conserved in eukaryotes, which maintains the folding and basepairing functions of the anticodon. However, the responsible enzymes METTL2 and METTL6 were identified only in recent years. The loss of human METTL6 (hMETTL6) affects the translational process and proteostasis in cells, while in mESCs cells, it leads to defective pluripotency potential. Despite its important functions, the catalytic mechanism of the C32 methylation by this enzyme is poorly understood. Here we present the 1.9 Å high-resolution crystal structure of hMETTL6 bound by SAH. The key residues interacting with the ligand were identified and their roles were confirmed by ITC. We generated a docking model for the hMETTL6-SAH-CMP ternary complex. Interestingly, the CMP molecule binds into a cavity in a positive patch with the base ring pointing to the inside, suggesting a flipped-base mechanism for methylation. We further generated a model for the quaternary complex with tRNASer as a component, which reasonably explained the biochemical behaviors of hMETTL6. Taken together, our crystallographic and biochemical studies provide important insight into the molecular recognition mechanism by METTL6 and may aid in the METTL-based rational drug design in the future.

Modifications of the human tRNA anticodon loop and their associations with genetic diseases

Transfer RNAs (tRNAs) harbor the most diverse posttranscriptional modifications. Among such modifications, those in the anticodon loop, either on nucleosides or base groups, compose over half of the identified posttranscriptional modifications. The derivatives of modified nucleotides and the crosstalk of different chemical modifications further add to the structural and functional complexity of tRNAs. These modifications play critical roles in maintaining anticodon loop conformation, wobble base pairing, efficient aminoacylation, and translation speed and fidelity as well as mediating various responses to different stress conditions. Posttranscriptional modifications of tRNA are catalyzed mainly by enzymes and/or cofactors encoded by nuclear genes, whose mutations are firmly connected with diverse human diseases involving genetic nervous system disorders and/or the onset of multisystem failure. In this review, we summarize recent studies about the mechanisms of tRNA modifications occurring at tRNA anticodon loops. In addition, the pathogenesis of related disease-causing mutations at these genes is briefly described.

Applications and potentials of nanopore sequencing in the (epi)genome and (epi)transcriptome era

The Human Genome Project opened an era of (epi)genomic research, and also provided a platform for the development of new sequencing technologies. During and after the project, several sequencing technologies continue to dominate nucleic acid sequencing markets. Currently, Illumina (short-read), PacBio (long-read), and Oxford Nanopore (long-read) are the most popular sequencing technologies. Unlike PacBio or the popular short-read sequencers before it, which, as examples of the second or so-called Next-Generation Sequencing platforms, need to synthesize when sequencing, nanopore technology directly sequences native DNA and RNA molecules. Nanopore sequencing, therefore, avoids converting mRNA into cDNA molecules, which not only allows for the sequencing of extremely long native DNA and full-length RNA molecules but also document modifications that have been made to those native DNA or RNA bases. In this review on direct DNA sequencing and direct RNA sequencing using Oxford Nanopore technology, we focus on their development and application achievements, discussing their challenges and future perspective. We also address the problems researchers may encounter applying these approaches in their research topics, and how to resolve them.

Artificial intelligence: A powerful paradigm for scientific research

Artificial intelligence (AI) coupled with promising machine learning (ML) techniques well known from computer science is broadly affecting many aspects of various fields including science and technology, industry, and even our day-to-day life. The ML techniques have been developed to analyze high-throughput data with a view to obtaining useful insights, categorizing, predicting, and making evidence-based decisions in novel ways, which will promote the growth of novel applications and fuel the sustainable booming of AI. This paper undertakes a comprehensive survey on the development and application of AI in different aspects of fundamental sciences, including information science, mathematics, medical science, materials science, geoscience, life science, physics, and chemistry. The challenges that each discipline of science meets, and the potentials of AI techniques to handle these challenges, are discussed in detail. Moreover, we shed light on new research trends entailing the integration of AI into each scientific discipline. The aim of this paper is to provide a broad research guideline on fundamental sciences with potential infusion of AI, to help motivate researchers to deeply understand the state-of-the-art applications of AI-based fundamental sciences, and thereby to help promote the continuous development of these fundamental sciences.

CCR5 is a prognostic biomarker and an immune regulator for triple negative breast cancer

This study aims to explore the clinical implications and potential mechanistic functions of CCR5 in triple negative breast cancer. Briefly, we demonstrated that CCR5 is overexpressed in TNBC and is associated with better prognosis of TNBC. CCR5 expression is positively correlated with tumor immune cell infiltration and tumor immune response related pathways. Multi-omics data analyses identified CCR5 associated genomic and proteomic changes. CCR5 overexpression was associated with better overall survival in TNBC patients with TP53 mutation. We also summarized the latest findings on ICB efficacy related genes and explored the association between CCR5 and those genes. These results indicated that CCR5 is a potential tumor suppressor gene and individualized therapeutic strategy could be established based on multi-omics background and expression pattern of ICB related genes. In conclusion, CCR5 is associated with better survival of TNBC patients with TP53 mutation, which may exert its roles through tumor immune environment.

Comprehensive Analysis of the Expression of TGF-β Signaling Regulators and Prognosis in Human Esophageal Cancer

More and more evidences show that TGF-β has a crucial role in tumor initiation and development. However, the mechanism of the TGF-β signal regulator in esophageal cancer (EC) is still unclear. Here, we use a variety of bioinformatics methods to analyze the expression and survival data of TGF-β signal regulators in patients with EC. We extracted the expression of the S-TGF-β signal regulator from The Cancer Genome Atlas (TCGA). The cBioPortal database was used to assess the frequency of genetic variation. The TGF-β signal regulator is expressed in EC and normal tissues. The objective is to use the Kaplan-Meier plotter database to investigate the prognostic value of TGF-β signal regulators in cancer patients. The DAVID and clusterProfiler software package were used for functional enrichment analysis. We found that patients with TGF-β signaling mutations have shorter overall survival, disease-free survival, disease-specific survival, platinum overall survival, and platinum-free progression survival. We found that compared with the noncancerous tissues of patients with EC, ZFYVE9, BMPR1B, TGFB3, TGFBRAP1, ACVRL1, TGFBR2, SMAD4, SMAD7, ACVR2A, BMPR1, and SMAD9 were significantly downregulated in tumor tissues, while ACVR1 and Smad1 were significantly upregulated in tumor samples. Univariate survival analysis showed that ACVR1, TGFBR3, TGFBRAP1, BMPR1A, SMAD4, and TGFBR2 were positively correlated with overall survival (OS) prolongation. In addition, TGF-β signal transduction regulators could be divided into two classes. Subclass 1 was involved in regulating cell adhesion, PI3K-Akt signaling, and Rap1 signaling. Subclass 2 was related to regulating angiogenesis and PI3K signaling. In short, all members of TGF-β signal regulators can be used as biomarkers to predict the prognosis of patients with EC.

Upregulation of miR-29c-3p Hinders Melanoma Progression by Inhibiting CDCA4 Expression

OBJECTIVE

To investigate the expression and regulation mechanism of miR-29c-3p and cell division cycle associated 4 (CDCA4) in melanoma (MM). Data and Methods. Fifty-nine patients with MM admitted to our hospital were enrolled as the MM group. They were followed up for 3 years to analyze the prognostic factors; meanwhile, 51 healthy subjects were allocated into a normal group. MM cell lines (M21 and C8161) were transfected with miR-29c-3p-mimics, miR-29c-3p-inhibitor, miR-NC, si-CDCA4, and sh-CDCA4. The expression of miR-29c-3p, CDCA4, Bax, Caspase3, Bcl-2, N-cadherin, vimentin, and E-cadherin was quantified, and cell proliferation, migration, invasion, and apoptosis, as well as epithelial-mesenchymal transition (EMT), were determined.

RESULTS

Serum miR-29c-3p was lowly expressed and CDCA4 was highly expressed in the MM group. The area under the curve (AUC) of both for diagnosing MM was greater than 0.9. miR-29c-3p and CDCA4 were related to regional lymph node staging (N staging), distant metastasis (M staging), tumor diameter, and pathological differentiation. Low miR-29c-3p and high CDCA4 were associated with poor prognosis of MM. Overexpression of miR-29c-3p and suppression of CDCA4 hindered cell proliferation, migration, invasion, and expression of Bax, Caspase3, N-cadherin, and vimentin, but cell apoptosis and expression of Bcl-2 and E-cadherin were enhanced. Dual-luciferase reporter (DLR) assay confirmed the targeted relationship between miR-29c-3p and CDCA4. After miR-29c-3p-mimics+sh-CDCA4 was transfected into M21 and C8161 cells, the proliferation, invasion, and apoptosis were not different from those in the miR-NC group transfected with unrelated sequences.

CONCLUSION

Overexpression of miR-29c-3p suppresses CDCA4 expression and decreases proliferation, migration, invasion, apoptosis, and EMT of MM cells, thus hindering MM progression.

CircITGA7 Suppresses Gastric Cancer Progression Through miR-1471/MTDH Axis

In recent years, there have been reports about the involvement of circular RNAs (circRNAs) in the pathogenesis of gastric cancer (GC), but the molecular mechanism in cell proliferation, invasion, and migration is still unclear. Based on The Cancer Genome Atlas (TCGA) database, we analyzed differentially expressed circRNAs between GC and non-tumor tissues. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analysis were used to clarify the functional role in GC. Here, we showed that circITGA7 was lowly expressed in GC tissues based on the TCGA database. In vitro, silencing the expression of circITGA7 increased cell proliferation and metastasis, whereas overexpression did the opposite. Mechanistically, miR-1471 has circITGA7 as a sponge, and miR-1471 has metadherin (MTDH) as a target gene. Consequently, functional analysis showed that the tumor suppressor effect of circITGA7 was the result of regulating the miR-1471/MTDH axis. Overall, the circITGA7/miR-1471/MTDH signaling pathway may play a crucial role in GC, providing a new potential mechanism involved in GC progression.

Mutually exclusive substrate selection strategy by human m3C RNA transferases METTL2A and METTL6

tRNAs harbor the most diverse posttranscriptional modifications. The 3-methylcytidine (m3C) is widely distributed at position C32 (m3C32) of eukaryotic tRNAThr and tRNASer species. m3C32 is decorated by the single methyltransferase Trm140 in budding yeasts; however, two (Trm140 and Trm141 in fission yeasts) or three enzymes (METTL2A, METTL2B and METTL6 in mammals) are involved in its biogenesis. The rationale for the existence of multiple m3C32 methyltransferases and their substrate discrimination mechanism is hitherto unknown. Here, we revealed that both METTL2A and METTL2B are expressed in vivo. We purified human METTL2A, METTL2B, and METTL6 to high homogeneity. We successfully reconstituted m3C32 modification activity for tRNAThr by METT2A and for tRNASer(GCU) by METTL6, assisted by seryl-tRNA synthetase (SerRS) in vitro. Compared with METTL2A, METTL2B exhibited dramatically lower activity in vitro. Both G35 and t6A at position 37 (t6A37) are necessary but insufficient prerequisites for tRNAThr m3C32 formation, while the anticodon loop and the long variable arm, but not t6A37, are key determinants for tRNASer(GCU) m3C32 biogenesis, likely being recognized synergistically by METTL6 and SerRS, respectively. Finally, we proposed a mutually exclusive substrate selection model to ensure correct discrimination among multiple tRNAs by multiple m3C32 methyltransferases.