Skip to main content

The polymorphism analysis and therapy vaccine target epitopes screening of HPV-35 E6 E7 among the threaten α-9 HPV in Sichuan area

Abstract

High-risk human papilloma virus (HR-HPV) persistent infection is closely associated with the development of cervical cancer and squamous intraepithelial lesion (SIL).The α-9 HPVs, which is predominantly composed of HR-HPV types, account for 75% of HR-HPV infection in Sichuan. The oncoproteins E6 and E7 of HPV play a crucial role in tumor initiation and progression. Notably, HPV-35 is the only HR-HPV type within the α-9 genus that is not included in the nine-valent HPV prophylactic vaccine. Cervical cell samples obtained from Sichuan were collected for HPV detection and genotyping. Among the 406 HPV-positive samples, 31 HPV-35 were detected, 24 HPV-35 E6 and 26 E7 were successfully amplified and sequenced, five nucleotide mutations in E6 and three in E7 were detected, T232C, T434G of E6 (W78R, I145R) and C67T, G84T of E7 (H23Y, L28F) were non-synonymy mutation. PAML 4.8 server was used to detect positive selection sites of HPV-35 E6, E7, and E6 is W78R. Phyre2 were used to predict and analyze protein structures, W78R made influences on protein structure. IEDB were used to screen epitopes vaccine target for HPV-35 affection therapy, and 5 HPV-35 E6 and 3 HPV-35 E7 most potential epitopes were obtained, the most potential peptides for therapy vaccine design were 79-91YRYSVYGETLEKQ, 45-60FACYDLCIVREGQPY, 124-135RFHNIGGRWTGR of E6; 3-19GEITTLQDYVLDLEPEA, 38-47TIDGPAGQAK, 70-88VQSTHIDIRKLEDLLMGTF of E7 and W78R mainly decreased the epitopes affinity.

Conclusions Amino acid substitution in the positive selection sites of HPV-35 E6 and E7 genes have been found to influence protein structure and to decrease the overall affinity of antigen epitopes. This observation aligns with the evolutionary significance of positive selection site, which may confer advantages to the virus by making infected cells more challenging for the immune system to detect, thereby enhancing HPV’s adaptability to the host environment. The polymorphism analysis of HPV-35 E6, E7 contributes to the enrichment of α-9 HPV data in Sichuan China, which is instrumental in improving the effectiveness of clinical detection. Furthermore, these findings provide a relevant theoretical foundation for the prevention and treatment of HPV-related diseases.

Introduction

Cervical cancer, alongside liver cancer, is a type of malignancy associated with viral infection with 99.7% of cases linked to persistent infection with HR-HPVs. The α-9 genus HPVs predominantly comprises mucosal HR-HPVs, representing approximately 75% of HR-HPV infection in Sichuan. This genus includes HPV types 16, 31, 33, 35, 52 and 58. Notably, HPV-35 is the only HR α-9 HPV that has been omitted from the nine-valent prophylactic vaccine and has not been extensively studied. The scarcity of data on HPV-35 hinders our understanding of its characteristics. Since HPV characteristic can vary by geography and ethnicity, it is particularly crucial to conduct research in Sichuan, to provide a theoretical basis for the prevention and treatment of HPV-35 infections.

HR-HPV are known to induce the occurrence and progression of cervical cancer and other malignant tumors primarily through the action of their E6 and E7 oncoproteins [1]. Protein structure determines its function, HPV E6 is consists of one N-terminal (residues 1–36), one C-terminal (residues 147–158) and two Zinc finger (residues 37–73 and 110–146). These Zinc finger motifs are characterized by a consensus CXXC-(29x)-CXXC amino acid sequence and are critical for binding zinc ions, which are essential for various cellular processes, including RNA transcription, transformation, cell immortalization, and protein–protein interactions with host cell proteins. Moreover, the two Zinc finger domains fold into a deep pocket, which, through their common alpha-helical structure, is involved in the ubiquitination and degradation of the tumor suppressor protein p53 [2,3,4,5,6,7]. Distinct from low-risk HPVs (LR-HPVs), the carboxyl terminal of HR-HPV E6 contains a highly conserved and unique amino acid sequence that serves as a PDZ-binding motif. This motif is known to interact with PDZ domain-containing proteins, which are involved in intracellular cell adhesion, cell junction formation and cell polarization. The binding of E6 to these specific proteins can confer a more aggressive phenotype on cells or maintain high levels of protein expression [8, 9]. The PDZ-binding motif (residues 145–149) is a key target for the E6 protein’s involvement in cellular transformation, and the carboxy terminal half, principally responsible for p53 binding, appears to be crucial for the induction of cellular immortality [10]. HPV E7, on the other hand, promotes cell immortalization and inhibits apoptosis in infected cells by binding to the retinoblastoma protein (pRB). The E7 protein is divided into three conserved regions (CR1:1–15, CR2:16–37, CR3:38–98). The CR2 region contains a strictly conserved sequence LXCXE, with amino acids at positions 21–29 comprising the pRB-binding regions [11]. The CR2 region exhibits a high affinity for the RB protein and considered the active region of the E7 protein. Additionally, the LXCXE sequence serves as a phosphorylation site for casein kinase (CK) II [12], further modulating the protein’s interactions and functions within the host cell. The CR3 region of HPV E7 contains two CXXC sequences separated by 29 and 30 amino acids, respectively, which constitute the zinc finger binding domain. The amino acids at positions 58, 61, 91 and 94 of the HPV E7 protein serve as the Zn-binding sites. The zinc finger binding domain is involved in the interaction with a variety of functional proteins [13]. Substitutions within this region have the potential to alter the conformation, and consequently, the function of the E7 protein to some degree.

