Identification of m6A-related genes and m6A RNA methylation regulators in pancreatic cancer and their association with survival

m6A RNA Methylation and Implications for Hepatic Lipid Metabolism

This review presents a summary of recent progress in research on the N6-methyladenosine (m6A) modification and regulatory roles in hepatic lipid metabolism. As the most abundant internal modification of eukaryotic RNA, the m6A modification is a dynamic and reversible process of the m6A enzyme system, which includes writers, erasers, and readers. m6A methylation depressed lipid synthesis and facilitated lipolysis in liver. The depletion of m6A methyltransferase Mettl14/Mettl3 raised fatty acid synthase (FAS), stearoyl-CoA desaturase-1 (SCD1), acetyl-CoA carboxylase (ACC), and elongase of very long chain fatty acids 6 (ELOVL6) in rodent liver, causing increases in liver weight, triglyceride (TG) production, and content in hepatocytes. FTO catalyzed m6A demethylation and the suppression m6A reader YTHDC2 promoted hepatocellular TG generation and hepatic steatosis in C57BL/6 mice through sterol regulatory element-binding protein 1c (SREBP-1c) signaling pathway, which upregulated the lipogenic genes FAS, SCD1, ACC, recombinant acetyl coenzyme a carboxylase alpha, and cell death-inducing DNA fragmentation factor-like effector C (CIDEC). Furthermore, FTO overexpression did not only enhance mitochondrial fusion to impair mitochondrial function and lipid oxidation but also promoted lipid peroxidation, accompanied by excessive TG in hepatocytes and rodent liver. Elevated m6A modification potently suppressed hepatic lipid accumulation, while the shrinkage of m6A modification arose hepatic lipid deposition. These findings have highlighted the beneficial role of m6A RNA methylation in hepatic lipid metabolism, potentially protecting liver from lipid metabolic disorders.

Comprehensive analysis of m6A reader YTHDF2 prognosis, immune infiltration, and related regulatory networks in hepatocellular carcinoma

Background

N6-Methyladenosine (m6A) RNA modification is the most prevalent internal modification pattern in eukaryotic mRNAs and plays critical roles in diverse physiological and pathological processes. However, the expression of m6A regulator YTHDF2, its prognostic value, its biological function, its correlation with tumor microenvironment (TME) immune infiltrates, and related regulatory networks in hepatocellular carcinoma (HCC) remain determined.

Methods

TCGA, GTEx, and GEO databases were used to investigate the expression profile of YTHDF2 in HCC. We performed differentially expressed genes (DEGs) analysis and constructed a PPI network to explore the biological processes of YTHDF2 in HCC. Kaplan-Meier curves and Cox regression analysis were used to assess the prognostic value of YTHDF2 and then a clinical prognostic nomogram was constructed. Additionally, ssGSEA was performed to assess the correlation between YTHDF2 and immune infiltration levels. The TISIDB database was applied to explore the expression of YTHDF2 in immune and molecular subtypes of HCC. GSEA identifies the YTHDF2-related signaling pathways. Finally, we utilized miRNet and starBase database to construct regulatory networks for HCC based on lncRNA-miRNA and miRNA-YTHDF2 interactions.

Results

YTHDF2 was significantly upregulated in HCC tumor tissues compared with the adjacent normal tissues. HCC patients in the high YTHDF2 expression group had poorer survival. Multivariate Cox analysis suggested that YTHDF2 may be a new independent prognostic indicator for HCC patients, with the prognostic nomogram exhibiting satisfactory results. YTHDF2 expression was significantly correlated with TME immune cell-infiltrating characteristics. Strong correlations were also shown in immune subtypes, molecular subtypes and immune checkpoints. Further analysis revealed that the combination of YTHDF2 expression and immune cell score was considerably associated with survival outcome in HCC patients. GESA analysis demonstrated that high YTHDF2 expression is associated with multiple biological processes and oncogenic pathways. Moreover, 14 possible regulatory networks were constructed, which are associated with HCC progression.

Conclusion

Our findings revealed that YTHDF2 may serve as a promising prognostic biomarker for HCC and may regulate the tumor immune microenvironment to provide effective therapeutic strategies.

Exploring potential genes and mechanisms linking erectile dysfunction and depression

Background

The clinical correlation between erectile dysfunction (ED) and depression has been revealed in cumulative studies. However, the evidence of shared mechanisms between them was insufficient. This study aimed to explore common transcriptomic alterations associated with ED and depression.

Materials and methods

The gene sets associated with ED and depression were collected from the Gene Expression Omnibus (GEO) database. Comparative analysis was conducted to obtain common genes. Using R software and other appropriate tools, we conducted a range of analyses, including function enrichment, interactive network creation, gene cluster analysis, and transcriptional and post-transcriptional signature profiling. Candidate hub crosslinks between ED and depression were selected after external validation and molecular experiments. Furthermore, subpopulation location and disease association of hub genes were explored.

Results

A total of 85 common genes were identified between ED and depression. These genes strongly correlate with cell adhesion, redox homeostasis, reactive oxygen species metabolic process, and neuronal cell body. An interactive network consisting of 80 proteins and 216 interactions was thereby developed. Analysis of the proteomic signature of common genes highlighted eight major shared genes: CLDN5, COL7A1, LDHA, MAP2K2, RETSAT, SEMA3A, TAGLN, and TBC1D1. These genes were involved in blood vessel morphogenesis and muscle cell activity. A subsequent transcription factor (TF)-miRNA network showed 47 TFs and 88 miRNAs relevant to shared genes. Finally, CLDN5 and TBC1D1 were well-validated and identified as the hub crosslinks between ED and depression. These genes had specific subpopulation locations in the corpus cavernosum and brain tissue, respectively.

Conclusion

Our study is the first to investigate common transcriptomic alterations and the shared biological roles of ED and depression. The findings of this study provide insights into the referential molecular mechanisms underlying the co-existence between depression and ED.

A glycosylation-related signature predicts survival in pancreatic cancer

BACKGROUND

Tumor initiation and progression are closely associated with glycosylation. However, glycosylated molecules have not been the subject of extensive studies as prognostic markers for pancreatic cancer. The objectives of this study were to identify glycosylation-related genes in pancreatic cancer and use them to construct reliable prognostic models.

MATERIALS AND METHODS

The Cancer Genome Atlas and Gene Expression Omnibus databases were used to assess the differential expression of glycosylation-related genes; four clusters were identified based on consistent clustering analysis. Kaplan-Meier analyses identified three glycosylation-related genes associated with overall survival. LASSO analysis was then performed on The Cancer Genome Atlas and International Cancer Genome Consortium databases to identify glycosylation-related signatures. We identified 12 GRGs differently expressed in pancreatic cancer and selected three genes (SEL1L, TUBA1C, and SDC1) to build a prognostic model. Thereafter, patients were divided into high and low-risk groups. Eventually, we performed Quantitative real-time PCR (qRT-PCR) to validate the signature.

RESULTS

Clinical outcomes were significantly poorer in the high-risk group than in the low-risk group. There were also significant correlations between the high-risk group and several risk factors, including no-smoking history, drinking history, radiotherapy history, and lower tumor grade. Furthermore, the high-risk group had a higher proportion of immune cells. Eventually, three glycosylation-related genes were validated in human PC cell lines.

CONCLUSION

This study identified the glycosylation-related signature for pancreatic cancer. It is an effective predictor of survival and can guide treatment decisions.

Hypoxia promotes immune escape of pancreatic cancer cells by lncRNA NNT-AS1/METTL3-HuR-mediated ITGB1 m6A modification

Pancreatic cancer (PC) cell immune escape is a crucial element in PC malignant development. Some previous studies have reported that LncRNA NNT-AS1 played a carcinogenic role in various tumors. However, the effect of lncRNA NNT-AS1 in PC cell immune escape remains unclear. To evaluate PC cell immune escape, PC cells were co-cultured with CD8+ T cells under a hypoxic condition. PC cell proliferation and migration were evaluated using the colony formation assay and transwell assay. CD8+ T cell proliferation and aoptosis were measured using the carboxy fluorescein diacetate succinimidyl ester (CFSE) assay and flow cytometry. The secretion of antitumor cytokines was assessed using enzyme-linked immunosorbent assay (ELISA). The molecular interactions were analyzed using chromatin immunoprecipitation (ChIP), RNA immunoprecipitation (RIP), or dual-luciferase reporter gene assays. A tumor xenograft model was established to evaluate the effects of lncRNA NNT-AS1 on PC in vivo. It was found that lncRNA NNT-AS1 was highly expressed in PC, and its silencing inhibited hypoxia-induced PC cell growth and immune escape in vivo and in vitro. Mechanically, HIF-1α transcriptionally activated NNT-AS1 expression and NNT-AS1 increased ITGB1 stability and expression in a METTL3-HuR dependent manner. ITGB1 overexpression reversed the inhibitory effects of NNT-AS1 knockdown on hypoxia-induced PC cell immune escape. In conclusion, Hypoxia promoted PC cell immune escape through lncRNA NNT-AS1/METTL3-HuR-mediated m6A modification to increase ITGB1 expression, which provided a theoretical foundation and a potential therapeutic target for PC.

Research progress on branched-chain amino acid aminotransferases

