Exploring the Papillomaviral Proteome to Identify Potential Candidates for a Chimeric Vaccine against Cervix Papilloma Using Immunomics and Computational Structural Vaccinology

The human papillomavirus (HPV) 58 is considered to be the second most predominant genotype in cervical cancer incidents in China. HPV type-restriction, non-targeted delivery, and the highcost of existing vaccines necessitate continuing research on the HPV vaccine. We aimed to explore the papillomaviral proteome in order to identify potential candidates for a chimeric vaccine against cervix papilloma using computational immunology and structural vaccinology approaches. Two overlapped epitope segments (23⁻36) and (29⁻42) from the N-terminal region of the HPV58 minor capsid protein L2 are selected as capable of inducing both cellular and humoral immunity. In total, 318 amino acid lengths of the vaccine construct SGD58 contain adjuvants (Flagellin and RS09), two Th epitopes, and linkers. SGD58 is a stable protein that is soluble, antigenic, and non-allergenic. Homology modeling and the structural refinement of the best models of SGD58 and TLR5 found 96.8% and 93.9% favored regions in Rampage, respectively. The docking results demonstrated a HADDOCK score of -62.5 ± 7.6, the binding energy (-30 kcal/mol) and 44 interacting amino acid residues between SGD58-TLR5 complex. The docked complex are stable in 100 ns of simulation. The coding sequences of SGD58 also show elevated gene expression in Escherichia coli with 1.0 codon adaptation index and 59.92% glycine-cysteine content. We conclude that SGD58 may prompt the creation a vaccine against cervix papilloma.

Ancient Evolution and Dispersion of Human Papillomavirus 58 Variants

Human papillomavirus 58 (HPV58) is found in 10 to 18% of cervical cancers in East Asia but is rather uncommon elsewhere. The distribution and oncogenic potential of HPV58 variants appear to be heterogeneous, since the E7 T20I/G63S variant is more prevalent in East Asia and confers a 7- to 9-fold-higher risk of cervical precancer and cancer. However, the underlying genomic mechanisms that explain the geographic and carcinogenic diversity of HPV58 variants are still poorly understood. In this study, we used a combination of phylogenetic analyses and bioinformatics to investigate the deep evolutionary history of HPV58 complete genome variants. The initial splitting of HPV58 variants was estimated to occur 478,600 years ago (95% highest posterior density [HPD], 391,000 to 569,600 years ago). This divergence time is well within the era of speciation between Homo sapiens and Neanderthals/Denisovans and around three times longer than the modern Homo sapiens divergence times. The expansion of present-day variants in Eurasia could be the consequence of viral transmission from Neanderthals/Denisovans to non-African modern human populations through gene flow. A whole-genome sequence signature analysis identified 3 amino acid changes, 16 synonymous nucleotide changes, and a 12-bp insertion strongly associated with the E7 T20I/G63S variant that represents the A3 sublineage and carries higher carcinogenetic potential. Compared with the capsid proteins, the oncogenes E7 and E6 had increased substitution rates indicative of higher selection pressure. These data provide a comprehensive evolutionary history and genomic basis of HPV58 variants to assist further investigation of carcinogenic association and the development of diagnostic and therapeutic strategies.IMPORTANCE Papillomaviruses (PVs) are an ancient and heterogeneous group of double-stranded DNA viruses that preferentially infect the cutaneous and mucocutaneous epithelia of vertebrates. Persistent infection by specific oncogenic human papillomaviruses (HPVs), including HPV58, has been established as the primary cause of cervical cancer. In this work, we reveal the complex evolutionary history of HPV58 variants that explains the heterogeneity of oncogenic potential and geographic distribution. Our data suggest that HPV58 variants may have coevolved with archaic hominins and dispersed across the planet through host interbreeding and gene flow. Certain genes and codons of HPV58 variants representing higher carcinogenic potential and/or that are under positive selection may have important implications for viral host specificity, pathogenesis, and disease prevention.

Designing of Epitope-Focused Vaccine by Targeting E6 and E7 Conserved Protein Sequences: An Immuno-Informatics Approach in Human Papillomavirus 58 Isolates

