Abstract
Human papillomavirus (HPV) E6 and E7 genes are biomarkers and drivers of the progression of cervical cancer (CxCa). The aim of this study was to investigate the relationship between HPV16 E6, E7 gene mutations and the occurrence and development of CxCa. Cervical exfoliated cells and clinical data of patients with cervical diseases were collected. Sample DNA was extracted, the E6 and E7 gene fragments were amplified by PCR, and the mutations were detected by Sanger sequencing and compared with standard sequences. Microarray was used to sequence the transcriptome of cells. Data of transcriptome analyzed and visualized using R software and its packages. Analysis of clinical characteristics demonstrated the association of HPV16 infection with CxCa (P < .05). Sanger sequencing results showed that the mutation sites of E6 gene included T178G/A, T350G, A131C, and T241G; among these, A131C and T241G were synonymous mutations. The mutation sites of E7 gene included A647G, T846C, G666A, T843C, and T760C, and all of them were synonymous mutations except A647G. There was no significant difference in the distribution of HPV16 E6, E7 mutations among CxCa, cervical intraepithelial neoplasia, and infection groups (P > .05). Compared with the non- CxCa group, gene ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analysis of differentially expressed genes (DEGs) showed more significant enrichment of DEGs in the biological processes, pathways, and diseases closely related to cancer. Compared with the non-mutation group, the DEGs in the E6, E7 gene mutation group were significantly enriched in the events related to infection and immunity. To summarize, HPV16 may be associated with the occurrence and development of CxCa, but HPV16 E6 and E7 gene mutations have little effect on the occurrence and development of CxCa. Individual differences may have a greater effect on the progression of CxCa.