Cervical Cancer Development: Implications of HPV16 E6E7-NFX1-123 Regulated Genes

SIMPLE SUMMARY: High-risk human papillomavirus (HPV) causes 4.5% of cancers and nearly all cervical cancers. HPV’s carcinogenic potential depends on its misappropriation of cellular proteins by HPV’s oncoproteins E6 and E7. High-risk HPV type 16 (HPV16) E6 binds directly to the cellular protein NFX1-123 and dysregulates proliferation, differentiation, and immunity genes. The effect of HPV16 E7 has not been studied in relation to HPV16 E6-NFX1-123-mediated dysregulation. As HPV expresses both oncogenes, and HPV carcinogenesis requires E6 and E7, it is valuable to investigate what dysregulations occur in this context. It is also important to understand their clinical and prognostic ramifications. This study’s goal was to define the gene expression profile regulated by HPV16 E6, E7, and NFX1-123 across cervical precancers and cancers, identify genes correlating with disease progression, assess patient survival, and validate findings in cell models. Finding correlates of survival and disease progression aids in biomarker identification and focuses future studies. ABSTRACT: High-risk human papillomavirus (HR HPV) causes nearly all cervical cancers, half of which are due to HPV type 16 (HPV16). HPV16 oncoprotein E6 (16E6) binds to NFX1-123, and dysregulates gene expression, but their clinical implications are unknown. Additionally, HPV16 E7’s role has not been studied in concert with NFX1-123 and 16E6. HR HPVs express both oncogenes, and transformation requires their expression, so we sought to investigate the effect of E7 on gene expression. This study’s goal was to define gene expression profiles across cervical precancer and cancer stages, identify genes correlating with disease progression, assess patient survival, and validate findings in cell models. We analyzed NCBI GEO datasets containing transcriptomic data linked with cervical cancer stage and utilized LASSO analysis to identify cancer-driving genes. Keratinocytes expressing 16E6 and 16E7 (16E6E7) and exogenous NFX1-123 were tested for LASSO-identified gene expression. Ten out of nineteen genes correlated with disease progression, including CEBPD, NOTCH1, and KRT16, and affected survival. 16E6E7 in keratinocytes increased CEBPD, KRT16, and SLPI, and decreased NOTCH1. Exogenous NFX1-123 in 16E6E7 keratinocytes resulted in significantly increased CEBPD and NOTCH1, and reduced SLPI. This work demonstrates the clinical relevance of CEBPD, NOTCH1, KRT16, and SLPI, and shows the regulatory effects of 16E6E7 and NFX1-123.