Skip to main content

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



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.


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.


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.


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.


Cervical cancer is a common female malignant tumor. According to the statistic results in 2020 all over the world [1], cervical cancer is the fourth most frequently diagnosed cancer and fourth leading cause of cancer death in women, with an estimated 604,172 new cases and 341,831 deaths worldwide in 2020. That account for 6.5% and 7.7% of new cancer cases and deaths worldwide, respectively. The burden of cervical cancer is higher in developing country, has a tremendous impact on the lifetime of millions of women [2].

Human papillomavirus (HPV) is a spherical and non-enveloped virus with a double-stranded circular DNA genome, showing strong squamous epithelium-like characteristics to infect damaged epithelial cells [3]. More than 200 types of HPV have been fully characterized [4]. HPVs have been divided into low-risk types (HPV-6,11, ect.) and high-risk types (HPV-16, 18, 52, 58, ect.) based on the risk of the virus to cause cancer in the uterine cervix [5].

It has been proved that long-term and persistent infection of HR-HPV is closely related to the occurrence and development of cervical cancer [6]. However, there are significant differences in HPV susceptible types in different regions. In China, HPV-16, HPV-52, and HPV-58 are the most common types, and the infection rate of HPV-58 is 15.31% [7], which is significantly higher than the global infection rate (4.7%) [8]. HPV-58 belongs to the Alpha-papillomavirus 9 (α9 HPV), and is relatively prevalent in China and other Asian countries [9]. HPV-58 genome can be divided into early region (E1, E2, E4, E5, E6, and E7), late region (L1 and L2) and long control region (LCR). Two HPV proteins, E6 and E7 (encoded by E6 gene and E7 gene, respectively), are the main drivers of cervical cancer development. The E6 protein binds to a ubiquitin ligase (E6AP) in the cell to form a complex, then binds specifically to p53 tumor suppressor protein. As a result, p53 is quickly degraded by proteases, resulting in decreasing the content of P53. E7 protein interacts with the retinoblastoma gene product (pRb), to inhibit the activity of the tumor suppressor protein through ubiquitination [10]. When HPV integrates into the host cell, E6 and E7 are invariably reserved and uncontrolled expressed, leading to the cell immortalization and malignant transformation [11,12,13].

Virus variations may cause differences in immune response and oncogenic potential [14]. There are some indications that the HPV-58 prototype and HPV-58 variants show a difference in their biological and biochemical characteristics to some extent. Several variations of HPV-58 have been reported in Latin America, Japan, southwest and northwest China [15,16,17,18], but few data were available on the variants of HPV-58 E6 and E7 genes in central China. In this study, we aim to analyze the gene polymorphisms and epidemic characteristics of HPV-58 E6 and E7 genes in Jingzhou area of central China and further evaluate the possible influence of E6-E7 sequence variations on immune reaction of antigens through describing the multiple putative B cell epitopes affected by the amino acid changes in the HPV-58 E6-E7. Our results could provide experimental data for the further study of therapeutic vaccines development, epidemiology, and prevention on HPV-58.

Materials and methods

Specimen collection

In this study, samples of cervical exfoliated cells were collected from the female patients with possible HR-HPV infection who underwent routine cervical screenings at the Department of Gynecology of Jingzhou Hospital from Spe. 2019 to Oct. 2021. All the samples were stored at 4 °C, and DNA was extracted within 1 day. Informed consent of patients had been obtained, and the study subjects’ privacy was protected before sample collection.

HPV DNA extraction and typing

Based on the method of magnetic beads, DNA was extracted according to the instruction of the nucleic acid extracting kit (Guangzhou Magen Biotechnology Co., Ltd.). 4 μL extraction products were used for HPV typing as a PCR template. According to the instruction of the HR-HPV typing test kit (Shanghai ZJ Bio-Tech Co., Ltd.), DNA extraction products were detected through the method of real-time quantitative PCR, including 15 types of HR-HPV (16,18,31,33,35,39,45,51,52,56,58,59,66,68,82), with advantages of excellent sensitivity, specificity and typing accurately (minimum detection limit: 104copies/mL; coincidence rate of positive reference: 100%; coincidence rate of negative reference: 100%; no cross-reaction with HPV-11, -6, -42, -43, -44, -73, -26, -53 types, Chlamydia trachomatis, Ureaplasma urealyticum, Neisseria gonorrhoeae, Mycoplasma hominis, Mycoplasma genitalium, Candida albicans, HIV-1, group B streptococcus). Finally, DNA extraction products of single positive samples of HPV-58 were selected and stored at − 80 °C for the follow-up experiments.