Human leukocyte antigens (HLAs) are molecules presen on the surface of all human cells that plays an crucial role in distinguishing immune responses by recognizing autoantigen epitope targets and either stimulating or modulating these responses. This mechanism enables the body to resist and eliminate HPV infections. Currently, there is no specific drug available to cure HPV infection, instead, the primary clearance mechanism is through the immune system, with T cell immunity being the key component. In the context of HPV infection-related immune responses, CD8+ cytotoxic T lymphocytes (CTL) and CD4+ helper T lymphocytes (Th) cells are the primary functional cells. The activation and proliferation of CTLs can directly kill tumor cells and secrete cytokines to combat tumors, making them the primary effectors against HPV infection and cervical cancer cells [14, 15]. Activated CD4+ Th cells can produce lymphocytokines that enhance the function of CTLs and natural killer cells (NK cells), activate antigen-presenting cells (APCs), and assist in the elimination of HPV-infected and cervical cancer cells [16]. The E6 and E7 proteins are early proteins, distinct from the targets of preventive vaccines such as L1 and L2, which are expressed in large quantities following infection to activate the immune response. These oncoproteins are considered as important targets for therapeutic vaccines. T-cell epitopes of E6 and E7 have been found to be located in the C-segment and zinc finger region, and HLA-A epitopes may also be present in the N-terminal.

The E6 and E7 oncoproteins of HR-HPVs are critical for the induction of cervical cancer. These proteins not only possess transforming and carcinogenic capabilities but also exhibit trans-activation of viral genes and cellular gene transcription. Upon HPV infection, the immune system recognizes, processes, and presents foreign proteins. Nucleotide non-synonymy mutations can lead to amino acid substitutions in the protein sequence, altering the physical, chemical properties, structure, and immunogenicity of the protein. These changes can affect the immune response to the virus and its pathogenicity. Positive selective mutation can result in a stable increase in gene frequency, enhancing the adaptability of the virus to the environment. This adaptability can lead to changes in infection rates, allowing those variants that are better suited to the natural environment to survive and thrive. Such as specific mutations L83V (L90V), D25E (D32E) in HPV-16 E6 [17]. This adaptability can lead to changes in infection rates, allowing those variants that are better suited to the natural environment to survive and thrive.

Epidemic data indicate that HPV-35 is one of the most prevalent genotypes in cervical cancer cases and is associated with a risk of cervical cancer that is second only to that of HPV-16, making it a significant HR-HPV type. While, there is relevant research data on HPV-35 in North America, Canada, and other places, but lack of relevant research on HPV-35 in China, especially in Sichuan. The types of HPV and the prevalence of specific mutations can vary between different regions and ethnic groups. Consequently, it is imperative to investigate the polymorphism of E6 and E7, positive selection sites, epitopes, and protein structures of HPV-35 in Sichuan. Such research will aid in identifying therapeutic vaccine targets and provide crucial data for the prevention and control of HPV-35 infection-related diseases within this specific region.

Materials and methods

Samples resource

Cervical cell specimens were randomly collected from January 2019 to December 2022 at Sichuan Provincial Maternity and Child Health Care Hospital, the Affiliated Hospital of the Sichuan Reproductive Health Research Center, Chengdu Jinsha hospital, the Infertility Hospital Affiliated to Chengdu Medical College, Chengdu Song zi niao Sterility Hospital. These specimens were stored in a -20℃ antiseptic buffer. Prior to sample collection, written informed consent was obtained from all patients or their legal guardians. The study was ethically approved by the Medical Ethics Committee of the Sichuan Provincial Maternity and Child Health Care Hospital, and stringent measures were implemented to ensure patient privacy and confidentiality.

Genomic DNA extraction and HPV typing

HPV DNA was extracted and genotyped by Human Papillomavirus Genotyping test Kit (Hailes Bio, Ningbo, China), in accordance with the manufacturer’s instructions.

PCR amplification

Primers specific to the HPV-35 E6, E7 gene regions were designed by PRIMER version 5.0 software and validated through the NCBI (National Center for Biotechnology Information) Primer-Blast tool, based on the reference sequences of HPV-35 (GenBank No: HQ537730.1). These primers were subsequently synthesized by TSINGKE (Chengdu, China). The primer sequences and the polymerase chain reaction (PCR) reaction mix composition are detailed in Supplementary Table S1. The PCR products were visualized by gel electrophoresis in a 2% agarose gel (Sangon Biotech Co., Ltd.), and purified PCR products were subjected to bidirectional DNA sequencing by TSINGKE at least twice (Chengdu, China) to ensure accuracy.

HPV-35 E6, E7 polymorphisms, evolutionary and positive selection sites analysis

