Skip to main content
Download PDF

Phylogenetic analysis and antigenic epitope prediction for E6 and E7 of Alpha-papillomavirus 9 in Taizhou, China

BMC Genomics volume 25, Article number: 507 (2024) Cite this article

Abstract

Background

Alpha-papillomavirus 9 (α-9) is a member of the human papillomavirus (HPV) α genus, causing 75% invasive cervical cancers worldwide. The purpose of this study was to provide data for effective treatment of HPV-induced cervical lesions in Taizhou by analysing the genetic variation and antigenic epitopes of α-9 HPV E6 and E7.

Methods

Cervical exfoliated cells were collected for HPV genotyping. Positive samples of the α-9 HPV single type were selected for E6 and E7 gene sequencing. The obtained nucleotide sequences were translated into amino acid sequences (protein primary structure) using MEGA X, and positive selection sites of the amino acid sequences were evaluated using PAML. The secondary and tertiary structures of the E6 and E7 proteins were predicted using PSIPred, SWISS-MODEL, and PyMol. Potential T/B-cell epitopes were predicted by Industrial Engineering Database (IEDB).

Results

From 2012 to 2023, α-9 HPV accounted for 75.0% (7815/10423) of high-risk HPV-positive samples in Taizhou, both alone and in combination with other types. Among these, single-type-positive samples of α-9 HPV were selected, and the entire E6 and E7 genes were sequenced, including 298 HPV16, 149 HPV31, 185 HPV33, 123 HPV35, 325 HPV52, and 199 HPV58 samples. Compared with reference sequences, 34, 12, 10, 2, 17, and 17 nonsynonymous nucleotide mutations were detected in HPV16, 31, 33, 35, 52, and 58, respectively. Among all nonsynonymous nucleotide mutations, 19 positive selection sites were selected, which may have evolutionary significance in rendering α-9 HPV adaptive to its environment. Immunoinformatics predicted 57 potential linear and 59 conformational B-cell epitopes, many of which are also predicted as CTL epitopes.

Conclusion

The present study provides almost comprehensive data on the genetic variations, phylogenetics, positive selection sites, and antigenic epitopes of α-9 HPV E6 and E7 in Taizhou, China, which will be helpful for local HPV therapeutic vaccine development.

Peer Review reports

Background

Cervical cancer remains the fourth most common cancer affecting women’s health worldwide, especially in China. Persistent high-risk HPV infection is considered the primary aetiological factor for cervical cancer [1]. HPV is a small, double-stranded, circular DNA virus that exclusively infects epithelial cells of the skin or mucosa. The genome size of HPV is approximately 8.0 kb, containing six early genes (E1, E2, E4, E5, E6, E7), two late genes (L1, L2), and long control region (LCR). The oncoproteins encoded by the E6 and E7 genes play a crucial role in HPV-driven viral carcinogenesis and cancers [2, 3]. E6 binds to the tumour suppressor protein p53 and prevents its translocation, and mediate the cellular transformation by inhibiting the ability of p53 to activate. E7 binds to the retinoblastoma protein (Rb) and induces cells to enter into premature S-phase by disrupting Rb-E2F complexes. These processes lead to impaired p53 and Rb functions, involving DNA repair, cell cycle, apoptosis, and ultimately result in immortalization of HPV-infected cells [4].

To date, over 200 HPVs have been identified, which are divided into five main genera referred to as alpha (α)-, beta (β)-, gamma (γ)-, Mu(μ)-, and Nu(ν)-papillomaviruses [5]. In addition, these HPVs are divided into high- and low-risk types based on their carcinogenic potential. Of which, at least 14 are classified as high-risk HPV types, such as HPV16, 18, 31, 33, 35, 39, 45, 51, 52, 56, 58, 59, 66 and 68. It is worth noting that almost all high-risk HPV types belong to the α genus, including α-5 (HPV51), α-6 (HPV56, 66), α-7 (HPV18, 39, 45, 59, 68), and α-9 (HPV16, 31, 33, 35, 52, 58). Of which, α-9 species causes approximately 75% invasive cervical cancers worldwide [6]. Our previous study suggested that women infected with HPV16, 18, 31, 33, 52, or 58 are more likely to have CIN2 + lesions [7]. Moreover, the incidence of CIN2 + caused by HPV33, 52, and 58 was greater than that caused by HPV 18 in Taizhou [7]. Different HPV types exhibit differences in their immunogenicity, adaptability, and carcinogenicity, which may be due to different genetic variations and selection pressures during the evolutionary process [8, 9]. Therefore, focus should be placed on these common carcinogenic HPV types, such as α-9 species in the Taizhou region.

