Discovery and validation of DNA methylation markers for overall survival prognosis in patients with thymic epithelial tumors

BACKGROUND: The current prognosis of thymic epithelial tumors (TETs) is according to the World Health Organization (WHO) histologic classification and the Masaoka staging system. These methods of prognosis have certain limitations in clinical application and there is a need to seek new method for determining the prognosis of patients with TETs. To date, there have been no studies done on the use of DNA methylation biomarkers for prognosis of TETs. The present study was therefore carried out to identify DNA methylation biomarkers that can determine the overall survival in patients with TETs. METHODS: Bioinformatic analysis of TCGA 450 K methylation array data, transcriptome sequencing data, WHO histologic classification and Masaoka staging system was performed to identify differentially expressed methylation sites between thymoma and thymic carcinoma as well as the different DNA methylation sites associated with the overall survival in patients with TETs. Using pyrosequencing, 4 different methylation sites (cg05784862, cg07154254, cg02543462, and cg06288355) were sequenced from tumor tissues of 100 Chinese patients with TETs. A prognostic model for TETs was constructed using these four methylation sites. RESULTS: The TCGA dataset showed 5155 and 6967 hyper- and hypomethylated CpG sites in type A–B3 group and type C group, respectively, of which 3600 were located within the gene promoter regions. One hundred thirty-four genes were silenced by promoter hypermethylation and 174 mRNAs were upregulated. Analysis of univariate and multivariate Cox regression showed significant association between the methylation levels of 187 sites and the overall survival in patients with TETs. cg05784862(KSR1), cg07154254(ELF3), cg02543462(ILRN), and cg06288355(RAG1) were identified as independent prognostic factors for overall survival in patients with TETs after adjusting for Masaoka staging in 100 Chinese patients. The prognostic model which consists of the four abovementioned genes had higher accuracy for predicting the 5-year overall survival in patients with TETs as compared to the Masaoka clinical staging. (Time-dependent ROC analysis AUC 1.000 vs 0.742, P = 2.7 × 10(−6)). CONCLUSIONS: The methylation levels of cg05784862(KSR1), cg07154254(ELF3), cg02543462(ILRN), and cg06288355(RAG1) sites are associated with the progression of TETs and may serve as new biomarkers for predicting the overall survival in patients with TETs. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (10.1186/s13148-019-0619-z