Immune Signature-Based Subtypes of Cervical Squamous Cell Carcinoma Tightly Associated with Human Papillomavirus Type 16 Expression, Molecular Features, and Clinical Outcome

Substantial heterogeneity exists within cervical cancer that is generally infected by human papillomavirus (HPV). However, the most common histological subtype of cervical cancer, cervical squamous cell carcinoma (CSCC), is poorly characterized regarding the association between its heterogeneity and HPV oncoprotein expression. We filtered out 138 CSCC samples with infection of HPV16 only as the first step; then we compressed HPV16 E6/E7 expression as HPVpca and correlated HPVpca with the immunological profiling of CSCC based on supervised clustering to discover subtypes and to characterize the differences between subgroups in terms of the HPVpca level, pathway activity, epigenetic dysregulation, somatic mutation frequencies, and likelihood of responding to chemo/immunotherapies. Supervised clustering of immune signatures revealed two HPV16 subtypes (namely, HPV16-IMM and HPV16-KRT) that correlated with HPVpca and clinical outcomes. HPV16-KRT is characterized by elevated expression of genes in keratinization, biological oxidation, and Wnt signaling, whereas HPV16-IMM has a strong immune response and mesenchymal features. HPV16-IMM exhibited much more epigenetic silencing and significant mutation at FBXW7, while MUC4 and PIK3CA were mutated frequently for HPV16-KRT. We also imputed that HPV16-IMM is much more sensitive to chemo/immunotherapy than is HPV16-KRT. Our characterization tightly links the expression of HPV16 E6/E7 with biological and clinical outcomes of CSCC, providing valuable molecular-level information that points to decoding heterogeneity. Together, these results shed light on stratifications of CSCC infected by HPV16 and shall help to guide personalized management and treatment of patients