Human papillomavirus (HPV) is a DNA virus that belongs to the papillomavirus family and is capable of infecting humans. Currently, few vaccines are available to prevent infection by HPV. However, they are not so much effective and provide little benefit to women who have already been infected with HPV. The aim of this study was to design epitope-based vaccines of HPV58 by targeting E6 and E7 proteins of HPV58. Proteomic sequences were retrieved from different isolates at different time periods and later analyzed by performing alignment of these sequences. To ensure the capacity of humoral and cell-mediated immunity, both B cell and T cell immunity were checked for the peptides. For E6 protein, the peptide sequence from 48 to 54 amino acids and one 9-m epitope ETSVHEIEL were the most potential B cell and T cell epitopes, respectively. This peptide could interact with as many as eight MHC-1 alleles and showed high population coverage up to 90.31 %. On the other hand, the peptide region for the E7 protein ranged from 27 to 33 amino acids and two 9-m epitopes QAQPATANY, SSDEDEIGL were found as the most potential B cell and T cell epitopes, respectively. The peptide sequences could interact with as many as seven MHC-1 alleles and showed population coverage up to 90.31 %. Furthermore, conservancy analysis was also performed using in silico tools and showed a conservancy of 100 % for all the selected epitopes. In addition to this, the allergenicity of the epitopes was also evaluated. Although the study requires further in vitro and in vivo screening, this epitope-focused peptide vaccine designing opens up a new skyline that holds a prospective future in HPV research.

Genetic variability and oncogenic risk association of human papillomavirus type 58 E6 and E7 genes in Taizhou area, China

It is well recognized that the association of human papillomavirus (HPV) and cervical carcinogenesis is based on the presence of HPV DNA sequence. The E6 and E7 oncoproteins encoded by high-risk HPV types play a key role in carcinogenesis. HPV58 type accounts for a larger share of cervical disease in China, whereas data on HPV58 genetic variability in China is limited. We aimed to evaluate the diversity of HPV58 genetic variants by sequencing the entire E6 and E7 genes. Phylogenetic trees were constructed by Maximum likelihood method by MEGA 5.05 software. In this study, the overall HPV infection rate was 22.6% (2891/12780) in Southeast China and the prevalence of HPV58 infection rate was 2.6% (335/12780). 26 nucleotides substitutions were observed in E6 and E7 genes with 10 novel substitutions and 17 non-synonymous substitutions. We obtained 25 distinct variation patterns which the accession GenBank numbers as MH348918-MH348942. All of HPV58 variants belong to lineage A, while no lineage B, C and D were detected in Taizhou area, Southeast China. The sublineage A1, A2, and A3 variants were found in 136 (68.3%), 39 (19.6%), and 24 (12.1%) of HPV58 isolates, respectively. The sublineage A3 variants with T20I/G63S substitutions at E7 oncoprotein carried a significantly higher risk for high-grade cervical intraepithelial neoplasia (CIN2 or worse, CIN2+) when compared with other HPV58 variants (odds ratio = 4.41, P < 0.05). Nevertheless, there was no association between HPV58 (sub) lineages and cervical lesions. These data provide the critical characteristics of HPV58 variants to assist further investigation of carcinogenic association and the development of next generation vaccines and diagnostic assays in China.

Genetic diversity of human papillomavirus types 35, 45 and 58 in cervical cancer in Brazil

In Brazil, most studies of intra-type variants of human papillomavirus (HPV) have focused on HPV16 and HPV18, but other high-risk HPV types have not been studied. Here, we report the prevalence of lineages and variants of HPV35, HPV45 and HPV58 in cervical cancers from the Amazonian and Southeast Brazilian regions. The most frequent sublineages were A1 for HPV35, B2 for HPV45, and A2 for HPV58. The Southeast region had a higher frequency of the B2 sublineage of HPV45, and for HPV35, the genetic and nucleotide sequence diversity were higher in the Southeast region, suggesting that regional factors are influencing the diversity and lineage prevalence.

Oncogenic comparison of human papillomavirus type 58 E7 variants