Branched-chain amino acid aminotransferases, widely present in natural organisms, catalyze bidirectional amino transfer between branched-chain amino acids and branched-chain α-ketoacids in cells. Branched-chain amino acid aminotransferases play an important role in the metabolism of branched-chain amino acids. In this paper, the interspecific evolution and biological characteristics of branched-chain amino acid aminotransferases are introduced, the related research of branched-chain amino acid aminotransferases in animals, plants, microorganisms and humans is summarized and the molecular mechanism of branched-chain amino acid aminotransferase is analyzed. It has been found that branched-chain amino acid metabolism disorders are closely related to various diseases in humans and animals and plants, such as diabetes, cardiovascular diseases, brain diseases, neurological diseases and cancer. In particular, branched-chain amino acid aminotransferases play an important role in the development of various tumors. Branched-chain amino acid aminotransferases have been used as potential targets for various cancers. This article reviews the research on branched-chain amino acid aminotransferases, aiming to provide a reference for clinical research on targeted therapy for various diseases and different cancers.

METTL14 inhibits malignant progression of oral squamous cell carcinoma by targeting the autophagy-related gene RB1CC1 in an m6A-IGF2BP2-dependent manner

N6-methyladenosine (m6A) plays crucial roles in tumorigenesis and autophagy. However, the underlying mechanisms mediated by m6A and autophagy in the malignant progression of oral squamous cell carcinoma (OSCC) remain unclear. In the present study, we revealed that down-regulated expression of METTL14 was correlated with advanced clinicopathological characteristics and poor prognosis in OSCC. METTL14 knockdown significantly inhibited autophagy and facilitated malignant progression in vitro, and promoted tumor growth and metastasis in vivo. A cell model of rapamycin-induced autophagy was established to identify RB1CC1 as a potential target gene involved in m6A-regulated autophagy in OSCC, through RNA sequencing and methylated RNA immunoprecipitation sequencing (meRIP-seq) analysis. Mechanistically, we confirmed that METTL14 posttranscriptionally enhanced RB1CC1 expression in an m6A-IGF2BP2-dependent manner, thereby affecting autophagy and progression in OSCC, through methylated RNA immunoprecipitation qRT-PCR (meRIP-qPCR), RNA stability assays, mutagenesis assays and dual-luciferase reporter. Collectively, our findings demonstrated that METTL14 serves as an OSCC suppressor by regulating the autophagy-related gene RB1CC1 through m6A modification, which may provide a new insight for the diagnosis and therapy of OSCC.

m6A eraser FTO impairs gemcitabine resistance in pancreatic cancer through influencing NEDD4 mRNA stability by regulating the PTEN/PI3K/AKT pathway

BACKGROUND

Gemcitabine resistance has brought great challenges to the treatment of pancreatic cancer. The N6-methyladenosine (m6A) mutation has been shown to have a significant regulatory role in chemosensitivity; however, it is not apparent whether gemcitabine resistance can be regulated by fat mass and obesity-associated protein (FTO).

METHODS

Cells with established gemcitabine resistance and tissues from pancreatic cancer patients were used to evaluate FTO expression. The biological mechanisms of the effects of FTO on gemcitabine resistant cells were investigated using CCK-8, colony formation assay, flow cytometry, and inhibitory concentration 50. Immunoprecipitation/mass spectrometry, MeRIP-seq, RNA sequencing and RIP assays, RNA stability, luciferase reporter, and RNA pull down assays were employed to examine the mechanism of FTO affecting gemcitabine resistant pancreatic cancer cells.

RESULTS

The results revealed that FTO was substantially expressed in cells and tissues that were resistant to gemcitabine. Functionally, the gemcitabine resistance of pancreatic cancer could be enhanced by FTO, while its depletion inhibited the growth of gemcitabine resistant tumor cells in vivo. Immunoprecipitation/mass spectrometry showed that the FTO protein can be bound to USP7 and deubiquitinated by USP7, leading to the upregulation of FTO. At the same time, FTO knockdown significantly decreased the expression level of NEDD4 in an m6A-dependent manner. RNA pull down and RNA immunoprecipitation verified YTHDF2 as the reader of NEDD4, which promoted the chemoresistance of gemcitabine resistant cells. FTO knockdown markedly increased the PTEN expression level in an NEDD4-dependent manner and influenced the chemosensitivity to gemcitabine through the PI3K/AKT pathway in pancreatic cancer cells.

CONCLUSION

In conclusion, we found that gemcitabine resistance in pancreatic cancer can be influenced by FTO that demethylates NEDD4 RNA in a m6A-dependent manner, which then influences the PTEN expression level and thereby affects the PI3K/AKT pathway. We also identified that the FTO level can be upregulated by USP7.

Comprehensive analysis of m6A regulators and relationship with tumor microenvironment, immunotherapy strategies in colorectal adenocarcinoma

BACKGROUND

The N6-methyladenosine (m6A) RNA modification is the most prevalent and abundant type found in eukaryotic cells. It plays a crucial role in the initiation and progression of cancers. In this study, we aimed to comprehensively investigate the landscape of m6A regulators and their association with tumor microenvironment (TME), immunotherapeutic strategies in colon adenocarcinoma (COAD).

RESULTS

The differential expression, mutation, CNV frequency and prognostic value of 27 m6A regulators were systematically analyzed in COAD. Patients were classified into two clusters based on m6A regulators through consistent clustering analysis, with cluster A showing significant survival benefits. Most of the m6A regulators were negatively correlated with immune cells, except for WTAP, IGF2BP3, FTO, ALKBH5, which showed a positive correlation. We developed an m6A scoring system to calculate the m6Ascore for each patient. Patients with a high-m6Ascore had a better outcome, with the AUC of 0.775. An independent cohort of 416 COAD patients acquired from GSE38832 database was used to validate the prognosis prediction ability of m6Ascore. Moreover, the m6Ascore was negatively correlated with infiltration of anti-tumor immune cells. Additionally, patients with a high-m6Ascore responded better to anti-PD1 and anti-CTLA4 therapies, and those with MSI-H had a higher m6Ascore. Finally, we investigated the value of m6Ascore in predicting the response of patients to 15 commonly used drugs.

CONCLUSIONS

We comprehensively analyzed m6A regulators in COAD, including RNA expression, CNV changes, mutations and their correlation with TME. Our results showed that the m6A scoring system had significant predictive power for the prognosis of COAD patients, potentially leading to new personalized immunotherapy strategies.

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

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

RNA demethylase ALKBH5 suppresses tumorigenesis via inhibiting proliferation and invasion and promoting CD8+ T cell infiltration in colorectal cancer

BACKGROUND

ALKBH5 belongs to the ALKB family consists of a Fe (II) and a-ketoglutarate-dependent dioxygenase. ALKBH5 directly catalyzes the oxidative demethylation of m6A-methylated adenosine. ALKBH5 involves in tumorigenesis and tumor progression, and is often dysregulated in a wide range of cancers, including colorectal cancer. Emerging evidence indicates that the expression of ALKBH5 is associated with the abundance of infiltrating immune cells in the microenvironment. However, how ALKBH5 affects immune cell infiltration in the microenvironment in colorectal cancer (CRC) has not been reported. The aim of this study was to identify how the expression of ALKBH5 affects the biological behaviors of CRC cell lines and regulates the effects on infiltrating CD8⁺ T cells in CRC microenvironment with its specific mechanism.

METHODS

Firstly, the transcriptional expression profiles of CRC were downloaded from TCGA database and integrated via R software (4.1.2). Between CRC and normal colorectal tissues, ALKBH5 mRNA expressions were compared (Wilcoxon rank-sum). We further identified the expression levels of ALKBH5 in CRC tissues and cell lines through quantitative PCR, western blot, and immunohistochemistry. Then, how ALKBH5 affects the biological behaviors of CRC cells were confirmed by gain- and loss-of-function analysis. Furthermore, the relationship between ALKBH5 level and 22 tumor-infiltrating immune cells was examined through CIBERSORT in R software. Furthermore, we explored the correlation between ALKBH5 expression and tumor-infiltrated CD8⁺, CD4⁺ and regulatory T cells by utilizing the TIMER database. Finally, the association between chemokines and CD8⁺ T cells infiltration in CRC was analyzed using GEPIA online database. qRT-PCR, WB and IHC were used to further determine the effect of ALKBH5 on NF-κB-CCL5 signaling axis and CD8⁺ T cells infiltration.

RESULTS

Clinically, ALKBH5 expression was downregulated in CRC and low levels of ALKBH5 expression were correlated with poor overall survival (OS). Functionally, overexpression of ALKBH5 reduced the proliferation, migration and invasion of CRC cells, and vice versa. Overexpression of ALKBH5 suppresses NF-κB pathway, thus reduces CCL5 expression and promotes CD8⁺ T cells infiltration in CRC microenvironment.

CONCLUSIONS

ALKBH5 is poorly expressed in CRC, and overexpression of ALKBH5 attenuates CRC malignant progression by inhibiting CRC cell proliferation, migration, invasion and promoting CD8⁺ T cells infiltration in the tumor microenvironment through NF-κB-CCL5 axis.

The role of demethylase AlkB homologs in cancer

