Genome-wide DNA methylation profile of thymomas and potential epigenetic regulation of thymoma subtypes

The aim of the present study was to examine the whole-genome DNA methylation status of thymomas and identify differences in thymoma DNA methylation profiles. DNA methylation profiles of tissues (n=12) were studied using the Infinium MethylationEPIC BeadChip microarray (850K) and analyzed in relation to gene expression data. Functional annotation analysis of DNA methylation between the different groups was performed using the online tool GeneCodis3. In order to assess the diagnostic value of candidate DNA methylation markers, receiver operation characteristic (ROC) analysis was performed using the pROC package. A total of 10,014 CpGs were found to be differentially methylated (Δβ>0.2) between two thymoma types (type A and B). Combination analysis showed that 36 genes had differentially methylated CpG sites in their promoter region. ‘Pathways in cancer’, ‘focal adhesion’ and ‘regulation of actin cytoskeleton’ were the most enriched KEGG pathways of differentially methylated genes between tumor and controls. Among the 29 genes that were hypomethylated with a high expression, zinc finger protein 396 and Fraser extracellular matrix complex subunit 1 had the largest area under the curve. The present results may provide useful insights into the tumorigenesis of thymomas and a strong basis for future research on the molecular subtyping of epigenetic regulation in thymomas.