Mutations on the FG surface loop of human papillomavirus type 16 major capsid protein affect recognition by both type-specific neutralizing antibodies and cross-reactive antibodies

Genetic diversity of HPV35 in Chad and the Central African Republic, two landlocked countries of Central Africa: A cross-sectional study

Human Papillomavirus (HPV)-35 accounts for up 10% of cervical cancers in Sub-Saharan Africa. We herein assessed the genetic diversity of HPV35 in HIV-negative women from Chad (identified as #CHAD) and HIV-infected men having sex with men (MSM) in the Central African Republic (CAR), identified as #CAR. Ten HPV35 DNA from self-collected genital secretions (n = 5) and anal margin samples (n = 5) obtained from women and MSM, respectively, were sequenced using the ABI PRISM® BigDye Sequencing technology. All but one HPV35 strains belonged to the A2 sublineage, and only #CAR5 belonged to A1. HPV35 from #CAR had higher L1 variability compared to #CHAD (mean number of mutations: 16 versus 6). L1 of #CAR5 showed a significant variability (2.29%), suggesting a possible intra-type divergence from HPV35H. Three (BC, DE, and EF) out of the 5 capsid loops domains remained totally conserved, while FG- and HI- loops of #CAR exhibited amino acid variations. #CAR5 also showed the highest LCR variability with a 16bp insertion at binding sites of the YY1. HPV35 from #CHAD exhibited the highest variability in E2 gene (P<0.05). E6 and E7 oncoproteins remained well conserved. There is a relative maintenance of a well conserved HPV35 A2 sublineage within heterosexual women in Chad and MSM with HIV in the Central African Republic.

A bacterially expressed triple-type chimeric vaccine against human papillomavirus types 51, 69, and 26

Persistent infection of high-risk human papillomavirus (HPV) is a leading cause of some cancers, including cervical cancer. However, with over 20 carcinogenic HPV types, it is difficult to design a multivalent vaccine that can offer complete protection. Here, we describe the design and optimization of a HPV51/69/26 triple-type chimeric virus-like particle (VLP) for vaccine development. Using E. coli and a serial N-terminal truncation strategy, we created double- and triple-type chimeric VLPs through loop-swapping at equivalent surface loops. The lead candidate, H69-51BC-26FG, conferred similar particulate properties as that of its parental VLPs and comparable immunogenicity against HPV51, -69 and -26. When produced in a GMP-like facility, these H69-51BC-26FG VLPs were verified to have excellent qualities for the development of a multivalent HPV vaccine. This study showcases an amenable way to create a single VLP using type-specific epitope clustering for the design of a triple-type vaccine.

Genetic variation within the human papillomavirus type 16 genome is associated with oropharyngeal cancer prognosis

PURPOSE

A significant barrier to adoption of de-escalated treatment protocols for human papillomavirus-driven oropharyngeal cancer (HPV-OPC) is that few predictors of poor prognosis exist. We conducted the first large whole-genome sequencing (WGS) study to characterize the genetic variation of the HPV type 16 (HPV16) genome and to evaluate its association with HPV-OPC patient survival.

PATIENTS AND METHODS

A total of 460 OPC tumor specimens from two large United States medical centers (1980-2017) underwent HPV16 whole-genome sequencing. Site-specific variable positions [single nucleotide polymorphisms (SNPs)] across the HPV16 genome were identified. Cox proportional hazards model estimated hazard ratios (HRs) and 95% confidence intervals (CIs) for overall survival by HPV16 SNPs. Harrell C-index and time-dependent positive predictive value (PPV) curves and areas under the PPV curves were used to evaluate the predictive accuracy of HPV16 SNPs for overall survival.

RESULTS

A total of 384 OPC tumor specimens (83.48%) passed quality control filters with sufficient depth and coverage of HPV16 genome sequencing to be analyzed. Some 284 HPV16 SNPs with a minor allele frequency ≥1% were identified. Eight HPV16 SNPs were significantly associated with worse survival after false discovery rate correction (individual prevalence: 1.0%-5.5%; combined prevalence: 15.10%); E1 gene position 1053 [HR for overall survival (HRos): 3.75, 95% CI 1.77-7.95; Pfdr = 0.0099]; L2 gene positions 4410 (HRos: 5.32, 95% CI 1.91-14.81; Pfdr = 0.0120), 4539 (HRos: 6.54, 95% CI 2.03-21.08; Pfdr = 0.0117); 5050 (HRos: 6.53, 95% CI 2.34-18.24; Pfdr = 0.0030), and 5254 (HRos: 7.76, 95% CI 2.41-24.98; Pfdr = 0.0030); and L1 gene positions 5962 (HRos: 4.40, 95% CI 1.88-10.31; Pfdr = 0.0110) and 6025 (HRos: 5.71, 95% CI 2.43-13.41; Pfdr = 0.0008) and position 7173 within the upstream regulatory region (HRos: 9.90, 95% CI 3.05-32.12; Pfdr = 0.0007). Median survival time for patients with ≥1 high-risk HPV16 SNPs was 3.96 years compared with 18.67 years for patients without a high-risk SNP; log-rank test P < 0.001. HPV16 SNPs significantly improved the predictive accuracy for overall survival above traditional factors (age, smoking, stage, treatment); increase in C-index was 0.069 (95% CI 0.019-0.119, P < 0.001); increase in area under the PPV curve for predicting 5-year survival was 0.068 (95% CI 0.015-0.111, P = 0.008).

CONCLUSIONS

HPV16 genetic variation is associated with HPV-OPC prognosis and can improve prognostic accuracy.

Structural basis for the shared neutralization mechanism of three classes of human papillomavirus type 58 antibodies with disparate modes of binding

Human papillomavirus type 58 (HPV58) is associated with cervical cancer and poses a significant health burden worldwide. Although the commercial 9-valent HPV vaccine covers HPV58, the structural and molecular-level neutralization sites of the HPV58 complete virion are not fully understood. Here, we report the high-resolution (∼3.5 Å) structure of the complete HPV58 pseudovirus (PsV58) using cryo-electron microscopy (cryo-EM). Three representative neutralizing monoclonal antibodies (nAbs 5G9, 2H3 and A4B4) were selected through clustering from a nAb panel against HPV58. Bypassing the steric hindrance and symmetry-mismatch in the HPV Fab-capsid immune-complex, we present three different neutralizing epitopes in the PsV58, and show that, despite differences in binding, these nAbs share a common neutralization mechanism. These results offer insight into HPV58 genotype specificity and broaden our understanding of HPV58 neutralization sites for antiviral research.IMPORTANCE Cervical cancer primarily results from persistent infection with high-risk types of human papillomavirus (HPV). HPV type 58 (HPV58) is an important causative agent, especially within Asia. Despite this, we still have limited data pertaining to the structural and neutralizing epitopes of HPV58, and this encumbers our in-depth understanding of the virus mode of infection. Here, we show that representative nAbs (5G9, 10B11, 2H3, 5H2 and A4B4) from three different groups share a common neutralization mechanism that appears to prohibit the virus from associating with the extracellular matrix and cell surface. Furthermore, we identify that the nAbs engage via three different binding patterns: top-center binding (5G9 and 10B11), top-fringe binding (2H3 and 5H2), and fringe binding (A4B4). Our work shows that, despite differences in the pattern in binding, nAbs against HPV58 share a common neutralization mechanism. These results provide new insight into the understanding of HPV58 infection.

In silico prediction of T-cell and B-cell epitopes of human papillomavirus type 16 L1 protein

Human papillomavirus type 16 (HPV-16) is one of the most important cause of developing cervical cancer. Therefore, effective epitope-based vaccine design for HPV-16 would be of major medical benefit. The aim of our study was to identify B- and T-cell epitopes of HPV-16 L1 protein. In this study, the HPV-16 L1 gene was isolated from HPV recovered from five vaginal swab samples using specific primers and finally sequenced. The ExPASy translate tool (http://web.expasy.org/translate/) was used to convert nucleotide sequence into amino acid sequence. Bioinformatic analysis was employed to predict suitable B- and T-cell epitopes and immunogenicity, allergenicity, and toxicity of predicted epitopes were then evaluated. Afterward, the selected T-cell epitopes were docked using Molegro Virtual Docker software. The two epitopes ₂₀₇ AMDFTTLQA₂₁₅ and ₂₀₀ MVDTGFGAM₂₀₈ have showed a very strong binding affinity to HLA-A0201 and HLA-B3501 molecules, respectively. Outcome of B-cell epitope prediction showed that epitope ₄₇₅ KAKPKFTLGKRK ATPTTSSTSTTAKRKK₅₀₂ contained overlapped epitope, which might be the epitope associated with the production of neutralizing antibody response. Based on this finding, the predicted B- and T-cell epitopes are promising targets for epitope-based vaccine development against HPV-16. Further in vivo and in vitro experiments are needed to confirm our findings.

Worldwide genetic variations in high-risk human papillomaviruses capsid L1 gene and their impact on vaccine efficiency

BACKGROUND

Human papillomavirus is the most common sexually transmitted infection. It is associated with different cancers, mainly cervical cancer, which remains the fourth most frequent cancer among women worldwide; it is also related to anogenital (anus, vulvar, vagina, and penis) and oropharyngeal cancers. Vaccination against HPV infection is the major way of prevention, and it has demonstrated impressive efficacy in reducing cervical cancer incidence. Nowadays, all the licensed HPV recombinant vaccines were designed based on HPV major capsid L1 protein. However, some variations in the HPV L1 gene sequence may induce structural changes within the L1 protein, which may alter the affinity and interaction of monoclonal antibodies (MAbs) with L1 protein epitopes, and influence host immune response and recognition. Hence, the importance of accuracy in delineating epitopes relevant to vaccine design and defining genetic variations within antigenic regions in the L1 gene to predict its impact on prophylactic vaccine efficiency. The present review reports the sequence variations in HR-HPV L1 gene isolates from different countries around the world, which may help to understand the effect of HPV L1 gene variations on vaccine efficiency.

METHODS

Research studies were retrieved from PubMed, Google Scholar, Science direct, and the National Center for Biotechnology Information (NCBI) database. A total of 31 articles describing genetic variations within the major capsid L1 gene and conducted in Africa, Europe, America and Asia were found. Only 26 studies conducted on HPV16, 18, 31, 33, 58, 45 and 52 which are the targets of HPV prophylactic vaccines, and which reported genetic variations within the L1 gene, were selected and evaluated in this review.

FINDINGS

We found a total of 87, 49, 11, 7, 22, 3, and 17 non-synonymous single nucleotide polymorphisms (SNPs) within HPV16, HPV18, HPV31, HPV58, HPV45, and HPV52 L1 gene, respectively. Four mutations were frequently observed in HPV16 L1 sequences: T353P in the HI loop, H228D in the EF loop, T266A in the FG loop, and T292A in the FG loop. Two mutations in HPV58 L1 sequences: T375N in the HI loop and L150F in the DE loop. Three mutations in HPV33 L1 sequences: T56N in the BC loop, G133S in the DE loop, T266K in the FG loop. Other mutations were found in HPV18, HPV45, and HPV52 L1 sequences. Some were found in different countries, and others were specific to a given population. Furthermore, some variations were located on peptide binding epitopes and lead to a modification of epitopes, which may influence MAbs interactions. Others need further investigations due to the lack of studies.

CONCLUSION

This study investigated the major capsid L1 genetic diversity of HPV16, 18, 31, 33, 58, 45, and 52 circulating in different populations around the world. Further investigations should be conducted to confirm their effect on immunogenicity and prophylactic vaccine efficiency.

Identification and Analysis of Unstructured, Linear B-Cell Epitopes in SARS-CoV-2 Virion Proteins for Vaccine Development

The efficacy of SARS-CoV-2 nucleic acid-based vaccines may be limited by proteolysis of the translated product due to anomalous protein folding. This may be the case for vaccines employing linear SARS-CoV-2 B-cell epitopes identified in previous studies since most of them participate in secondary structure formation. In contrast, we have employed a consensus of predictors for epitopic zones plus a structural filter for identifying 20 unstructured B-cell epitope-containing loops (uBCELs) in S, M, and N proteins. Phylogenetic comparison suggests epitope switching with respect to SARS-CoV in some of the identified uBCELs. Such events may be associated with the reported lack of serum cross-protection between the 2003 and 2019 pandemic strains. Incipient variability within a sample of 1639 SARS-CoV-2 isolates was also detected for 10 uBCELs which could cause vaccine failure. Intermediate stages of the putative epitope switch events were observed in bat coronaviruses in which additive mutational processes possibly facilitating evasion of the bat immune system appear to have taken place prior to transfer to humans. While there was some overlap between uBCELs and previously validated SARS-CoV B-cell epitopes, multiple uBCELs had not been identified in prior studies. Overall, these uBCELs may facilitate the development of biomedical products for SARS-CoV-2.

The DE and FG loops of the HPV major capsid protein contribute to the epitopes of vaccine-induced cross-neutralising antibodies

The human papillomavirus (HPV) vaccines consist of major capsid protein (L1) virus-like particles (VLP) and are highly efficacious against the development of cervical cancer precursors attributable to oncogenic genotypes, HPV16 and HPV18. A degree of vaccine-induced cross-protection has also been demonstrated against genetically-related genotypes in the Alpha-7 (HPV18-like) and Alpha-9 (HPV16-like) species groups which is coincident with the detection of L1 cross-neutralising antibodies. In this study the L1 domains recognised by inter-genotype cross-neutralising antibodies were delineated. L1 crystallographic homology models predicted a degree of structural diversity between the L1 loops of HPV16 and the non-vaccine Alpha-9 genotypes. These structural predictions informed the design of chimeric pseudovirions with inter-genotype loop swaps which demonstrated that the L1 domains recognised by inter-genotype cross-neutralising antibodies comprise residues within the DE loop and the late region of the FG loop. These data contribute to our understanding of the L1 domains recognised by vaccine-induced cross-neutralising antibodies. Such specificities may play a critical role in vaccine-induced cross-protection.

Lessons learned from successful human vaccines: Delineating key epitopes by dissecting the capsid proteins

Recombinant VLP-based vaccines have been successfully used against 3 diseases caused by viral infections: Hepatitis B, cervical cancer and hepatitis E. The VLP approach is attracting increasing attention in vaccine design and development for human and veterinary use. This review summarizes the clinically relevant epitopes on the VLP antigens in successful human vaccines. These virion-like epitopes, which can be delineated with molecular biology, cryo-electron microscopy and x-ray crystallographic methods, are the prerequisites for these efficacious vaccines to elicit functional antibodies. The critical epitopes and key factors influencing these epitopes are discussed for the HEV, HPV and HBV vaccines. A pentamer (for HPV) or a dimer (for HEV and HBV), rather than a monomer, is the basic building block harboring critical epitopes for the assembly of VLP antigen. The processing and formulation of VLP-based vaccines need to be developed to promote the formation and stabilization of these epitopes in the recombinant antigens. Delineating the critical epitopes is essential for antigen design in the early phase of vaccine development and for critical quality attribute analysis in the commercial phase of vaccine manufacturing.

Evidences of Changes in Surface Electrostatic Charge Distribution during Stabilization of HPV16 Virus-Like Particles

The stabilization of human papillomavirus type 16 virus-like particles has been examined by means of different techniques including dynamic and static light scattering, transmission electron microscopy and electrophoretic mobility. All these techniques provide different and often complementary perspectives about the aggregation process and generation of stabilized virus-like particles after a period of time of 48 hours at a temperature of 298 K. Interestingly, static light scattering results point towards a clear colloidal instability in the initial systems, as suggested by a negative value of the second virial coefficient. This is likely related to small repulsive electrostatic interactions among the particles, and in agreement with relatively small absolute values of the electrophoretic mobility and, hence, of the net surface charges. At this initial stage the small repulsive interactions are not able to compensate binding interactions, which tend to aggregate the particles. As time proceeds, an increase of the size of the particles is accompanied by strong increases, in absolute values, of the electrophoretic mobility and net surface charge, suggesting enhanced repulsive electrostatic interactions and, consequently, a stabilized colloidal system. These results show that electrophoretic mobility is a useful methodology that can be applied to screen the stabilization factors for virus-like particles during vaccine development.

Cross-protection of the Bivalent Human Papillomavirus (HPV) Vaccine Against Variants of Genetically Related High-Risk HPV Infections

BACKGROUND

Results from the Costa Rica Vaccine Trial (CVT) demonstrated partial cross-protection by the bivalent human papillomavirus (HPV) vaccine, which targets HPV-16 and HPV-18, against HPV-31, -33, and -45 infection and an increased incidence of HPV-51 infection.

METHODS

A study nested within the CVT intention-to-treat cohort was designed to assess high-risk HPV variant lineage-specific vaccine efficacy (VE). The 2 main end points were (1) long-term incident infections persisting for ≥2 years and/or progression to high-grade squamous intraepithelial lesions (ie, cervical intraepithelial neoplasia grade 2/3 [CIN 2/3]) and (2) incident transient infections lasting for <2 years. For efficiency, incident infections due to HPV-16, -18, -31, -33, -35, -45, and -51 resulting in persistent infection and/or CIN 2/3 were matched (ratio, 1:2) to the more-frequent transient viral infections, by HPV type. Variant lineages were determined by sequencing the upstream regulatory region and/or E6 region.

RESULTS

VEs against persistent or transient infections with HPV-16, -18, -33, -35, -45, and -51 did not differ significantly by variant lineage. As the possible exception, VEs against persistent infection and/or CIN 2/3 due to HPV-31 A/B and HPV-31C variants were -7.1% (95% confidence interval [CI], -33.9% to 0%) and 86.4% (95% CI, 65.1%-97.1%), respectively (P = .02 for test of equal VE). No difference in VE was observed by variant among transient HPV-31 infections (P = .68).

CONCLUSIONS

Overall, sequence variation at the variant level does not appear to explain partial cross-protection by the bivalent HPV vaccine.

Naturally Occurring Capsid Protein Variants of Human Papillomavirus Genotype 31 Represent a Single L1 Serotype

We investigated naturally occurring variation within the major (L1) and minor (L2) capsid proteins of oncogenic human papillomavirus (HPV) genotype 31 (HPV31) to determine the impact on capsid antigenicity. L1L2 pseudoviruses (PsVs) representing the three HPV31 variant lineages, variant lineages A, B, and C, exhibited comparable particle-to-infectivity ratios and morphologies. Lineage-specific L1L2 PsVs demonstrated subtle differences in susceptibility to neutralization by antibodies elicited following vaccination or preclinical L1 virus-like particle (VLP) immunization or by monoclonal antibodies; however, these differences were generally of a low magnitude. These data indicate that the diagnostic lineage-specific single nucleotide polymorphisms within the HPV31 capsid genes have a limited effect on L1 antibody-mediated neutralization and that the three HPV31 variant lineages belong to a single L1 serotype. These data contribute to our understanding of HPV L1 variant antigenicity.

IMPORTANCE

The virus coat (capsid) of the human papillomavirus contains major (L1) and minor (L2) capsid proteins. These proteins facilitate host cell attachment and viral infectivity and are the targets for antibodies which interfere with these events. In this study, we investigated the impact of naturally occurring variation within these proteins upon susceptibility to viral neutralization by antibodies induced by L1 VLP immunization. We demonstrate that HPV31 L1 and L2 variants exhibit similar susceptibility to antibody-mediated neutralization and that for the purposes of L1 VLP-based vaccines, these variant lineages represent a single serotype.

Identification of the neutralizing epitopes of Merkel cell polyomavirus major capsid protein within the BC and EF surface loops

Merkel cell polyomavirus (MCPyV) is the first polyomavirus clearly associated with a human cancer, i.e. the Merkel cell carcinoma (MCC). Polyomaviruses are small naked DNA viruses that induce a robust polyclonal antibody response against the major capsid protein (VP1). However, the polyomavirus VP1 capsid protein epitopes have not been identified to date. The aim of this study was to identify the neutralizing epitopes of the MCPyV capsid. For this goal, four VP1 mutants were generated by insertional mutagenesis in the BC, DE, EF and HI loops between amino acids 88-89, 150-151, 189-190, and 296-297, respectively. The reactivity of these mutants and wild-type VLPs was then investigated with anti-VP1 monoclonal antibodies and anti-MCPyV positive human sera. The findings together suggest that immunodominant conformational neutralizing epitopes are present at the surface of the MCPyV VLPs and are clustered within BC and EF loops.

Amino acid sequence diversity of the major human papillomavirus capsid protein: implications for current and next generation vaccines

•

We evaluated amino acid diversity of the major capsid protein of HPV.

•

Residues displaying high entropy were found within surface-exposed domains.

•

We discuss the implications of this diversity on the current and next generation HPV vaccines.

Despite the fidelity of host cell polymerases, the human papillomavirus (HPV) displays a degree of genomic polymorphism resulting in distinct genotypes and intra-type variants. The current HPV vaccines target the most prevalent genotypes associated with cervical cancer (HPV16/18) and genital warts (HPV6/11). Although these vaccines confer some measure of cross-protection, a multivalent HPV vaccine is in the pipeline that aims to broaden vaccine protection against other cervical cancer-associated genotypes including HPV31, HPV33, HPV45, HPV52 and HPV58. Both current and next generation vaccines comprise virus-like particles, based upon the major capsid protein, L1, and vaccine-induced, type-specific protection is likely mediated by neutralizing antibodies targeting L1 surface-exposed domains.

The aim of this study was to perform an in silico analysis of existing full length L1 sequences representing vaccine-relevant HPV genotypes in order to address the degree of naturally-occurring, intra-type polymorphisms. In total, 1281 sequences from the Americas, Africa, Asia and Europe were assembled. Intra-type entropy was low and/or limited to non-surface-exposed residues for HPV6, HPV11 and HPV52 suggesting a minimal effect on vaccine antibodies for these genotypes. For HPV16, intra-type entropy was high but the present analysis did not reveal any significant polymorphisms not previously identified. For HPV31, HPV33, HPV58, however, intra-type entropy was high, mostly mapped to surface-exposed domains and in some cases within known neutralizing antibody epitopes. For HPV18 and HPV45 there were too few sequences for a definitive analysis, but HPV45 displayed some degree of surface-exposed residue diversity. In most cases, the reference sequence for each genotype represented a minority variant and the consensus L1 sequences for HPV18, HPV31, HPV45 and HPV58 did not reflect the L1 sequence of the currently available HPV pseudoviruses. These data highlight a number of variant amino acid residues that warrant further investigation for vaccine and natural history studies of HPV.

Identification of neutralizing conformational epitopes on the human papillomavirus type 31 major capsid protein and functional implications

The aim of this study was to characterize the conformational neutralizing epitopes of the major capsid protein of human papillomavirus type 31. Analysis of the epitopes was performed by competitive epitope mapping using 15 anti-HPV31 and by reactivity analysis using a HPV31 mutant with an insertion of a seven-amino acid motif within the FG loop of the capsid protein. Fine mapping of neutralizing conformational epitopes on HPV L1 was analyzed by a new approach using a system displaying a combinatorial library of constrained peptides exposed on E. coli flagella. The findings demonstrate that the HPV31 FG loop is dense in neutralizing epitopes and suggest that HPV31 MAbs bind to overlapping but distinct epitopes on the central part of the FG loop, in agreement with the exposure of the FG loop on the surface of HPV VLPs, and thus confirming that neutralizing antibodies are mainly located on the tip of capsomeres. In addition, we identified a crossreacting and partially crossneutralizing conformational epitope on the relatively well conserved N-terminal part of the FG loop. Moreover, our findings support the hypothesis that there is no correlation between neutralization and the ability of MAbs to inhibit VLP binding to heparan sulfate, and confirm that the blocking of virus attachment to the extracellular matrix is an important mechanism of neutralization.

Type-specific and cross-reactive antibodies induced by human papillomavirus 31 L1/L2 virus-like particles

The aim of this study was to determine whether antibodies induced by human papillomavirus (HPV) type 31 L1/L2 virus-like particles (VLPs) could cross-react with VLPs of the closely related HPV-16 and distantly related HPV-11, and to investigate the potential role of the L2 protein in L1/L2 VLPs in inducing cross-neutralizing antibodies. Antisera were prepared from rabbits immunized with intact or denatured HPV-31 L1/L2 VLPs. Cross-reaction and cross-neutralization were analysed by Western blotting and ELISA, and by haemagglutination inhibition, respectively. Western blotting results showed that H31 L1/L2 (D) antiserum (antiserum from a rabbit immunized with denatured HPV-31 L1/L2 VLPs) could cross-react with the L1 protein of HPV-11 and -16. HPV-31 L1/L2 VLP antiserum showed strong cross-reaction with and cross-neutralization of HPV-16 VLPs, but this was significantly less with HPV-11 VLPs. In addition, the cross-neutralizing activity against HPV-16 L1/L2 VLPs was slightly higher than that against HPV-16 L1 VLPs, although the difference was not statistically significant. Epitope-blocking ELISA showed that mAb H16.V5 could partially inhibit the cross-reaction of HPV-31 L1/L2 VLP antiserum with HPV-16 L1/L2 VLPs. These results suggested that (i) H31 L1/L2 (D) antiserum could cross-react with L1 protein from both closely related and distantly related HPV types, but HPV-31 L1/L2 VLP antiserum could only cross-neutralize closely related HPV types, (ii) surface-exposed epitopes of the L2 protein in L1/L2 VLPs may induce only low titres of cross-neutralizing antibodies and (iii) certain epitopes that cross-reacted with HPV-31 L1/L2 VLP antiserum are located close to the epitopes recognized by mAb H16.V5. These findings may provide suggestions for the design of multivalent VLP vaccines.