Identification of human papillomavirus-16 E6 variation in cervical cancer and their impact on T and B cell epitopes

A Comprehensive in Silico Analysis for Identification of Immunotherapeutic Epitopes of HPV-18

Human papillomavirus (HPV) remains the major cause of cervical cancer, globally. High risk HPV (Hr-HPV) 16 and 18 together account for more than 70% of cervical cancer cases, whereas the hr-HPV-18 is the second most prevalent hr-HPV type, causing about 5.2% of all cancers worldwide. Considering the high prevalence and mortality rate, cervical cancer remains a noteworthy health problem among women. As of now, no registered immunotherapies are available after the HPV infection. Thus, developing an immunotherapeutic candidate against hr-HPV would be of major clinical benefit. Nowadays, the T-and B-cell peptide based targeted vaccines have been considered as the best candidate for vaccine development against viral infections. In this study, both prophylactic and therapeutic vaccine candidates against hr-HPV-18 were predicted. To achieve this, the prediction of T-and B-cell epitopes of major histocompatibility complex (MHC) were accomplished, that can be used for HPV immunotherapy. For MHC-I, a maximum number (20) of potent peptides were found, against HLA-B*51:01 (L1 = 9, L2 = 6, E2 = 4, and E4 = 1) having percentile value < 1 and, immunogenicity scores higher than 0.5, followed by HLA-A*11:01 (L1 = 8, E2 = 7 L2 = 2, and E6 = 1, E7 = 1); 19 epitopes. For MHC-II, the highest number of peptides found, against the HLA-DRB1*04:01 (L2 = 10, E5 = 7, and E4 = 4), HLA-DRB1*04:05 (E5 = 7, E2 = 5, E4 = 5, and L1 = 4) HLA-DPA1*01:03/DPB1*04:01 (E7 = 7, E6 = 5, L2 = 5, and E2 = 2), HLA-DRB5*01:01(E6 = 6, L1 = 6, and L2 = 6); peptides 21, 21, 19 and 18 respectively. For B-cell, total 94, 16 amino acid long B-cell epitopes were predicted. In conclusion, these predicted epitopes can be valuable candidates for in vitro or in vivo therapeutic vaccine studies against hr-HPV-18 associated cancer.

The genetic variability, phylogeny and functional significance of E6, E7 and LCR in human papillomavirus type 52 isolates in Sichuan, China

Background

Variations in human papillomavirus (HPV) E6 and E7 have been shown to be closely related to the persistence of the virus and the occurrence and development of cervical cancer. Long control region (LCR) of HPV has been shown multiple functions on regulating viral transcription. In recent years, there have been reports on E6/E7/LCR of HPV-16 and HPV-58, but there are few studies on HPV-52, especially for LCR. In this study, we focused on gene polymorphism of the HPV-52 E6/E7/LCR sequences, assessed the effects of variations on the immune recognition of viral E6 and E7 antigens, predicted the effect of LCR variations on transcription factor binding sites and provided more basic date for further study of E6/E7/LCR in Chengdu, China.

Methods

LCR/E6/E7 of the HPV-52 were amplified and sequenced to do polymorphic and phylogenetic analysis. Sequences were aligned with the reference sequence by MEGA 7.0 to identify SNP. A neighbor-joining phylogenetic tree was constructed by MEGA 7.0, followed by the secondary structure prediction of the related proteins using PSIPRED 4.0. The selection pressure of E6 and E7 coding regions were estimated by Bayes empirical Bayes analysis of PAML 4.9. The HLA class-I and II binding peptides were predicted by the Immune Epitope Database server. The B cell epitopes were predicted by ABCpred server. Transcription factor binding sites in LCR were predicted by JASPAR database.

Results

50 SNP sites (6 in E6, 10 in E7, 34 in LCR) were found. From the most variable to the least variable, the nucleotide variations were LCR > E7 > E6. Two deletions were found between the nucleotide sites 7387–7391 (TTATG) and 7698–7700 (CTT) in all samples. A deletion was found between the nucleotide sites 7287–7288 (TG) in 97.56% (40/41) of the samples. The combinations of all the SNP sites and deletions resulted in 12 unique sequences. As shown in the neighbor-joining phylogenetic tree, except for one belonging to sub-lineage C2, others sequences clustered into sub-lineage B2. No positive selection was observed in E6 and E7. 8 non-synonymous amino acid substitutions (including E3Q and K93R in the E6, and T37I, S52D, Y59D, H61Y, D64N and L99R in the E7) were potential affecting multiple putative epitopes for both CD4+ and CD8+ T-cells and B-cells. A7168G was the most variable site (100%) and the binding sites for transcription factor VAX1 in LCR. In addition, the prediction results showed that LCR had the high probability binding sites for transcription factors SOX9, FOS, RAX, HOXA5, VAX1 and SRY.