The AlkB family (ALKBH1-8 and FTO), a member of the Fe (II)- and α-ketoglutarate-dependent dioxygenase superfamily, has shown the ability to catalyze the demethylation of a variety of substrates, including DNA, RNA, and histones. Methylation is one of the natural organisms’ most prevalent forms of epigenetic modifications. Methylation and demethylation processes on genetic material regulate gene transcription and expression. A wide variety of enzymes are involved in these processes. The methylation levels of DNA, RNA, and histones are highly conserved. Stable methylation levels at different stages can coordinate the regulation of gene expression, DNA repair, and DNA replication. Dynamic methylation changes are essential for the abilities of cell growth, differentiation, and division. In some malignancies, the methylation of DNA, RNA, and histones is frequently altered. To date, nine AlkB homologs as demethylases have been identified in numerous cancers’ biological processes. In this review, we summarize the latest advances in the research of the structures, enzymatic activities, and substrates of the AlkB homologs and the role of these nine homologs as demethylases in cancer genesis, progression, metastasis, and invasion. We provide some new directions for the AlkB homologs in cancer research. In addition, the AlkB family is expected to be a new target for tumor diagnosis and treatment.

The function and clinical implication of YTHDF1 in the human system development and cancer

YTHDF1 is a well-characterized m6A reader protein that is essential for protein translation, stem cell self-renewal, and embryonic development. YTHDF1 regulates target gene expression by diverse molecular mechanisms, such as promoting protein translation or modulating the stability of mRNA. The cellular levels of YTHDF1 are precisely regulated by a complicated transcriptional, post-transcriptional, and post-translational network. Very solid evidence supports the pivotal role of YTHDF1 in embryonic development and human cancer progression. In this review, we discuss how YTHDF1 influences both the physiological and pathological biology of the central nervous, reproductive and immune systems. Therefore we focus on some relevant aspects of the regulatory role played by YTHDF1 as gene expression, complex cell networking: stem cell self-renewal, embryonic development, and human cancers progression. We propose that YTHDF1 is a promising future cancer biomarker for detection, progression, and prognosis. Targeting YTHDF1 holds therapeutic potential, as the overexpression of YTHDF1 is associated with tumor resistance to chemotherapy and immunotherapy.

Expression patterns of platinum resistance-related genes in lung adenocarcinoma and related clinical value models

The purpose of this study was to explore platinum resistance-related biomarkers and mechanisms in lung adenocarcinoma. Through the analysis of gene expression data of lung adenocarcinoma patients and normal patients from The Cancer Genome Atlas, Gene Expression Omnibus database, and A database of genes related to platinum resistance, platinum resistance genes in lung adenocarcinoma and platinum resistance-related differentially expressed genes were obtained. After screening by a statistical significance threshold, a total of 252 genes were defined as platinum resistance genes with significant differential expression, of which 161 were up-regulated and 91 were down-regulated. The enrichment results of up-regulated gene Gene Ontology (GO) showed that TOP3 entries related to biological processes (BP) were double-strand break repair, DNA recombination, DNA replication, the down-regulated gene GO enriches the TOP3 items about biological processes (BP) as a response to lipopolysaccharide, muscle cell proliferation, response to molecule of bacterial origin. Gene Set Enrichment Analysis showed that the top three were e2f targets, g2m checkpoint, and rgf beta signaling. A prognostic model based on non-negative matrix factorization classification showed the characteristics of high- and low-risk groups. The prognostic model established by least absolute shrinkage and selection operator regression and risk factor analysis showed that genes such as HOXB7, NT5E, and KRT18 were positively correlated with risk score. By analyzing the differences in m6A regulatory factors between high- and low-risk groups, it was found that FTO, GPM6A, METTL3, and YTHDC2 were higher in the low-risk group, while HNRNPA2B1, HNRNPC, TGF2BP1, IGF2BP2, IGF2BP3, and RBM15B were higher in the high-risk group. Immune infiltration and drug sensitivity analysis also showed the gene characteristics of the platinum-resistant population in lung adenocarcinoma. ceRNA analysis showed that has-miR-374a-5p and RP6-24A23.7 were lower in the tumor expression group, and that the survival of the low expression group was worse than that of the high expression group. In conclusion, the results of this study show that platinum resistance-related differentially expressed genes in lung adenocarcinoma are mainly concentrated in biological processes such as DNA recombination and response to lipopolysaccharide. The validation set proved that the high-risk group of our prognostic model had poor survival. M6A regulatory factor analysis, immune infiltration, and drug sensitivity analysis all showed differences between high and low-risk groups. ceRNA analysis showed that has-miR-374a-5p and RP6-24A23.7 could be protective factors. Further exploration of the potential impact of these genes on the risk and prognosis of drug-resistant patients with lung adenocarcinoma would provide theoretical support for future research.

Integrated bioinformatic analysis identified a novel prognostic pan-programmed cell death signature for bladder cancer

Programmed cell death (PCD) refers to a molecularly regulated form of cell death that functions as an essential anticancer defense mechanism and serves as a target of anticancer therapies. Multiple types of PCD comprehensively regulate tumorigenesis and tumor progression and metastasis. However, a systemic exploration of the multiple types of PCD in cancers, especially bladder cancer, is lacking. In this study, we evaluated the expression pattern of genes associated with multiple types of PCD in bladder cancer using the "ssGSEA" method and conceptualized the multiple types of PCD as being collectively involved in "Pan-PCD". Based on the differentially expressed genes related to Pan-PCD, we developed a Pan-PCD-related prognostic signature (PPRPS) to predict patient prognosis via univariate and multivariate Cox regression analysis. The PPRPS is an independent prognostic factor, and the AUC (Area Under Curve) for 3-year overall survival was 0.748. Combined with age and stage, PPRPS displayed excellent predictive ability. Based on the PPRPS, higher levels of immune cell infiltration, tumor microenvironment, and immune checkpoint molecules were observed in the high-PPRPS group. Furthermore, PPRPS enabled accurate risk prediction for metastatic urothelial carcinoma after anti-PD-L1 monoclonal antibody treatment. Patients in the high-PPRPS group had poor prognoses. Docetaxel, staurosporine, and luminespib were identified as potentially effective drugs for high-PPRPS bladder cancer patients. In summary, we developed the Pan-PCD signature to improve the accuracy of bladder cancer prognostic predictions and to provide a novel classification method to guide treatment selection.

RNA Modifications Meet Tumors

RNA modifications occur through the whole process of gene expression regulation, including transcription, translation, and post-translational processes. They are closely associated with gene expression, RNA stability, and cell cycle. RNA modifications in tumor cells play a vital role in tumor development and metastasis, changes in the tumor microenvironment, drug resistance in tumors, construction of tumor cell-cell "internet", etc. Several types of RNA modifications have been identified to date and have various effects on the biological characteristics of different tumors. In this review, we discussed the function of RNA modifications, including N 6-methyladenine (m6A), 5-methylcytosine (m5C), N 7-methyladenosine (m7G), N 1-methyladenosine (m1A), pseudouridine (Ψ), and adenosine-to-inosine (A-to-I), in the microenvironment and therapy of solid and liquid tumors.

Linc-UROD stabilizes ENO1 and PKM to strengthen glycolysis, proliferation and migration of pancreatic cancer cells

Pancreatic cancer (PC) is a fatal malignancy, threatening human health in worldwide. Long non-coding RNAs (lncRNAs) have been acknowledged to be essential regulators in various biological processes of human cancers. However, the role of some novel lncRNAs in PC remain to be explored. In this study, we focused on the function and molecular mechanism of a novel lncRNA linc-UROD (also named TCONS_00002016 or XLOC_000166) in PC. The expression of linc-UROD was found to be upregulated in PC cells. The results of loss-of-function assays demonstrated that linc-UROD knockdown suppressed cell proliferation and migration, induced cell cycle G0/G1 arrest, and accelerated apoptosis of PC cells. Through mechanistic experiments, we found that IGF2BP3 stabilized linc-UROD through METTL3-mediated m6A modification. In addition, linc-UROD enhances the stability of ENO1 and PKM through interacting with them to inhibit ubiquitination. Detection on glucose consumption, pyruvate kinase activity and lactate production indicated that linc-UROD accelerated glycolysis of PC cells through PKM/ENO1-mediated pathway. To summarize, linc-UROD stabilized by IGF2BP3/METTL3 contributes to glycolysis and malignant phenotype of PC cells by stabilizing ENO1 and PKM. The findings suggest that linc-UROD may be a novel therapeutic target for PC patients.

A Risk Signature Consisting of Eight m6A Methylation Regulators Predicts the Prognosis of Glioma

Glioma progression seriously correlates to the epigenetic context. This study aims to identify glioma subtypes by clustering analysis of patients using the multi-omics data of N6-methyladenosine (m⁶A) methylation regulators and to construct a risk signature for investigating the role of m⁶A methylation regulators in the prognosis of glioma. Multi-omics data of glioma and normal control tissues were obtained through The Cancer Genome Atlas (TCGA) database. The clustering analysis of multi-omics data of patients was conducted using the R package iClusterPlus software. The risk model was constructed by univariate and multivariate Cox analysis, and the glioma expression data and related clinical data were obtained by Chinese Glioma Genome Atlas (CGGA) datasets to verify the risk model. By analyzing the glioma data in TCGA, we found that the risk signature could be constructed according to the eight genes with m⁶A methylation modification function, including ALKBH5, HNRNPA2B1, IGF2BP2, IGF2BP3, RBM15, WTAP, YTHDF1, and YTHDF2. Meanwhile, we found that IGF2BP2 and IGF2BP3 were highly expressed in glioma subtypes with high-risk scores and closely related to the prognosis of glioma patients. m⁶A methylation regulators, especially IGF2BP2 and IGF2BP3, play important roles in the malignant progression of glioma. The risk signature constructed by eight m⁶A methylation regulators can predict the prognosis of glioma. IGF2BP2 and IGF2BP3 may be the key regulatory factors of m⁶A methylation regulators involved in the occurrence and development of glioma, and can serve as molecular markers for the prognosis of glioma.

