Expression of the E5 Oncoprotein of HPV16 Impacts on the Molecular Profiles of EMT-Related and Differentiation Genes in Ectocervical Low-Grade Lesions

Infection with human papillomavirus type 16 (HPV16) is one of the major risk factors for the development of cervical cancer. Our previous studies have demonstrated the involvement of the early oncoprotein E5 of HPV16 (16E5) in the altered isoform switch of fibroblast growth factor receptor 2 (FGFR2) and the consequent expression in human keratinocytes of the mesenchymal FGFR2c isoform, whose aberrant signaling leads to EMT, invasiveness, and dysregulated differentiation. Here, we aimed to establish the possible direct link between these pathological features or the appearance of FGFR2c and the expression of 16E5 in low-grade squamous intraepithelial lesions (LSILs). Molecular analysis showed that the FGFR2c expression displayed a statistically significant positive correlation with that of the viral oncoprotein, whereas the expression values of the epithelial FGR2b variant, as well as those of the differentiation markers keratin 10 (K10), loricrin (LOR) and involucrin (INV), were inversely linked to the 16E5 expression. In contrast, the expression of EMT-related transcription factors Snail1 and ZEB1 overlapped with that of 16E5, becoming a statistically significant positive correlation in the case of Snail2. Parallel analysis performed in human cervical LSIL-derived W12 cells, containing episomal HPV16, revealed that the depletion of 16E5 by siRNA was able to counteract these molecular events, proving to represent an effective strategy to identify the specific role of this viral oncoprotein in determining LSIL oncogenic and more aggressive profiles. Overall, coupling in vitro approaches to the molecular transcript analysis in ectocervical early lesions could significantly contribute to the characterization of specific gene expression profiles prognostic for those LSILs with a greater probability of direct neoplastic progression.