Conclusion

This study provides basic data for understanding the relation among E6/E7/LCR mutations, lineages and carcinogenesis. Furthermore, it provides an insight into the intrinsic geographical relatedness and biological differences of the HPV-52 variants, and contributes to further research on the HPV-52 therapeutic vaccine development.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12985-021-01565-5.

Identification of variants and therapeutic epitopes in HPV-33/HPV-58 E6 and E7 in Southwest China

BACKGROUND

Human papillomavirus (HPV) E6 and E7 oncoproteins play a crucial role in HPV-related diseases, such as cervical cancer, and can be used as ideal targets for therapeutic vaccines. Human leukocyte antigen (HLA) participates in the immune response to block HPV infection and invasion by its target/recognition function. HPV-33 and HPV-58 are highly prevalent among Chinese women. Therefore, it is of great significance to study the E6 and E7 region-specific gene polymorphisms of HPV-33 and HPV-58 in Southwest China and to identify ideal epitopes for vaccine design. Both HPV-33 and HPV-58 belong to α-9 genus HPV and are highly homologous, so their correlations are included in our research.

METHODS

To study the E6 and E7 variations and polymorphisms of HPV-33 and HPV-58 in Southwest China, we collected samples, extracted and sequenced DNA, and identified variants. Nucleotide sequences were translated into amino acids by Mega 6.0 software. The physical/chemical properties, amino acid-conserved sequences and secondary structure of protein sequences were analysed by the Protparam server, ConSurf server and PSIPRED software. The T and B cell epitopes of the E6/E7 reference and variant sequences in HPV-33 and HPV-58 were predicted by the Immune Epitope Database (IEDB) analysis server and the ABCpred server, respectively.

RESULTS

Five and seven optimal HLA-I restricted T cell epitopes were selected from HPV-33 and HPV-58 E6, respectively, and these optimal epitopes are mainly located in _41-58EVYDFAFADLTVVYREGN of HPV-33 E6 and _40-60SEVYDFVFADLRIVYRDGNPF of HPV-58 E6. Six optimal HLA-I-restricted T cell epitopes were selected from HPV-33 and HPV-58 E7, and these epitopes are mainly located in _77-90RTIQQLLMGTVNIV of HPV-33 E7 and _78-91RTLQQLLMGTCTIV of HPV-58 E7.

CONCLUSIONS

HPV-33/HPV-58 E6/E7 gene polymorphisms and T/B cell epitopes of their reference and variant sequences were studied, and candidate epitopes were selected by bioinformatics techniques for therapeutic vaccine design for people in Southwest China. This study was the first to investigate the correlation of epitopes between HPV-33 and HPV-58. After experimental validation, these selected epitopes will be employed to induce a wide range of immune responses in heterogeneous HLA populations.

Structural and Epitope Analysis (T- and B-Cell Epitopes) of Hepatitis C Virus (HCV) Glycoproteins: An in silico Approach

BACKGROUND

Chronic infection with Hepatitis C Virus (HCV) poses a major risk for liver disease like cirrhosis, liver failure and hepatocellular carcinoma. In terms of percentage, the prevalence of HCV in India was found to be low to moderate (1-1.5%), but in terms of sheer numbers, India has a significant number of global HCV patients. Presently, HCV can be treated with direct acting-antibody drugs but there is no prophylactic or therapeutic vaccine available against it. In HCV infection, T- and B-cell immunity is important for clearing the virus. In the present study immunoinformatics was used to identify potent vaccine target for HCV vaccine development.

METHODS

Sequence of HCV was retrieved from NCBI and their structural analysis was done by using Protpram, PSIPRED, iTASSER and PDBsum servers. T-cell and B-cell epitopes were predicted by Immune Epitope Database and ACBPRED servers.

RESULTS

