- Research article
- Open Access
m1A methylation modification patterns and metabolic characteristics in hepatocellular carcinoma
BMC Gastroenterology volume 22, Article number: 93 (2022)
The dysregulation of RNA methylation has been demonstrated to contribute to tumorigenicity and progression in recent years. However, the alteration of N1-methyladenosine (m1A) methylation and its role in hepatocellular carcinoma (HCC) remain unclear.
We systematically investigated the modification patterns of 10 m1A regulators in HCC samples and evaluated the metabolic characteristics of each pattern. A scoring system named the m1Ascore was developed using principal component analysis. The clinical value of the m1Ascore in risk stratification and drug screening was further explored.
Three m1A modification patterns with distinct metabolic characteristics were identified, corresponding to the metabolism-high, metabolism-intermediate and metabolism-excluded phenotypes. Patients were divided into high- or low-m1Ascore groups, and a significant survival difference was observed. External validation confirmed the prognostic value of the m1Ascore. A nomogram incorporating the m1Ascore and other clinicopathological factors was constructed and had good performance for predicting survival. Two agents, mitoxantrone and doxorubicin, were determined to be potential therapeutic drugs for the high-risk group.
This study provided novel insights into m1A modification and metabolic heterogeneity in cancer, promoted risk stratification in the clinic from the perspective of m1A modification, and further guided individual treatment strategies.
N1-methyladenosine (m1A), an RNA methylation pattern, is a dynamic RNA modification regulated by methyltransferases, demethylases and binding proteins . The formation of m1A, which is catalyzed by methyltransferases, including TRMT10C, TRMT61B, TRMT6 and TRMT61A, is reversible by the ALKBH1 and ALKBH3 demethylases. YTHDF1, YTHDF2, YTHDF3 and YTHDC1 act as binding proteins that specifically recognize m1A sites and induce downstream effects . The pioneering works by Dan et al.  and Li et al.  provided transcriptome-wide mapping of m1A and uncovered that m1A modification is enriched in the vicinity of the start codon, highlighting the roles of m1A in mRNA stability and translation. Today, it is recognized that the dynamic regulation of m1A in response to physiological stress and the dysregulated expression of m1A regulators are correlated with tumorigenesis and cancer recurrence [5, 6].
Metabolic reprogramming, as a hallmark of cancer, exists in nearly all cancer cells . To support uncontrolled cell growth and proliferation, the metabolic manners of cancer cells are distinct from those of normal cells . A classic metabolic reprogramming method is aerobic glycolysis, and cancer cells secrete lactate despite the aerobic environment present . In addition, with a diversity of nutrients, including lipids, amid acids and nucleotides, found to take part in cancer metabolic reprogramming in addition to glucose , a consensus emerged based on the results from molecular-level clinical trials that showed different metabolic dependencies across cancer types , such as fatty acid metabolism in prostate cancer  and serine synthesis in breast cancer . This difference, which is also recognized as metabolic heterogeneity, is mainly associated with the heterogeneity of stomatic alterations among cancer cells . Drivers of oncogenesis vary in genomic alterations and subsequently affect different downstream pathways, finally reflecting heterogeneity in cancer metabolism phenotypes . Studies reported by Haider et al.  and Sander et al.  provided sufficient evidence in the molecular dimension, which showed widespread transcriptional differences in metabolic-related genes between normal tissue and tumor tissue. Moreover, the identification of cancer metabolism phenotypes contributed to the reasonable stratification of patients in clinical practice with regard to the application of metabolic inhibitors and the individual selection of target therapies .
Although the effective classification of cancer metabolic characteristics remains a controversial issue to date, several pioneering works that focus on the expression pattern of metabolic genes have been performed in this frontier field and have identified distinct metabolic phenotypes with different sensitivities to metabolic inhibitors [19, 20]. Here, we explored m1A modification and its relationship with metabolic characteristics in hepatocellular carcinoma (HCC). Three m1A modification patterns with distinct metabolic characteristics were identified, corresponding to the metabolism-excluded, metabolism-high and metabolism-intermediate phenotypes, respectively. Clinical characteristics and prognosis showed significant differences among these phenotypes. Further, an m1A-related scoring system was developed and showed reliable performance in distinguishing drug sensitivities and stratifying risk in the clinic.
Two RNA sequencing (RNA-seq) datasets were used in the current study: The Cancer Genome Atlas (TCGA, liver hepatocellular carcinoma (LIHC), https://www.cancer.gov) and International Cancer Genome Consortium (ICGC) (LIRI-JP, www.icgc.org). HCC patients with RNA-seq data (fragments per kilobase per million reads, FPKM) and corresponding clinical information were enrolled in March 2021.
R 4.0.1 was applied to carry out all statistical tests. All tests were two-sided, and a P value < 0.05 was considered statistically significant.
Additional detailed protocols are provided in the Supplementary Methods.
Landscape of genetic variations in m1A regulators in HCC
There was a prevalent alteration of the copy number variation (CNV) in m1A regulators (Fig. 1a). Four regulators, YTHDF2, TRMT61A, ALKBH1 and YTHDC1, showed a majority of loss alterations, and three regulators, YTHDF3, YTHDF1 and TRMT61B, mainly had gain alterations. The locations of the CNV alterations of each m1A regulator on chromosomes are shown in Fig. 1b. Somatic mutations were also investigated in the TCGA cohort, and only a few samples showed mutations in m1A regulators, with a frequency of 3.33% (Fig. 1c). The mRNA expression levels of regulators were compared between normal and HCC samples, and a remarkable upregulation of each regulator was observed in HCC tissues (Fig. 1d, Additional file 1: Fig. S1a). Similar results were observed when comparatively estimating the expression of regulators in paired HCC and adjacent normal tissues (Additional file 1: Fig. S1b). The above results revealed significant heterogeneity at the genetic and epigenetic levels of m1A regulators between normal and HCC tissues.
Identification of m1A methylation modification patterns
A network plot was used to visualize the connections and prognostic value of the regulators (Fig. 2a). A significantly positive correlation in expression was observed, and this correlation was prevalent despite the functional category of the m1A regulators. We found that TRMT10C, TRMT6, ALKBH1, YTHDF2, and YTHDF1 were positively correlated with the expression of other m1A regulators, and all of them were determined to be risk factors in univariate Cox analysis. In addition, none of the 10 m1A regulators was identified as a protective factor, revealing a similar prognostic effect among them.
Unsupervised cluster analysis was carried out using the expression data of m1A regulators, and three distinct modification patterns, which were named subtypes 1–3, were determined (Additional file 1: Fig. S2a-c). A total of 128 patients had subtype 1, 82 patients had subtype 2, and 120 patients had subtype 3. Subtype 2 had the highest m1A regulator expression levels compared with the other subtypes, while subtype 3 had the lowest expression levels (Fig. 2b). In accordance with the expression levels of m1A regulators in each subtype, a significant difference in the enrichment scores of mRNA methylation levels was observed among the three subtypes (Fig. 2c). Survival analysis showed that subtype 2 had a remarkably poorer prognosis, while subtype 1 and subtype 3 had better prognoses (Fig. 2d). These results revealed the modification heterogeneity of m1A regulators and distinct prognosis of HCC patients.
Metabolic characteristics of the three m1A modification subtypes
First, gene set variation analysis (GSVA) was performed in terms of the HALLMARK gene sets. As shown in Fig. 3a, b, the hallmarks related to metabolism, including fatty acid metabolism, xenobiotic metabolism and bile acid metabolism, were markedly downregulated in subtype 2. These results revealed low metabolic activity characterizing subtype 2, suggesting metabolic heterogeneity among the three subtypes. Thus, we further investigated the metabolic characteristics of each subtype. By estimating the GSVA enrichment scores of 49 metabolism-associated pathways, subtype 2 and subtype 1 had the lowest and highest enrichment scores, respectively, across the four major metabolic categories, respectively, while subtype 3 had modest values, strongly indicating that different m1A modification patterns were characterized by distinct metabolic activities (Fig. 3c–f, Additional file 1: Fig. S2d). Hence, subtype 1 was classified as a metabolism-high phenotype, corresponding to the hyperactive state in metabolism; subtype 2 was classified as a metabolism-excluded phenotype, corresponding to the hypoactive state; and subtype 3 was classified as a metabolism-intermediate phenotype.
We next explored the correlation between each m1A regulator and each metabolic pathway. Tumors with high activation of purine metabolism, pyrimidine metabolism, selenoamino acid metabolism, glucose metabolism, other glycan degradation and glycosaminoglycan metabolism showed a positive correlation with all regulators, while those with high activation of arginine and proline metabolism and linoleic acid metabolism were negatively correlated with all regulators. These results revealed that whether tumors with high expression of an m1A regulator showed high activation of a metabolic pathway was actually determined by the specific m1A regulator and metabolic pathway (Fig. 3g).
Development of an m1A scoring system
270 overlapping m1A subtype-related genes were extracted from the different expression genes (DEGs) (Additional file 1: Fig. S3a). Unsupervised cluster analysis was performed again based on these 270 DEGs to classify patients into three genomic patterns (m1A gene clusters 1–3) (Additional file 1: Fig. S3b-c).
The expression level of m1A regulators and mRNA methylation in the three gene clusters were estimated, and both showed significant differences (Fig. 4a, b). Survival analysis revealed that patients in gene cluster 1 had a better prognosis than patients in the other gene clusters (Fig. 4c). The metabolic characteristics of the three gene clusters were further explored, and a significant difference was observed. Gene clusters 1 and 3 had the highest and lowest enrichment scores, respectively, across the four major metabolic categories, which suggested that clusters 1 and 3 could be identified as the metabolic-high phenotype and metabolism-excluded phenotype, respectively, and cluster 2 could be identified as the metabolic-intermediate phenotype. (Fig. 4d).
We further established a scoring system termed the m1Ascore. The correlation between the known classifications in this study and the m1Ascore is presented as an alluvial diagram (Fig. 5a). Patients with subtype 2 had the highest median m1Ascore compared with patients with other subtypes (Fig. 5b), indicating that a high m1Ascore could be associated with the hypoactive metabolic state. Gene cluster 3 had the highest median m1Ascore, while cluster 1 had the lowest median score (Fig. 5c). A significant correlation between the m1Ascore and metabolic pathways was directly observed (Fig. 3g).
Risk stratification of HCC patients
We classified 330 patients in the TCGA cohort into high- or low-m1Ascore groups, and Kaplan–Meier (KM) curves showed that a high m1Ascore was associated with a short survival time (Fig. 6a). These results suggested the clinical value of the m1Ascore in risk stratification; thus, we named the high- and low-m1Ascore groups the high-risk group and low-risk group, respectively.
The clinicopathological information and molecular characteristics of the high-risk and low-risk groups were further assessed. Patients in the high-risk group were significantly associated with advanced stage, poor differentiation and higher Child–Pugh scores (Fig. 6b). The expression level of m1A regulators in the high-risk group was significantly higher than that in the low-risk group (Additional file 1: Fig. S4a). This result is consistent with the difference in mRNA methylation levels between the two groups (Additional file 1: Fig. S4b). A positive correlation between mRNA methylation and the m1Ascore was also observed (Additional file 1: Fig. S4c). According to GSVA, hallmarks related to cell signaling and proliferation were significantly activated in the high-risk group compared with the low-risk group, including G2M checkpoint, TGF-β signaling, Wnt-β catenin signaling and PI3K-AKT-mTOR signaling (Fig. 6c). The metabolic characteristics were investigated, and the results indicated that the low-risk group was associated with a hyperactive metabolic state, while the high-risk group was associated with a hypoactive metabolic state (Additional file 1: Fig. S4d-g). The somatic mutation landscapes are shown in Fig. 6d. TP53 showed the highest mutation frequency in the high-risk group (37%), while CTNNB1 and TTN both had the highest mutation frequencies in the low-risk group (24%). No significant differences were detected between the two groups in terms of tumor mutation burden, microsatellite instability or RNAss (Additional file 1: Fig. S4h-j).
Identification of potential therapeutic agents
We successfully identified HCC patients with poor prognosis according to the m1Ascore; thus, this score would be valuable in the clinic to select suitable drugs for the high-risk group in the context of personalized treatment. For this purpose, 1770 compounds were retrieved for drug screening in the CTRP and PRISM databases after excluding 160 duplicated compounds. We next eliminated compounds lacking information in more than 20% of the samples. Finally, 1645 compounds were included in the subsequent analysis.
The clinical samples were classified by deciles according to the m1Ascore, and the estimated area under curve (AUC) value of each compound in each decile was obtained. Compounds with lower estimated AUC values in the top decile group were identified when compared with the bottom decile (log2FC > 0.2). Moreover, compounds that showed a negative correlation between the m1Ascore and estimated AUC values by Spearman correlation analysis were recorded (Spearman’s r < -0.40 for PRISM and CTRP). We retained 22 candidate agents after screening, including 6 agents from CTRP and 16 agents from PRISM (Fig. 6e–f). These compounds were all characterized by lower estimated AUC values and a negative correlation with the m1Ascore, representing better drug sensitivity in the high-risk HCC group. However, such features are not reliable to support the therapeutic effect of these compounds in the clinic. Hence, we subsequently performed CMap analysis, and two compounds with CMap scores < -95 were finally identified: mitoxantrone and doxorubicin. In addition, the immunotherapeutic response was also evaluated; however, no significant response difference was detected between the high-risk and low-risk groups (Additional file 1: Fig. S5a-b).
External validation and construction of prognostic model
We subsequently performed multivariate Cox regression analysis, and the results demonstrated that a high m1Ascore was independently associated with poor survival (Fig. 7a). In the external validation LIRI-JP cohort, the expression levels of m1A regulators between normal and tumor tissues were analyzed, and the results confirmed the overexpression of m1A regulators in HCC patients (Additional file 1: Fig. S6a). The survival time of high-m1Ascore patients was significantly shorter than that of low-m1Ascore patients (Additional file 1: Fig. S6b). Multivariate analysis demonstrated that a high m1Ascore was an independent risk factor for survival (Additional file 1: Fig. S6c). By incorporating the other independent clinicopathological factors in the TCGA cohort, we built a nomogram for use in clinical practice (Fig. 7b). No remarkable deviation of predicted survival from actual survival was observed according to the calibration plots (Fig. 7c). tROC analysis showed that the accuracy of the nomogram in predicting survival was more satisfactory at 1, 2, and 3 years than AJCC stage or the m1Ascore alone (Fig. 7d). The AUC values of the nomogram for predicting 1-, 2-, and 3-year overall survival (OS) were 0.75, 0.69, and 0.73, respectively. In the internal validation, the adjusted C-index of the nomogram was 0.698. The decision curve analysis (DCA) results indicated that compared with the m1Ascore or AJCC stage, the utilization of the nomogram in the clinic had greater net benefits and net reduction (Fig. 7e–f). These results suggested that the nomogram performed well.
We further chosen data on 18 cancer types that had more than 5 normal samples from the TCGA database to comprehensively investigate the expression patterns of the m1A regulators . The expression distributions of the m1A regulators across all cancers are shown in Additional file 1: Fig. S7a, and the expression landscapes are shown in Fig. 8a. The expression level of each regulator was significantly upregulated in cholangiocarcinoma (CHOL) compared with the control, while downregulation was observed in kidney chromophobe (KICH) and thyroid carcinoma (THCA). In kidney renal clear cell carcinoma (KIRC), compared with normal tissues, we noticed that the expression levels of TRMT10C and TRMT6 were significantly higher, while those of ALKBH1, ALKBH3 and YTHDF1 were significantly lower (Additional file 1: Fig. S7b). The expression correlations among m1A regulators across cancers are shown in Fig. 8b, in which ALKBH1 and YTHDC1 had the highest association (r = 0.77, P < 0.01). These results indicated a prevalent m1A regulator expression difference between normal and tumor tissues, and the expression level of each regulator showed remarkable intratumor heterogeneity. Univariate Cox regression analysis was performed on 33 cancer types to explore the impact of m1A regulators on survival, and the results showed that the prognostic value of each regulator varied across cancer types (Fig. 8c). For example, overexpression of YTHDF2 and ALKBH1 is associated with poor survival in KICH and HCC but with better survival in KIRC. Finally, to better understanding the influence of m1A regulators in cancer biology, we provided a schematic to summarize the current knowledge of the roles of m1A regulators in metabolic reprogramming of hepatocellular carcinoma and pointed the potential therapeutic targets (Fig. 9).
m1A, which has been reported to play a critical role in a variety of functional activities, has drawn increasing interest in recent years. In RNA metabolism, it was found that the recruitment of the m1A reader YTHDF2 induces the degradation of m1A-modified RNAs . In translation processes, it was reported that m1A modification was enriched in the region near the start codon, and these sites tended to form secondary structures with higher translation efficiency . However, the impact of m1A in cancer biology has been poorly elucidated to date.
Here, we first reported three distinct m1A modification patterns in HCC patients and found that each pattern had distinct metabolic characteristics. Subtype 1 displayed the highest metabolic activity across the four major metabolic categories, characterized by a metabolism-high phenotype. Subtype 2 displayed the lowest metabolic activity, characterized by a metabolism-excluded phenotype. Subtype 3 displayed moderate metabolic activity, characterized by a metabolism-intermediate phenotype. Survival analysis showed that patients with subtype 2 had the worst prognosis compared with patients with the other subtypes. Yang C et al. reported three HCC subclasses (subclasses 1–3) with distinct metabolic characteristics based on the expression profiles of 2752 metabolism-related genes . Subclass 3, which displayed a hypoactive metabolic state, had a poorer prognosis than subclass 1. Shen et al. reported three clusters in HCC patients with distinct metabolic characteristics using the expression data of m6A regulators . Cluster 3 showed the worst OS and the lowest metabolic activity. The findings of these studies are consistent with our results. We also observed that among seven major processes of cancer, hallmarks related to metabolism, including fatty acid metabolism, xenobiotic metabolism and bile acid metabolism, were significantly downregulated in subtype 2 compared with the other subtypes. Combined with the metabolic characteristics of each subtype, these results demonstrated the hypoactive metabolism of subtype 2, favoring the credibility of our classification of metabolic phenotypes for different m1A modification patterns. We further elucidated the association between m1A modification patterns and metabolic characteristics by analyzing the m1A subtype-related DEGs. The three genomic patterns were also found to be related to distinct metabolic states and prognoses, which again confirmed the metabolic heterogeneity of each m1A modification pattern.
We next developed a scoring system, the m1Ascore, and found that subtype 2 and gene cluster 3 had the highest median scores among the three subtypes and gene clusters. As expected, patients belonging to the high-m1Ascore group had highly expressed m1A regulators with a hypoactive metabolic state. These results suggested that the m1Ascore could be used as a tool to estimate the m1A modification pattern in individual patients, further reflecting the metabolic state. Further analysis confirmed that the m1Ascore was an independent risk factor for HCC patient prognosis, and a nomogram incorporating the m1Ascore and other clinicopathological factors was successfully constructed and showed good discriminatory capacity and accuracy. These results showed that m1A regulators contribute to the prognosis of patients with HCC. The roles of m1A regulators in HCC also had been investigated in previous studies. The expression levels of readers, including YTHDF1, YTHDF2, and YTHDF3, are reported to upregulate in HCC tissues compared with adjacent tissues, and a positive correlation was detected between immune cell infiltration and the expression levels of YTHDF1, YTHDF2, and YTHDF3, revealing the functional roles of m1A regulators in modulating the infiltration of immune cells in HCC tissues . As for erasers, consistent with our results, Wang et al. reported an elevated ALKBH3 expression in HCC compared with adjuvant non-tumorous tissues, and found that patients with high ALKBH3 expression level displayed poor prognosis and the knockdown of ALKBH3 inhibits tumor cells proliferation, indicating the functional role of m1A modification in promoting cell cycle . A recent study reported that the expression levels of writers TRMT6 and TRMT61A in HCC tissues were significantly increased and showed that the increased TRMT6 and TRMT61A expression levels are negatively associated with patients’ prognosis, and found that TRMT6/TRMT61A-mediated m1A methylation can initiate liver tumorigenesis by regulating lipid metabolism . The modification effect of ALKBH3 on m1A methylation is opposite to that of TRMT6 and TRMT61A, however, as revealed by the above studies and our results, a poor prognosis was consistently observed in HCC patients with an elevated TRMT6/ TRMT61A or ALKBH3 expression level, suggesting the crucial role of m1A methylation dysfunction in HCC. It has been reported that m1A modulates the PI3K/AKT/mTOR pathway in gastrointestinal cancers , and the dysregulation of this pathway has been identified as a key factor contributing to poor survival in several cancers [29,30,31]. Thus, the mechanism underlying the impact of m1A on survival is probably associated with the dysregulation of the PI3K/AKT/mTOR pathway. Moreover, the PI3K/AKT/mTOR pathway is known as an important process in regulating glucose intake  and is deemed a crosstalk center between epigenetics and metabolic heterogeneity for the ability to drive the Warburg effect . Further research on the PI3K/AKT/mTOR pathway is needed to enhance our understanding of the metabolic heterogeneity between each m1A modification pattern. In addition, consistent with the results reported by Shi et al. that CNV events commonly occurred in m1A regulators in HCC , our study found a prevalence of CNVs in all m1A regulators, indicating that the alternations to m1A regulators in HCC are prevalent. It had been reported that recessive mutations in TRMT10C affect aerobic respiration processing and result in lactic acidosis . However, very few studies have focused on the effects of the genetic alternations of m1A regulators on HCC biology to date. TP53 mutation, which has been demonstrated as an important tumorigenesis procedure  and to enhance aerobic glycolysis of tumor cells , was found to positively correlate with the genetic changes of m1A regulators in HCC, and mutation frequency of TP53 in low- and high-m1Ascore groups were also found significantly different in our study (22% vs. 37%, P = 0.003). These results promoted us that the genetic variants of m1A regulators might cooperate with the mutations of cancer-causing genes in the oncogenesis and metabolic reprogramming of HCC.
Identifying promising treatment strategies based on molecular characteristics to maximize the therapeutic effect is also one of the main purposes of this study. We subsequently identified two potential agents for the treatment of high-risk HCC patients, mitoxantrone and doxorubicin. Mitoxantrone is an antineoplastic drug that is widely used in treating acute myeloid leukemia (AML)  and breast cancer . The efficacy of mitoxantrone for HCC patients has already been investigated [40, 41], and the results suggested that patients with smaller tumor masses and good liver reserves may benefit from a therapeutic regimen consisting of mitoxantrone, 5-fluorouracil and cisplatin. Our study provided new insights to identify patients suitable for treatment with mitoxantrone. Doxorubicin is one of the first-line chemotherapeutic agents used in the management of hematological tumors  and has also been approved for the treatment of HCC patients . Cardiotoxicity is a common side effect in the application of doxorubicin . It has been reported that 48% of patients using 700 mg/m2 doxorubicin will develop heart failure . Our study indicated that the clinical use of doxorubicin in patients with a high m1Ascore might be more reasonable compared with that in patients with a low m1Ascore. This result is helpful for reducing the risk of cardiotoxicity in patients and avoiding ineffective overtreatment.
There are some limitations in the current study. First, we estimated the metabolic characteristics of each m1A modification pattern using the GSVA score of metabolic pathways rather than directly analyzing the metabolite profiling or expression patterns of metabolic genes. Some metabolic pathways related to fatty acid metabolism, which shared common metabolic-related genes, might have similar scores. Thus, we have to acknowledge that there might be discrepancies to a certain extent between the estimated metabolic state and actual metabolic state of each pattern. Second, the impacts of each m1A regulator on key metabolic enzymes were not discussed here, and the molecular mechanism underlying the correlation between m1A modification and metabolic heterogeneity remains unknown. Finally, the results of the current study are all based on bioinformatics analyses. The lack of experimental verification might undermine the persuasiveness of our conclusion; hence, further in vivo studies are needed to promote the application of our findings in clinical practice.
Our findings highlighted the role of m1A methylation modification in the crosstalk between epigenetics and metabolic heterogeneity in cancer. The estimation of m1A modification in HCC patients will promote our understanding of metabolic characteristics and be beneficial for survival stratification, further guiding personalized clinical treatment.
Hepatocellular carcinoma (HCC);
The Cancer Genome Atlas
International Cancer Genome Consortium
Copy number variation
Gene set variation analysis
Different expression genes
Area under curve
Decision curve analysis
Xiong X, Li X, Yi C. N1-methyladenosine methylome in messenger RNA and non-coding RNA. Curr Opin Chem Biol. 2018;45:179–86.
Shi H, Wei J, He C. Where, when, and how: context-dependent functions of RNA methylation writers, readers, and erasers. Mol Cell. 2019;74:640–50.
Dominissini D, Nachtergaele S, Moshitch-Moshkovitz S, Peer E, Kol N, Ben-Haim MS, et al. The dynamic N(1)-methyladenosine methylome in eukaryotic messenger RNA. Nature. 2016;530:441–6.
Li X, Xiong X, Wang K, Wang L, Shu X, Ma S, et al. Transcriptome-wide mapping reveals reversible and dynamic N(1)-methyladenosine methylome. Nat Chem Biol. 2016;12:311–6.
Zhang C, Jia G. Reversible RNA modification N1-methyladenosine (m1A) in mRNA and tRNA. Genomics Proteomics Bioinformatics. 2018;16:155–61.
Han X, Wang M, Zhao Y-L, Yang Y, Yang Y-G. RNA methylations in human cancers. Semin Cancer Biol. 2020.
Faubert B, Solmonson A, DeBerardinis RJ. Metabolic reprogramming and cancer progression. Science. 2020;368.
Vernieri C, Casola S, Foiani M, Pietrantonio F, de Braud F, Longo V. Targeting cancer metabolism: dietary and pharmacologic interventions. Cancer Discov. 2016;6:1315–33.
Vander Heiden MG, Cantley LC, Thompson CB. Understanding the Warburg effect: the metabolic requirements of cell proliferation. Science. 2009;324:1029–33.
DeBerardinis RJ, Chandel NS. Fundamentals of cancer metabolism. Sci Adv. 2016;2:e1600200.
Vander Heiden MG, DeBerardinis RJ. Understanding the intersections between metabolism and cancer biology. Cell. 2017;168:657–69.
Berquin IM, Edwards IJ, Kridel SJ, Chen YQ. Polyunsaturated fatty acid metabolism in prostate cancer. Cancer Metastasis Rev. 2011;30:295–309.
Possemato R, Marks KM, Shaul YD, Pacold ME, Kim D, Birsoy K, et al. Functional genomics reveal that the serine synthesis pathway is essential in breast cancer. Nature. 2011;476:346–50.
Kim J, DeBerardinis RJ. Mechanisms and implications of metabolic heterogeneity in cancer. Cell Metab. 2019;30:434–46.
Peng X, Chen Z, Farshidfar F, Xu X, Lorenzi PL, Wang Y, et al. Molecular characterization and clinical relevance of metabolic expression subtypes in human cancers. Cell Rep. 2018;23:255-269.e4.
Haider S, McIntyre A, van Stiphout RGPM, Winchester LM, Wigfield S, Harris AL, et al. Genomic alterations underlie a pan-cancer metabolic shift associated with tumour hypoxia. Genome Biol. 2016;17:140.
Reznik E, Sander C. Extensive decoupling of metabolic genes in cancer. PLoS Comput Biol. 2015;11:e1004176.
Mukherjee S, Mukherjee N, Gayen P, Roy P, Babu SPS. Metabolic inhibitors as antiparasitic drugs: pharmacological, biochemical and molecular perspectives. Curr Drug Metab. 2016;17:937–70.
Gong Y, Ji P, Yang Y-S, Xie S, Yu T-J, Xiao Y, et al. Metabolic-pathway-based subtyping of triple-negative breast cancer reveals potential therapeutic targets. Cell Metab. 2021;33:51-64.e9.
Daemen A, Peterson D, Sahu N, McCord R, Du X, Liu B, et al. Metabolite profiling stratifies pancreatic ductal adenocarcinomas into subtypes with distinct sensitivities to metabolic inhibitors. Proc Natl Acad Sci U S A. 2015;112:E4410-4417.
Zhao Y, Xu L, Wang X, Niu S, Chen H, Li C. A novel prognostic mRNA/miRNA signature for esophageal cancer and its immune landscape in cancer progression. Mol Oncol. 2021;15:1088–109.
Seo KW, Kleiner RE. YTHDF2 recognition of N1-Methyladenosine (m1A)-Modified RNA Is associated with transcript destabilization. ACS Chem Biol. 2020;15:132–9.
Yang C, Huang X, Liu Z, Qin W, Wang C. Metabolism-associated molecular classification of hepatocellular carcinoma. Mol Oncol. 2020;14:896–913.
Shen X, Hu B, Xu J, Qin W, Fu Y, Wang S, et al. The m6A methylation landscape stratifies hepatocellular carcinoma into 3 subtypes with distinct metabolic characteristics. Cancer Biol Med. 2020;17:937–52.
Liu M, Zhao Z, Cai Y, Bi P, Liang Q, Yan Y, et al. YTH domain family: potential prognostic targets and immune-associated biomarkers in hepatocellular carcinoma. Aging. 2021;13:24205–18.
Wang Q, Wang G, Wang Y, Liu C, He X. Association of AlkB homolog 3 expression with tumor recurrence and unfavorable prognosis in hepatocellular carcinoma. J Gastroenterol Hepatol. 2018. https://doi.org/10.1111/jgh.14117.
Wang Y, Wang J, Li X, Xiong X, Wang J, Zhou Z, et al. N1-methyladenosine methylation in tRNA drives liver tumourigenesis by regulating cholesterol metabolism. Nat Commun. 2021;12:6314.
Zhao Y, Zhao Q, Kaboli PJ, Shen J, Li M, Wu X, et al. m1A Regulated genes modulate PI3K/AKT/mTOR and ErbB Pathways in gastrointestinal cancer. Transl Oncol. 2019;12:1323–33.
Cai J, Sun H, Chen L, Xie M, Zhuang J, Gao L, et al. NAT1 is a critical prognostic biomarker and inhibits proliferation of colorectal cancer through modulation of PI3K/Akt/mTOR. Future Oncol Lond Engl. 2021.
Cai M, Li H, Chen R, Zhou X. MRPL13 promotes tumor cell proliferation, migration and EMT process in breast cancer through the PI3K-AKT-mTOR pathway. Cancer Manag Res. 2021;13:2009–24.
W W, X G, H D. α2A-Adrenergic Receptor Inhibits the Progression of Cervical Cancer Through Blocking PI3K/AKT/mTOR Pathway. OncoTargets Ther. 2020; doi:https://doi.org/10.2147/OTT.S264409.
Elstrom RL, Bauer DE, Buzzai M, Karnauskas R, Harris MH, Plas DR, et al. Akt stimulates aerobic glycolysis in cancer cells. Cancer Res. 2004;64:3892–9.
Thakur C, Chen F. Connections between metabolism and epigenetics in cancers. Semin Cancer Biol. 2019;57:52–8.
Shi Q, Xue C, Yuan X, He Y, Yu Z. Gene signatures and prognostic values of m1A-related regulatory genes in hepatocellular carcinoma. Sci Rep. 2020;10:15083.
Metodiev MD, Thompson K, Alston CL, Morris AAM, He L, Assouline Z, et al. Recessive mutations in TRMT10C cause defects in mitochondrial RNA processing and multiple respiratory chain deficiencies. Am J Hum Genet. 2016;98:993–1000.
Timmerman DM, Remmers TL, Hillenius S, Looijenga LHJ. Mechanisms of TP53 pathway inactivation in embryonic and somatic cells-relevance for understanding (Germ Cell) tumorigenesis. Int J Mol Sci. 2021;22:5377.
Eriksson M, Ambroise G, Ouchida AT, Lima Queiroz A, Smith D, Gimenez-Cassina A, et al. Effect of mutant p53 proteins on glycolysis and mitochondrial metabolism. Mol Cell Biol. 2017;37:e00328-e417.
Deng L, Zhang C, Ying S, Cai B, Zhou F. Effect of dose ratio on mitoxantrone and daunorubicin in acute myeloid leukemia: a systematic review and meta-analysis of randomized controlled trials. Clin Lymphoma Myeloma Leukemia. 2021;21(1):e10–20.
Schmid P, Flath B, Akrivakis K, Heilmann V, Mergenthaler H-G, Sezer O, et al. Gemcitabine and mitoxantrone in metastatic breast cancer: a phase-I-study. Invest New Drugs. 2005;23:349–56.
Yang T-S, Chang H-K, Chen J-S, Lin Y-C, Liau C-T, Chang W-C. Chemotherapy using 5-fluorouracil, mitoxantrone, and cisplatin for patients with advanced hepatocellular carcinoma: an analysis of 63 cases. J Gastroenterol. 2004;39:362–9.
Ikeda M, Okusaka T, Sato Y, Furuse J, Mitsunaga S, Ueno H, et al. A Phase I/II trial of continuous hepatic intra-arterial infusion of 5-fluorouracil, mitoxantrone and cisplatin for advanced hepatocellular carcinoma. Jpn J Clin Oncol. 2017;47:512–9.
Dimitrakis P, Romay-Ogando M-I, Timolati F, Suter TM, Zuppinger C. Effects of doxorubicin cancer therapy on autophagy and the ubiquitin-proteasome system in long-term cultured adult rat cardiomyocytes. Cell Tissue Res. 2012;350:361–72.
Llovet JM. Updated treatment approach to hepatocellular carcinoma. J Gastroenterol. 2005;40:225–35.
Shabalala S, Muller CJF, Louw J, Johnson R. Polyphenols, autophagy and doxorubicin-induced cardiotoxicity. Life Sci. 2017;180:160–70.
Swain SM, Whaley FS, Ewer MS. Congestive heart failure in patients treated with doxorubicin: a retrospective analysis of three trials. Cancer. 2003;97:2869–79.
We acknowledge the TCGA database, ICGC database and Servier Medical Art for free use.
There was no funding for this study regarding the design of the study as well as collection, analysis, interpretation of data, or writing the manuscript. In general, work in the group of WW and CYH was founded from the Wannan Medical College (WK2017F20) and Yijishan Hospital of Wannan Medical College (YR201608) to WW, and the Anhui Province Key Research and Development Program Project (201904a07020028) to CYH.
Ethics approval and consent to participate
Consent for publication
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Clinical Characteristics of the HCC patients in TCGA and ICGC; Figure S1. Expression levels of m1A regulators in HCC and normal tissues; Figure S2. Unsupervised clustering of m1A regulators and metabolic characteristics in each cluster; Figure S3. m1A subtype-related overlap genes and unsupervised clustering of the overlap genes; Figure S4. Characteristics in high and low m1Ascore groups; Figure S5. Difference in immunotherapy response between high- and low-m1Ascore groups; Figure S6. External validation of m1Ascore in ICGC-JP cohort; Figure S7. Expression levels of m1A regulators in TCGA cohort across pan-cancer analysis.
About this article
Cite this article
Tong, C., Wang, W. & He, C. m1A methylation modification patterns and metabolic characteristics in hepatocellular carcinoma. BMC Gastroenterol 22, 93 (2022). https://doi.org/10.1186/s12876-022-02160-w
- Hepatocellular carcinoma