Human papillomavirus 58 (HPV58) ranks the second or third in East Asian cervical cancers. Current studies on HPV58 are scarce and focus on the prototype. Previously, we identified the three most common circulating HPV58 E7 strains contained amino acid alterations: G41R/G63D (51%), T20I/G63S (22%) and T74A/D76E (14%) respectively. Among them, the T20I/G63S variant (V1) had a stronger epidemiological association with cervical cancer. We therefore suggested that V1 possessed stronger oncogenicity than the other two variants. Here, we performed phenotypic assays to characterize and compare their oncogenicities with HPV58 E7 prototype. Our results showed that overexpression of V1 conferred a higher colony-forming ability to primary murine epithelial cells than prototype (P < 0.05) and other variants, implicating its higher immortalising potential. Further experiments showed that both V1 and prototype enhanced the anchorage-independent growth of NIH/3T3 cells (P < 0.001), implicating their stronger transforming power than the two other variants. Moreover, they possessed an increased ability to degrade pRb (P < 0.001), which is a major effector pathway of E7-driven oncogenesis. Our work represents the first study to compare the oncogenicities of HPV58 E7 prototype and variants. These findings deepened our understanding of HPV58 and might inform clinical screening and follow-up strategy.

Prognosis and related factors of HPV infections in postmenopausal Uyghur women

With the aim to explore the characteristics of persistent HPV infections in postmenopausal Uyghur women and analyse the possible related risk factors, from September 2012 to September 2013; postmenopausal Uyghur women with HPV positive and pathologically diagnosed as non-cervical intraepithelial neoplasia (CIN) lesions and non-cervical cancer were recruited. Their clinical course was closely followed up for 24-36 months, and the risk factors were analysed by a logistic regression model. One hundred and sixteen positive women were followed for 36 months. The total persistent HPV infection rate was 67.9%, and the type-specific persistent infection rate was 73.7% at 36 months. Nine (32.1%) women were naturally cleared of their HPV infection at 36 months. We found that an HPV16 infection and an HPV58 infection, and time since menopause over 2 years were closely related with a persistent HPV infection. More attention should be paid to the women above 2 years of menopause who were infected with HPV16 and HPV58 in their further cervical carcinoma screening. Impact statement What is already known on this subject? Previous study revealed that menopause was a risk factor for a persistent HPV infection in Uyghur women. What do the results of this study add? The present study presented the characteristics of HPV persistent infection and the risk factors in Uyghur postmenopausal women. More attention should be paid to the women above 2 of years of menopause who are infected with HPV16 and HPV58. What are the implications of these findings for clinical practice and/or further research? This study would offer a theoretical basis for a better screening design, especially the women above 2 years' menopause who have been infected with HPV16 and HPV58 in the Xinjiang region.

Combined truncations at both N- and C-terminus of human papillomavirus type 58 L1 enhanced the yield of virus-like particles produced in a baculovirus system

Human papillomavirus (HPV) major capsid protein L1 virus-like particles (VLPs) produced in baculovirus system are highly immunogenic, but the relatively high production cost limits its application in the development of broad-spectrum vaccines. Here we report a novel method for enhancing VLP production in this system. We incorporated respectively 4, 8 or 13 residues truncation mutations in the N-terminus of L1ΔC, a C-terminal 25-residue-deleted L1 of HPV58, to construct three mutants. After expression in Sf9 cells, L1ΔN4C exhibited 2.3-fold higher protein production, 2.0-fold mRNA expression and lower rate of mRNA decay, compared to L1ΔC. More importantly, L1ΔN4C protein was easily purified by two-step chromatography with a VLP yield of up to 60 mg/L (purity > 99 %), 5-fold that of L1ΔC, whereas L1ΔN8C and L1ΔN13C behaved similarly to L1ΔC either in protein or mRNA expression. Moreover, L1ΔN4C VLPs showed similar binding activities with six HPV58 neutralizing monoclonal antibodies and induced comparable level of neutralizing antibody in mice to that of L1ΔC VLPs. Our results demonstrate that certain N- and C-terminal truncations of HPV58 L1 can enhance VLP yield. This method may be used to reduce production costs of other L1VLPs or chimeric VLPs to developing pan-HPV vaccines using baculovirus system.

Viral load and integration status of HPV58 associated with cervical lesion severity in women of Northeast China

Objective

We aimed to evaluate the viral load and integration status of human papillomavirus 58 in women with different grades of cervical lesions to determine whether viral load and integration status are related to malignant transformation in HPV58-infected women.