m6A Methyltransferase METTL3 Reduces Hippocampal Neuron Apoptosis in a Mouse Model of Autism Through the MALAT1/SFRP2/Wnt/β-catenin Axis

OBJECTIVE

Hippocampal neuron apoptosis contributes to autism, while METTL3 has been documented to possess great potentials in neuron apoptosis. Our study probed into the role of METTL3 in neuron apoptosis in autism and to determine the underlying mechanism.

METHODS

Bioinformatics analysis was used to analyze expressed genes in autism samples. Institute of Cancer Research mice were treated with valproic acid to develop autism models. The function of METTL3 in autism-like symptoms in mice was analyzed with behavioral tests and histological examination of their hippocampal tissues. Primary mouse hippocampal neurons were extracted for in vitro studies. Downstream factors of METTL3 were explored and validated.

RESULTS

METTL3, MALAT1, and Wnt/β-catenin signaling were downregulated, while SFRP2 was upregulated in the hippocampal tissues of a mouse model of autism. METTL3 stabilized MALAT1 expression by promoting m6A modification of MALAT1. MALAT1 promoted SFRP2 methylation and led to reduced SFRP2 expression by recruiting DNMT1, DNMT3A, and DNMT3B to the promoter region of SFRP2. Furthermore, SFRP2 facilitated activation of the Wnt/β-catenin signaling. By this mechanism, METTL3 suppressed autism-like symptoms and hippocampal neuron apoptosis.

CONCLUSION

This research suggests that METTL3 can reduce autism-like symptoms and hippocampal neuron apoptosis by regulating the MALAT1/SFRP2/Wnt/β-catenin axis.

The role and regulatory mechanism of m6A methylation in the nervous system

N6-methyladenosine (m⁶A) modification regulates RNA translation, splicing, transport, localization, and stability at the post-transcriptional level. The m⁶A modification has been reported to have a wide range of effects on the nervous system, including neurogenesis, cerebellar development, learning, cognition, and memory, as well as the occurrence and development of neurological disorders. In this review, we aim to summarize the findings on the role and regulatory mechanism of m⁶A modification in the nervous system, to reveal the molecular mechanisms of neurodevelopmental processes, and to promote targeted therapy for nervous system-related diseases.

m6A Methylation Regulators Are Predictive Biomarkers for Tumour Metastasis in Prostate Cancer

Simple Summary

Recurrence and metastatic progression always lead to dismal outcomes in prostate cancer (PCa). There is no reliable biomarker for the prediction of recurrence and metastasis other than the Prostate Cancer Antigen (PCA). N6-methyladenosine (m⁶A) is the most common post-transcriptional mRNA modification and is regulated by m⁶A regulators dynamically. Since m⁶A modification is associated with cancer development and outgrowth, we performed a consensus clustering on PCa with regard to the gene expression of all m⁶A regulators. We identified three subtypes of Pca with distinct m⁶A expression patterns and enriched biological pathways. We also established an m⁶A score for metastasis prediction based on our clustering, which is potentially a predictive biomarker for Pca metastasis.

Abstract

Prostate cancer (PCa) is one of the most common cancers in men. Usually, most PCas at initial diagnosis are localized and hormone-dependent, and grow slowly. Patients with localized PCas have a nearly 100% 5-year survival rate; however, the 5-year survival rate of metastatic or progressive PCa is still dismal. N6-methyladenosine (m⁶A) is the most common post-transcriptional mRNA modification and is dynamically regulated by m⁶A regulators. A few studies have shown that the abnormal expression of m⁶A regulators is significantly associated with cancer progression and immune cell infiltration, but the roles of these regulators in PCa remain unclear. Here, we examined the expression profiles and methylation levels of 21 m⁶A regulators across the Cancer Genome Atlas (TCGA), 495 PCas by consensus clustering, and correlated the expression of m⁶A regulators with PCa progression and immune cell infiltration. Consensus clustering was applied for subtyping Pca samples into clusters based on the expression profiles of m⁶A regulators. Each subtype’s signature genes were obtained by a pairwise differential expression analysis. Featured pathways of m⁶A subtypes were predicted by Gene Ontology. The m⁶A score was developed to predict m⁶A activation. The association of the m⁶A score with patients’ survival, metastasis and immune cell infiltration was also investigated. We identified three distinct clusters in PCa based on the expression profiles of 21 m⁶A regulators by consensus clustering. The differential expression and pathway analyses on the three clusters uncovered the m⁶A regulators involved in metabolic processes and immune responses in PCa. Moreover, we developed an m⁶A score to evaluate the m⁶A regulator activation for PCa. The m⁶A score is significantly associated with Gleason scores and metastasis in PCa. The predictive capacity of the m⁶A score on PCa metastasis was also validated in another independent cohort with an area under the curve of 89.5%. Hence, our study revealed the critical role of m⁶A regulators in PCa progression and the m⁶A score is a promising predictive biomarker for PCa metastasis.

The role, mechanism, and application of RNA methyltransferase METTL14 in gastrointestinal cancer

Gastrointestinal cancer is the most common human malignancy characterized by high lethality and poor prognosis. Emerging evidences indicate that N6-methyladenosine (m6A), the most abundant post-transcriptional modification in eukaryotes, exerts important roles in regulating mRNA metabolism including stability, decay, splicing, transport, and translation. As the key component of the m6A methyltransferase complex, methyltransferase-like 14 (METTL14) catalyzes m6A methylation on mRNA or non-coding RNA to regulate gene expression and cell phenotypes. Dysregulation of METTL14 was deemed to be involved in various aspects of gastrointestinal cancer, such as tumorigenesis, progression, chemoresistance, and metastasis. Plenty of findings have opened up new avenues for exploring the therapeutic potential of gastrointestinal cancer targeting METTL14. In this review, we systematically summarize the recent advances regarding the biological functions of METTL14 in gastrointestinal cancer, discuss its potential clinical applications and propose the research forecast.

N6-Methyladenosine RNA-Binding Protein YTHDF1 in Gastrointestinal Cancers: Function, Molecular Mechanism and Clinical Implication

Simple Summary

N⁶-methyladenosine (m⁶A) is the most abundant internal modification in eukaryotic mRNA and plays a crucial role in the occurrence and development of diseases. YTHDF1 is the most powerful and abundant m⁶A-encoded RNA reader. In this review, we summarize the evidence of the involvement of YTHDF1 in gastrointestinal cancers, its molecular mechanisms of action, and therapeutic implications.

Abstract

N⁶-methyladenosine (m⁶A) is the most abundant internal modification in eukaryotic cell mRNA, and this modification plays a key role in regulating mRNA translation, splicing, and stability. Emerging evidence implicates aberrant m⁶A as a crucial player in the occurrence and development of diseases, especially GI cancers. Among m⁶A regulators, YTHDF1 is the most abundant m⁶A reader that functionally connects m⁶A-modified mRNA to its eventual fate, mostly notably protein translation. Here, we summarized the function, molecular mechanisms, and clinical implications of YTHDF1 in GI cancers. YTHDF1 is largely upregulated in multiple GI cancer and its high expression predicts poor patient survival. In vitro and in vivo experimental evidence largely supports the role of YTDHF1 in promoting cancer initiation, progression, and metastasis, which suggests the oncogenic function of YTHDF1 in GI cancers. Besides, YTHDF1 overexpression is associated with changes in the tumor microenvironment that are favorable to tumorigenesis. Mechanistically, YTHDF1 regulates the expression of target genes by promoting translation, thereby participating in cancer-related signaling pathways. Targeting YTHDF1 holds therapeutic potential, as the overexpression of YTHDF1 is associated with tumor resistance to chemotherapy and immunotherapy. In summary, YTHDF1-mediated regulation of m⁶A modified mRNA is an actionable target and a prognostic factor for GI cancers.

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.

N(6)-methyladenosine-mediated miR-380-3p maturation and upregulation promotes cancer aggressiveness in pancreatic cancer

