Global gene methylation profiling of common warts caused by human papillomaviruses infection

Infection with the human papillomaviruses (HPV) often involves the epigenetic modification of the host genome. Despite its prevalence among the population, host genome methylation in HPV-induced warts is not clearly understood. In this study, genome-wide methylation profiling was carried out on paired healthy skin and wart samples in order to investigate the effects that benign HPV infection has on gene methylation status. To overcome this gap in knowledge, paired wart (n = 12) and normal skin (n = 12) samples were obtained from Arab males in order to perform DNA extraction and subsequent genome-wide methylation profiling on the Infinium Methylation EPIC Bead Chip microarray. Analysis of differential methylation revealed a clear pattern of discrimination between the wart and normal skin samples. In warts, the most differentially methylated (DM) genes included long non-coding RNAs (AC005884, AL049646.2, AC126121.2, AP001790.1, and AC107959.3), microRNAs (MIR374B, MIR596, MIR1255B1, MIR26B, and MIR196A2),snoRNAs (SNORD114-22, SNORD70, and SNORD114-31), pseudogenes (AC069366.1, RNU4ATAC11P, AC120057.1, NANOGP3, AC106038.2, TPT1P2, SDC4P, PKMP3, and VN2R3P), and protein-coding genes (AREG, GJB2, C12orf71, AC020909.2, S100A8, ZBED2, FABP7, and CYSLTR1). In addition, pathway analysis revealed that, among the most differentially methylated genes, STAT5A, RARA, MEF2D, MAP3K8, and THRA were the common regulators. It can be observed that HPV-induced warts involve a clear and unique epigenetic alteration to the host genome.