Abstract
HPV infection is one of the most commonly transmitted diseases among the global population. While it can be asymptomatic, non-genital HPV infection often gives rise to cutaneous warts, which are benign growths arising from the epidermal layer of the skin. This study aimed to produce a global analysis of the ways in which cutaneous wart formation affected the CpG island methylome. The Infinium MethylationEPIC BeadChip microarray was utilized in order to quantitatively interrogate CpG island methylation in genomic DNA extracted from 24 paired wart and normal skin samples. Differential methylation analysis was carried out by means of assigning a combined rank score using RnBeads. The 1000 top-ranking CpG islands were then subject to Locus Overlap Analysis (LOLA) for enrichment of genomic ranges, while signaling pathway analysis was carried out on the top 100 differentially methylated CpG islands. Differential methylation analysis illustrated that the most differentially methylated CpG islands in warts lay within the ITGB5, DTNB, RBFOX3, SLC6A9, and C2orf27A genes. In addition, the most enriched genomic region sets in warts were Sheffield's tissue-clustered DNase hypersensitive sites, ENCODE's segmentation and transcription factor binding sites, codex sites, and the epigenome sites from cistrome. Lastly, signaling pathway analysis showed that the GRB2, GNB1, NTRK1, AXIN1, and SKI genes were the most common regulators of the genes associated with the top 100 most differentially methylated CpG islands in warts. Our study shows that HPV-induced cutaneous warts have a clear CpG island methylation profile that sets them apart from normal skin. Such a finding could account for the temporary nature of warts and the capacity for individuals to undergo clinical remission.