Thyroid carcinoma is the most prevalent endocrine cancer globally and the primary cause of cancer-related mortality. Epigenetic modifications are progressively being linked to metastasis. This study aimed to examine whole-genome DNA methylation patterns and the gene expression profiles in thyroid cancer tissue samples using a MethylationEPIC BeadChip (850K), RNA sequencing, and a targeted bisulfite sequencing assay. The results of the Illumina Infinium human methylation kit (850K) analyses identified differentially methylated CpG locations (DMPs) and differentially methylated CpG regions (DMRs) encompassing nearly the entire genome with high resolution and depth. Gene ontology and KEGG pathway analyses revealed that the genes associated with DMRs belonged to various domain-specific ontologies, including cell adhesion, molecule binding, and proliferation. The RNA-Seq study found 1627 differentially expressed genes, 1174 of which that were up-regulated and 453 of which that were down-regulated. The targeted bisulfite sequencing assay revealed that CHST2, DPP4, DUSP6, ITGA2, SLC1A5, TIAM1, TNIK, and ABTB2 methylation levels were dramatically lowered in thyroid cancer patients when compared to the controls, but GALNTL6, HTR7, SPOCD1, and GRM5 methylation levels were significantly raised. Our study revealed that the whole-genome DNA methylation patterns and gene expression profiles in thyroid cancer shed new light on the tumorigenesis of thyroid cancer.
Keywords: DNA methylation; MethylationEPIC BeadChip (850K); RNA-Seq; targeted bisulfite sequencing assay; thyroid carcinoma.