N(6)-methyladenosine (m6A)-modified microRNAs (miRNAs) are relevant to cancer progression. Also, although the involvement of miR-380-3p in regulating cancer progression in bladder cancer and neuroblastoma has been preliminarily explored, its role in other types of cancer, such as pancreatic cancer (PC), has not been studied. Thus, this study aimed to investigate the role of miR-380-3p in regulating PC progression. Here, through performing Real-Time qPCR, we evidenced that miR-380-3p was significantly upregulated in the clinical pancreatic cancer tissues and cells compared to their normal counterparts. Interestingly, miR-380-3p was enriched with m6A modifications, and elimination of m6A modifications by deleting METTL3 and METTL14 synergistically suppressed miR-380-3p expressions in PC cells. Next, the gain and loss-of-function experiments verified that knockdown of miR-380-3p suppressed cell proliferation, epithelial-mesenchymal transition (EMT), and tumorigenesis in PC cells in vitro and in vivo, whereas miR-380-3p overexpression had opposite effects. Furthermore, the underlying mechanisms were uncovered, and our data suggested that miR-380-3p targeted the 3' untranslated regions (3'UTRs) of PTEN for its inhibition and degradation, resulting in the activation of the downstream Akt signal pathway. Moreover, the rescuing experiments validated that both PTEN overexpression and Akt pathway inhibitor LY294002 abrogated the promoting effects of miR-380-3p overexpression on cancer aggressiveness in PC cells. Collectively, this study firstly investigated the role of the m6A-associated miR-380-3p/PTEN/Akt pathway in regulating PC progression, which provided novel therapeutic and diagnostic biomarkers for this cancer.

Transcriptome-Wide N6-Methyladenosine Methylome Alteration in the Rat Spinal Cord After Acute Traumatic Spinal Cord Injury

Recent studies showed that RNA N6-methyladenosine (m6A) plays an important role in neurological diseases. We used methylated RNA immunoprecipitation sequencing (MeRIP-Seq) technology to generate the m6A modification map after traumatic spinal cord injury (TSCI). A total of 2,609 differential m6A peaks were identified after TSCI. Our RNA sequencing results after TSCI showed 4,206 genes with significantly altered expression. Cross-link analysis of m6A sequencing results and RNA sequencing results showed that 141 hyper-methylated genes were upregulated, 53 hyper-methylated genes were downregulated, 57 hypo-methylated genes were upregulated, and 197 hypo-methylated genes were downregulated. Among these, the important inflammatory response factor Tlr4 and the important member of the neurotrophin family Ngf were both upregulated and hyper-methylated after TSCI. This study provides that in the future, the epigenetic modifications of the genes could be used as an indicator of TSCI.

The Potential Role of m6A RNA Methylation in the Aging Process and Aging-Associated Diseases

N6-methyladenosine (m⁶A) is the most common and conserved internal eukaryotic mRNA modification. m⁶A modification is a dynamic and reversible post-transcriptional regulatory modification, initiated by methylase and removed by RNA demethylase. m⁶A-binding proteins recognise the m⁶A modification to regulate gene expression. Recent studies have shown that altered m⁶A levels and abnormal regulator expression are crucial in the ageing process and the occurrence of age-related diseases. In this review, we summarise some key findings in the field of m⁶A modification in the ageing process and age-related diseases, including cell senescence, autophagy, inflammation, oxidative stress, DNA damage, tumours, neurodegenerative diseases, diabetes, and cardiovascular diseases (CVDs). We focused on the biological function and potential molecular mechanisms of m⁶A RNA methylation in ageing and age-related disease progression. We believe that m⁶A modification may provide a new target for anti-ageing therapies.

Signature of m5C-Related lncRNA for Prognostic Prediction and Immune Responses in Pancreatic Cancer

BACKGROUND

Pancreatic cancer (PC) has a high mortality and dismal prognosis, predicting to be the second most lethal malignancy. 5-Methylcytosine (m5C) and long noncoding RNAs (lncRNAs) are both crucial in the prognostic outcome and immunotherapeutic effect for PC patients. Therefore, we aimed to create an m5C-related lncRNA signature (m5C-LS) for PC patients' prognosis and treatment.

METHODS

Clinicopathological information and RNAseq data were acquired from The Cancer Genome Atlas (TCGA) database. Pearson's correlation analysis was used to extract m5C-related lncRNAs in PC. Univariate, least absolute shrinkage and selection operator (LASSO), and multivariate Cox analyses were adopted to build an m5C-LS. Kaplan-Meier (K-M), principal component analysis (PCA), and nomogram were utilized to assess model accuracy. In addition, we explored the model's possible immunotherapeutic responses and drug sensitivity targets.

RESULTS

Three m5C-related lncRNAs were finally established to construct the risk signature, which has a good and independent predictive ability for PC patients. Based on the m5C-LS, patients were classified into the low- and high-m5C-LS group, with the latter having a worse prognosis. Furthermore, the m5C-LS allowed us to better discriminate the immunotherapeutic responses of PC patients in different subgroups.

CONCLUSIONS

Our study constructed an m5C-LS and established a nomogram model that accurately predicted the prognosis of PC patients, as well as provides promising immunotherapeutic strategies in the future.

RNA demethylase ALKBH5 in cancer: from mechanisms to therapeutic potential

RNA demethylase ALKBH5 takes part in the modulation of N6-methyladenosine (m6A) modification and controls various cell processes. ALKBH5-mediated m6A demethylation regulates gene expression by affecting multiple events in RNA metabolism, e.g., pre-mRNA processing, mRNA decay and translation. Mounting evidence shows that ALKBH5 plays critical roles in a variety of human malignancies, mostly via post-transcriptional regulation of oncogenes or tumor suppressors in an m6A-dependent manner. Meanwhile, increasing non-coding RNAs are recognized as functional targets of ALKBH5 in cancers. Here we reviewed up-to-date findings about the pathological roles of ALKBH5 in cancer, the molecular mechanisms by which it exerts its functions, as well as the underlying mechanism of its dysregulation. We also discussed the therapeutic implications of targeting ALKBH5 in cancer and potential ALKBH5-targeting strategies.

Comprehensive Analysis of m6A RNA Methylation Regulators and the Immune Microenvironment to Aid Immunotherapy in Pancreatic Cancer

Pancreatic cancer (PAAD) is one of the most malignant cancers and immune microenvironment has been proved to be involved in pathogenesis of PAAD. m6A modification, related to the expression of m6A regulators, participates in the development of multiple cancers. However, the correlation between m6A regulators and immune microenvironment was largely unknown in PAAD. And because of the small sample size of pancreatic cancer in the TCGA database, it is not enough to draw a convincing conclusion. In the present study, we downloaded seven pancreatic cancer datasets with survival data and removed batch effects among these datasets to be used as the PAAD cohort to analyze the immune landscape of PAAD and the expression pattern of m6A regulators and divided the integrated dataset into cluster 1 and cluster 2 by consensus clustering for m6A regulators. Lower m6A regulators were found to be related to higher immune cell infiltration and a better survival. Moreover, we identified six m6A regulators and constructed the prognostic signature of m6A regulators. Patients with low-risk score had a higher response to immune checkpoint inhibitor and a longer overall survival. To figure out the underlying mechanism, we analyzed the cancer immunity cycle, most altered genes, gene set enrichment analysis (GSEA) and gene set variation analysis (GSVA) in risk subtypes. In summary, the present study proved m6A regulators modulated the PAAD immune microenvironment. And risk scores served as predictive indicator for immunotherapy and played a prognostic role for PAAD patients. Our study provided novel therapeutic targets to improve immunotherapy efficacy.

YTH domain family: potential prognostic targets and immune-associated biomarkers in hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is the most common high malignancy with insidious onset, invasive fast-growing, high recurrence rate and fatality. YTH domain family plays essential roles in development of HCC. However, the biological function of YTH domain family in HCC have not been clarified. Here, through evaluating the expression profiles of YTH domain family, we found that upregulated YTHDF1 might be more significant and valuable in development and progression of HCC. There was a strong correlation between YTHDC1, YTHDF1 and YTHDF2 and pathological stage of HCC patients. Kaplan-Meier plotter revealed that HCC patients with high level of YTHDF1 and YTHDF2 were highly related to a shorter overall survival time, and low level of YTHDF1 (p = 0.0017) has an important association with a longer progression-free survival time. Genetic alterations using cBioPortal revealed that the alteration rates of YTHDF3 were the highest. We also found that the functions of YTH domain family were linked to several cancer-associated pathways, including peptidyl-serine modification, peptidyl-tyrosine modification and negative regulation of cellular component movement. TIMER database indicated that the YTH domain family had a strong relationship with the infiltration of six types of immune cells (macrophages, neutrophils, CD8+ T-cells, B-cells, CD4+ T-cells and dendritic cells). Next, Ualcan databases revealed that the global methylation levels of YTHDC1 was higher in HCC patients, while YTHDF2 was lower in HCC patients. In conclusion, our findings will enhance the understanding of YTH domain family in HCC pathology, and provide novel insights into YTH-targeted therapy for HCC patients.

Classification and immune invasion analysis of breast cancer based on m6A genes

Background

Breast cancer (BRCA) shows genetic, epigenetic, and phenotypic diversity. Methylation of N6-methyladenosine (m6A) affects the occurrence, development, and therapeutic efficacy of BRCA. However, the characteristics and prognostic value of m6A in BRCA remain unclear. We aimed to classify and construct a scoring system for the m6A regulatory gene in BRCA, and to explore its potential mechanisms.

Methods

In this study, we selected 23 m6A regulatory genes and analyzed their genetic variation in BRCA, including copy number variation (CNV) data, expression differences, mutations, gene types, and correlations between genes. Survival curves were drawn by the Kaplan-Meier method, and a log-rank P<0.05 was considered statistically significant. The partitioning around medoids (PAM) algorithm was used for molecular subtype analysis of m6A, single-sample Gene Set Enrichment Analysis (ssGSEA) algorithm was used to quantify the relative infiltration levels of various immune cell subgroups, and a scoring system was built based on principal component analysis (PCA).