Nucleotide mutation and the corresponding amino-substitution of sequenced HPV-35 E6 and E7 in Sichuan were detected and analyzed by NCBI BLAST, Primer Premier 5 and DNAMAN5.2.2 (Table 1), difference in protein structure were predicted by PSIPred (http://bioinf.cs.ucl.ac.uk/psipred/).

Table 1 The polymorphism and protein function effect analysis of HPV-35 E6 and E7

Evolutionary analysis of HPV-35 E6, E7 genes were conducted by MEGA 6.0, based on the reference sequence of HPV-35 E6 and E7 (GenBank No.: HQ537730.1-HQ537708.1) [18]. The positive selection sites of HPV-35 E6 and E7 were conducted by Maximum Likelihood 4.8 (PAML 4.8, http://abacus.gene.ucl.ac.uk/software/paml.html) through ù = dN/dS ratio calculation (dN: non-synonymous mutation rate, dS: synonymous mutation rate) [19].

The effect of HPV-35 E6, E7 polymorphism on RNA, protein structure and function

The impact of non-synonymous mutation on protein function was assessed by the Polymorphism Phenotyping v2 (PolyPhen-2, http://genetics.bwh.harvard.edu/pph2) and the Protein Variation Effect Analyzer (PROVEAN, http://provean.jcvi.org) Online software [20, 21].

The impact of non-synonymous mutation on RNA splicing pattern and structure were investigated by Alternative Splice Site Predictor (ASSP) Online software (http://wangcomputing.com/assp/) and RNAsnp software (http://rth.dk/resources/rnasnp/) respectively [22, 23].

The influence of non-synonymous mutation on protein secondary structure was predicted by Phyre2 (http://www.sbg.bio.ic.ac.uk/~phyre2/html/page.cgi?id=index) [24]. Additionally, the impact of these non-synonymous mutation on protein hydrophobicity and stability was evaluated by ExPASy-ProtScale (https://web.expasy.org/protscale) and I-Mutant 2.0 (http://folding.biofold.org/i-mutant/i-mutant2.0.html) respectively [25, 26].

HPV-35 E6, E7 optimal T-cell antigen epitopes selection

Based on the major histocompatibility complex database (dbMHC), HLA alleles with an average frequency of over 5% in the Chinese population were designated as high frequency allele for the selection of T-cell epitopes in Sichuan. This selection was performed by the Immunology Database and Analysis Project (IEDB) resources (http://www.iedb.org/) (Supplementary Table S1). For potential vaccine design targets, epitopes with lower peptides percentile rank (PR) were considered to have better affinity. Therefore, HLA-I epitopes with PR < 1.0 and HLA-II epitopes with PR < 5.0 were chosen for potential vaccine design target [27].

Results

HPV-35 prevalence and E6, E7 polymorphisms in Sichuan

Among the 406 HPV positive cell specimens, 31 HPV-35 positive samples were detected (7.63%), among that 24 E6 and 26 E7 of HPV-35 were succeed sequenced. Compare with HPV-35 reference sequence (HQ533730.1), all E6 and 24 E7 were variants, 5 nucleotide mutations were detected in E6 sequence and 3 in E7, details were shown in Table 1. In E6, T232C, T434G non-synonymous mutations made W78R and I145R, and in E7, C67T, G84T made H23Y and L28F respectively.

Evolutionary and positive selection sites analysis of HPV-35 E6, E7

The phylogenetic tree constructed based on the HPV-35 sample sequences, reference sequences and known sequences is presented in Fig. 1. In this study, no novel mutation emerged, and all sequenced mutant strains were identified. Notably, the sample sequences QV31639 and QV29782 exhibited the highest homology to the reference sequence. Similarly, the sequences HPV35E6/E7 01 and HPV-35 E6/E7 10 demonstrated the highest homology to the reference sequence, suggesting that these mutant strains may have originated from reference sequences. The HPV-35 A1 and A2 line variants previously identified by Chen Zigui et al. were not included in this analysis. However, the HPV E610 and HPV E620 sample were sequenced and are numbered in Fig. 1. These sequences have been deposited in Genbank under the names: HPV-35 E6 01- HPV-35 E6 24, HPV-35 E7 01- HPV-35 E7 24.

Fig.1
figure 1

Evolutionary tree of HPV-35 E6, E7

Calculated by Codeml software using Naive NEB and Bayes Empirical Bayes models, the positive selection site of HPV-35 E6 was W78R.

Non-synonymous mutation effects on alternative splice site, structure of HPV-35 E6, E7 RNA

In the HPV-35 E6 gene, the nucleotide mutation T232C (W78R) made a different in Intron GC, Alt./Cryptic, Constitutive, Confidence of sequence splice site in 220 bp, Score, Alt./Cryptic, Constitutive, Confidence of sequence splice site in 233 bp, 244 bp. T434G (I145R) made a different in Alt./Cryptic, Constitutive, Confidence of sequence splice site in 220 bp, Constitutive, Confidence of sequence splice site in 233 bp, T232C and T434G indicated that c.T272G mutation influences the RNA splice site (Table 2). In HPV-35 E7, G84T (L28F) made 2 new sequence splice sites in 93 bp and 98 bp.

Table 2 Effect of non-synonymous mutation on RNA splicing pattern

By RNAsnp web server, the wild-type, mutant T232C (distance: 0.0527, p-value = 0.4029, 190 to 340 nt highlighted), mutant T434G (distance: 0.0128, p-value = 0.6848, 398 to 447 nt highlighted) are difference in HPV E6 RNA structure, c.T232C (p.W78R) decreased the minimum free energy of RNA from -80.20 kcal/mol kcal/mol to -82.1 kcal/mol; c.U434G (p.W78R) increased the minimum free energy of RNA from -40.20 kcal/mol to -38.50 kcal/mol, T/C and T/G transition results in the RNA secondary structure change (Fig. 1)

In the HPV E7 gene, the wild-type, mutant c.C67T (distance: 0.0406, p-value = 0.5424, 55 to 180 nt highlighted), mutant G84U (distance: 0.0811, p-value = 0.3055, 59 to 108 nt highlighted) are difference in RNA structure, c.C67T (p.H23Y) increased the minimum free energy of RNA from -58.40 kcal/mol to -56.60 kcal/mol; c.G84T (p.L28F) decreased the minimum free energy of RNA from -40.20 kcal/mol to -38.50 kcal/ mol, the T/C and G/T transition results in the RNA secondary structure change (Fig. 2).

Fig. 2
figure 2

The effect of HPV-35 E6, E7 polymorphism on RNA structure

The effects of non-synonymous mutation on HPV-35 E6, E7 protein structure, structural stability, and hydrophilicity

In the E6 protein, the mutations W78R, I145R were not situated within the alpha-helix or beta-sheet. The W78R resulted in an increase in the lowest free energy of the mutant protein increased from -80.20 to -82.10 kcal/mol, indicating enhanced hydrophilicity, but decreased protein structural stability. Similarly, the mutation I145R led to an increase in the lowest free energy from -40.20 to -38.50 kcal/mol, again suggesting enhanced hydrophilicity and reduced protein structural stability (Fig. 3 and Table 1).

Fig. 3
figure 3

The effect of non-synonymous mutation on HPV-35 E6, E7 hydrophobicity

In the E7 protein, the mutation H23Y caused the disappearance of the alpha-helix located between the 15th and 89th amino acids, which is remote from the site of substitution. This mutation increased the lowest free energy of the mutant protein from -58.40 to -56.60 kcal/mol, indicating decreased hydrophilicity and increased protein structural stability. Conversely, the mutation L28F resulted in a decrease in the lowest free energy of the mutant protein from -70.00 to -71.10 kcal/mol, suggesting enhanced hydrophilicity and decreased protein structural stability (Fig. 3 and Table 1).

T-cell antigen epitopes of HPV-35 E6, E7

In HPV-35 E6 reference sequence, 92 HLA-I and 84 HLA-II epitopes were predicted. In W78R variant, 90 HLA-I and 89 HLA-II epitopes were predicted, W78R made 5 HLA-I, 19 HLA-II new epitopes appear, 7 HLA-I, 14 HLA-II epitopes disappear, and 5 HLA-I, 41 HLA-II epitopes affinity increase; amino-acid substitution made HLA-Iepitopes affinity, number decrease and HLA-II epitopes affinity, number increase. In I145R variant, 93 HLA-I and 84 HLA-II epitopes were predicted, I145R made a new HLA-I epitopes appear. Details were shown in Tables 3 and 4.

Table 3 The effect of HPV-35 E6, E7 polymorphism on HLA-I epitopes
Table 4 The effect of HPV-35 E6, E7 polymorphism on HLA-II epitopes

In HPV-35 E7 reference sequence, 27 HLA-I and 147 HLA-II epitopes were predicted. In H23Y variant, 28 HLA-I and 147 HLA-II epitopes were predicted, H23Y made a new HLA-I epitopes. In L28F variant, 27 HLA-I and 147 HLA-II epitopes were predicted, L28F have no influence on HLA-I and HLA-II epitopes (Tables 3 and 4) (Fig. 4).

Fig. 4
figure 4

The effect of non-synonymous mutation on HPV-35 E6, E7 protein structure

For the design of therapeutic vaccine, epitopes devoid of mutation sites were concatenated into continuous polypeptide sequences. This approach aims to preserve the native structure and immune recognition properties of the epitopes, thereby enhancing the potential efficacy of the vaccine. In HPV-35 E6, HLA-I predicted epitopes occurred in 1–19, 17–32, 35–45, 41–55, 52–64, 67–77, 79–96, 124–140 segment, and HLA-II were in 45–62, 79–91, 120–135; 81-88YSVYGETL and 52-60IVYREGQPY were the most frequency E6 HLA-I epitopes, as well as 15 potential HLA-I epitopes located in HLA-I and HLA-II common segment, details were shown in Tables 5 and 6.

Table 5 The HLA-I, HLA-II distribution analysis of HPV-35 E6 and E7
Table 6 The T-cell optimal epitopes of HPV-35 E6, E7

In HPV-35 E7, HLA-I predicted epitopes occurred in 3–19, 38–55, 70–80, 79–88, 90–98 segment, and HLA-II were in 1–22, 30–49, 70–91; 45-53QAKPDTSNY was the lowest PR (0.1) and most frequently E7 HLA-I epitopes, as well as 12 HLA-I epitopes located in HLA-I and HLA-II common segment, details were shown in Tables 5 and 6.

The HLA-II epitope binding core is a key important region for CD4+ T-cell recognition, binding, and activation, residing in the peptide-binding groove of HLA-II molecules. Both CD4+ and CD8+ T-cell participate in HPV-related immune response, CD8+ T-cell are considered as the primary effectors in eliminating HPV-infected and cervical cancer cells. Consequently, optimal epitopes that are the potential MHC-I epitopes that contain at least one core region of HLA-II-binding peptide is valuable. These epitopes can potentially activate both CD4+ and CD8+ T-cells, thereby enhancing the breadth and effectiveness of the immune response against HPV.

Finally, the optimal epitopes were selected as fellow step: 1. Selected the epitopes do not contain mutation sites. 2. Selected epitopes located in the common region of HLA-I and HLA-II predicted epitopes. 3. Selected epitopes contain at least one HLA-II binding core region peptide. 4. Epitopes that bind more HLA-I alleles and lower the percentile rank is better. 5. The integrated HLA-I epitope peptides are also excellent vaccine targets peptides in HLA-I and HLA-II predicted epitopes common region and called optimal peptides. In HPV-35 E6, 5 optimal epitopes were selected, they were 52-60IVYREGQPY, 79-88YRYSVYGETL, 79-90YRYSVYGETLEK, 80-90RYSVYGETLEK, 81-90YSVYGETLEK; the optimal peptides were 79-91YRYSVYGETLEKQ, 45-60FACYDLCIVYREGQPY, 124-135RFHNIGGRWTGR. In HPV-35 E7, optimal epitopes were 3-13GEITTLQDYVL, 5-13ITTLQDYVL, 11-19YVLDLEPEA; the optimal peptides: 3-19GEITTLQDYVLDLEPEA, 38-47TIDGPAGQA, 70-88VQSTHIDRKLEDLLMGT.

Discussion

The paucity of data on HPV-35 hinders our understanding of this virus, highlighting the significance of conducting research in Sichuan to establish a theoretical framework for the prevention and treatment of HPV35. The HPV characteristics vary across different geographical and ethnic groups, highlighting the significance of conducting research in Sichuan to establish a theoretical framework for the prevention and treatment of HPV35. HPV-35 is the only HR α-9 genus HPV that is understudied, largely due to its historically low prevalence. International research has identified HPV-35 as a significant risk factor for cervical cancer, ranking it second only to HPV16, in terms of pathogenicity, with a threat to human health that is comparable to that of HPV16. The limited data on HPV35, coupled with its low infection rate, can lead to a tendency to overlook its risks, yet latent HPV-35 may pose an even greater threat than HPV16. In our study, only 31 HPV35-positive samples were detected among 406 positive samples (7.67%), which aligns with the 8.1% positive detection data reported by Basto Diogo Lisboa et al. The low infection rate of HPV-35 reflects the scarcity of genomic data, which may be attributed to the highly conservative nature of HPV-35 mutations, rendering them less susceptible to genomic mutations and inheritance within host cells. Additionally, HPV-35 may be less likely to cross-infect with other HPV types or variants, or it may have a longer incubation period in host cells than other HPV types, making it more difficult to detect [28]. The decades-long incubation period of HPV-35 compounds our challenges, and the lack of comprehensive data on HPV-35 impedes our awareness of its implications. Given the distinct characteristics of HPV across different geographical and ethnic groups, it is imperative to conduct systematic research in Sichuan to provide a solid theoretical foundation for HPV-35 prevention and treatment.

HPV-35 isolates are genetically related to HPV31 and HPV16, and they shared a recent common ancestor with HPV33 and HPV58 isolates [28]. In our study, 5 nucleotide-acid mutations in E6 (2 non-synonymous mutations: T232C (W78R) and T434G (I146R), 3 synonymous mutations: T18C, T27C, and G262C) were detected and 3 in E7 (2 non-synonymous mutations: C67T (H23Y) and G84T (L28F), 1 synonymous mutation: C114T), and due to the limitation of sample, some mutation sites may not have been identified. The variants of HPV-35 E6 and E7 in Sichuan, China, are distinct and do not belong to the well-known A1 and A2 series. They exhibit strong regional characteristics. The homology between these variants and the reference sequences is notably higher for the E6 variant compared to the E7 variant. Furthermore, the E7 variant appears to have a more distant evolutionary lineage. HPV 35E6/E7 01, HPV 35 E6/E7 10 was found to have the highest degree of homology with the reference sequence among the HPV-35 E6/E7 sample sequence. This isolate in question likely originated from a reference sequence, with HPV-35 E6/E7 01, 35 E6/E7 10 and the A1 line belonging to the same clade, while the remaining sequences also grouped within the same clade. As observed by Chen Zigui and Basto Diogo Lisboa et al., all variants of HPV-35 are highly conserved, displaying minimal intratype genetic diversity, and failing to meet the criteria for classification into multiple lineages. The highest infection rate for HPV-35 was found in the 56-years age bracket (42.9%) [29]. These finding suggest the prevalence pattern may be influenced by the existence of another variant lineage of HPV-35 that could be present in isolated and/or unsampled populations. Intratype mutation within host cell may be infrequent, and the virus may have a longer incubation period, making it more challenging to detect. The mutant strains in the Sichuan area of China exhibit distinct regional characteristics, which should be taken into account, when designing vaccines and probes specific to the Sichuan population. The mutation rate in the E6 gene was higher than that of E7 (78.94%/52.63%). Consequently, the more stable E7 should be prioritized, when considering the design of probe for HPV35.

HR-HPV E6, E7 are critical oncoprotein that differ from those of LR-HPV. They are capable of interacting with RB, p53 tumor suppressor proteins and several cycle regulatory proteins. E6, E7 non-synonymous mutations can result in alterations to the amino acid composition, protein structure and function. For instance, the mutations W78R, I145R in E6 do not reside within the alpha-helix or beta-sheet region. However, I145R is located in the highly conserved carboxyl terminal, which is involved in binding to PDZ domain-containing proteins that play roles in intracellular cell adhesion, cell linkage and cell polarization. W78R is situated near the Znic-figure and within the groove formed by two Znic-blinding figures, which are crucial for p53 binding. In E7, the mutations H23Y, L28F are not located within the alpha-helix or beta-sheet regions either. The region spanning amino acids 16–37 of E7 exhibits a high affinity for RB protein and is considered the active region of E7 protein. The residues 21–29 positions constitute the pRB binding regions, with the LXCXE motif being strictly conserved as the phosphorylation site of casein kinase (CK) II. Both H23Y and L28F are positioned within this functional region of E7. Additionally, H23Y induces a change from an alpha-helix to a coil conformation at positions 15 and 89, which are distant from the site of substitution. The structural change at the 89th amino acid caused by H23Y is in close proximity to amino acids 58, 61, 91, and 94 of HPV E7, which are zinc-binding sites involved in the binding of various functional proteins. The Znic-binding figures, carboxyl terminal and PDZ domain binding motif are essential structural and functional domains of E6, E7. Substitutions in these regions can lead to differences in the binding abilities of these oncoproteins to host proteins such as p53, pRB, and others, which is considered key factors in their carcinogenic potential.

HPV E6, E7, which are early protein, are considered prime candidates for prophylactic vaccine targets. When these proteins are digested and presented to the body's immune system, they can active immune response that eliminate HPV infection and reduces the risk of HPV-related diseases. E6 and E7 T-cell epitopes have been found to locate in the C-terminal segment and Zinc-finger region, with HLA-A epitopes potentially also be in the N-terminal. Research has indicated that mutations can impact the effectiveness of HPV vaccine, particularly those designed to target specific epitopes, as the efficacy may be compromised if mutations are not considered. Non-synonymous mutation can alter epitope characteristics. For instance, the W78R mutation in E6 has been shown to increase the affinity and number of HLA-I epitopes, potentially leading to the emergence of new and improved epitopes. Similarly, the I145R mutation enhances HLA-I epitope affinity and number. H23Y in E7 increases HLA-I epitope affinity as well. The E7 protein exhibits a greater number and affinity of epitopes compared to E6, suggesting that E7 may not be as readily recognized by host cells, and its carcinogenic potential may exceed that of E6. These hypotheses require confirmation in subsequent studies. Positive selection sites play a crucial role in enhancing a species’ adaptability to its environment by increasing the frequency of mutated genes [30]. Analysis using the software package paml has revealed that the HPV-35 E6 positive selection site W78R (mutation risk:22/24) is a frequent non-synonymous mutation, indicating that W78R contributes to the adaptation of HPV-35 E6 and has become widely prevalent. The conformational and functional changes caused by this substitution may influence the virus’s pathogenicity [31]. Specific intramotypic mutations in HPV E6 can result in differences in virus infection ability and pathogenicity. For example, positive selection site D32E of HPV16, has been confirmed to be associated with the development of cervical cancer [17]; Similarly, W78R in HPV-35 E6, located in an active region, affects protein conformation, function and may reduce the immunogenicity of the peptide containing this sites. Substitutions at positive selection sites can impact protein structure, decrease epitope affinity, make HPV-infected cells more difficult to detect by the immune system, enhance HPV’s adaptability to the environment, improve the efficiency of virus infection, and potentially promote the development of cervical cancer.

This study represents the first investigation into the interplay between the protein structure, positive selection sites, epitopes, and pathogenicity of HPV-35 E6 and E7 in Sichuan Province. It aims to elucidate the relationship between HPV E6, E7 and the varying abilities of HPV to infect, with the ultimate goal of enhancing the development of therapeutic vaccine targeting HPV-35 E6 and E7 in Sichuan Province.

Conclusion

HPV-35 variants exhibit high conservation and can remain latent in host cells for extended periods. Consequently, the design of detection probes for HPV-35 should be tailored to enhance the detection rate. Early detection and treatment of HPV-35 infections carry significant clinical and social importance. The enrichment of HPV-35 polymorphism data will provide a theoretical foundation for the development of effective prevention and treatment strategies.

Among the 406 HPV-positive cell specimens, 31 HPV-35 positive samples were detected (7.63%), 2 non-synonymous mutations were detected in E6, E7 respectively,. These mutations are in active region, and have the potential to affect protein conformation, function and immunogenicity of those peptide containing these sites. The potential epitopes of HPV-35 E6 were: 52-60IVYREGQPY, 79-88YRYSVY GETL, 79-90YRYSVYGETLEK, 80-90RYSVYGETLEK, 81-90YSVYGETLEK, optimal peptides were 79–91 YRYSVYGETLEKQ, 45-60FACYDLCIVYREGQPY, 124-135RFHNIGGRWTGR; potential epitopes of HPV-35 E7 were: 3-13GEITTL QDYVL, 5-13ITTLQDYVL, 11-19YVLDLEPEA, optimal peptides were 3-19GEITT LQDYVLDLEPEA, 38-47TIDGPAGQA, 70-88VQSTHIDRKLEDLLMGT.

Availability of data and materials

No datasets were generated or analysed during the current study.

References

  1. Chen Z, Schiffman M, Herrero R, DeSalle R, Anastos K, Segondy M, et al. Evolution and taxonomic classification of Alphapapillomavirus 7 complete genomes: HPV18, HPV39, HPV45, HPV59, HPV68 and HPV70. PLoS ONE. 2013;8:1–16.

    CAS  Google Scholar 

  2. Nominé Y, Charbonnier S, Ristriani T, Stier G, Masson M, Cavusoglu N, et al. Domain substructure of HPV E6 oncoprotein: Biophysical characterization of the E6 C-terminal DNA-binding domain. Biochemistry. 2003;42:4909–17.

    Article  PubMed  Google Scholar 

  3. Martinez-Zapien D, Ruiz FX, Poirson J, Mitschler A, Ramirez J, Forster A, et al. Structure of the E6/E6AP/p53 complex required for HPV-mediated degradation of p53 Nature. Nature. 2016;529:541–5 Publishing Group.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Chemes LB, Camporeale G, Sánchez IE, De Prat-Gay G, Alonso LG. Cysteine-rich positions outside the structural zinc motif of human papillomavirus E7 provide conformational modulation and suggest functional redox roles. Biochemistry. 2014;53:1680–96.

    Article  CAS  PubMed  Google Scholar 

  5. Li YL, Ma ZL, Zhao Y, Zang J. Immunization with mutant HPV16 E7 protein inhibits the growth of TC-1 cells in tumor-bearing mice. Oncol Lett. 2015;9:1851–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Huibregtse JM, Beaudenon SL. Mechanism of HPV E6 proteins in cellular transformation. Semin Cancer Biol. 1996;7:317–26.

    Article  CAS  PubMed  Google Scholar 

  7. Alvarez-Salas L M , Cullinan A E , Siwkowski A ,et al.inhibition of hpv-16 e6e7 immortalization of normal keratinocytes by hairpin ribozymes[J]. 2018. https://doi.org/10.1073/pnas.95.3.1189.

  8. Ristriani T, Masson M, Nominé Y, Laurent C, Lefèvre JF, Weiss É, et al. HPV oncoprotein E6 is a structure-dependent DNA-binding protein that recognizes four-way junctions. J Mol Biol. 2000;296:1189–203.

    Article  CAS  PubMed  Google Scholar 

  9. Ganti K, Broniarczyk J, Manoubi W, Massimi P, Mittal S, Pim D, et al. The human Papillomavirus E6 PDZ binding motif: from life cycle to malignancy. Viruses. 2015;7:3530–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Yoshimatsu Y, Nakahara T, Tanaka K, Inagawa Y, Narisawa-saito M, Yugawa T, et al. Roles of the PDZ-binding motif of HPV 16 E6 protein in oncogenic transformation of human cervical keratinocytes. Cancer Sci. 2017;108:7–13.

    Article  Google Scholar 

  11. Robles-Rodríguez OA, Pérez-Trujillo JJ, Barrón-Cantú JA, Torres-Cerda A, Gutiérrez-Puente Y, García-García A, et al. Antitumor effect of adenoviruses expressing mutant non-oncogenic E7 versions from HPV-16 fused to calreticulin. J BUON. 2020;25:543–8.

    Google Scholar 

  12. Jing Y, Wang T, Chen Z, Ding X, Xu J, Mu X, et al. Phylogeny and polymorphism in the long control regions E6, E7, and L1 of HPV Type 56 in women from southwest China. Mol Med Rep. 2018;1–11. Available from: http://www.spandidos-publications.com/https://doi.org/10.3892/mmr.2018.8743

  13. Ding W, Hu Z, Zhu D, Jiang X, Yu L, Wang X, et al. Zinc finger nucleases targeting the human Papillomavirus E7 oncogene induce E7 disruption and a transformed phenotype in HPV16 / 18-positive cervical cancer cells. Clin Cancer Cells. 2014;20:6495–503.

    Article  CAS  Google Scholar 

  14. Singh AK, Rath SK, Misra K. Identification of epitopes in Indian human papilloma virus 16 E6: A bioinformatics approach. J Virol Methods. 2011;177:26–30. https://doi.org/10.1016/j.jviromet.2011.06.006. Elsevier B.V.

    Article  CAS  PubMed  Google Scholar 

  15. Roden RBS, Ling M, Wu TC. Vaccination to prevent and treat cervical cancer. Hum Pathol. 2004;35:971–82.

    Article  PubMed  Google Scholar 

  16. Ferns DM, Heeren AM, Samuels S, Bleeker MCG, de Gruijl TD, Kenter GG, et al. Classical and non-classical HLA class I aberrations in primary cervical squamous- and adenocarcinomas and paired lymph node metastases. J Immunother Cancer. 2016;4:78.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Simas B, Vinicius M, Batista DA, Crovella S, Pavla A, Diniz A, et al. Infection, genetics and evolution novel E6 and E7 oncogenes variants of human papillomavirus type 31 in Brazilian women with abnormal cervical cytology. Infect Genet Evol. 2013;16:13–8.

    Article  Google Scholar 

  18. Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: Molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30:2725–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Yang Z. PAML : Phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2014;24:1586–91.

    Article  Google Scholar 

  20. Adzhubei I, Jordan DM, Sunyaev SR. Predicting functional effect of human missense mutations using PolyPhen-2. Curr Protoc Hum Genet. 2013;Chapter 7:Unit 7.20.

    Google Scholar 

  21. Choi Y, Chan AP. PROVEAN web server: A tool to predict the functional effect of amino acid substitutions and indels. Bioinformatics. 2015;31:2745–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Wang M, Marín A. Characterization and prediction of alternative splice sites. Gene. 2006;366:219–27.

    Article  CAS  PubMed  Google Scholar 

  23. Sabarinathan R, Tafer H, Seemann SE, Hofacker IL, Stadler PF, Gorodkin J. RNAsnp: efficient detection of local RNA secondary structure changes induced by SNPs. Hum Mutat. 2013;34:546–56.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Kelley LA, Mezulis S, Yates CM, Wass MN, Sternberg MJE. The Phyre2 web portal for protein modeling, prediction and analysis. Nat Protoc. 2015;10:845–58 Nature Publishing Group.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Wilkins MR, Gasteiger E, Bairoch A, Sanchez JC, Williams KL, Appel RD, et al. Protein identification and analysis tools in the ExPASy server. Methods Mol Biol. 1999;112:31.

    Google Scholar 

  26. Capriotti E, Fariselli P, Casadio R. I-Mutant2.0: Predicting stability changes upon mutation from the protein sequence or structure. Nucleic Acids Res. 2005;33:306–10.

    Article  Google Scholar 

  27. Bui HH, Sidney J, Li W, Fusseder N, Sette A. Development of an epitope conservancy analysis tool to facilitate the design of epitope-based diagnostics and vaccines. BMC Bioinformatics. 2007;8:1–6.

    Article  Google Scholar 

  28. 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 

  29. Basto D L , Vidal J P ,Pontes, Valéria Barbosa,et al.Genetic diversity of human papillomavirus types 35, 45 and 58 in cervical cancer in Brazil[J].Archives of Virology, 2017. https://doi.org/10.1007/s00705-017-3439-5.

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

    Article  CAS  PubMed  Google Scholar 

  31. Chen AA, Heideman DAM, Boon D, Chen Z, Burk RD, De Vuyst H, et al. Human papillomavirus 33 worldwide genetic variation and associated risk of cervical cancer. Virology. 2014;448:356.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

The authors would like to thank the following hospitals for the sample collection: Sichuan Provincial Maternity and Child Health Care Hospital, Sichuan Reproductive Health Research Center Affiliated Hospital, Chengdu Song ziniao Sterility Hospital, Infertility Hospital Affiliated to Chengdu Medical College, Angel Women's and Children's Hospital, and Chengdu Jinsha hospital.

Funding

No funding.

Author information

Authors and Affiliations

Authors

Contributions

J.H., T.J., conceived and designed the study; C.C., N.L., D.L. and P.G. collected the samples; J.H., C.C., N.L., R.L., D.L. and P.G., performed the experiments. J.H., T.J., and X.D. wrote the manuscript. All authors reviewed the manuscript.

Corresponding author

Correspondence to Xianping Ding.

Ethics declarations

Ethics approval and consent to participate

The present study was approved by Sichuan Provincial Maternity and Child Health Care Hospital Medical Ethics Committee (Chengdu, China; approval number: 20201015–84). Before sample collection, a written informed consent was obtained from the patients or their guardians, and the privacy of patient/study was protected carefully.

Consent to publication

All patients participating in this study have consent for publication.

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

Additional file 1:

Table S1. The selected HLA-I and HLA-IIalleles with average frequency over 5% in Chinese. Table S2. The HLA-I predicted epitopes of HPV-35 E6. Table S3. The HLA-I predicted epitopes of HPV-35 E7. Table S4. The HLA-II predicted epitopes of HPV-35 E6. Table S5. The HLA-II predicted epitopes of HPV-35 E7.

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/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

He, J., Li, T., Cheng, C. et al. The polymorphism analysis and therapy vaccine target epitopes screening of HPV-35 E6 E7 among the threaten α-9 HPV in Sichuan area. Virol J 21, 213 (2024). https://doi.org/10.1186/s12985-024-02357-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12985-024-02357-3

Keywords