Persistent infection of human epithelial cells by HPV leads to integration of the viral DNA into the host genome, usually disrupting the E1 and/or E2 genes [10]. HPV integration is a key event for cervical carcinogenesis, leading to structural aberrations in the host genome or abnormal gene expression of target genes. Of which, the expression of E6 and E7 oncoproteins cause cellular immortalization and neoplastic transformation. The most frequent integration sites include SHKBP1, ERBB3, CASP8, HLA-A, HLA-B, TGFBR2, PIK3CA, EP300, FBXW7, PTEN, NFE2L2, ARID1A, KRAS, MAPK1, etc. [11, 12]. The transformed cervical cells show expression of HPV E6 or E7, with antigenic epitopes on their protein surface that stimulate B lymphocytes to produce antibodies [13, 14]. Viral antigen peptides are presented at the cell surface through human leukocyte antigen (HLA) and are recognized by CD8 + cytotoxic T lymphocytes (CTLs) [15,16,17]. Therefore, E6 or E7 molecules are considered ideal targets for HPV therapeutic vaccines, inducing cell-mediated immunity by stimulating CTLs in immune response strategies or eliciting humoral immunity by activating B lymphocytes to produce specific antibodies [18, 19].

Unfortunately, as a major issue for cervical cancer, local data available on HPV therapeutic vaccines are still limited in China, and almost no consideration is given to E6 or E7 gene mutations. Hence, there is an urgent need to further study the genetic variation, positive selection site, protein structure, and antigen epitopes of α-9 HPV E6 or E7 to provide data to explore effective treatment of HPV-induced cervical lesions in Taizhou, China.

Materials and methods

Study population

This study was ethically approved by the Institutional Medical Ethics Review Board of Taizhou Hospital, China. Cervical exfoliated cells were collected from women who underwent cervical cancer screening at Taizhou Hospital. The specimens were collected by cervical scraping and stored in 2.5 ml of cell preservation buffer at -20 °C. Before specimen collection, written informed consent was obtained from all participants, and the participants’ privacy was strictly protected.

HPV genotyping

HPV types were identified using the HPV Genotyping Kit (Tellgen Corporation, China), which was approved by the China’s FDA (Certified No. (2014): 3,400,847). The protocol for HPV genotyping has been described in detail in our previous studies [7, 20]. Briefly, the kit uses a set of biotinylated amplimers and multiplex HPV genotyping methods with bead-based Luminex suspension array technology, which is able to simultaneously identify 14 high-risk HPV types including 16, 18, 31, 33, 35, 39, 45, 51, 52, 56, 58, 59, 66, 68 and β-globin gene (internal control).

PCR amplification and sequencing