Methods

A total of 212 cervical specimens were collected from women in Northeast China who had undergone human papillomavirus genotyping and were HPV58-positive. The HPV58 viral load was determined using real-time polymerase chain reaction, and the integration status was discriminated using the ratio of HPV58 E2 gene copy number to E6 gene copy number.

Results

The median HPV58 viral load in women with normal cervix or cervicitis, low-grade squamous cell intraepithelial lesion, high-grade squamous cell intraepithelial lesion and cervical cancer was 352.12, 864.21, 1199.75 and 693.04 copies/genome, respectively. High significance was obtained when comparing the viral load of infected women presenting normal/cervicitis with that of the women either with precancerous cervical lesions or cervical cancer (P < 0.05). The HPV58 genome was in the episomal form in 35 samples (16.5%), mixed episomal and integrated forms in 165 (77.8%) samples, and completely integrated into the host genome in 12 (5.7%) samples. The HPV58 E2/E6 copy number ratio in the cervical cancer group was significantly lower than that in the other groups (P < 0.01).

Conclusions

The HPV58 viral load in patients with precancerous cervical lesions or cervical cancer increases significantly with disease progression. The HPV58 E2/E6 copy number ratio in patients with cervical cancer is lower than that for less severe cervical lesions, suggesting a high degree of viral integration may be a considerable risk factor for cervical cancer.

A Comparative Study on Delivery of Externally Attached DNA by Papillomavirus VLPs and Pseudoviruses

Human papillomavirus (HPV) 16 capsids have been chosen as a DNA delivery vehicle in many studies. Our preliminary studies suggest that HPV58 capsids could be better vehicles than HPV16 capsids to deliver encapsidated DNA in vitro and in vivo. In the current study, we compared HPV16, HPV58, and the cottontail rabbit papillomavirus (CRPV) capsids either as L1/L2 VLPs or pseudoviruses (PSVs) to deliver externally attached GFP-expressing DNA. Both rabbit and human cells were used to test whether there was a species-specific effect. DNA delivery efficiency was determined by quantifying either GFP-expressing cell populations or mean fluorescent intensities (MFI) by flow cytometry. Interestingly, CRPV and 58-VLPs and PSVs were significantly more efficient at delivering attached DNA when compared to 16-VLPs and PSVs. A capsid/DNA ratio of 2:1 showed the highest efficiency for delivering external DNA. The PSVs with papillomavirus DNA genomes also showed higher efficiency than those with irrelevant plasmid DNA. HPV16L1/58L2 hybrid VLPs displayed increased efficiency compared to HPV58L1/16L2 VLPs, suggesting that L2 may play a critical role in the delivery of attached DNA. Additionally, we demonstrated that VLPs increased in vivo infectivity of CRPV DNA in rabbits. We conclude that choosing CRPV or 58 capsids to deliver external DNA could improve DNA uptake in in vitro and in vivo models.

Development of HPV58 type-specific antibodies and detection kit

Currently, human papillomavirus (HPV) vaccines are in short supply, so the development of HPV vaccines has a broad market prospect. The 3-, 9-, and 15-valent HPV vaccines developed by ourselves all contain HPV58-derived antigen components. It is important to detect HPV 58 during vaccine production. Here, we introduced a development process of HPV58 type-specific antibodies and a detection kit. Briefly, HPV58 L1-Virus Like Particles (VLPs) were used as antigens to immunize mice, followed by extraction of the ascites to prepare hybridoma cells. After culturing, the supernatants containing secreted antibodies were harvested, purified, and screened to obtain monoclonal antibodies (mAbs). In the pool of attained monoclonal antibodies, we selected 2F7 and 2G7 to evaluate their subtypes, specificity, neutralizing activity, serum competition, binding affinity and gene sequencing. Finally, an enzyme-linked immunosorbent assay (ELISA) detection kit was assembled with 2F7 and 2G7 mAbs which possessed high specificity to HPV58 L1-VLPs. The detection kit developed by 2F7 and 2G7 could be adopted to specifically detect HPV58 L1 protein with good linearity and detection range, which could be widely used in clinical testing and quality control in the production of HPV vaccines.Abbreviations: BSA: Bovine serum albumin; CDRs: Complementarity-determining regions; CV: Coefficient of variation; DTT: Dithiothreitol; ELISA: Enzyme-linked immunosorbent assay; HAT: Hypoxanthine-aminopterin-thymidine; HPV: Human Papillomavirus; IC₅₀: 50% inhibition rate; IC₉₀: 90% inhibition rate; mAbs: Monoclonal antibodies; VLP: Virus-like particle.