Results

In BRCA, m6A regulatory gene mutation frequency is not high, while that of CNV mutation is high, which is related to gene expression and closely related to prognosis. In this study, we identified 3 different m6A subtypes, which are closely related to the level of immune cell infiltration. We further constructed an m6A score system, in which lower scores were correlated with low tumor mutation burden (TMB), later clinical staging, programmed cell death 1 ligand 1 (PD-L1) expression, and triple-negative breast cancer (TNBC).

Conclusions

This study highlights the diversity and complexity of the role of m6A in BRCA. The classification of BRCA based on the m6A regulatory gene can help us understand the characteristics of BRCA and help develop individualized immunotherapy regimens.

Radiogenomics Map Reveals the Landscape of m6A Methylation Modification Pattern in Bladder Cancer

We aimed to develop a noninvasive radiomics approach to reveal the m6A methylation status and predict survival outcomes and therapeutic responses in patients. A total of 25 m6A regulators were selected for further analysis, we confirmed that expression level and genomic mutations rate of m6A regulators were significantly different between cancer and normal tissues. Besides, we constructed methylation modification models and explored the immune infiltration and biological pathway alteration among different models. The m6A subtypes identified in this study can effectively predict the clinical outcome of bladder cancer (including m6AClusters, geneClusters, and m6Ascore models). In addition, we observed that immune response markers such as PD1 and CTLA4 were significantly corelated with the m6Ascore. Subsequently, a total of 98 obtained digital images were processed to capture the image signature and construct image prediction models based on the m6Ascore classification using a radiomics algorithm. We constructed seven signature radiogenomics models to reveal the m6A methylation status, and the model achieved an area under curve (AUC) degree of 0.887 and 0.762 for the training and test datasets, respectively. The presented radiogenomics models, a noninvasive prediction approach that combined the radiomics signatures and genomics characteristics, displayed satisfactory effective performance for predicting survival outcomes and therapeutic responses of patients. In the future, more interdisciplinary fields concerning the combination of medicine and electronics remains to be explored.

Dynamic m6A-ncRNAs association and their impact on cancer pathogenesis, immune regulation and therapeutic response

Several types of modifications have been proven to participate in the metabolism and processing of different RNA types, including non-coding RNAs (ncRNAs). N-6-methyladenosine (m6A) is a dynamic and reversible RNA modification that is closely involved in the ncRNA homeostasis, and serves as a crucial regulator for multiple cancer-associated signaling pathways. The ncRNAs usually regulate the epigenetic modification, mRNA transcription and other biological processes, displaying enormous roles in human cancers. In this review, we summarized the significant implications of m6A-ncRNA interaction in various types of cancers. In particular, the interplay between m6A and ncRNAs in cancer pathogenesis and therapeutic resistance are being widely recognized. We also discussed the relevance of m6A-ncRNA interaction in immune regulation, followed by the interference on cancer immunotherapeutic procedures. In addition, we briefly highlighted the computation tools that could identify the accurate features of m6A methylome among ncRNAs. In summary, this review would pave the way for a better understanding of the biological functions of m6A-ncRNA crosstalk in cancer research and treatment.

Identification of N6-Methyladenosine-Associated Long Non-coding RNAs for Immunotherapeutic Response and Prognosis in Patients With Pancreatic Cancer

Pancreatic cancer is a highly aggressive disease with poor prognosis. N6-methyladenosine (m6A) is critical for post-transcriptional modification of messenger RNA (mRNA) and long non-coding RNA (lncRNA). However, the m6A-associated lncRNAs (m6A-lncRNA) and their values in predicting clinical outcomes and immune microenvironmental status in pancreatic cancer patients remain largely unexplored. This study aimed to evaluate the importance of m6A-lncRNA and established a m6A-lncRNA signature for predicting immunotherapeutic response and prognosis of pancreatic cancer. The m6A-lncRNA co-expression networks were constructed using data from the TCGA and GTEx database. Based on the least absolute shrinkage and selection operator (LASSO) analysis, we constructed an 8 m6A-lncRNA signature risk model, and selection operator (LASSO) analysis, and stratified patients into the high- and low-risk groups with significant difference in overall survival (OS) (HR = 2.68, 95% CI = 1.74-4.14, P < 0.0001). Patients in the high-risk group showed significantly reduced OS compared to patients in the low-risk group (P < 0.001). The clinical characteristics and m6A-lncRNA risk scores were used to construct a nomogram which accurately predicted the OS in pancreatic cancer. TIMER 2.0 were used to investigate tumor immune infiltrating cells and its relationship with pancreatic cancer. CIBERSORT analysis revealed increased higher infiltration proportions of M0 and M2 macrophages, and lower infiltration of naive B cell, CD8⁺ T cell and Treg cells in the high-risk group. Compared to the low-risk group, functional annotation using ssGSEA showed that T cell infiltration and the differential immune-related check-point genes are expressed at low level in the high-risk group (P < 0.05). In summary, our study constructed a novel m6A-associated lncRNAs signature to predict immunotherapeutic responses and provided a novel nomogram for the prognosis prediction of pancreatic cancer.

Integrative bioinformatics and experimental analysis revealed down-regulated CDC42EP3 as a novel prognostic target for ovarian cancer and its roles in immune infiltration

Ovarian cancer is a significant clinical challenge as no effective treatments are available to enhance patient survival. Recently, N6-methyladenosine (m⁶A) RNA modification has been demonstrated to play a pivotal role in tumorigenesis and progression. However, the roles of m⁶A target genes in ovarian cancer haven’t been clearly illustrated. In this study, we presented a comprehensive bioinformatics and in vitro analysis to evaluate the roles of m⁶A target genes. Cell division cycle 42 effector protein 3 (CDC42EP3), one probable m6A target gene, was identified to be down-regulated in ovarian cancer tissues and cells. Meanwhile, quantitative PCR (qPCR) and western blot were used to confirm the down-regulated CDC42EP3 in ovarian cancer cells A2780 and TOV112D. The biological function of CDC42EP3 in ovarian cancer was further validated with several algorithms, such as PrognoScan, K-M plotter, LinkedOmics and TISIDB. These findings indicated that lower expression of CDC42EP3 was correlated with poor prognosis in patients with ovarian cancer. In addition, CDC42EP3 expression was significantly associated with a diverse range of tumor-infiltrating immune cells, including natural killer cells (NK), T central memory cells (Tcm), T gamma delta cells (Tgd), etc. Taken together, this study uncovered the potential roles of m⁶A target gene CDC42EP3 in the regulation of immune microenvironment in the ovarian cancer, and identified CDC42EP3 as a novel prognostic target.

Characterization of Modification Patterns, Biological Function, Clinical Implication, and Immune Microenvironment Association of m6A Regulators in Pancreatic Cancer

Objective: N⁶-methyladenosine (m⁶A) modification may modulate various biological processes. Nonetheless, clinical implications of m⁶A modification in pancreatic cancer are undefined. Herein, this study comprehensively characterized the m⁶A modification patterns in pancreatic cancer based on m⁶A regulators. Methods: Genetic mutation and expression pattern of 21 m⁶A regulators and their correlations were assessed in pancreatic cancer from TCGA dataset. m⁶A modification patterns were clustered using unsupervised clustering analysis in TCGA and ICGC datasets. Differences in survival, biological functions and immune cell infiltrations were assessed between modification patterns. A m⁶A scoring system was developed by principal component analysis. Genetic mutations and TIDE scores were compared between high and low m⁶A score groups. Results: ZC3H13 (11%), RBM15B (9%), YTHDF1 (8%), and YTHDC1 (6%) frequently occurred mutations among m⁶A regulators. Also, most of regulators were distinctly dysregulated in pancreatic cancer. There were tight crosslinks between regulators. Two m⁶A modification patterns were constructed, with distinct prognoses, immune cell infiltration and biological functions. Furthermore, we quantified m⁶A score in each sample. High m⁶A scores indicated undesirable clinical outcomes. There were more frequent mutations in high m⁶A score samples. Lower TIDE score was found in high m⁶A score group, with AUC = 0.61, indicating that m⁶A scores might be used for predicting the response to immunotherapy. Conclusion: Collectively, these data demonstrated that m⁶A modification participates pancreatic cancer progress and ornaments immune microenvironment, providing an insight into pancreatic cancer pathogenesis and facilitating precision medicine development.

The bioinformatics and experimental analysis of AlkB family for prognosis and immune cell infiltration in hepatocellular carcinoma

Background

Serving as N6-methyladenosine demethylases, the AlkB family is involved in the tumorigenesis of hepatocellular carcinoma (HCC). However, the molecular profiles and clinical values of the AlkB family in HCC are not well known.

Methods

Several bioinformatics tools and in vitro experiments were used to identify the immune-related profiles and prognostic values of AlkB family in HCC.

Results

