Genetic variability of human papillomavirus type 18 based on E6, E7 and L1 genes in central China

Abstract Background High-risk human papillomavirus (HR-HPV) infection is an important factor for the development of cervical cancer. HPV18 is the second most common HR-HPV after HPV16. Methods In this study, MEGA11 software was used to analyze the variation and phylogenetic tree of HPV18 E6-E7 and L1 genes. The selective pressure to E6, E7 and L1 genes was estimated using pamlX. In addition, the B cell epitopes of L1 amino acid sequences and T cell epitopes of E6-E7 amino acid sequences in HPV18 were predicted by ABCpred server and IEDB website, respectively. Results A total of 9 single nucleotide variants were found in E6-E7 sequences, of which 2 were nonsynonymous variants and 7 were synonymous variants. Twenty single nucleotide variants were identified in L1 sequence, including 11 nonsynonymous variants and 9 synonymous variants. Phylogenetic analysis showed that E6-E7 and L1 sequences were all distributed in A lineage. In HPV18 E6, E7 and L1 sequences, no positively selected site was found. The nonconservative substitution R545C in L1 affected hypothetical B cell epitope. Two nonconservative substitutions, S82A in E6, and R53Q in E7, impacted multiple hypothetical T cell epitopes. Conclusion The sequence variation data of HPV18 may lay a foundation for the virus diagnosis, further study of cervical cancer and vaccine design in central China.