Genetic variability of E6 and E7 genes of human papillomavirus type 58 in Jingzhou, Hubei Province of central China

Background Cervical cancer is a common malignant tumor in women, with a high mortality rate, has great harm to women's health. Long-term and persistent infection of high-risk human papillomavirus (HR-HPV) is the main reason of the occurrence and development of cervical cancer. Methods The infection rate of HPV-58 is higher in the Jingzhou area. In this study, 172 complete HPV-58 E6-E7 sequences were amplified by polymerase chain reaction (PCR), the amplified products were sequenced, and the gene variations of HPV-58 E6-E7 were analyzed. A Neighbor-Joining phylogenetic tree was constructed by MEGA 11. The secondary structure of E6 and E7 protein was investigated. PAML X was used to analyze the selective pressure. The B cell epitopes of E6 and E7 proteins in HPV-58 were predicted by ABCpred server. Results In E6 sequences, 10 single nucleotide variants were observed, including 7 synonymous and 3 non-synonymous variants. In E7 sequences, 12 single nucleotide variants were found, including 3 synonymous variants and 9 non-synonymous variants. There are 5 novel variants. The phylogenetic analysis showed that all the E6-E7 sequences were distributed in A lineage. No positively selected site was found in E6 sequence, but G63 in E7 sequences was identified as positively selected site. Some amino acid substitutions affected multiple B cell epitopes. Conclusion Various E6 and E7 mutational data may prove useful for development of better diagnostic and vaccines for the region of Jingzhou, Hubei province of central China.