PCR amplification and sequencing

The E6 and E7 genes were amplified from HPV-58 single positive samples. According to the reference sequence provided by GenBank (D90400), software Primer Primer 6.0 was used to design sequence-specific primers of HPV-58 E6 and E7 genes following primer design principles. All primers were synthesized by Sangon Biotech (Sangon Biotech (Shanghai) Co., Ltd.). As followed, E6 F:5′- TTGGGTCACATTGTTCATG-3′, E6 R:5′-TCATAGCAGAATAGGTCAGT-3′; E7 F:5′- ATTTCGGGTCGTTGGACAGG-3′, E7 R:5′- TCTGTACCACTATCGTCTGCTGTT-3′. 25 μL reaction volume included 2 μL HPV DNA extraction, 15.875 μL ddH2O, 2.5 μL 10 × Buffer (Mg2+ plus) (Takara), 2 μL dNTPs (Takara), 0.125 μL Taq polymerase (Takara), and 1.25 μL forward and reverse primers (20 μM) (Sangon Biotech). The amplification conditions of E6 were initial denaturation at 95 °C for 3 min, followed by 35 cycles of 94 °C for 45 s, 58.2 °C for 45 s, 72 °C for 1 min, and the 72 °C for 10 min. The E7 amplification conditions were initial denaturation at 95 °C for 3 min, followed by 35 cycles of 94 °C for 45 s, 62.1 °C for 45 s, 72 °C for 1 min, and the 72 °C for 10 min. The PCR products were analyzed by electrophoresis in 2.5% agarose gel stained with GelRed nucleic acid dye (Biotium), and then sent to Sangon Biotech for sequencing.

Analysis of E6 and E7 genes

The nucleotide variants and amino acid substitutions were analyzed by comparing all sequencing results with HPV-58 reference sequence (D90400) via the DNAstar and NCBI blast.

Phylogenetic analysis

The Neighbor-Joining tree of HPV-58 E6 and E7 sequences was constructed by MEGA 11 software, Kimura 2-parameter model was used, and the Bootstrap replications was set at 1000. The following reference sequences collected from GeneBank sequences database were used to construct the phylogenetic branches: D90400 (A1), HQ537752 (A2), HQ537758 (A3), HQ537762 (B1), HQ537764 (B2), HQ537774 (C), HQ537768 (D1), HQ537770 (D2) [19].

Secondary structure analysis and selective pressure analysis