Based on the reference sequences of α-9 HPV types in GenBank, specific primer pairs for the entirety of the E6 and E7 regions were designed using the Primer-BLAST tool (ncbi.nlm.nih.gov/tools/primer-blast). The primers, PCR conditions, amplicon size, and reference sequences are listed in Table S1. Genomic DNA was extracted using a DNA Extraction Kit (#GK0122, GENEray, China). PCR products were purified and sequenced at BGI, and all data were confirmed by repeating PCR and sequencing reactions at least twice. In this study, genetic variant data for HPV31 and 35 were combined with data from our previous studies on HPV16, 33, 52 and 58 [21,22,23,24].

Phylogenetic analysis and homology comparison

All successfully acquired nucleotide sequences were aligned by BioEdit. Then, a phylogenetic tree of α-9 HPV E6 and E7 variation patterns was constructed by the maximum-likelihood method with one thousand bootstrap replicates using MEGA X. The phylogenetic tree was constructed using 201 complete α-9 HPV E6 and E7 sequences, including 64 HPV16, 16 HPV31, 15 HPV33, 5 HPV35, 27 HPV52, 25 HPV58, and 49 reference sequences downloaded from NCBI.

Selective pressure analysis

The CodeML program in PAML (abacus.gene.ucl.ac.uk/software/paml.html) was used to calculate the nonsynonymous (dN)/synonymous (dS) ratio (ω) for selective pressure analysis. If nonsynonymous mutations are favoured by Darwinian selection, they will be fixed at a higher rate than synonymous mutations, resulting in dN > dS, ω = dN/dS > 1 [25].

Protein structure analysis

The obtained nucleotide sequences were translated into amino acid sequences (protein primary structure) using MEGA X. Then, the secondary and tertiary structures of the E6 and E7 proteins of α-9 HPV were predicted using PSIPred (bioinf.cs.ucl.ac.uk/psipred) and Swiss-model (swissmodel.ExPASy.org), respectively.

Prediction of linear and conformational B-cell epitopes

Linear B-cell epitopes were predicted from the primary sequence of the E6 or E7 protein using sequence-based methods (Kolaskar and Tongaonkar’s antigenicity, tools.immuneepitope.org/bcell/). Conformational B-cell epitopes were predicted from the tertiary structure of E6 or E7 using ElliPro, which identifies protrusions in antigen surfaces (tools.immuneepitope.org/ellipro).

Prediction of cytotoxic T-cell epitopes

CTL epitopes were predicted from the E6 or E7 protein of α-9 HPV types using the NetCTL server, which accepts the FASTA format (tools.immuneepitope.org/netchop). NetCTL gives results for 9-mer peptides together with their predicted MHC binding affinity, binding affinity rescale value, C-terminal cleavage affinity, and TAP transport efficiency. C-terminal cleavage weights and TAP transport efficiency were calculated using the default values 0.15 and 0.05, respectively; 9-mer peptides with a prediction score > 0.75 were considered to be potential CTL epitopes.

Results

The prevalence of α-9 HPV types in Taizhou

A total of 60,259 women underwent the first round of HPV screening in Taizhou Hospital from December 2012 to June 2023. The overall rate of high-risk HPV infection was 17.3% (10,423/60,259), including infection with multiple HPV types (4.5%). All high-risk HPV-positive samples belonged to the α genus, including α-5 (782, 7.5%), α-6 (1466, 14.1%), α-7 (3521, 33.8%), and α-9 (7815, 75.0%) (Fig. 1, Table S2). α-9 HPV accounted for 75.0% of the high-risk HPV-positive samples in Taizhou, both alone and in combination with other types. In this study, we selected the α-9 HPV single-type-positive samples (HPV 16, 31, 33, 35, 52, 58) for subsequent analyses, including 298 HPV16, 149 HPV31, 185 HPV33, 123 HPV35, 325 HPV52, and 199 HPV58 samples. The clinicopathological information of the study population are shown in Table S3.

Fig. 1
figure 1

Distribution of different high-risk HPV genotypes in Taizhou from 2012 to 2023

Full size image

Nucleotide variations of α-9 HPV types

Compared with the reference sequences from GenBank, 53 nucleotide mutations were detected in HPV16, with 34 (34/53, 64.2%) nonsynonymous mutations (24 in E6 and 10 in E7). In addition, 12 (12/24, 50.0%), 10 (10/20, 50.0%), 2 (2/7, 28.6%), 17 (17/35, 48.6%), and 17 (17/26, 65.4%) nonsynonymous nucleotide mutations were detected in the HPV31, 33, 35, 52, and 58 E6 and E7 genes, respectively. A summary of nucleotide variations and amino acid substitutions in the E6 and E7 regions of α-9 HPV are shown in Tables 1 and 2.

Table 1 Nucleotide variations and amino acid substitutions in α-9 HPV E6 region
Full size table
Table 2 Nucleotide variations and amino acid substitutions in α-9 HPV E7 region
Full size table

Of all genetic variations in our study, 31 newly reported nonsynonymous mutations have never been reported: T105G(M8R), C110G(Q10E), G126A(R15Q), A134C(K18Q), T137G(L19V), A152G(T24A), A296C(K72Q), T310G(F76L), T434G(C118G) in HPV16 E6, and A619T(T20S), A647C(N29T), G676C(D39H), T730C(F57L), G791T(R77L), G823T(G88*) in HPV16 E7; A299G(T64A), A446C(E113D), C537G(P143G) in HPV31 E6, and C788T(R77C) in HPV31 E7; T196G (C30G), G458T (R117L) in HPV33 E6; G108C(E3Q), G267A(D56N), T292G(I64S) in HPV52 E6; A194C(K29Q), G225A(R39Q), C228T(S40F), A390C(K94T), G432A(R108K), and C534A(P142H) in HPV58 E6 and T678A(D35E) in HPV58 E7.

Phylogenetic construction

All nucleotide sequences of this study were submitted to GenBank, and accession numbers were obtained (HPV16E6E7: MT681266-MT681329, HPV31E6E7: OR540563-OR540578, HPV33E6E7: OQ672665-OQ672679, HPV35E6E7: OR540579-OR540583, HPV52E6E7: ON529577-ON529603, HPV58E6E7: MH348918-MH348942). A phylogenetic tree was constructed from 201 complete α-9 HPV E6 and E7 nucleotide sequences (152 obtained from our study and 49 downloaded from GenBank) (Fig. 2). Furthermore, we constructed phylogenetic trees based on the E6 and E7 amino acid sequences (Figures S1- S2).

Fig. 2
figure 2

Phylogenetic tree of α-9 HPV variants. Maximum-likelihood analysis (with MEGA X) of E6 and E7 nucleotide sequences was inferred from 201 complete α-9 HPV E6/E7 sequences, including 64 HPV16, 16 HPV31, 15 HPV33, 5 HPV35, 27 HPV52, 25 HPV58, and 49 reference sequences

Full size image

Selective pressure analysis

The results of dN/dS analysis for α-9 HPV E6 and E7 genes showed a high nonsynonymous mutation rate for positively selected sites. Nineteen positive selection sites were detected, including 32E of HPV16 E6, 93K of HPV33 E6, 93R of HPV52 E6, 29S*, 77R*, 88-** of HPV16 E7, 45A, 97Q of HPV33 E7, 63E of HPV35 E7, 52S, 64D of HPV52 E7, 9R, 20T, 35D, 41G, 61T, 63G, 64T, and 77V of HPV58 E7. However, no reliable HPV31 E6, HPV35 E6, HPV58 E6, or HPV31 E7 positive selection sites were detected (Table 3).

Table 3 Positive section sites of α-9 HPV E6 and E7 sequences
Full size table

Protein structure analysis and homology modelling

Nucleotide nonsynonymous substitution changes the amino acid composition, which affects the structure and function of proteins. Our analysis showed that α-9 E6 and E7 are composed of residues 148–158 and 97–99, respectively. The template-target pairwise sequence alignment for α-9 HPV E6 and E7. More details are shown in Table S1 and Figures S3-S4. All amino acid substitutions in E6 and E7 are shown in Fig. 3. As depicted in Figure S5, six E6 or E7 proteins of α-9 genus HPV are highly homologous, so their correlations are included in our research. As shown in Fig. 3, the majority of amino acid substitutions are located on the outer edge of E6 or E7 proteins and near the zinc granule, which is situated in the active site. Interestingly, we found that the 93rd residue is not only a common nonsynonymous mutation in the E6 region but also a positive selection site.

Fig. 3
figure 3

3D structural models of α-9 HPV E6 or E7 proteins

Full size image

We selected the HPV variant with the highest infection rate in each type for homology modelling. The template coverage for the E6 protein of α-9 HPV included most of the protein sequence, no involving only a short N-terminal stretch (see Figure S5, A-F for the homology models and Figure S3 for the template-target pairwise sequence alignments). However, the template coverage for the E7 protein was low, excluding approximately 50 residues at the N-terminal stretch (see Figure S5, G-L for the homology models and Figure S4 for the template-target pairwise sequence alignments).

Prediction of linear and conformational B-cell epitopes

The predicted linear B-cell epitopes of α-9 HPV E6 or E7 proteins are presented in Table 4. Conformational B-cell epitopes were predicted from protein tertiary structure models using ElliPro. Amino acid residues, the number of residues, the sequence location and the PI score of the predicted conformational epitopes are given in Table 5, and the graphical depiction of these epitopes are provided in Figure S6. Immunoinformatics predicted 57 potential linear and 59 conformational B-cell epitopes, many of which are also predicted as CTL epitopes (Tables 4, 5 and 6).

Table 4 Predicted linear B-cell epitopes of E6 and E7 proteins of α-9 genus HPV
Full size table
Table 5 ElliPro predicted conformational B-cell epitopes of E6 and E7 proteins of α-9 genus HPV
Full size table
Table 6 NetCTL predicted CTL epitopes of E6 and E7 proteins of α-9 genus HPV
Full size table

Prediction of cytotoxic T-cell epitopes

Based on the literature, we selected the HLA-A*02 and HLA-B*62 supertypes for NetCTL epitope prediction, as well as CTL epitopes that overlap with the sites for the predicted B-cell epitopes (linear and/or conformational). The predicted CTL epitope 9-mer peptides with their predicted MHC-I binding affinity, rescaled binding affinity, proteasomal C-terminal cleavage affinity, and TAP transport efficiency are indicated in Table 6, with an overall prediction score threshold of 0.75.

Discussion

Cervical cancer remains a major challenge for women’s health, with approximately 600,000 new cases and 340,000 deaths worldwide every year [1, 26]. High-risk HPV genotypes are the main cause of cervical cancer, which can be largely prevented through HPV vaccination and cervical screening. However, the cervical cancer burden for women remains heavy in China, with only 3% of women aged 9–45 years receiving complete HPV vaccination [27, 28]. Although prophylactic vaccines are effective at protecting against approximately 90% of HPV infections, their benefits in eliminating preexisting infection are limited. Therefore, potential therapeutic vaccines need to be developed for treatment of persistent HPV infection or cervical lesions, with the goal of activating adaptive immune responses by presenting viral antigen peptides to the immune system. However, because the genetic variations of HPV genotypes show some degree of geographical differences, it is recommended that the ideal therapeutic vaccine be based on a local type.

The Taizhou area is located along the central coast of Zhejiang Province in China, with high prevalence and pathogenicity of α-9 HPV (HPV 16, 31, 33, 35, 52, 58), especially HPV52 and 58 [20]. According to our previous findings, the odds ratio (OR) for CIN2 + in women infected with α-9 HPV is 3.2 when compared to women infected with α-5, α-6, or α-7 [7]. The genetic variation of E6 or E7 genes may correlate highly with cancer risk in Taizhou [21,22,23]. Therefore, E6 or E7 molecules might be regarded as ideal targets for HPV therapeutic vaccine development and cervical cancer treatment. The purpose of this study was to provide data for effective prevention and treatment of HPV-induced cervical lesions in Taizhou by analysing the phylogenetic tree and epitope prediction of α-9 HPV E6 or E7.

In Taizhou, we found that the most common nucleotide mutations observed in the E6 gene were T178G (D32E, 192/298, positive selection) of HPV16; C520T (I138V, 63/149) of HPV31; A387C (K93N, 45/185, positive selection) of HPV33; T341C (W78R, 122/123) of HPV35; A379G (K93R, 313/325, positive selection) of HPV52; and A388C (K93N, 87/199, positive selection) of HPV58. Most high-frequency nonsynonymous mutations are positive selection sites, which contribute to HPV adaptation. In addition, their corresponding genetic variants have been widespread in Taizhou, including MT681266 (38.9%, 116/298, HPV16), OR540577 (24.8%, 37/149, HPV31), OQ672675 (20%, 37/185, HPV33), OR540579 (62.6%, 77/123, HPV35), ON529581 (60.9%, 198/325, HPV52), MH348927 (22.1%, 44/199, HPV58). Interestingly, we found that the 93rd residue is not only a common nonsynonymous mutation in the E6 region but also a positive selection site, which may have evolutionary significance in rendering HPV33, 52, and 58 adaptive to their environment. Notably, three amino acid substitutions at the 93rd residue were detected in HPV52 E6: K93R (96.3%), K93G (2.5%), and K93Q (2.2%). In addition, the K93N substitution was observed in both HPV33 and HPV58. These 93rd-amino acid substitutions may lead to differences in the ability of these viruses to bind to p53, thus affecting their carcinogenicity [29]. The D32E substitution in HPV16 E6 has been confirmed to be associated with cancer risk, possibly because the amino acid substitution affects the E6-p53 interaction [30].

In the E7 gene, the most common nucleotide mutation observed was A647G (N29S, 195/298, positive selection) of HPV16. Similarly, there were three amino acid substitutions at the 29th residue in HPV16 E7: N29S (65.44%), N29H (9.40%), and N29T (0.34%). In addition, all HPV58 E7 with G63S carries T20I, and all HPV58 E7 samples with G63D carry G41R. In our previous study, it was reported that HPV58 E7 T20I/G63S substitutions increase risk of HPV carcinogenesis [24]. Boon et al. [31] reported that T20I/G63S substitutions possess greater ability to degrade Rb, immortalize, and transform primary cells.

Consistently with our selective pressure analysis, positive selection sites for 32E in HPV16 E6, 29S in HPV16 E7 were found in Kunming, Southwest of China [32]. Positive selection sites for 93K in HPV33 E6, 45A, 97Q in HPV33 E7 were found in Sichuan province, Southwest of China [33]. Notably, positive selection site for 63G in HPV58 E7, which has been widely reported in other regions of China, including Sichuan province [33] and Yunnan province [34, 35] in Southwest China, Hubei province in central China [36], as well as the present study (Southeast China). However, there were no reports of positive selection sites for HPV58 E6 in China [33,34,35,36]. No positive selection sites for HPV16 E6 and E7 genes have been reported in central China [37]. No positive selection sites for HPV52 E6 and E7 genes have been reported in central and Southwest China [38, 39]. In addition, selective pressure analysis from non-Asian population showed different results, with 10R, 14H, 83V in HPV E6 and 85G in HPV16 E7 under positive selective pressure [40, 41]. Therefore, these results indicate that genetic variations among HPV types may lead to biological advantages through fixed mutations in their genomes and that even small variations might lead to minor adaptive improvements. Furthermore, the genetic variations of HPV may differ in terms of infectivity, carcinogenic potential, and host immune response. Therefore, the data provided in this study may have significant implications for understanding the biological differences of HPVs in Taizhou, as well as for developing local therapeutic vaccines.

HLA participates in the local immune response of viral infection through its target recognition function, blocking HPV infection or preventing tumour cell invasion and metastasis. However, a minority of infected cells can escape host immune surveillance, causing persistent HPV infection. Immunoinformatics provides new strategies for identifying ideal epitopes for HPV therapeutic vaccine targets. Our predicted T- and B-cell epitopes may be used for development of vaccines targeting specific HPV variants, and our results suggest that amino acid substitution may influence these epitopes. For example, the prediction score of the HPV16 E6 CTL epitope 25-33ELQTTIHDI was 1.0904, and because of the D32E substitution the score of the epitope 25-33ELQTTIHEI increased to 1.3272; the HPV16 E6 predicted epitope 31-39HEIILECVY became a new epitope because of the mutation D32E. Therefore, the substitution of the positive selection site D32E in HPV16 E6 influences its antigenic epitopes, which may make it more difficult to detect by the immune system, thereby enhancing adaptability of HPV to the environment. In addition, Li et al. [37] suggested that non-conservative substitutions of amino acids should be fully considered when developing therapeutic vaccines, such as H31Y, D32N, D32E, I34M, L35V, E36Q, L45P, N65S, and K75T in HPV16 E6. Chen et al. [33] suggested that K93N in HPV33 E6, Q97L in HPV33 E7, R145K in HPV58 E6, and T20I in HPV58 E7 belong to ideal B-cell epitopes.

In addition, our results suggest that the CTL epitope of HPV16 E6 is 75-83KFYSKISEY; homologous epitopes were identified in HPV31 E6 68-76RFYSKVSEF, HPV33 E6 68-76RFLSKISEY, HPV35 E6 68-76KFYSKISEY, HPV52 E6 68-76RFLSKISEY, and HPV58 E6 68-76RLLSKISEY. Homologous epitopes were also identified in HPV16 E7 11-19YMLDLQPET, HPV31 E7 11-19YVLDLQPEA, HPV35 E7 11-19YVLDLEPEA, and HPV52 E7 11-19YILDLQPET. These E7 epitopes are restricted to HLA-A*02, which is the first noteworthy allele of the HPV-restricted epitope [42,43,44]. Riemer et al. [45] found that CTLs can recognize and lyse target cells presenting HPV16 E7 11-19YMLDLQPET. However, further experiments are required to validate the potential HPV therapeutic vaccines predicted through immunoinformatics.

Conclusions

In summary, this is the first almost comprehensive study to explore the genetic variations, phylogenetics, positive selection sites, and antigenic epitopes of α-9 HPV E6 and E7 molecules in Taizhou, China, and the results will be helpful for local HPV therapeutic vaccine development.

Availability of data and materials

All data generated during this study are included in this published article. All sequences of α-9 HPV E6 and E7 genes are available in the GenBank database with the accession codes of MT681266 to MT681329, OR540563 to OR540578, OQ672665 to OQ672679, OR540579 to OR540583, ON529577 to ON529603, and MH348918 to MH348942. The links are https://www.ncbi.nlm.nih.gov/nuccore/MT681266 ~ https://www.ncbi.nlm.nih.gov/nuccore/MH348942.

References

  1. Cohen PA, Jhingran A, Oaknin A, Denny L. Cervical cancer. Lancet. 2019;393:169–82.

    Article  PubMed  Google Scholar 

  2. Alejo M, Alemany L, Clavero O, Quiros B, Vighi S, Seoud M, et al. Contribution of Human papillomavirus in neuroendocrine tumors from a series of 10,575 invasive cervical cancer cases. Papillomavirus Res. 2018;5:134–42.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Matsumoto K, Oki A, Furuta R, Maeda H, Yasugi T, Takatsuka N, et al. Predicting the progression of cervical precursor lesions by human papillomavirus genotyping: a prospective cohort study. Int J Cancer. 2011;128:2898–910.

    Article  CAS  PubMed  Google Scholar 

  4. Ferenczy A, Franco E. Persistent human papillomavirus infection and cervical neoplasia. Lancet Oncol. 2002;3:11–6.

    Article  CAS  PubMed  Google Scholar 

  5. Arroyo Mühr LS, Lagheden C, Hassan SS, Eklund C, Dillner J. The International Human Papillomavirus Reference Center: Standardization, collaboration, and quality assurance in HPV research and diagnostics. J Med Virol. 2023;95:e29332.

    Article  PubMed  Google Scholar 

  6. Chen Z, Schiffman M, Herrero R, Desalle R, Anastos K, Segondy M, et al. Evolution and taxonomic classification of human papillomavirus 16 (HPV16)-related variant genomes: HPV31, HPV33, HPV35, HPV52, HPV58 and HPV67. PLoS One. 2011;6:e20183.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Xu H, Lin A, Shao X, Shi W, Zhang Y, Yan W. Diagnostic accuracy of high-risk HPV genotyping in women with high-grade cervical lesions: evidence for improving the cervical cancer screening strategy in China. Oncotarget. 2016;7:83775–83.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Tommasino M. The human papillomavirus family and its role in carcinogenesis. Semin Cancer Biol. 2014;26:13–21.

    Article  CAS  PubMed  Google Scholar 

  9. Fehrmann F, Laimins LA. Human papillomaviruses: targeting differentiating epithelial cells for malignant transformation. Oncogene. 2003;22:5201–7.

    Article  CAS  PubMed  Google Scholar 

  10. Tsakogiannis D, Gortsilas P, Kyriakopoulou Z, Ruether IG, Dimitriou TG, Orfanoudakis G, et al. Sites of disruption within E1 and E2 genes of HPV16 and association with cervical dysplasia. J Med Virol. 2015;87:1973–80.

    Article  CAS  PubMed  Google Scholar 

  11. Hu Z, Zhu D, Wang W, Li W, Jia W, Zeng X, et al. Genome-wide profiling of HPV integration in cervical cancer identifies clustered genomic hot spots and a potential microhomology-mediated integration mechanism. Nat Genet. 2015;47:158–63.

    Article  CAS  PubMed  Google Scholar 

  12. Cancer Genome Atlas Research Network; Albert Einstein College of Medicine; Analytical Biological Services; Barretos Cancer Hospital; Baylor College of Medicine; Beckman Research Institute of City of Hope; Buck Institute for Research on Aging, et al. Integrated genomic and molecular characterization of cervical cancer. Nature. 2017;543:378–84.

    Article  Google Scholar 

  13. Yang A, Farmer E, Wu TC, Hung CF. Perspectives for therapeutic HPV vaccine development. J Biomed Sci. 2016;23:75.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Gomes D, Silvestre S, Duarte AP, Venuti A, Soares CP, Passarinha L, et al. In silico approaches: a way to unveil novel therapeutic drugs for cervical cancer management. Pharmaceuticals (Basel). 2021;14:741.

    Article  CAS  PubMed  Google Scholar 

  15. Mukherjee AG, Wanjari UR, Gopalakrishnan AV, Kannampuzha S, Murali R, Namachivayam A, et al. Exploring the Molecular Pathogenesis, Pathogen Association, and Therapeutic Strategies against HPV Infection. Pathogens. 2022;12:25.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Jhunjhunwala S, Hammer C, Delamarre L. Antigen presentation in cancer: insights into tumour immunogenicity and immune evasion. Nat Rev Cancer. 2021;21:298–312.

    Article  CAS  PubMed  Google Scholar 

  17. Markowitz LE, Unger ER. Human papillomavirus vaccination. N Engl J Med. 2023;388:1790–8.

    Article  PubMed  Google Scholar 

  18. Mo Y, Ma J, Zhang H, Shen J, Chen J, Hong J, et al. Prophylactic and therapeutic HPV vaccines: current scenario and perspectives. Front Cell Infect Microbiol. 2022;12:909223.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Tian Y, Hu D, Li Y, Yang L. Development of therapeutic vaccines for the treatment of diseases. Mol Biomed. 2022;3:40.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Xu HH, Lin A, Chen YH, Dong SS, Shi WW, Yu JZ, et al. Prevalence characteristics of cervical human papillomavirus (HPV) genotypes in the Taizhou area, China: a cross-sectional study of 37 967 women from the general population. BMJ Open. 2017;7:e014135.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Dai MZ, Qiu Y, Di XH, Shi WW, Xu HH. Association of cervical carcinogenesis risk with HPV16 E6 and E7 variants in the Taizhou area. China BMC Cancer. 2021;21:769.

    Article  CAS  PubMed  Google Scholar 

  22. Yan ZY, Di XH, Qiu Y, Ying YY, Gan J, Xu HH. Cervical carcinogenesis risk association of HPV33 E6 and E7 genetic variations in Taizhou, Southeast China. Virol J. 2023;20:156.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Yang Z, He ZH, Zhang Y, Di XH, Zheng DF, Xu HH. Genetic variability in the E6 and E7 oncogenes of HPV52 and its prevalence in the Taizhou area. China Virol J. 2022;19:194.

    Article  CAS  PubMed  Google Scholar 

  24. Yu JH, Shi WW, Zhou MY, Liu JM, Han QY, Xu HH. Genetic variability and oncogenic risk association of human papillomavirus type 58 E6 and E7 genes in Taizhou area. China Gene. 2019;686:171–6.

    Article  CAS  PubMed  Google Scholar 

  25. Álvarez-Carretero S, Kapli P, Yang Z. Beginner’s guide on the use of PAML to detect positive selection. Mol Biol Evol. 2023;40:msad041.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, 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:209–49.

    Article  PubMed  Google Scholar 

  27. Wang H, Jiang Y, Wang Q, Lai Y, Holloway A. The status and challenges of HPV vaccine programme in China: an exploration of the related policy obstacles. BMJ Glob Health. 2023;8:e012554.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Rahangdale L, Mungo C, O’Connor S, Chibwesha CJ, Brewer NT. Human papillomavirus vaccination and cervical cancer risk. BMJ. 2022;379:e070115.

    Article  PubMed  Google Scholar 

  29. Ainsworth J, Thomas M, Banks L, Coutlee F, Matlashewski G. Comparison of p53 and the PDZ domain containing protein MAGI-3 regulation by the E6 protein from high-risk human papillomaviruses. Virol J. 2008;5:67.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Ai W, Wu C, Jia L, Xiao X, Xu X, Ren M, et al. Deep sequencing of HPV16 E6 Region reveals unique mutation pattern of HPV16 and predicts cervical cancer. Microbiol Spectr. 2022;10:e0140122.

    Article  PubMed  Google Scholar 

  31. Boon SS, Xia C, Lim JY, Chen Z, Law PTY, Yeung ACM, et al. Human papillomavirus 58 E7 T20I/G63S variant isolated from an East Asian population possesses high Oncogenicity. J Virol. 2020;94:e00090-e120.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Yang L, Yang H, Wu K, Shi X, Ma S, Sun Q. Prevalence of HPV and variation of HPV 16/HPV 18 E6/E7 genes in cervical cancer in women in South West China. J Med Virol. 2014;86:1926–36.

    Article  CAS  PubMed  Google Scholar 

  33. Chen Z, Jing Y, Wen Q, Ding X, Wang T, Mu X, et al. E6 and E7 Gene Polymorphisms in Human Papillomavirus Types-58 and 33 Identified in Southwest China. PLoS One. 2017;12(1):e0171140.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Xi J, Chen J, Xu M, Yang H, Wen S, Pan Y, et al. The polymorphisms of LCR, E6, and E7 of HPV-58 isolates in Yunnan, Southwest China. Virol J. 2018;15(1):76.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Yang L, Yang H, Chen J, Huang X, Pan Y, Li D, et al. Genetic variability of HPV-58 E6 and E7 genes in Southwest China. Infect Genet Evol. 2014;21:395–400.

    Article  PubMed  Google Scholar 

  36. Yang Z, Zhang C, Luo P, Ye M, Gong Q, Mei B. Genetic variability of E6 and E7 genes of human papillomavirus type 58 in Jingzhou, Hubei Province of central China. Virol J. 2022;19(1):71.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Li T, Yang Z, Zhang C, Wang S, Mei B. Genetic variation of E6 and E7 genes of human papillomavirus type 16 from central China. Virol J. 2023;20(1):217.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Li S, Ye M, Chen Y, Gong Q, Mei B. Genetic variation of E6 and E7 genes of human papillomavirus 52 from Central China. J Med Virol. 2021;93(6):3849–56.

    Article  CAS  PubMed  Google Scholar 

  39. Song Z, Cui Y, Li Q, Deng J, Ding X, He J, et al. The genetic variability, phylogeny and functional significance of E6, E7 and LCR in human papillomavirus type 52 isolates in Sichuan, China. Virol J. 2021;18(1):94.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Carvajal-Rodríguez A. Detecting recombination and diversifying selection in human alpha-papillomavirus. Infect Genet Evol. 2008;8(5):689–92.

    Article  PubMed  Google Scholar 

  41. Bletsa G, Zagouri F, Amoutzias GD, Nikolaidis M, Zografos E, Markoulatos P, et al. Genetic variability of the HPV16 early genes and LCR. Present and future perspectives. Expert Rev Mol Med. 2021;23:e19.

    Article  CAS  PubMed  Google Scholar 

  42. He Y, Li J, Mao W, Zhang D, Liu M, Shan X, et al. HLA common and well-documented alleles in China. HLA. 2018;92:199–205.

    Article  CAS  PubMed  Google Scholar 

  43. Yao Y, Huang W, Yang X, Sun W, Liu X, Cun W, et al. HPV-16 E6 and E7 protein T cell epitopes prediction analysis based on distributions of HLA-A loci across populations: an in silico approach. Vaccine. 2013;31:2289–94.

    Article  CAS  PubMed  Google Scholar 

  44. Chenzhang Y, Wen Q, Ding X, Cao M, Chen Z, Mu X, et al. Identification of the impact on T- and B-cell epitopes of human papillomavirus type-16 E6 and E7 variant in Southwest China. Immunol Lett. 2017;181:26–30.

    Article  CAS  PubMed  Google Scholar 

  45. Riemer AB, Keskin DB, Zhang G, Handley M, Anderson KS, Brusic V, et al. A conserved E7-derived cytotoxic T lymphocyte epitope expressed on human papillomavirus 16-transformed HLA-A2+ epithelial cancers. J Biol Chem. 2010;285:29608–22.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Funding

This work was supported by grants from National Natural Science Foundation of China (No.81901625), by Science and Technology Bureau of Taizhou (No.23ywa01), and by Shaoxing University (No.Y20230324).

Author information

Author notes

  1. Haobo Yuan, Ziyi Yan and Jun Gan contributed equally to this work.

Authors and Affiliations

  1. School of Medicine, Shaoxing University, Shaoxing, Zhejiang, 312000, P. R. China

    Haobo Yuan

  2. Medical Research Center, Taizhou Hospital of Zhejiang Province, Wenzhou Medical University, Linhai, Zhejiang, 317000, P. R. China

    Ziyi Yan, Jun Gan, Xinghong Di, Yi Qiu & Huihui Xu

  3. Key Laboratory of Minimally Invasive Techniques & Rapid Rehabilitation of Digestive System Tumor of Zhejiang Province, Linhai, Zhejiang, 317000, P. R. China

    Jun Gan & Huihui Xu

  4. Scientific Research Department, Taizhou Hospital of Zhejiang Province, Wenzhou Medical University, Linhai, Zhejiang, 317000, P. R. China

    Huihui Xu

Authors
  1. Haobo Yuan
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Ziyi Yan
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Jun Gan
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Xinghong Di
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Yi Qiu
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Huihui Xu
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

H.H.X. designed the experiments, performed analysis and drafted the manuscript. X.H.D. and Y.Q. performed HPV genotyping. Z.Y.Y., H.B.Y., and J.G. carried out the sample collection and PCR amplification. All authors have read and approved the final manuscript.

Corresponding author

Correspondence to Huihui Xu.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the Institutional Medical Ethics Review Board of Taizhou Hospital of Zhejiang Province (approval # MERB-2017–020). All patients provided written informed consent before specimen collection, and the patients’ privacy is strictly protected.

Consent for publication

Not applicable.

competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

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

Supplementary Information

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 http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) 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

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Yuan, H., Yan, Z., Gan, J. et al. Phylogenetic analysis and antigenic epitope prediction for E6 and E7 of Alpha-papillomavirus 9 in Taizhou, China. BMC Genomics 25, 507 (2024). https://doi.org/10.1186/s12864-024-10411-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-024-10411-1

Keywords

BMC Genomics

ISSN: 1471-2164

Contact us