On epitope prediction, 25 and 55 potent MHC-I epitopes and 7 and 13 potent B-cell epitopes were predicted for E1 and E2 protein respectively. Their antigenicity score was also calculated. The most potent MHC-I epitopes were MMMNWSPAV and MAWDMMMNW for HLA-A*02:01 and HLA-B*53:01 and most potent B-cell epitope was TGHRMAWDMMMNWSPA for E1 protein. For E2, four MHC-I epitopes having the lowest binding energy and most potent B-cell epitope was DRPYCWHYAPRPCDTI.

CONCLUSION

In the present study, most potent epitopes for HCV was determined on the basis of their antigenicity along with 3D modeling and docking. Identified B- and T-cell epitopes can be used for the development of potent vaccine against most prevalent HCV type in India to limit its infection.

An unusually high substitution rate in transplant-associated BK polyomavirus in vivo is further concentrated in HLA-C-bound viral peptides

Infection with human BK polyomavirus, a small double-stranded DNA virus, potentially results in severe complications in immunocompromised patients. Here, we describe the in vivo variability and evolution of the BK polyomavirus by deep sequencing. Our data reveal the highest genomic evolutionary rate described in double-stranded DNA viruses, i.e., 10⁻³–10⁻⁵ substitutions per nucleotide site per year. High mutation rates in viruses allow their escape from immune surveillance and adaptation to new hosts. By combining mutational landscapes across viral genomes with in silico prediction of viral peptides, we demonstrate the presence of significantly more coding substitutions within predicted cognate HLA-C-bound viral peptides than outside. This finding suggests a role for HLA-C in antiviral immunity, perhaps through the action of killer cell immunoglobulin-like receptors. The present study provides a comprehensive view of viral evolution and immune escape in a DNA virus.

Little is known about the mechanisms of evolution and viral immune escape in double-stranded DNA (dsDNA) viruses. Here, we study the evolution of BK polyomavirus and observe the highest genomic evolutionary rate described so far for a dsDNA virus, in the range of RNA viruses, which usually evolve rapidly. Furthermore, the prediction of viral peptides to determine immune escape suggests a specific role of HLA-C in antiviral immunity. These findings are helpful for future advances in antiviral therapies and provide a step forward in our understanding of in vivo viral evolution in humans.

The current state of therapeutic and T cell-based vaccines against human papillomaviruses

Human papillomavirus (HPV) is known to be a necessary factor for many gynecologic malignancies and is also associated with a subset of head and neck malignancies. This knowledge has created the opportunity to control these HPV-associated cancers through vaccination. However, despite the availability of prophylactic HPV vaccines, HPV infections remain extremely common worldwide. In addition, while prophylactic HPV vaccines have been effective in preventing infection, they are ineffective at clearing pre-existing HPV infections. Thus, there is an urgent need for therapeutic and T cell-based vaccines to treat existing HPV infections and HPV-associated lesions and cancers. Unlike prophylactic vaccines, which generate neutralizing antibodies, therapeutic, and T cell-based vaccines enhance cell-mediated immunity against HPV antigens. Our review will cover various therapeutic and T cell-based vaccines in development for the treatment of HPV-associated diseases. Furthermore, we review the strategies to enhance the efficacy of therapeutic vaccines and the latest clinical trials on therapeutic and T cell-based HPV vaccines.

Perspectives for therapeutic HPV vaccine development

Background

Human papillomavirus (HPV) infections and associated diseases remain a serious burden worldwide. It is now clear that HPV serves as the etiological factor and biologic carcinogen for HPV-associated lesions and cancers. Although preventative HPV vaccines are available, these vaccines do not induce strong therapeutic effects against established HPV infections and lesions. These concerns create a critical need for the development of therapeutic strategies, such as vaccines, to treat these existing infections and diseases.

Main Body

Unlike preventative vaccines, therapeutic vaccines aim to generate cell-mediated immunity. HPV oncoproteins E6 and E7 are responsible for the malignant progression of HPV-associated diseases and are consistently expressed in HPV-associated diseases and cancer lesions; therefore, they serve as ideal targets for the development of therapeutic HPV vaccines. In this review we revisit therapeutic HPV vaccines that utilize this knowledge to treat HPV-associated lesions and cancers, with a focus on the findings of recent therapeutic HPV vaccine clinical trials.

Conclusion

Great progress has been made to develop and improve novel therapeutic HPV vaccines to treat existing HPV infections and diseases; however, there is still much work to be done. We believe that therapeutic HPV vaccines have the potential to become a widely available and successful therapy to treat HPV and HPV-associated diseases in the near future.

