Identification of cross-talk between m6A and 5mC regulators associated with onco-immunogenic features and prognosis across 33 cancer types
Journal of Hematology & Oncology volume 13, Article number: 22 (2020)
Methylation of RNA and DNA, notably in the forms of N6-methyladenosine (m6A) and 5-methylcytosine (5mC) respectively, plays crucial roles in diverse biological processes. Currently, there is a lack of knowledge regarding the cross-talk between m6A and 5mC regulators. Thus, we systematically performed a pan-cancer genomic analysis by depicting the molecular correlations between m6A and 5mC regulators across ~ 11,000 subjects representing 33 cancer types. For the first time, we identified cross-talk between m6A and 5mC methylation at the multiomic level. Then, we further established m6A/5mC epigenetic module eigengenes by combining hub m6A/5mC regulators and informed a comprehensive epigenetic state. The model reflected status of the tumor-immune-stromal microenvironment and was able to predict patient survival in the majority of cancer types. Our results lay a solid foundation for epigenetic regulation in human cancer and pave a new road for related therapeutic targets.
To the Editor,
Nucleotide methylation, notably in the forms of 5-methylcytosine (5mC) in DNA and N6-methyladenosine (m6A) in mRNA, carries important information for gene regulation . Recent research advances highlight the biological importance of m6A methylation as a dynamic and reversible post-transcriptional modification . 5mC DNA methylation, a conserved epigenetic modification along with m6A RNA modification, also plays critical roles in fundamental biological processes [3, 4]. In addition, recent studies have identified 5mC methylation as a modulator of alternative mRNA splicing at the post-transcriptional level [5, 6]. Although Zhou and colleagues  established a molecular link between 5mC DNA methylation and m6A mRNA methylation during fruit ripening, the potential cross-talk still remains uncharacterized in human cancers.
To address this issue, we curated a catalog of 20 and 21 genes that function mainly as regulators of RNA and DNA methylation, respectively (Fig. 1a). The genome-wide omics data comprising of 11,080 human samples across 33 cancer types from the The Cancer Genome Atlas (TCGA) were obtained for analyses (please see Methods and Table S1). First, most of the m6A and 5mC regulators were found to exhibit comparable expression levels across the 33 cancer types (Supplementary Fig. S1). Basing on the Gene Set Cancer Analysis (GSCA) web server , we further assessed the gene set differential expression profiles among 14 cancer types with available paired tumor-normal tissue expression data. Across multiple cancer types, the differentially expressed genes (upregulated or downregulated) included both m6A and 5mC regulators (Supplementary Fig. S2). Then, we investigated the mutation frequencies of the m6A and 5mC regulators. Intriguingly, m6A and 5mC regulators exhibited comparable levels of mutation frequency, and significant co-occurrences of genetic alterations were observed between the two regulators (Fig. 1b). Our results showed correlated expression patterns for genes within the same regulator class and even high correlations between the expression of m6A and 5mC regulators (Fig. 1b). Moreover, these m6A and 5mC regulators interacted with one another frequently in protein-protein interaction networks (Fig. 1d).
To identify the hub regulators involved in RNA and DNA methylation, we then applied weighted gene coexpression network analysis (WGCNA) to determine the hub genes in m6A and 5mC regulators (Fig. 2a). Strikingly, the number of hub m6A regulators was highly correlated with that of hub 5mC regulators in different cancer types (R = 0.84; Fig. 2b), which may be explained by the cross-talk. We then combined the hub m6A/5mC genes to develop an epigenetic module eigengene (EME), which may reflect both the pre- and post-transcriptional modification statuses. Next, we examined the correlation between EMEs and the activity of hallmark oncogenic pathways (Fig. 2c). Interestingly, our results indicate that high expression of the EME may reflect a highly proliferative and aggressive status in the majority of tumors. In addition, we applied GSCA  to analyze the effect (activation or inhibition) of m6A/5mC regulators on cancer-related pathways and confirmed that the m6A and 5mC regulators may be functionally related (Supplementary Fig. S3).
In addition to the tumor compartments, we further investigated the associations between the EME and immuno-stromal signatures representing different statuses of the immune and stromal cells (Table S2) across cancer types. In general, relatively low expression of inflammatory markers and low infiltration of immune cells were observed in the EMEhigh versus EMElow subgroups across cancer types (Fig. 2d). Interestingly, the high enrichment of stromal-related signatures was observed in the EMEhigh subgroups in almost all cancer types, indicating that hub m6A/5mC regulators may generally be involved in stroma activation (Fig. 2d).
Finally, we assessed the prognostic value of the EME in various types of cancers. We found that the EME showed oncogenic features in most cancer types, with overall survival (OS) hazard ratios larger than one (Supplementary Fig. S4a–c). Of these, high expression of the EME was significantly associated with unfavorable OS in cancer types such as KICH, ACC, and LGG (Supplementary Fig. S4a, b). Among HNSC, KIRC, and READ, improved survival was observed in the EMEhigh versus EMElow groups (Supplementary Fig. S4a, c).
In summary, to our best knowledge, this is the first study suggesting potential cross-talk between m6A and 5mC regulators in human cancers. This study provides essential insights into epigenetic regulation in cancer and paves new ways for related therapeutic targets.
Availability of data and materials
The datasets used in this study are publicly available. All other relevant data and R and other custom scripts are available upon request.
Cervical squamous cell carcinoma and endocervical adenocarcinoma
Copy number variation
Epigenetic module eigengene
Gene set enrichment analysis
Head and neck squamous carcinoma
Kidney renal clear cell carcinoma
Kidney renal papillary cell carcinoma
Brain low-grade glioma
Liver hepatocellular carcinoma
Methylation of N6 adenosine
Pheochromocytoma and paraganglioma
Skin cutaneous melanoma
The Cancer Genome Atlas
Uterine corpus endometrial carcinoma
Weighted gene coexpression network analysis
Chen K, Zhao BS, He C. Nucleic acid modifications in regulation of gene expression. Cell Chem Biol. 2016;23(1):74–85.
Yu J, Shen L, Liu Y, Ming H, Zhu X, Chu M, et al. The m6A methyltransferase METTL3 cooperates with demethylase ALKBH5 to regulate osteogenic differentiation through NF-kappaB signaling. Mol Cell Biochem. 2020;463(1-2):203–210.
Schubeler D. Function and information content of DNA methylation. Nature. 2015;517(7534):321–6.
Dor Y, Cedar H. Principles of DNA methylation and their implications for biology and medicine. Lancet. 2018;392(10149):777–86.
Lev Maor G, Yearim A, Ast G. The alternative role of DNA methylation in splicing regulation. Trends Genet. 2015;31(5):274–80.
Yearim A, Gelfman S, Shayevitch R, Melcer S, Glaich O, Mallm JP, et al. HP1 is involved in regulating the global impact of DNA methylation on alternative splicing. Cell Rep. 2015;10(7):1122–34.
Zhou L, Tian S, Qin G. RNA methylomes reveal the m(6)A-mediated regulation of DNA demethylase gene SlDML2 in tomato fruit ripening. Genome Biol. 2019;20(1):156.
Liu CJ, Hu FF, Xia MX, Han L, Zhang Q, Guo AY. GSCALite: a web server for gene set cancer analysis. Bioinformatics. 2018;34(21):3771–2.
We would like to thank the staff members involved in the TCGA Research Network.
This work was supported by grants from the National Natural Science Foundation of China (81802707, 81702682) and the Fundamental Research Funds for the Central Universities (19ykpy186). YPC also received support from the Sun Yat-sen University Cancer Center Promotion Program for Talented Youth of the National Natural Science Foundation of China.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Gene expression profile of m6A/5mC regulators across 33 cancer types. Pan-cancer normalized RNA-Seq by Expectation-Maximization (RSEM) data were used. For a given m6A/5mC regulators in a given cancer type, the non-scaled median expression level is presented.
Gene set differential expression profile of m6A/5mC regulators among 14 cancer types with available paired tumor-normal tissue expression data calculated by Gene Set Cancer Analysis (GSCA).
Heatmap showing percentage of cancers in which a pathway may be activated (red) or inhibited (blue) by the m6A/5mC regulators calculated by Gene Set Cancer Analysis (GSCA). Reverse phase protein array (RPPA) data of 32 cancer types from The Cancer Proteome Atlas (TCPA) are used for the calculation; acute myeloid leukemia (LAML) is not included. A total of 10 cancer related pathways are included (i.e., TSC/mTOR, RTK, RAS/MAPK, PI3K/AKT, Hormone ER, Hormone AR, EMT, DNA Damage Response, Cell Cycle, and Apoptosis pathways), and only m6A/5mC regulators that have function (activate or inhibit) in at least five cancer types are shown by GSCA.
Clinical relevance of the EMEs across 33 cancer types. a Forest plots showing the hazard ratios (HRs; squares) and 95% confidence intervals (CIs; horizontal ranges) of overall survival (OS) across 33 cancer types. Significant results are indicated by red (unfavorable prognosticators) or blue (favorable prognosticators) squares. b Kaplan-Meier plots showing unfavorable OS in the EMEhigh versus EMElow groups for KICH, ACC, LGG, CESC, SARC, LIHC, BRCA, and LUAD. P values for the two-sided log-rank test are shown. c Kaplan-Meier plots showing improved OS in the EMEhigh versus EMElow groups for HNSC, KIRC, and READ. P values for the two-sided log-rank test are shown.
Materials and Methods.
Details of the 33 cancer types from the TCGA.
Immuno-stromal signatures used in the current study.
About this article
Cite this article
Chen, YT., Shen, JY., Chen, DP. et al. Identification of cross-talk between m6A and 5mC regulators associated with onco-immunogenic features and prognosis across 33 cancer types. J Hematol Oncol 13, 22 (2020). https://doi.org/10.1186/s13045-020-00854-w
- m6A regulators
- 5mC regulators
- Pan-cancer analyses
- Genomic alterations
- Tumor microenvironment