Pan-Cancer Genome-Wide DNA Methylation Analyses Revealed That Hypermethylation Influences 3D Architecture and Gene Expression Dysregulation in HOXA Locus During Carcinogenesis of Cancers

DNA methylation dysregulation during carcinogenesis has been widely discussed in recent years. However, the pan-cancer DNA methylation biomarkers and corresponding biological mechanisms were seldom investigated. We identified differentially methylated sites and regions from 5,056 The Cancer Genome Atlas (TCGA) samples across 10 cancer types and then validated the findings using 48 manually annotated datasets consisting of 3,394 samples across nine cancer types from Gene Expression Omnibus (GEO). All samples’ DNA methylation profile was evaluated with Illumina 450K microarray to narrow down the batch effect. Nine regions were identified as commonly differentially methylated regions across cancers in TCGA and GEO cohorts. Among these regions, a DNA fragment consisting of ∼1,400 bp detected inside the HOXA locus instead of the boundary may relate to the co-expression attenuation of genes inside the locus during carcinogenesis. We further analyzed the 3D DNA interaction profile by the publicly accessible Hi-C database. Consistently, the HOXA locus in normal cell lines compromised isolated topological domains while merging to the domain nearby in cancer cell lines. In conclusion, the dysregulation of the HOXA locus provides a novel insight into pan-cancer carcinogenesis.