Genome-wide DNA methylation profiling of HPV-negative leukoplakia and gingivobuccal complex cancers

Gingivobuccal complex oral squamous cell carcinoma (GBC-OSCC) is an aggressive malignancy with high mortality often preceded by premalignant lesions, including leukoplakia. Previous studies have reported genomic drivers in OSCC, but much remains to be elucidated about DNA methylation patterns across different stages of oral carcinogenesis. There is a serious lack of biomarkers and clinical application of biomarkers for early detection and prognosis of gingivobuccal complex cancers. Hence, in search of novel biomarkers, we measured genome-wide DNA methylation in 22 normal oral tissues, 22 leukoplakia, and 74 GBC-OSCC tissue samples. Both leukoplakia and GBC-OSCC had distinct methylation profiles as compared to normal oral tissue samples. Aberrant DNA methylation increases during the different stages of oral carcinogenesis, from premalignant lesions to carcinoma. We identified 846 and 5111 differentially methylated promoters in leukoplakia and GBC-OSCC, respectively, with a sizable fraction shared between the two sets. Further, we identified potential biomarkers from integrative analysis in gingivobuccal complex cancers and validated them in an independent cohort. Integration of genome, epigenome, and transcriptome data revealed candidate genes with gene expression synergistically regulated by copy number and DNA methylation changes. Regularised Cox regression identified 32 genes associated with patient survival. In an independent set of samples, we validated eight genes (FAT1, GLDC, HOXB13, CST7, CYB5A, MLLT11, GHR, LY75) from the integrative analysis and 30 genes from previously published reports. Bisulfite pyrosequencing validated GLDC (P = 0.036), HOXB13 (P