The online server PsiPred ( was used to predict the secondary structure of E6 and E7 amino acid sequences. The selective pressure was analyzed by PAML X software. The non-synonymous and synonymous nucleotide divergence (dN/dS) was calculated by Codeml program. Using the Bayes Empirical Bayes (BEB) analysis, the sites with the posterior probability ≥ 0.95 were identified as the positively selected sites.

B cell epitope prediction

We predicted linear B cell epitopes of E6 and E7 proteins by online software ABCpred (ABCpred submission page ( according to the default parameters. The higher the predicted score, the better the affinity of the epitope.


Epidemic characteristics of HR-HPV in Jingzhou, Hubei province

A total of 18,266 samples were collected in Jingzhou Hospital, Yangtze University from Sep. 2019 to Oct. 2021. The overall prevalence of HR-HPV infection was 18.58% (3394/18266), with 14.62% (2671/18266) of participants being positive for a single HPV type, and the age group with highest infection rate was 58 to 67 years (32.43%). The detailed results are shown in Table 1. The most frequently detected HR-HPV types were type 52 (1201/3394, 35.39%), 58 (577/3394, 17.00%), 16 (465/3394, 13.70%) and 51 (313/3394, 9.22%), the results are shown in Fig. 1. In 577 HPV-58 positive samples, 319 cases (319/577, 55.29%) were single infection, 258 cases (258/577, 44.71%) were multiple infection. Among the HPV-58 co-infection, mixed infection with HPV-68 was the most common (77/258, 29.84%), followed by HPV-52 (76/258, 29.46%). The peak age of HPV-58 infection mainly at older than 67 years old and 57 to 67 years old, with the positive rate of 8.93% and 8.68%, respectively (Table 1).

Table 1 Multiple infections and positive for high-risk HPV infection status in different age groups [(n, %)]
Fig. 1
figure 1

Prevalence of HPV in Jingzhou area. The horizontal axis shows different HPV types, the vertical axis shows infection rate. The numeric value above the histogram represents the HPV infection rate

Variation of E6 and E7 genes

There were 319 single positive samples of HPV-58, and some samples were failed amplified due to the low viral load or degradation of viral nucleic acid. Finally, 172 entire E6-E7 genes were successfully amplified and further analyzed. Compared with the reference sequence D90400, the prototype accounted for 18.02% (31/172) of the total sequencing samples, and the remaining 141 sequences were divided into 21 different variant groups named 58HB01-58HB21, which were submitted to the GenBank database and got the accession numbers were OL989036 to OL989056.

A total of 22 single nucleotide polymorphisms (SNPs) were observed in E6-E7 sequences. In E6 sequences, 10 SNPs were found, including 7 synonymous variants and 3 non-synonymous variants. The prevalent variants were A388C (K93N) (81/172) and C307T (57/172). In E7 sequences, 12 SNPs were observed, including 3 synonymous variants and 9 non-synonymous variants. T744G (104/172), G761A (G63D) (51/172), G694A (G41R) (47/172) and T803C (V77A) (46/172) were common variants. G261A (R51K), A310G, A340G and T407G in E6 sequence and A597G in E7 sequence were novel variants. Two non-synonymous variants (R51K and L100V) occurred in E6 sequences encoding the alpha helix, two non-synonymous variants (R9K and V77A) were found in E7 sequences encoding the strand. No insertion or deletion was detected in any sequence. The detailed results are shown in Table 2.

Table 2 The variations of HPV-58 E6-E7 gene

Phylogenetic tree

The phylogenetic tree was constructed with 21 unique HPV-58 E6-E7 sequences and 8 sub-lineage reference sequences by MEGA 11 software. The phylogenetic tree is shown in Fig. 2. All the variants were distributed in A lineage, with A1 sub-lineage was most abundant (115/172, 66.86%), followed by A2 (47/172, 27.33%) and A3 (10/172, 5.81%).

Fig. 2
figure 2

Neighbor-joining phylogenetic tree of HPV-58 E6–E7. The strains with black dot represented the reference sequences of sub-lineages, and the others were variant sequences. Only the bootstrap value > 70% was displayed

Selective pressure analysis

The variable dN/dS rate ratios were tested among the various lineages using the PAML X. In HPV-58 E6 sequence, there was no positively selected site. In E7, codon site G63 was identified as the positively selected site. All the results are shown in Table 3.

Table 3 Positively selected sites of HPV-58 E6-E7 gene

B cell epitope prediction of HPV-58 E6 and E7 proteins

To evaluate the influence of HPV-58 E6 and E7 sequence polymorphisms on the host immune response and recognition, online server ABCpred was used to predict the potential B cell epitopes of E6 and E7 proteins. The epitope prediction scores should be greater than 0.85. The results are shown in Table 4. The occurrence of non-synonymous variants affected the B cell epitopes. For E6, the A to C substitution at the nucleotide site 388 led to the change of the amino acid site 93 (K to N), where B cell epitope YSLYGDTLEQTLKKCL (E6 81–96) was located, and the epitope scores decreased from 0.90 to 0.88. The T407G (L100V) variant located in the B cell epitope TLEQTLKKCLNEIL (E6 87–100), made the epitope scores decrease to 0.86. For E7, the non-synonymous variants R9K affected the B cell epitope RGNNPTLREYILDLHPEP (E7 2–19). G41R affected the B cell epitope DEDEIGLDGPDGQAQP (E7 33–48). T61N/G63D/T64A affected the B cell epitope YYIVTCCYCGTTV (E7 53–66), and G63S also affected this epitope. T74A and V77A affected the B cell epitope RLCINSTTTDVRLCINS (E7 67–82). Except G41R variant, which made the epitope prediction scores increase from 0.86 to 0.91, the scores of all the other variants dropped below 0.85.

Table 4 The putative B cell epitopes of HPV-58 E6 and E7 proteins


Persistent infection of HR-HPV is the main factor for the occurrence and development of cervical cancer. The E6 and E7 proteins encode by E6 and E7 genes are the main carcinogenic proteins and play an important role in the cancer progression. HPV variants may be linked to the evolution of human groups to get an intrinsic geographical difference in prevalence and infections. The variants analysis of HPV E6 and E7 gene from clinical samples is of great significant for further study on vaccine and cervical cancer therapy in local area.

In this present study, we found that HPV-52, 58, and 16 types were prevalent types in Jingzhou area. The infection rate of HPV-58 was 17.00%, ranked second. This unusual high prevalence has also been reported in Chinese populations in Xi’an, Guangxi, and Chongqing [20,21,22].The prevalence of HPV-58 reached a peak when the patients were older than 67 years old in this study, suggesting that more attention should be paid to HPV screening in middle-aged women, which provided a basis for prevention, screening, and treatment of cervical cancer in Jingzhou area.

The results of E6-E7 sequencing showed that the sequence variability of E7 sequence (4.0%) was higher than that of E6 sequence (2.2%), suggesting that E6 may be a more suitable target for HPV therapeutic vaccine than E7. In E6 gene, C307T and A388C (K93N) were common variants. In E7 gene, T744G was most prevalent synonymous variant, G761A (G63D), G694A (G41R), and T803C (V77A) were common non-synonymous variants. There were four non-synonymous variants located in the sequences encoding the secondary structure of E6 and E7 proteins, may cause changes in oncoprotein folding and function, led to the differences in their ability to interact with tumor suppressor proteins, thus affecting the pathogenicity of HPV-58 [23]. Among the 21 variant groups, 58HB04, with the C307T/A694G/T744G/G761A variant, was most prevalent (43/172, 25.00%), which is consistent with a previous study in Colombian [24]. At present, there have been some researches on HPV-58 E6-E7 gene polymorphisms, and the common variants of E6 and E7 also have a high frequency in other areas [15, 16, 25,26,27]. Some studies have examined the HPV-58 E6 and E7 genetic diversity and their risk association with the development of cervical lesion. The A388C (K93N) substitution in E6 is more prevalent in women with normal cervix or low-grade lesions, shows a statistically significant negative trend of association with the severity of neoplasia in Hong Kong [28], which is consistent with the results of the study on HPV-58 in Shanghai area that A388C (K93N) substitution in E6 significantly reduces the risk of high grade squamous intraepithelial lesions (HSIL) [29]. In E7, the occurrence of C632T (T20I) and G760A (G63S) variants are positively correlated with the severity of neoplasia in Hong Kong [28], which is also found in a study in Zhejiang by Ding et al. [30]. The exact mechanism for the increased oncogenicity of the variants needs further investigation. G261A (R51K), A310G, A340G, T407G (L100V) of E6 gene and A597G of E7 gene were newly identified variants in this study, which have not been reported in previous studies, may provide real and effective data for the development of regional therapeutic vaccines. Our study can provide basic data for further studies of the relationship between E6-E7 gene variations and the severity of the cervix lesion.

HPV-58 type has four lineages of A, B, C and D. Lineage A is significantly more frequent in Asia than in America, lineage B is mainly distributed in America, and lineage C and D seem to be more frequent in Africa [31]. In this study, lineage A was most prevalent in Jingzhou area, especially for sub-lineage A1, which was also prevalent in Yunnan, Sichuan, Zhejiang etc. [16, 25, 32]. The higher oncogenicity and infection rate of sub-lineage A1 may explain the reason why HPV-58 has a high incidence of cervical cancer in China.

At present, all HPV vaccines on the market (Cervarix, Gardasil, and Gardasil-9) are prophylactic vaccines based on L1 protein immunogenicity, which play an effective role in preventing HPV infection by producing specific neutralizing antibodies. However, the vaccines do not prevent cancer development in individuals who are already infected [33]. Therefore, the development of therapeutic vaccines has broad application prospects in the management and control of HPV infection, and HPV E6 and E7 genes are ideal targets for the design of therapeutic vaccines against HPV [34, 35]. ABCpred server was used to predict the B cell epitopes of HPV-58 E6 and E7 proteins and analyzed the effect of amino acid changes caused by gene variation on the B cell epitope. In this study, only G694A (G41R) nucleotide substitution increase the epitope score of DEDEIGLDGPDGQAQP (E7 33–48), while the other variations all led to the decrease of epitope score. B cells play a significant role in HPV-related tumor immunotherapy and response to cervical lesion and cancers caused by HPV [36]. By analyzing the impacts of the occurrence of non-synonymous variants on the affinity of B cell epitopes, we can further understand their influence on the immune function of host cells. The B cell epitopes predicted in this study can be considered as potential candidates for therapeutic vaccines against HPV-58 development in the Chinese population, and the occurrence of amino acid substitution must be considered in the development of vaccine. These predicted epitopes need to be verified in vitro and in vivo experiments in subsequent studies.

This study analyzed the positive selection of HPV-58 E6 and E7 sequences in Jingzhou area. The main feature of positive selection is that it causes an unusual repaid increase in the frequency of alleles, which allows a species to adapt quickly to environmental changes [37]. The positively selected sites G63 found in E7 gene may have evolutionary significance in HPV-58 adaption, meaning that this variation is beneficial for the virus to adapt to the changes in the environment. However, relevant research is minimal, and more studies are needed to confirm this hypothesis.


In this study, we analyzed the gene polymorphism of HPV-58 E6-E7 sequence and amino acid substitution of E6-E7 proteins in Jingzhou area, constructed a phylogenetic tree to analyze the lineage characteristics, and predicted the B cell epitope to clarify the regional specificity of HPV-58 gene variation. Some new variations and B cell epitopes were observed in this study, and we also got a general knowledge of the gene polymorphism of HPV-58 in Jingzhou area. This study can supplement the knowledge of gene diversity of HPV-58 in central China, and provide an experimental basis for subsequent studies on epidemiology, evolution, pathogenicity and therapeutic targets of HPV-58.

Availability of data and materials

The data generated during the current study are available in the NCBI repository (Home—Nucleotide—NCBI ( The sequence data were submitted to GenBank with accession numbers OL989036-OL989056.



Human papillomavirus type 58


High-risk human papillomavirus


Polymerase chain reaction


Single nucleotide polymorphisms


Long control region


High grade squamous intraepithelial lesion

The amino acid acronyms R, K, N, L, V, T, I, G, S, D and A:

Arginine, lysine, asparagine, leucine, valine, threonine, isoleucine, glycine, serine, aspartic acid, and alanine


  1. Sung H, et al. Global Cancer Statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2021;71(3):209–49.

    Article  Google Scholar 

  2. Huang J, et al. Global distribution, risk factors, and recent trends for cervical cancer: a worldwide country-level analysis. Gynecol Oncol. 2022;164(1):85–92.

    Article  Google Scholar 

  3. Roden RBS, Stern PL. Opportunities and challenges for human papillomavirus vaccination in cancer. Nat Rev Cancer. 2018;18(4):240–54.

    CAS  Article  Google Scholar 

  4. Bzhalava D, Eklund C, Dillner J. International standardization and classification of human papillomavirus types. Virology. 2015;476:341–4.

    CAS  Article  Google Scholar 

  5. Van Doorslaer K, et al. Papillomaviruses: evolution, Linnaean taxonomy and current nomenclature. Trends Microbiol. 2011;19(2):49–50 (author reply 50-1).

    Article  Google Scholar 

  6. Burd EM. Human papillomavirus laboratory testing: the changing paradigm. Clin Microbiol Rev. 2016;29(2):291–319.

    CAS  Article  Google Scholar 

  7. Zhang J, Cheng K, Wang Z. Prevalence and distribution of human papillomavirus genotypes in cervical intraepithelial neoplasia in China: a meta-analysis. Arch Gynecol Obstet. 2020;302(6):1329–37.

    CAS  Article  Google Scholar 

  8. Li N, et al. Human papillomavirus type distribution in 30,848 invasive cervical cancers worldwide: variation by geographical region, histological type and year of publication. Int J Cancer. 2011;128(4):927–35.

    CAS  Article  Google Scholar 

  9. Chan PK. Human papillomavirus type 58: the unique role in cervical cancers in East Asia. Cell Biosci. 2012;2(1):17.

    Article  Google Scholar 

  10. Narisawa-Saito M, Kiyono T. Basic mechanisms of high-risk human papillomavirus-induced carcinogenesis: roles of E6 and E7 proteins. Cancer Sci. 2007;98(10):1505–11.

    CAS  Article  Google Scholar 

  11. Munger K, Howley PM. Human papillomavirus immortalization and transformation functions. Virus Res. 2002;89(2):213–28.

    CAS  Article  Google Scholar 

  12. Doorbar J, et al. The biology and life-cycle of human papillomaviruses. Vaccine. 2012;30(Suppl 5):F55-70.

    CAS  Article  Google Scholar 

  13. zur Hausen H. Papillomaviruses and cancer: from basic studies to clinical application. Nat Rev Cancer. 2002;2(5):342–50.

    CAS  Article  Google Scholar 

  14. Chen Z, et al. Correction: E6 and E7 gene polymorphisms in human papillomavirus types-58 and 33 identified in Southwest China. PLoS ONE. 2017;12(7):e0181475.

    Article  Google Scholar 

  15. Liu JH, et al. Variations of human papillomavirus type 58 E6, E7, L1 genes and long control region in strains from women with cervical lesions in Liaoning province, China. Infect Genet Evol. 2012;12(7):1466–72.

    CAS  Article  Google Scholar 

  16. Xi J, et al. The polymorphisms of LCR, E6, and E7 of HPV-58 isolates in Yunnan, Southwest China. Virol J. 2018;15(1):76.

    Article  Google Scholar 

  17. Xin CY, et al. Analysis of E6 variants of human papillomavirus type 33, 52 and 58 in Japanese women with cervical intraepithelial neoplasia/cervical cancer in relation to their oncogenic potential. Cancer Lett. 2001;170(1):19–24.

    CAS  Article  Google Scholar 

  18. Conde-Ferraez L, et al. Genetic variability in E6 and E7 oncogenes from human papillomavirus type 58 in Mexican women. Intervirology. 2017;60(6):235–46.

    CAS  Article  Google Scholar 

  19. Burk RD, Harari A, Chen Z. Human papillomavirus genome variants. Virology. 2013;445(1–2):232–43.

    CAS  Article  Google Scholar 

  20. Han X, et al. Prevalence and genotype distribution of human papillomavirus infection among women aged 30–65 years in Xi’an, China: a population-based study of 14,655 women. Hum Vaccin Immunother. 2021;17(12):5439–46.

    CAS  Article  Google Scholar 

  21. Wei X, Lu Q, Wang S. Prevalence characteristics of cervical human papillomavirus genotypes in Nanning, China: a 10-year survey of 77,756 women from one medical center. J Med Virol. 2021;94:2787–95.

    Article  Google Scholar 

  22. Luo Q, et al. Prevalence and genotype distribution of high-risk human papillomavirus infection among women with cervical cytological abnormalities in Chongqing, China, 2014–2020. Diagn Cytopathol. 2021;49(12):1237–43.

    Article  Google Scholar 

  23. He J, et al. The polymorphism analysis and epitope predicted of Alphapapillomavirus 9 E6 in Sichuan, China. Virol J. 2022;19(1):14.

    CAS  Article  Google Scholar 

  24. Molano M, et al. Follow-up study of HPV58 variants in women with incident HPV58 infection from a Colombian cohort. Oncol Lett. 2018;16(2):2511–6.

    PubMed  PubMed Central  Google Scholar 

  25. Yu JH, et al. Genetic variability and oncogenic risk association of human papillomavirus type 58 E6 and E7 genes in Taizhou area, China. Gene. 2019;686:171–6.

    CAS  Article  Google Scholar 

  26. Yang L, et al. Genetic variability of HPV-58 E6 and E7 genes in Southwest China. Infect Genet Evol. 2014;21:395–400.

    Article  Google Scholar 

  27. Wu EQ, et al. Profile of physical status and gene variation of human papillomavirus 58 genome in cervical cancer. J Gen Virol. 2009;90(Pt 5):1229–37.

    CAS  Article  Google Scholar 

  28. Chan PK, et al. Association of human papillomavirus type 58 variant with the risk of cervical cancer. J Natl Cancer Inst. 2002;94(16):1249–53.

    CAS  Article  Google Scholar 

  29. Zhao J, et al. Phylogeny and polymorphism in the E6 and E7 of human papillomavirus: alpha-9 (HPV16, 31, 33, 52, 58), alpha-5 (HPV51), alpha-6 (HPV53, 66), alpha-7 (HPV18, 39, 59, 68) and alpha-10 (HPV6, 44) in women from Shanghai. Infect Agent Cancer. 2019;14:38.

    Article  Google Scholar 

  30. Ding T, et al. Distribution of human papillomavirus 58 and 52 E6/E7 variants in cervical neoplasia in Chinese women. Gynecol Oncol. 2010;119(3):436–43.

    Article  Google Scholar 

  31. Chan PK, et al. Identification of human papillomavirus type 58 lineages and the distribution worldwide. J Infect Dis. 2011;203(11):1565–73.

    Article  Google Scholar 

  32. Wu H, et al. Lineage distribution and E2 sequence variation of high-risk human papillomavirus types isolated from patients with cervical cancer in Sichuan province, China. Arch Virol. 2015;160(11):2845–55.

    CAS  Article  Google Scholar 

  33. Clark KT, Trimble CL. Current status of therapeutic HPV vaccines. Gynecol Oncol. 2020;156(2):503–10.

    CAS  Article  Google Scholar 

  34. Yang A, et al. Perspectives for therapeutic HPV vaccine development. J Biomed Sci. 2016;23(1):75.

    Article  Google Scholar 

  35. Hoppe-Seyler K, et al. The HPV E6/E7 oncogenes: key factors for viral carcinogenesis and therapeutic targets. Trends Microbiol. 2018;26(2):158–68.

    CAS  Article  Google Scholar 

  36. He J, et al. Identification of variants and therapeutic epitopes in HPV-33/HPV-58 E6 and E7 in Southwest China. Virol J. 2019;16(1):72.

    Article  Google Scholar 

  37. Iyengar VK, Reeve HK, Eisner T. Paternal inheritance of a female moth’s mating preference. Nature. 2002;419(6909):830–2.

    CAS  Article  Google Scholar 

Download references


Not application.


This work was supported by Central Funds Guiding the Local Science and Technology Development of Hubei Province, Grant/Award Number: 2019ZYYD066.

Author information




ZPY and BM conceived and designed the study. ZPY, CLZ, MXY, and PL collected the samples. ZPY performed the experiments and drafted the manuscript. BM and QG revised the paper. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Bing Mei.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the Ethics Committee of Jingzhou Hospital, Yangtze University, and all the work followed the ethics guidelines of the hospital. Informed consent of patients had been obtained and the privacy of the study subjects were protected before sample collection.

Consent for publication

Not application.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Yang, Z., Zhang, C., Luo, P. et al. Genetic variability of E6 and E7 genes of human papillomavirus type 58 in Jingzhou, Hubei Province of central China. Virol J 19, 71 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Human papillomavirus type 58 (HPV-58)
  • Variants
  • Epitope
  • Phylogenetic analysis
  • Selection pressure