In this study expression levels of ALKBH1/2/3/4/7 were all remarkably increased in HCC tissues when compared with normal tissues. Quantitative PCR (qPCR) and immunohistochemistry were used to validate the expression of AlkB family members in HCC tissues and normal liver tissues. In addition, high expression levels of ALKBH4 were negatively correlated with overall survival (OS) and disease-free survival (DFS) in patients with HCC. Increased ALKBH4 was also associated with pathological stage in HCC patients. The molecular profiles of AlkB family in HCC were mainly associated with peptidyl-serine modification, peptidyl-tyrosine modification, regulation of metal ion transport, etc. Furthermore, tumor-infiltrating immune cell analysis indicated that ALKBH1/2/3/4/5/6/7/8 and FTO were related to the infiltration of different immune cell, such as CD8+ T cells, macrophages, neutrophils, dendritic cells and CD4+ T cells. We also discovered that the methylation levels of ALKBH1/2/4/5/6/8 and FTO were remarkably reduced in HCC tissues.

Conclusions

Collectively, our findings may deepen the understanding of specific molecular profiles of the AlkB family in HCC pathology. In particular, ALKBH4 could serve as a promising prognostic candidate for treating HCC, and these results might potentiate the development of more reliable therapeutic strategies for patients with HCC.

Identification of an N6-methyladenosine (m6A)-related signature associated with clinical prognosis, immune response, and chemotherapy in primary glioblastomas

Background

N6-methyladenosine (m6A) RNA methylation regulators play crucial role in tumorigenicity and progression. However, their biological significance in primary glioblastomas (GBM) has not been fully elucidated.

Methods

In the present study, we evaluated the 22 m6A RNA regulators using the integrated data of primary GBM samples from The Cancer Genome Atlas and Chinese Glioma Genome Atlas databases. The different m6A modification patterns and m6A-related gene signature in primary GBM were distinguished by using principal component analysis. Single-sample gene set enrichment analysis was introduced to assess the relative level of immune infiltration. Gene set variation analysis was performed to calculate the enrichment score of the signaling pathways for different clusters. An m6A scoring scheme was established to evaluate the m6A modification pattern in individual tumors in order to predict prognosis and evaluate tumor microenvironment (TME) cell infiltration, immune response, and chemotherapy effect in primary GBM.

Results

Two distinct m6A modification subgroups associated with different clinical features and biological pathways were identified among the 371 primary GBM. Based on 132 prognostic m6A phenotype-related differentially expressed genes (DEGs) between 2 m6A cluster subgroups, an m6A scoring model was constructed to assess the m6A modification pattern in individual tumors. The high-m6A score group was associated with better prognosis and immune response and worse chemotherapy effect.

Conclusions

The findings of the present study indicate the potential role of m6A modification in primary GBM, which will help enhance our understanding of TME characteristics, predict clinical prognosis, and provide important insight into effective immunotherapy and chemotherapy.

Emerging Roles of N6-Methyladenosine Demethylases and Its Interaction with Environmental Toxicants in Digestive System Cancers

Digestive system cancers are common cancers with high cancer deaths worldwide. They have become a major threat to public health and economic burden. As one of the most universal RNA modifications in eukaryotes, the N6-methyladenosine (m6A) modification is involved in the occurrence, development, prognosis, and treatment response of various cancers, including digestive system cancers. M6A demethylases shape the m6A landscape dynamically, playing important roles in cancers. In addition, accumulating evidence reveal that many environmental toxicants are the established risk factors for digestive system cancers and associated with m6A modification. In this review, we summarize the multiple functions of M6A demethylases (fat mass and obesity-associated protein (FTO), AlkB homolog 5 (ALKBH5) and AlkB homolog 3 (ALKBH3)) in digestive system cancers, which are aberrantly expressed and affect cancer progression. We also discuss the potential roles of m6A demethylases in the assessment of environmental exposure, the signature for prevention and diagnosis of digestive system cancers.

The tumor-suppressive effects of alpha-ketoglutarate-dependent dioxygenase FTO via N6-methyladenosine RNA methylation on bladder cancer patients

N6-methyladenosine (m6A) methylation participates in the progression of bladder cancer (BCa). Nevertheless, the regulatory mechanism of alpha-ketoglutarate-dependent dioxygenase FTO influencing the BCa progression has still remained elusive. In this study, to investigate the tumor-suppressive effects of FTO via m6A RNA methylation on BCa patients, a total of 15 cancer tissues and adjacent normal tissues (ANTs) were collected from BCa patients who received tumor resection in our hospital from September 2015 to December 2019. We found that the FTO expression was significantly reduced in cancer tissues compared with that in ANTs, which indicated a lower malignant potential and a higher overall survival rate. It was revealed that overexpression of FTO in two human urinary BCa cell lines (HT-1197 and HT-1376) significantly decreased the cell proliferation and invasion abilities compared with the negative controls, whereas the cell apoptosis was markedly enhanced. In addition, we noted that the changes in m6A methylation level mainly appeared at 5ʹ untranslated region (5ʹ UTR) of MALAT1 and NOTCH1 transcripts, and at 3ʹ UTR of CSNK2A2 and ITGA6 transcripts, responding to the overexpression of FTO. Mechanistically, we found that the splicing factor, proline- and glutamine-rich (SFPQ) could influence the FTO-mediated m6A RNA demethylation, eventually affecting the gene expression. This study provided a new insight into the relationship between the FTO expression and the m6A RNA methylation, assisting scholars to better understand the pathogenesis of BCa.

Expression of N6-methyladenosine (m6A) regulators correlates with immune microenvironment characteristics and predicts prognosis in diffuse large cell lymphoma (DLBCL)

This study conducted a comprehensive analysis of the clinical significance of N⁶-methyladenosine (m⁶A) regulators and their relationship with immune microenvironment characteristics in diffuse large cell lymphoma (DLBCL). Consensus clustering was performed to molecularly discriminate DLBCL subtypesbased on m⁶A regulators’ expression. Using the Cox and Lasso regression algorithm, survival-associated m⁶A regulators were identified, and a m⁶A-based prognostic signature was established. The influence of m⁶A risk on immune cell infiltration, immune checkpoint genes, cancer immunity cycle, and immunotherapeutic response was evaluated. Potential molecular pathways related to m⁶A risk were investigated using gene set enrichment analysis. The m⁶A regulators showed satisfactory performance in distinguishing DLBCL subgroups with distinct clinical traits and outcomes. A six m⁶A regulator-based prognostic signature was established and validated as an independent predictor, which separated patients into low- and high-risk groups. High-risk m⁶A indicated worse survival. The B cells naïve, T cells gamma delta, and NK cells resting were the three most affected immune cells by m⁶A risk. Up-regulated (PDCD1 and KIR3DL1) and down-regulated (TIGIT, IDO1, and BTLA) immune checkpoint genes in the high-risk group were identified. The m⁶A risk was found to influence several steps in the cancer immunity cycle. Patients with high-risk m⁶A were more likely to benefit from immunotherapy. Biological function enrichment analysis revealed that high-risk m⁶A to be tended related to malignant tumor characteristics, while low-risk m⁶A showed trend to be related to defensive response processes. Collectively, the m⁶A-based prognostic signature could be a practical prognostic predictor for DLBCL and immune microenvironment characteristics affected by m⁶A may be part of the mechanism.

Function and clinical significance of N6-methyladenosine in digestive system tumours

RNA modification, like DNA methylation, histone modification, non-coding RNA modification and chromatin rearrangement, plays an important role in tumours. N6-methyladenosine (m6A) is the most abundant RNA modification in cells, and it regulates RNA transcription, processing, splicing, degradation, and translation. m6A-associated proteins have been used as new biomarkers and therapeutic targets for tumour prediction and monitoring. There are three main types of proteins involved in m6A methylation: methyltransferases (METTL3, METTL14, WTAP, RBM15, ZC3H13 and KIAA1429), demethylases (FTO, ALKBH5 and ALKBH3) and RNA-binding proteins (YTHDF1-3, YTHDC1-2, IGF2BPs and HNRNPs). This article reviews the origins, characteristics and functions of m6A and its relationship with digestive system tumours based on recent research. The expression of m6A regulators can be used as an evaluation indicator of tumour growth and progression and as a prognostic indicator. In-depth research on m6A methylation in digestive system tumours may provide new directions for clinical prediction and further treatment.

Development and validation of a novel N6-methyladenosine (m6A)-related multi- long non-coding RNA (lncRNA) prognostic signature in pancreatic adenocarcinoma

Accumulating evidence has unveiled the pivotal roles of N6-methyladenosine (m6A) in pancreatic adenocarcinoma (PAAD). However, there are not many researches to predict the prognosis of PAAD using m6A-related long non-coding RNAs (lncRNAs). Raw data from The Cancer Genome Atlas (TCGA), International Cancer Genome Consortium (ICGC), and the Genotype-Tissue Expression project (GTEx) were utilized to comprehensively analyze the expression and prognostic performances of 145 m6A-related lncRNAs in PAAD and to develop and validate a novel m6A-related multi-lncRNA prognostic signature (m6A-LPS) for PAAD patients. In total, 57 differentially expressed m6A-related lncRNAs with prognostic values were identified. Based on LASSO-Cox regression analysis, m6A-LPS was constructed and verified by using five-lncRNA expression profiles for TCGA and ICGC cohorts. PAAD patients were then divided into high- and low-risKBIE_A_1933868k subgroups with different clinical outcomes according to the median risk score; this was further verified by time-dependent receiver operating characteristic curves. Risk scores were significantly associated with clinical parameters such as histological grade and cancer status among PAAD patients. A nomogram consisting of risk score, grade, and cancer status was generated to predict the survival probability of PAAD patients, as also demonstrated by calibration curves. Discrepancies in cellular processes, signaling pathways, and immune status between the high- and low-risk subgroups were investigated by functional and single-sample gene set enrichment analyses. In conclusion, the novel m6A-LPS for PAAD patients was developed and validated, which might provide new insight into clinical decision-making and precision medicine.