Distribution of HPV types among women with HPV-related diseases and exploration of lineages and variants of HPV 52 and 58 among HPV-infected patients in China: A systematic literature review

To summarize the distribution of types of human papillomavirus (HPV) associated with HPV-related diseases and investigate the potential causes of high prevalence of HPV 52 and 58 by summarizing the prevalence of lineages, sub-lineages, and mutations among Chinese women. We searched PubMed, EMBASE, CNKI, and WanFang from January, 2012 to June, 2023 to identify all the eligible studies. We excluded patients who had received HPV vaccinations. Data were summarized in tables and cloud/rain maps. A total of 102 studies reporting HPV distribution and 15 studies reporting HPV52/HPV58 variants were extracted. Among Chinese women, the top five prevalent HPV types associated with cervical cancer (CC) were HPV16, 18, 58, 52, and 33. In patients with vaginal cancers and precancerous lesions, the most common HPV types were 16 and 52 followed by 58. For women with condyloma acuminatum (CA), the most common HPV types were 11 and 6. In Chinese women with HPV infection, lineage B was the most prominently identified for HPV52, and lineage A was the most common for HPV58. In addition to HPV types 16, which is prevalent worldwide, our findings revealed the unique high prevalence of HPV 52/58 among Chinese women with HPV-related diseases. HPV 52 variants were predominantly biased toward lineage B and sub-lineage B2, and HPV 58 variants were strongly biased toward lineage A and sub-lineage A1. Further investigations on the association between the high prevalent lineage and sub-lineage in HPV 52/58 and the risk of cancer risk are needed. Our findings underscore the importance of vaccination with the nine-valent HPV vaccine in China.

An E7-retinoblastoma protein pathway mechanism may account for the higher carcinogenic ability of HPV16 over HPV58 in cervical cancer

Background

Among human papillomavirus (HPV) type, HPV16 displays the strongest carcinogenic capacity for cervical cancer, but the mechanism underlying this phenomenon remains unclear. We investigated the effect and the underlying mechanism of HPV16 on higher carcinogenic capacity than HPV58.

Methods

We collected 4,030 cervical exfoliated cell samples for genotyping HPV using HybriBio's proprietary flow-through hybridization technique, liquid-based cytology (LBC), colposcopy, and biopsies if indicated. Four plasmids containing E6 and E7 of HPV16 and 58 were constructed and transfected into 293T and U2OS cells. We detected the cell phenotype using Cell Counting Kit 8 (CCK8) assay, Transwell assay, flow cytometry, and apoptosis assay; the expression of retinoblastoma protein (Rb) and phosphorylated Rb (pRb) was determined via Western blot; and the cell activity was determined via a zebrafish model treated with or without roscovitine.

Results

The positive rates of HPV16 and 58 were, respectively, 18.9% and 19.7% in the ≤ low-grade squamous intraepithelial lesion (LSIL) group, 49.5% and 19.6% (P<0.001) in the high-grade squamous intraepithelial lesion (HSIL) group, 65.3% and 9.0% (P<0.001) in the cancer group. In vitro, both 293T and U2OS cells with overexpressed HPV16 E6 and E7 displayed significantly higher cell proliferation, faster cell invasion, decreased cell apoptosis, and accelerated cell cycle from G1 phase to S phase compared to those with overexpressed HPV58 E6 and E7 (all P values <0.05). Rb loss of function was observed in cells with HPV16 E7 overexpression, while a greater level of phosphorylated Rb was observed in cells with HPV58 E7 overexpression. Roscovitine restored Rb expression and decreased the cell activity in zebrafish.

Conclusions

HPV16 possesses a stronger carcinogenic ability than does HPV 58, and the mechanism underlying this effect may be the impairment of the E7-Rb pathway.