Abstract
Papillary thyroid carcinoma (PTC) is the most common type of thyroid cancer, and its incidence rate is rapidly growing. It is necessary to understand the pathogenesis of PTC to develop effective diagnosis methods. Promoter methylation has been recognized to contribute to the alterations in gene expression observed in tumorigenesis. Our RNA-seq data identified 1191 differentially expressed mRNAs and 147 differentially expressed lncRNAs in PTC. Next, promoter methylation of these genes was detected by reduced representation bisulfite sequencing (RRBS) technology and comprehensively analyzed to identify differential methylation. In total, 14 genes (13 mRNAs and 1 lncRNA), in which methylation was intimately involved in regulating gene expression, were proposed as novel diagnostic biomarkers. To gain insights into the relationships among these 14 genes, a core co-function network was constructed based on co-expression, co-function and co-methylation data. Notably, CXCL12 was identified as an essential gene in the network that was closely connected with the other genes. These data suggested that CXCL12 down-regulation in PTC may be caused by promoter hypermethylation. Our study was the first to perform an RRBS analysis for PTC and suggested that CXCL12 may contribute to PTC development by methylation-mediated epigenetic regulation of gene expression.