Crosstalk between 5-methylcytosine and N6-methyladenosine machinery defines disease progression, therapeutic response and pharmacogenomic landscape in hepatocellular carcinoma

Abstract Background Accumulated evidence highlights the significance of the crosstalk between epigenetic and epitranscriptomic mechanisms, notably 5-methylcytosine (5mC) and N6-methyladenosine (m6A). Herein, we conducted a widespread analysis regarding the crosstalk between 5mC and m6A regulators in hepatocellular carcinoma (HCC). Methods Pan-cancer genomic analysis of the crosstalk between 5mC and m6A regulators was presented at transcriptomic, genomic, epigenetic, and other multi-omics levels. Hub 5mC and m6A regulators were summarized to define an epigenetic and epitranscriptomic module eigengene (EME), which reflected both the pre- and post-transcriptional modifications. Results 5mC and m6A regulators interacted with one another at the multi-omic levels across pan-cancer, including HCC. The EME scoring system enabled to greatly optimize risk stratification and accurately predict HCC patients’ clinical outcomes and progression. Additionally, the EME accurately predicted the responses to mainstream therapies (TACE and sorafenib) and immunotherapy as well as hyper-progression. In vitro, 5mC and m6A regulators cooperatively weakened apoptosis and facilitated proliferation, DNA damage repair, G2/M arrest, migration, invasion and epithelial-to-mesenchymal transition (EMT) in HCC cells. The EME scoring system was remarkably linked to potential extrinsic and intrinsic immune escape mechanisms, and the high EME might contribute to a reduced copy number gain/loss frequency. Finally, we determined potential therapeutic compounds and druggable targets (TUBB1 and P2RY4) for HCC patients with high EME. Conclusions Our findings suggest that HCC may result from a unique synergistic combination of 5mC-epigenetic mechanism mixed with m6A-epitranscriptomic mechanism, and their crosstalk defines therapeutic response and pharmacogenomic landscape.