Alterations of 63 hub genes during lingual carcinogenesis in C57BL/6J mice

To identify potential biomarkers of lingual cancer, 75 female C57BL/6J mice were subjected to 16-week oral delivery of 4-nitroquinoline-1-oxide (4NQO; 50 mg/L), with 10 mice used as controls. Lingual mucosa samples representative of normal tissue (week 0) and early (week 12) and advanced (week 28) tumorigenesis were harvested for microarray and methylated DNA immunoprecipitation sequencing (MeDIP-Seq). Combined analysis with Short Time-series Expression Miner (STEM), the Cytoscape plugin cytoHubba, and screening of differentially expressed genes enabled identification of 63 hub genes predominantly altered in the early stage rather than the advanced stage. Validation of microarray results was carried out using qRT-PCR. Of 63 human orthologous genes, 35 correlated with human oral squamous cell carcinoma. KEGG analysis showed “pathways in cancer”, involving 13 hub genes, as the leading KEGG term. Significant alterations in promoter methylation were confirmed at Tbp, Smad1, Smad4, Pdpk1, Camk2, Atxn3, and Cdh2. HDAC2, TBP, and EP300 scored ≥10 on Maximal Clique Centrality (MCC) in STEM profile 11 and were overexpressed in human tongue cancer samples. However, expression did not correlate with smoking status, tumor differentiation, or overall survival. These results highlight potentially useful candidate biomarkers for lingual cancer prevention, diagnosis, and treatment.