Identification of the impact on T- and B- cell epitopes of human papillomavirus type-16 E6 and E7 variant in Southwest China

Cervical cancers almost are infected by human papillmavirus (HPV), encoding E6 and E7 oncoproteins which are regard as ideal targets on the mechanism of this disease and development of vaccines. HLA (human leukocyte antigen) participates in the local immune response to prevent tumor invasion and progression. But due to highly polymorphism of HLA, prediction shows its importance in this study. More effective immunoinformatics was used for predicting epitopes from HPV-16 E6 and E7, including T- and B-cell epitopes. Eight substitutions are detected. Specifically speaking, for HLA-I, HLA-A*33:03 (26), HLA-B*13:01 (14), HLA-C*03:02 (5) for E6 and HLA-A*02:01 (6), HLA-B*40:01 (5), HLA-C*03:04 (4) for E7 are most frequency. Epitope _41-48EVYDFAFR for HLA-A*33:03 (0.1) for E6 has best binding affinity, as well as HLA*02:01 and HLA-B*40:01 (0.2) for E7. The mutations of D25E and L83V of E6 and N29S of E7 produce new epitopes, and the percentile values change with them. For HLA-II, seventeen epitopes in the reference at percentile value from 0.22 to 4.76, while in variant from 0.22 to 4.96. For the B-cell epitopes, three most potent epitopes for E6 were listed, and N29S lead the growth of score from 0.81 to 0.83. In summary, E6_40-55REVYDFAFRDLCIVYR and E7_11-22YMLDLQPETTDL are the important regions, containing the majority of predicted epitopes. E6 72-83 for HLA-A*02:01 and E6 74-84 for HLA-B*15:02 maybe are the new direct for therapeutic vaccine aimed at L83V variants. HLA-DRB1*15:02 is better binder with T cell in our HLA class II. It is a systematic, detail recognition for T- and B-cell epitopes of HPV-16 E6 and E7 from Southwest China, which may be helpful to design vaccines specifically for women in Southwest China and testing methods specifically for this region. The results of our study may contribute to future researches on vaccines improvement, or screening methods for a particular population.

Genetic variability and lineage phylogeny of human papillomavirus type-16 and -53 based on the E6, E7, and L1 genes in Southwest China

Human papillomaviruses (HPVs) are circular double-stranded DNA viruses that are highly prevalent in the general population, and account for the cervical cancer burden in women worldwide. In this study, we analyzed HPV-16, the most prevalent type worldwide, and HPV-53, a possible high-risk type from infected women in Southwest China. To characterize mutations, intratypic variants, and genetic variability in the E6, E7, and L1 genes of HPV-16 (n=97) and HPV-53 (n=15), these genes were sequenced and submitted to GenBank. Phylogenetic trees were constructed using Bayesian trees, followed by secondary structure analysis and B-cell epitope prediction. Moreover, the selection pressures of the E6, E7, and L1 genes were estimated. In total, 27 novel variants of HPV-16 and 11 novel variants of HPV-53 were identified. In the HPV-16 E6-E7-L1 sequences, 73 nucleotide changes were observed with 40/73 being non-synonymous mutations (two in the alpha helix and five in the beta sheet) and 33/73 being synonymous. In the HPV-53 E6-E7-L1 sequences, 64 nucleotide changes were observed with 26/64 being non-synonymous mutations (three in the alpha helix and one in the beta sheet) and 38/64 being synonymous. Selective pressure analysis showed that most of these mutations did not reflect positive selection. The maximal divergence between any two variants within each gene of these two HPV types ranging from 0.94%(HPV-16 L1 gene)to 2.80%(HPV-53 E6 gene). Identifying new variants of HPV-16 and -53 from women in Southwest China may be helpful to design vaccines specifically for women in Southwest China and testing methods specifically for this region. The results of our study may contribute to future researches in diagnostic probes, vaccines improvement, or screening methods for a particular population.

The double-edge role of B cells in mediating antitumor T-cell immunity: Pharmacological strategies for cancer immunotherapy