N6-methyladenosine (m6A) in pancreatic cancer: Regulatory mechanisms and future direction

N6-methyladenosine (m6A), the most abundant RNA modification in eukaryotes, plays a pivotal role in regulating many cellular and biological processes. Aberrant m6A modification has recently been involved in carcinogenesis in various cancers, including pancreatic cancer. Pancreatic cancer is one of the deadliest cancers. It is a heterogeneous malignant disease characterized by a plethora of diverse genetic and epigenetic events. Increasing evidence suggests that dysregulation of m6A regulatory factors, such as methyltransferases, demethylases, and m6A-binding proteins, profoundly affects the development and progression of pancreatic cancer. In addition, m6A regulators and m6A target transcripts may be promising early diagnostic and prognostic cancer biomarkers, as well as therapeutic targets. In this review, we highlight the biological functions and mechanisms of m6A in pancreatic cancer and discuss the potential of m6A modification in clinical applications.

The Role of RNA Modifications and RNA-modifying Proteins in Cancer Therapy and Drug Resistance

The advent of new genome-wide sequencing technologies has uncovered abnormal RNA modifications and RNA editing in a variety of human cancers. The discovery of reversible RNA N6-methyladenosine (RNA: m⁶A) by fat mass and obesity-associated protein (FTO) demethylase has led to exponential publications on the pathophysiological functions of m⁶A and its corresponding RNA modifying proteins (RMPs) in the past decade. Some excellent reviews have summarized the recent progress in this field. Compared to the extent of research into RNA: m⁶A and DNA 5-methylcytosine (DNA: m⁵C), much less is known about other RNA modifications and their associated RMPs, such as the role of RNA: m⁵C and its RNA cytosine methyltransferases (RCMTs) in cancer therapy and drug resistance. In this review, we will summarize the recent progress surrounding the function, intramolecular distribution and subcellular localization of several major RNA modifications, including 5' cap N7-methylguanosine (m7G) and 2'-O-methylation (Nm), m⁶A, m⁵C, A-to-I editing, and the associated RMPs. We will then discuss dysregulation of those RNA modifications and RMPs in cancer and their role in cancer therapy and drug resistance.

Identification of m6A methyltransferase-related lncRNA signature for predicting immunotherapy and prognosis in patients with hepatocellular carcinoma

N6-methyladenosine (m6A) methyltransferase has been shown to be an oncogene in a variety of cancers. Nevertheless, the relationship between the long non-coding RNAs (lncRNAs) and hepatocellular carcinoma (HCC) remains elusive. We integrated the gene expression data of 371 HCC and 50 normal tissues from The Cancer Genome Atlas (TCGA) database. Differentially expressed protein-coding genes (DE-PCGs)/lncRNAs (DE-lncRs) analysis and univariate regression and Kaplan–Meier (K–M) analysis were performed to identify m6A methyltransferase-related lncRNAs. Three prognostic lncRNAs were selected by univariate and LASSO Cox regression analyses to construct the m6A methyltransferase-related lncRNA signature. Multivariate Cox regression analyses illustrated that this signature was an independent prognostic factor for overall survival (OS) prediction. The Gene Set Enrichment Analysis (GSEA) suggested that the m6A methyltransferase-related lncRNAs were involved in the immune-related biological processes (BPs) and pathways. Besides, we discovered that the lncRNAs signature was correlated with the tumor microenvironment (TME) and the expression of critical immune checkpoints. Tumor Immune Dysfunction and Exclusion (TIDE) analysis revealed that the lncRNAs could predict the clinical response to immunotherapy. Our study had originated a prognostic signature for HCC based on the potential prognostic m6A methyltransferase-related lncRNAs. The present study had deepened the understanding of the TME status of HCC patients and laid a theoretical foundation for the choice of immunotherapy.

The role of M6A modification in the regulation of tumor-related lncRNAs

N⁶-methyladenosine (m⁶A) is the most abundant modification in eukaryotic cells, and it regulates RNA transcription, processing, splicing, degradation, and translation. Long non-coding RNAs (lncRNAs), as transcriptional products with no or limited protein coding ability more than 200 nt in length, play an important role in epigenetic modification, mRNA transcription, splicing, stability, translation, and other biological functions. Extensive studies have shown that both m⁶A modification and lncRNAs are involved in the pathogenesis of various diseases, such as kinds of cancers, heart failure, Alzheimer’s disease, periodontitis, human abdominal aortic aneurysm, and obesity. To date, m⁶A modification has been identified as an important biological function in enrichment and regulation of lncRNAs. In this review, we summarize the role of m⁶A modification in the regulation and function of tumor-related lncRNAs. Moreover, we discuss the potential applications and possible future directions in the field.

Genomic analysis of pancreatic cancer reveals 3 molecular subtypes with different clinical outcomes

ABSTRACT

Pancreatic cancer has a very high mortality with a 5-year survival of <5%. The purpose of this study was to classify specific molecular subtypes associated with prognosis of pancreatic cancer using The Cancer Genome Atlas (TCGA) multiplatform genomic data.Multiplatform genomic data (N = 178), including gene expression, copy number alteration, and somatic mutation data, were obtained from cancer browser (https://genome-cancer.ucsc.edu, cohort: TCGA Pancreatic Cancer). Clinical data including survival results were analyzed. We also used validation cohort (GSE50827) to confirm the robustness of these molecular subtypes in pancreatic cancer.When we performed unsupervised clustering using TCGA gene expression data, we found three distinct molecular subtypes associated with different survival results. Copy number alteration and somatic mutation data showed different genomic patterns for these three subtypes. Ingenuity pathway analysis revealed that each subtype showed differentially altered pathways. Using each subtype-specific genes (200 were selected), we could predict molecular subtype in another cohort, confirming the robustness of these molecular subtypes of pancreatic cancer. Cox regression analysis revealed that molecular subtype is the only significant prognostic factor for pancreatic cancer (P = .042, 95% confidence interval 0.523-0.98).Genomic analysis of pancreatic cancer revealed 3 distinct molecular subtypes associated with different survival results. Using these subtype-specific genes and prediction model, we could predict molecular subtype associated with prognosis of pancreatic cancer.

Analysis of m6A RNA Methylation-Related Genes in Liver Hepatocellular Carcinoma and Their Correlation with Survival

N6-methyladenosine (m6A) modification on RNA plays an important role in tumorigenesis and metastasis, which could change gene expression and even function at multiple levels such as RNA splicing, stability, translocation, and translation. In this study, we aim to conduct a comprehensive analysis on m6A RNA methylation-related genes, including m6A RNA methylation regulators and m6A RNA methylation-modified genes, in liver hepatocellular carcinoma, and their relationship with survival and clinical features. Data, which consist of the expression of widely reported m6A RNA methylation-related genes in liver hepatocellular carcinoma from The Cancer Genome Atlas (TCGA), were analyzed by one-way ANOVA, Univariate Cox regression, a protein–protein interaction network, gene enrichment analysis, feature screening, a risk prognostic model, correlation analysis, and consensus clustering analysis. In total, 405 of the m6A RNA methylation-related genes were found based on one-way ANOVA. Among them, DNA topoisomerase 2-alpha (TOP2A), exodeoxyribonuclease 1 (EXO1), ser-ine/threonine-protein kinase Nek2 (NEK2), baculoviral IAP repeat-containing protein 5 (BIRC5), hyaluronan mediated motility receptor (HMMR), structural maintenance of chromosomes protein 4 (SMC4), bloom syndrome protein (BLM), ca-sein kinase I isoform epsilon (CSNK1E), cytoskeleton-associated protein 5 (CKAP5), and inner centromere protein (INCENP), which were m6A RNA methylation-modified genes, were recognized as the hub genes based on the protein–protein interaction analysis. The risk prognostic model showed that gender, AJCC stage, grade, T, and N were significantly different between the subgroup with the high and low risk groups. The AUC, the evaluation parameter of the prediction model which was built by RandomForest, was 0.7. Furthermore, two subgroups were divided by consensus clustering analysis, in which stage, grade, and T differed. We identified the important genes expressed significantly among two clusters, including uridine-cytidine kinase 2 (UCK2), filensin (BFSP1), tubulin-specific chaperone D (TBCD), histone-lysine N-methyltransferase PRDM16 (PRDM16), phosphorylase b ki-nase regulatory subunit alpha (PHKA2), serine/threonine-protein kinase BRSK2 (BRSK2), Arf-GAP with coiled-coil (ACAP3), general transcription factor 3C polypep-tide 2 (GTF3C2), and guanine nucleotide exchange factor MSS4 (RABIF). In our study, the m6A RNA methylation-related genes in liver hepatocellular carcinoma were analyzed systematically, including the expression, interaction, function, and prognostic values, which provided an important theoretical basis for m6A RNA methylation in liver cancer. The nine important m6A-related genes could be prognostic markers in the survival time of patients.