Prognostic and therapeutic prediction by screening signature combinations from transcriptome–methylome interactions in oral squamous cell carcinoma

DNA methylation pattern in oral squamous cell carcinoma (OSCC) remains poorly described. This study aimed to perform a genome-wide integrated analysis of the transcriptome and methylome and assess the efficacy of their prognostic signature model in patients with OSCC. We analyzed transcriptome and methylome data from 391 OSCC samples and 41 adjacent normal samples. A total of 8074 differentially expressed genes (DEGs) and 10,084 differentially expressed CpGs (DMCpGs) were identified. Then 241 DEGs with DMCpGs were identified. According to the prognostic analysis, the prognostic signature of methylation-related differentially expressed genes (mrDEGPS) was established. mrDEGPS consisted of seven prognostic methylation-related genes, including ESRRG, CCNA1, SLC20A1, COL6A6, FCGBP, CDKN2A, and ZNF43. mrDEGPS was a significant stratification factor of survival (P < 0.00001) irrespective of the clinical stage. The immune effector components, including B cells, CD4(+) T cells, and CD8(+) T cells, were decreased in the tumor environment of patients with high mrDEGPS. Immune checkpoint expressions, including CTLA-4, PD-1, LAG3, LGALS9, HAVCR2, and TIGHT, were comprehensively elevated (P < 0.001). The estimated half-maximal inhibitory concentration difference between low- and high-risk patients was inconsistent among chemotherapeutic drugs. In conclusion, the transcriptome–methylome interaction pattern in OSCC is complex. mrDEGPS can predict patient survival and responses to immunotherapy and chemotherapy and facilitate clinical decision-making in patients with OSCC.