Emerging evidence reveals the controversial role of B cells in antitumor immunity, but the underlying mechanisms have to be explored. Three latest articles published in the issue 521 of Nature in 2015 reconfirmed the puzzling topic and put forward some explanations of how B cells regulate antitumor T-cell responses both positively and negatively. This paper attempts to demonstrate that different B-cell subpopulations have distinct immunological properties and that they are involved in either antitumor responses or immunosuppression. Recent studies supporting the positive and negative roles of B cells in tumor development were summarized comprehensively. Several specific B-cell subpopulations, such as IgG(+), IgA(+), IL-10(+), and regulatory B cells, were described in detail. The mechanisms underlying the controversial B-cell effects were mainly attributed to different B-cell subpopulations, different B-cell-derived cytokines, direct B cell-T cell interaction, different cancer categories, and different malignant stages, and the immunological interaction between B cells and T cells is mediated by dendritic cells. Promising B-cell-based antitumor strategies were proposed and novel B-cell regulators were summarized to present interesting therapeutic targets. Future investigations are needed to make sure that B-cell-based pharmacological strategies benefit cancer immunotherapy substantially.

Current state in the development of candidate therapeutic HPV vaccines

The identification of human papillomavirus (HPV) as an etiological factor for HPV-associated malignancies creates the opportunity to control these cancers through vaccination. Currently, available preventive HPV vaccines have not yet demonstrated strong evidences for therapeutic effects against established HPV infections and lesions. Furthermore, HPV infections remain extremely common. Thus, there is urgent need for therapeutic vaccines to treat existing HPV infections and HPV-associated diseases. Therapeutic vaccines differ from preventive vaccines in that they are aimed at generating cell-mediated immunity rather than neutralizing antibodies. The HPV-encoded early proteins, especially oncoproteins E6 and E7, form ideal targets for therapeutic HPV vaccines since they are consistently expressed in HPV-associated malignancies and precancerous lesions, playing crucial roles in the generation and maintenance of HPV-associated disease. Our review will cover various therapeutic vaccines in development for the treatment of HPV-associated lesions and cancers. Furthermore, we review strategies to enhance vaccine efficacy and the latest clinical trials on therapeutic HPV vaccines.

Identification and validation of immunogenic potential of India specific HPV-16 variant constructs: In-silico &in-vivo insight to vaccine development

Cervical cancer is one of the most common gynecological cancers in the world but in India, it is the top most cancer among women. Persistent infection with high-risk human papillomaviruses (HR-HPVs) is the most important risk factor. The sequence variation(s) in the most common HR-HPV i.e. HPV type 16 leads to altered biological functions with possible clinical significance in the different geographical locations. Sixteen major variants (V1-V16) in full length L1 gene of HPV-16 were identified following analysis of 250 prospectively collected cervical cancer tissue biopsies and their effect on immunogenicity was studied. The effect of these major variations on the epitopes were predicted by in silico methods and the immunogenicity of variants and respective reference DNA vaccine constructs were evaluated by administration of prepared DNA vaccine constructs in female BALB/c mice to evaluate antibody titer. In the present study, L500F (V16) variation showed a significant ~2.7 fold (p < 0.002) increase in antibody titer, whereas T379P (V8) showed ~0.4 fold (p < 0.328) decrease after final injection. These results showed a promising roadmap for the development of DNA based vaccine and for the generation of effective response, though there is a need to study more prevalent variants of HPV in the Indian population.

Identification of immunotherapeutic epitope of E5 protein of human papillomavirus-16: An in silico approach

Cervical cancer is one of the most common gynaecological cancer in India and contributes 1/3rd of global burden. High risk-human papillomavirus (HR-HPV) is the major etiological factor for development of cervical cancer. Two available HPV vaccines provide protection against HPV induced cervical malignancy. However, vaccines having therapeutic values are of utmost priority. Till date, most of HPV therapeutic vaccines are focused on two major HPV oncoproteins (E6/E7). HPV-E5 which acts by altering the activity of cellular proteins, mainly growth factor pathways emerges as a new therapeutic target. In present study, we predicted the candidate B-cell and T-cell epitopes of HPV16-E5, which can be used for HPV immunotherapy. We identified that epitope SAFRCFIVYIIFVY as most potent peptide for HLA-A*11:01 having percentile value of 0.5 and immunogenicity score of 0.69558. For MHC-II, epitopes IPLFLIHTHARFLIT for HLA-DRB1*14:01 alleles have the lowest IC50 value (18.13 nM). The identification of structural feature and immunogenic epitopes provides the best information for development of drugs or vaccine. In conclusion, the expression of E5 protein was detected in the early phase of the HPV infection, which gives an opportunity to target HPV-E5 that would help in the prevention and progression of the precancerous lesion to invasive carcinomas.