Skip to main content

Association of naturally acquired type-specific HPV antibodies and subsequent HPV re-detection: systematic review and meta-analysis

Abstract

Background

Understanding the role of naturally acquired (i.e., infection-induced) human papillomavirus (HPV) antibodies against reinfection is important given the high incidence of this sexually transmitted infection. However, the protective effect of naturally acquired antibodies in terms of the level of protection, duration, and differential effect by sex remains incompletely understood. We conducted a systematic review and a meta-analysis to (1) strengthen the evidence on the association between HPV antibodies acquired through past infection and subsequent type-specific HPV detection, (2) investigate the potential influence of type-specific HPV antibody levels, and (3) assess differential effects by HIV status.

Methods

We searched Embase and Medline databases to identify studies which prospectively assessed the risk of type-specific HPV detection by baseline homologous HPV serostatus among unvaccinated individuals. Random-effect models were used to pool the measures of association of naturally acquired HPV antibodies against subsequent incident detection and persistent HPV positivity. Sources of heterogeneity for each type were assessed through subgroup analyses stratified by sex, anatomical site of infection, male sexual orientation, age group, and length of follow-up period. Evidence of a dose-response relationship of the association between levels of baseline HPV antibodies and type-specific HPV detection was assessed. Finally, we pooled estimates from publications reporting associations between HPV serostatus and type-specific HPV detection by baseline HIV status.

Results

We identified 26 publications (16 independent studies, with 62,363 participants) reporting associations between baseline HPV serostatus and incident HPV detection, mainly for HPV-16 and HPV-18, the most detected HPV type. We found evidence of protective effects of baseline HPV seropositivity and subsequent detection of HPV DNA (0.70, 95% CI 0.61–0.80, NE = 11) and persistent HPV positivity (0.65, 95% CI 0.42–1.01, NE = 5) mainly for HPV-16 among females, but not among males, nor for HPV-18. Estimates from 8 studies suggested a negative dose–response relationship between HPV antibody level and subsequent detection among females. Finally, we did not observe any differential effect by baseline HIV status due to the limited number of studies available.

Conclusion

We did not find evidence that naturally acquired HPV antibodies protect against subsequent HPV positivity in males and provide only modest protection among females for HPV-16. One potential limitation to the interpretation of these findings is potential misclassification biases due to different causes.

Introduction

Human papillomavirus (HPV) is one of the most common sexually transmitted viruses, affecting 50% of sexually active individuals at least once in their lifetime [1]. Infection with HPV is associated with 7–8% of all human malignancies and accountable for 96% of cervical cancer, 93% of anal cancers, 64% of vaginal cancers, 51% of vulvar cancer, 36% of penile cancers, and 63% of oropharyngeal carcinomas [2]. Among all high-risk HPV oncogenic types (HR-HPV), HPV-16 and HPV-18 account for 60–80% of all cervical cancers, one of the most commonly diagnosed cancers and a leading cause of cancer death in women in many low-income countries [3,4,5,6,7].

Most HPV infections clear within 1–2 year by cell-mediated immune response and generation of serum neutralizing antibodies (IgG) against the capsid L1 protein of HPV [1, 8]. However, while most of infections are cleared, some persist for years, which can lead to cell abnormalities and potentially to cancer, if not promptly diagnosed and adequately treated [9]. Studies have shown the potential protective effect of neutralizing antibodies against subsequent infections when the immune response is initiated by HPV vaccines [10]. Currently, several prophylactic vaccines are available, which include: bivalent vaccines that offer protection against HPV-16/18, quadrivalent vaccines against HPV-6/11/16/18, and nonavalent vaccines against HPV-6/11/16/18/31/33/45/52/58 [11]. The high protection (> 90%) against vaccine-targeted HPV-type incident infections, and HPV-related abnormalities and precancerous lesions, have been demonstrated by clinical studies [10, 12]. However, the extent to which naturally acquired, or infection-induced, antibodies can help prevent HPV reinfection remains poorly understood, especially with respect to the level and the duration of protection. One among many challenges in investigating the effect of naturally acquired antibodies against subsequent HPV infection is the role of latency and deposition, and consequently, the difficulties around determining whether the incident HPV DNA detection truly represents true re-infection, or rather a deposition from sexual partners and/or from cross-site contamination, or reactivation of the virus [13].

The evidence base regarding risk of subsequent infections with HR-HPV is mixed. Several studies found no effect of naturally-acquired antibodies in reducing the risk of subsequent HPV detection in men [14,15,16,17] and in women [18,19,20] while others found some levels of protection among women [21,22,23,24,25,26,27,28]. These mixed findings could be attributed to varying study designs, heterogeneous study populations, and different laboratory protocols used to confirm HPV detection, while small sample sizes might have increased the uncertainty of estimates. The most recent meta-analysis of naturally acquired HPV antibodies against subsequent genital HPV infection, published in 2016, has estimated a 30–35% reduction in the risk of subsequent HPV-16/18 genital infection among women, but not among men [29]. However, the study was not able to explore sources of heterogeneity aside from sex and it did not assess the influence of HPV antibody levels or HIV status on type-specific re-infection rates.

Current evidence points to multiple interactions between HIV and HPV infections, due to similar risk factors, and biological and immunological factors [30]. Studies have suggested that the risk of acquiring HPV infection, persistent infection, and disease progression is increased among people living with HIV (PLHIV). However, the impact of HIV infection on the natural immune response to HPV following natural infection is less understood [31, 32]. Acquiring quantitative estimates of the risk of subsequent HPV infection among PLHIV is particularly important for assessing the impact of HPV vaccination programs in settings where HIV prevalence is highest and informing parameterization of mathematical models to understand the epidemiology of HPV infection and cervical cancer in high-risk populations.

The objectives of this study were threefold. First, we updated the existing evidence on the association between naturally acquired immunity and subsequent type-specific HPV detection through a systematic review and meta-analysis. We also performed several subgroup analyses to further investigate sources of heterogeneity in the association between natural immunity and risk of subsequent HPV detection (stratified by anatomic site of detection, sex, sexual orientation, age group). Second, we investigated the influence of HPV antibody levels to type-specific HPV detection risk. Third, we examined if HPV detection risk differ by HIV status. We performed analyses on all HPV types where possible but reported results for HPV-16 and HPV-18 in the main results.

Methods

Search strategy and selection criteria

All results of this study are reported according to the PRISMA guidelines (Additional file 1: Table S1). We searched the Embase and MEDLINE databases for articles in English, French, and Japanese, published up to May 2022. We used key search terms related to four domains—HPV, study design, antibodies, and viral DNA—to capture prospective studies assessing an association between HPV serostatus at baseline (history of previous HPV infection) and subsequent type-specific HPV detection among HPV-unvaccinated individuals (HPV DNA detection at follow-up among DNA-negative participants at baseline; see Additional file 1: Table S2).

Two reviewers (KY and KG) independently screened titles and abstracts of potentially eligible articles after removing duplicate records and resolving any discrepancies in the selection. Full texts of all potentially relevant articles were then screened. We excluded publications that did not report estimates of association between baseline HPV serostatus and type-specific HPV incident detection at follow-up, or that did not provide sufficient data to derive them. Finally, the reference lists of all included publications were screened to identify additional relevant studies [29].

Data extraction

Data were retrieved by two reviewers (KY and KG) and any discrepancies were resolved by consensus. We extracted the reported crude and adjusted estimates, or the data to derive it, of any measures of association between type-specific serostatus of any reported HPV-type and incident detection of HPV DNA at follow-ups. These include the incidence rate ratio (IRR), cumulative incidence ratio or risk ratio (RR), hazards ratio (HR), or the odds ratio (OR) and their 95% confidence intervals (95%CI) or the data to derive unadjusted RR or IRR (herein referred to as self-calculated) from reported counts and/or incidence rates (Additional file 1: Text 2). We also extracted information on participants characteristics and HPV detection (i.e., sex, age, population type, HIV status, HPV type, infection site), study characteristics (i.e., country, study design, sample size, follow-up duration), and quality indicators (e.g., type of tests used for HPV serology and HPV DNA detection, variables adjusted for). For the third objective, we also extracted information on HIV disease stage, ART or other treatment status, and CD4+ cell counts of individuals living with HIV where possible.

We calculated the pooled measure of association (on the relative risk scale) and 95%CI using inverse variance weights and the DerSimonian-Laird random-effect method on the logarithmic scale (results are then presented on the original scale) [33]. To maximize the number of estimates, we pooled all measures of association (IRR, RR, HR, OR) in our main meta-analysis but assessed their influence in subgroup analyses. The I2 statistic was used to assess heterogeneity between study estimates [34]. To avoid including duplicated estimates on same participants, we only selected one estimate from studies reporting on the same study population (e.g., HIM cohort study) for each HPV type and infection site. Finally, estimates that were adjusted for potential confounders were chosen over unadjusted estimates, if available, for the pooled analysis. However, we also compared crude and adjusted estimates in a sensitivity analysis described below.

For the first objective, we used two different outcomes: (1) incident type-specific HPV detection, defined as the first detection of homologous HPV DNA at follow-up among HPV DNA-negative individuals at baseline (as defined by included studies), and (2) persistent HPV positivity, defined as the detection of homologous HPV DNA at two or more consecutive follow-up visits within a 6- or 12-month interval among HPV DNA-negative individuals at baseline. For both outcomes, we focused on HPV-16 and HPV-18 in our main meta-analyses, since these HPV types were the most commonly reported [35]. We presented results on other HPV types for the first outcome in Additional file 1 for completion. For the outcome of incident HPV detection, we also performed subgroup analyses to investigate whether heterogeneity between study estimates could be explained by sex, HPV detection site (cervical/cervicovaginal, penile/scrotal, anal, oral), sexual orientation for male (men who have sex with men [MSM], men who exclusively have sex with women [MSW]), and age group (≤ 30, >30 years). Subgroup analyses were restricted to HPV-16 and HPV-18 only, based on the availability of extracted data from studies. Differences between subgroups were tested using the Q-test or subgroup differences (two-sided, 5% alpha level). We did not perform subgroup analyses for the outcome of persistent infection due to low number of estimates available for stratification. We also performed sensitivity analyses to assess study quality and potential biases due to heterogeneity in study designs (described in the next section).

To assess evidence of a dose–response relationship between baseline type-specific HPV antibodies and same type HPV detection, we plotted the measures of association by levels of HPV antibodies (i.e., terciles, quartiles, low versus high) from publication reporting this data for HPV-16 and HPV-18. We also performed meta regression using the estimates from the publications above as individual data point. Given heterogeneity in the categorization of antibody titres, we compared and pooled the risk of detection in the highest category (i.e., “high”, “highest quartile”, “highest tercile”) to the lowest one (i.e., “low”, “lowest quartile”, “lowest tercile”).

Finally, for the third objective, we pooled estimates reporting associations between HPV serostatus and type-specific HPV detection by baseline HIV status. Differences between HIV positive and HIV negative individuals were assessed using the Q statistic. Again, these analyses were restricted to HPV-16 and HPV-18 only, based on the availability of extracted data from studies.

Sensitivity analyses

We tested the influence of each individual study estimate on the overall pooled estimates using leave-one-out sensitivity analyses. We also performed several sensitivity analyses to assess whether study estimates varied due to potential biases. Specifically, we compared (1) crude versus adjusted estimates; (2) reported versus estimates calculated from raw data reported in publications; (3) types of measures of associations (IRR, HR, RR, OR); (4) methods used for serotyping (neutralizing versus non-neutralizing assay); (5) primers set used for HPV-DNA detection (by decreasing test sensitivity: PGMY09/11, SPF10, MY09/11) [36]; and (6) lengths of follow-up periods. For the first sensitivity analysis, we used two different approaches: we compared crude versus adjusted estimates from all studies, as well as crude versus adjusted estimates in the subset of studies reporting both.

Study quality was assessed using a modified version of the Newcastle–Ottawa-Scale (NOS) for cohort studies (Additional file 1: Text 3) [37]. Scores were assigned on a continuous scale from 0 to 12 stars based on 8 criteria including the strength of study design for an inference, methods to assess the exposure/outcome, etc. Two independent reviewers assigned the scores independently, and any discrepancies were discussed to reach an agreement. We then categorized the studies into low, medium, and high quality for scores ranging from 0 to 4, 5 to 8, and 9 to 12, respectively. Finally, the risk of publication bias was further examined by funnel plots and Egger’s tests of symmetry [38]. We performed all analyses with R (version 4.0.3), using the “meta” package [39].

Results

Search results

The search strategy yielded 6474 records and an additional 129 publications were identified through screening of the reference lists of included publications (Fig. 1). After removal of duplicates, a total of 5453 publications were assessed for eligibility through title and abstract screening, of which 45,224 were excluded for non-relevance. In total, 26 publications reporting on 16 independent studies (n = 62,363 participants) met our eligibility criteria and were included in our study, which adds 12 more publications to the most recent meta-analysis. These publications reported a total of 92 different estimates of the associations between HPV serostatus at baseline and incident type-specific HPV detection for any HPV types reported.

Fig. 1
figure 1

Flow diagram of the study selection for the systematic review and meta-analysis of the association between baseline HPV serostatus and type-specific HPV re-detection

Study characteristics

Of the 16 included independent studies, 9 were prospective observational cohort studies and 7 were randomized clinical trials on HPV vaccines (Table 1). HPV-16 (number of publications [Np] = 26) and HPV-18 (Np = 14) were the two most common HR-HPV types documented in these studies (for other types, see Table 1). Eighteen publications reported estimates of incident detections and 4 publications on 6-month and 12-month persistent detections. Most publications (Np = 19) focused on females and HPV detection at cervical or cervicovaginal sites, while 6 publications focused on males, and 1 publication included both sexes as participants. Most publications reported on multi-country studies but only one study included cohorts from sub-Saharan Africa (Burkina Faso, South Africa). The most common infection sites assessed among males were penile and/or scrotal sites (Np = 4), followed by oral (Np = 3) and anal (Np = 2) sites. A total of 8 publications also provided estimates by baseline levels of HPV antibodies. One publication was restricted to HIV-positive females, and another four reported estimates stratified by HIV status among females or both sexes. SPF10 was the most common DNA primer set used for HPV DNA detection (Np = 14), and most of the studies used non-neutralizing serologic assay for the detection of serum HPV antibodies (Np = 24).

Table 1 Summary of study characteristics of publications included in the meta-analysis of the association between baseline HPV serostatus and type-specific HPV detection (Np = 26)

Four publications were classified as having “high” study quality, and 22 as “medium” study quality, based on the criteria regarding the assessment of exposure, comparability, and outcome assessment (Additional file 1: Table S3). There was no single unique criterion that determined the “high” quality of studies; however, one of the key criteria was the proportion of participants lost to follow-up.

Association between HPV serostatus and type-specific incident HPV detection

As the direction of association varied greatly by sex, the overall pooled estimates are difficult to interpret (Tables 2, 3) and we only presented overall estimates for completeness. For females, the risk of incident detection among baseline-HPV-16 seropositive was reduced compared to those baseline-seronegative females (0.70, 95% CI 0.61–0.80, NE = 11). In reverse, for males, the risk of incident detection among baseline-HPV seropositive individuals was increased compared to those baseline-seronegative (HPV-16: 1.43, 95% CI 1.12–1.83, NE = 5) (Table 2, Fig. 2). The results for HPV-18 were suggestive of similar direction of association but with wider confidence intervals (Females: 0.97, 95% CI 0.83–1.13, NE = 7; Males: 1.15, 95% CI 0.72–1.83, NE = 3) (Table 3, Fig. 2. For other types, see Additional file 1: Figs. S1, S2).

Table 2 Subgroup analyses of the association between baseline HPV-16 serostatus and HPV-16 incident detection
Table 3 Subgroup analyses of the association between baseline HPV-18 serostatus and HPV-18 incident detection
Fig. 2
figure 2

Forest plot of the association between baseline HPV serostatus and type-specific HPV incident detection by sex for A HPV-16, and B HPV-18. *Indicates estimates that are self-calculated using the provided data. **Indicates estimates that are adjusted for confounders

Consistent with observed sex differences, we found a negative association for HPV-16 incident detection when restricting to cervical/cervicovaginal sites (0.70, 95% CI 0.44–1.12, NE = 11). Among males, the associations remained positive when restricting to specific sites: penile/scrotal (1.36; 95% CI 1.01–1.83, NE = 2), anal (1.60, 95% CI 0.99–2.56, NE = 2), and oral sites (1.30, 95% CI 0.60–2.84, NE = 2) (Table 2, Additional file 1: Fig. S4). However, differences by site among males were inconclusive for HPV-18 and with large confidence intervals (Table 3).

We did not find any significant differences in the direction and/or magnitude of association by sexual orientation for males, nor by age groups across study estimates of the associations between baseline HPV serostatus and type-specific incident detection for both HPV-16 and HPV-18 (Tables 2, 3, Additional file 1: Fig. S5).

Association between HPV serostatus and type-specific persistent HPV positivity

Among females, we found some evidence of a protective effect of baseline antibody on HPV-16 persistent detection (6 months: 0.58, 95% CI 0.44–0.77, NE = 4; 12 months: 0.68, 95% CI: 0.56–0.83, NE = 5) (Tables 2, 3). Only one study estimate on persistent infection was available for males for HPV-16 (1.39, 95% CI 0.80–2.41) and HPV-18 (0.19, 95% CI 0.03–1.28) for the 6-month interval.

Association between HPV antibody level and incident HPV detection

A total of eight publications reported estimates of association by antibody levels for HPV-16 and HPV-18 separately for males and females. Compared to the lowest antibody titres categories, the highest ones were associated with reduced likelihood of HPV detection among females for HPV-16 (0.65, 95% CI 0.55–0.76, NE = 5) and for HPV-18 (0.55, 95% CI 0.41–0.75, NE = 4). Out of nine reported estimates of association among female participants for incident detection at cervical/cervicovaginal site, 4 estimates for HPV-16 and 4 estimates for HPV-18 showed negative associations between increasing baseline antibody levels and the likelihood of HPV detection. We found signs of a positive association between the highest antibody threshold categories and HPV detection risk among males for HPV-16 (1.93, 95% CI 0.67–5.57, NE = 2), however, the results were inconclusive. There was only one study reporting males’ estimates by tertiles for HPV-18 (Fig. 3, Additional file 1: Fig. S6). Finally, one estimate of association among both male and female participants combined, measured at oral sites for HPV-16 suggested a negative trend (Fig. 3).

Fig. 3
figure 3

Trends in estimates of association between baseline HPV serostatus and type-specific HPV incident detection by HPV antibody concentration level, stratified by sex (females, males, or both). Seven publications measured antibody levels in terms of ELISA units per mL (EU/mL) [14, 19, 21, 22, 24, 25, 41], and one, in Luminex units per mL (LU/mL) [15]. Five studies reported estimates by tertiles of antibody levels [14, 15, 22, 24, 41], two by quartiles [21, 24], and one dichotomized antibody levels (low versus high concentration) [19]. T1, T2, T3 refer to the lowest, middle, and upper tertiles of antibody levels used by each publication, respectively; Q1, Q2, Q3, Q4 refer to the first, second, third, and fourth quartiles, respectively; Low and High refer to low and high antibody concentrations, respectively. Horizontal dashed lines indicate the null value of the estimate of the association (1 on a relative scale). The y-axis indicates the risk of incident HPV detection on the relative scale. For instance, a blue dot with the highest y-value in the top right panel indicates a risk ratio of incident detection between individuals with no HPV antibodies versus individuals with naturally acquired HPV antibodies in the lowest tertile. The lower end of antibody level for the lowest tercile/quartile used by publications was between 7–8 EU/mL or 0.2–0.28 LU/mL, and the higher end of between 40–64 EU/mL or 0.46–0.55 LU/mL

Impact of HIV status on HPV serostatus and type-specific incident HPV detection

Results from publications assessing the associations between HPV serostatus and incident type-specific HPV detection by baseline HIV status were inconclusive. For HPV-16, the estimate of association between baseline serostatus and subsequent incident HPV detection was 1.68 (95% CI 0.63–4.46, NE = 1) among female individuals without HIV and 0.97 (95% CI 0.64–1.48, NE = 3) among female individuals living with HIV. For HPV-18, the results were also not statistically significant (female individuals living without HIV: 0.15, 95% CI 0.02–1.16, NE = 1; women living with HIV: 0.41, 95% CI 0.08–2.03, NE = 2) (Fig. 4, Tables 2, 3). Only one study reported estimates for male individuals living without HIV (HPV-16: 1.70, 95% CI 0.79–3.66; HPV-18: 1.10, 95% CI 0.40–3.06) and for male individuals living with HIV (HPV-16: 1.00, 95% CI 0.44–2.29; HPV-18: 2.30, 95% CI 0.73–7.22) (Table 2). Included studies did not provide enough information about the HIV disease stage or treatment status to perform subgroup analyses.

Fig. 4
figure 4

Forest plot of the association between baseline HPV serostatus and type-specific HPV incident detection by HIV status for A HPV-16 and B HPV-18. †Indicates estimates that are adjusted for confounders; NA = Not available. Participants included in the publication by Kelly et al. [40] were restricted to women living with HIV. The other publications reported separate estimates of measures of association by HIV status. CD4 cell counts of HIV positive individuals were reported in varying format by publications (threshold, median). All 5 publications adjusted for the baseline CD4 cell counts in the calculation of estimates in either main analysis or sensitivity analysis. Only two publications reported information regarding treatment status of HIV positive individuals at baseline and/or throughout the study period

Sensitivity analyses

For both HPV-16 and HPV-18 incident detections, we did not observe large differences between crude and adjusted estimates, methods used for serotyping, and by study quality, after stratifying by sex (Tables 2, 3, Additional file 1: Figs. S7, S8). The only difference was found in the primer sets used for HPV DNA detection among females for HPV-16 (PGMY09/11: 0.61, 95% CI 0.47–0.79, NE = 2; SPF10: 0.70, 95% CI 0.61–0.81, NE = 5; MY09/11: 0.91, 95% CI 0.67–1.22, NE = 3).

Leave-one-out analyses did not identify any individual study estimate that strongly impacted the pooled estimates of the association or the measured heterogeneity for both HPV-16 and HPV-18 (Additional file 1: Fig. S9). Finally, funnel plots for all HPV types as well as Egger’s test for HPV-16 and HPV-18 showed no evidence of publication bias (HPV-16 Egger P = 0.17; HPV-18 Egger P = 0.78; Additional file 1: Fig. S10).

Discussion

Pooling data on 62,363 participants, this systematic review found evidence of protective effects of baseline HPV seropositivity and subsequent detection of incident and persistent HPV, mainly for HPV-16 among females, but not among males, nor for HPV-18 or other HPV types. Moreover, our results suggested a protective effect of higher antibody levels against subsequent HPV-16 detection at cervical/cervicovaginal sites among female subjects, but not among males. Based on the 5 published studies available, there was no evidence of a differential impacts of naturally acquired antibodies on subsequent HPV-16/18 detection by HIV status.

Even with increasing uptake and coverage of HPV vaccines, improving our understanding of the heterogenous effects of naturally acquired antibodies on reinfection and persistent infection is important, especially to improve structural assumptions and parameterization of mathematical models of HPV, for instance. Our results suggest a partially protective effect of natural antibodies among females, but this level of protection is much lower than the one conferred by vaccines. For example, trial data show that a single dose of HPV 16/18 vaccine results in antibody geometric mean titres 5–9 times larger than the protection acquired from natural infection, even several years after vaccination [42].

Consistent with previous studies, the current meta-analysis found important sex differences [29]. The finding that, among males, HPV seropositivity increases risks of subsequent HPV detection needs to be interpreted cautiously. First, it is possible that the estimates are confounded by partially or unmeasured variables, such as sexual behaviours, as some of our included estimates were unadjusted in the original studies. Further, differences in immune response in extracervical sites, and its potential in explaining the sex differentials, requires further research. Several studies have noted differences in the level of immune response introduced after an infection at a mucosal epithelium of female genitals compared to the immune response after an infection at keratinized epithelium of male genitals (i.e., penis), suggesting a higher viral antigen level and a stronger antibody response at mucosa [60, 61]. MSM could be more frequently exposed to HPV through anal mucosa than men who have sex with women. However, our site-specific analyses do not show a reduced risk of subsequent HPV detection at anal sites in stratified analyses among MSM.

Besides a thorough examination of incident detections, we examined the impact of naturally acquired antibody on persistent HPV positivity. Similar to the association between baseline HPV serostatus and a subsequent incident detection (one time detection of HPV-DNA), there was a reduction in risk of persistent HPV positivity (two or more detections of HPV-DNA during consecutive follow-ups). Current knowledge points to multiple interactions between HPV infection and HIV status [30, 31, 40, 43, 44]. Our results remained inconclusive with wide confidence intervals in terms of the association between HPV antibodies and subsequent HPV detection by HIV status. Results stratified by HIV status were subject to limited precision due to a small number of estimates available for pooling as well as sources of unmeasured heterogeneity, such as antiretroviral treatment status and CD4 cell counts among PLHIV, that remained unadjusted due to low number of studies that measured these variables, which might have contributed to residual confounding.

Assessing the association between naturally acquired HPV antibodies and subsequent re-infection, as well as synthesizing the evidence of this association, are subject to challenges and limitations. First, it is challenging to differentiate between true re-infection and detection of an activated latent HPV or deposition when HPV DNA is detected [45,46,47]. Therefore, our results from pooled analysis require cautious interpretations when inferring about the role of natural history of HPV. Moreover, there is always a risk of misclassification of exposure and/or outcomes, depending on the type of serologic assay (i.e., neutralizing versus non-neutralizing) and DNA primer sets used for individual studies. Variability in types of serological tests or DNA primer sets used and difference in sensitivity by HPV type may have influenced the risk of misclassification in individual study results and introduced some degree of heterogeneity in the magnitude of the pooled estimates (rather than the direction of the association). A careful interpretation of the effect sizes while contextualizing the heterogeneity in study designs, population, and sensitivity of methodology used is required. Second, our analysis was restricted to type-specific risk of an incident detection and did not assess the effect of cross-protection of antibodies across different HPV types. Studies have found that infections with multiple HPV types was associated with persistence of HPV infections, while others did not find any differences [48, 49]. Third, some of our analyses had limited statistical power due to a small number of publications reporting estimates among men, MSM, by HIV serostatus, and age. Lastly, our study could not rule out the potential information bias and selection bias, as some publications did not always report characteristics of study populations such as mean/median age in a standardized format, nor participation rate/loss to follow up. Our estimates are also subject to confounding, as we pooled adjusted and unadjusted estimates provided by individual studies to increase the number of effect sizes pooled.

This review synthesises current knowledge on naturally acquired immunity against HPV. It substantially increases the precision of effect sizes estimates, especially for HPV-16. Compared to the previous review, we added 12 additional publications and nearly 35,000 new participants to the pooled analysis. Our results suggest differences in the effect of HPV antibodies between males and females and a potential negative dose–response relationship between antibody titre levels and subsequent detection for HPV-16 and HPV-18 in females. Although our results were inconclusive regarding differential impacts of naturally acquired immunity by HIV status, PLHIV could be at greater risk of reinfection, and consequently, at higher risk of developing HPV-attributable cancers.

Conclusion

We did not find evidence that naturally acquired HPV antibodies protect against subsequent HPV positivity in males and provide only modest protection among females for HPV-16. Being the first study to evaluate the dose–response relationship between antibody titre levels and subsequent detection of HPV-DNA, we found some evidence of a potential negative dose–response. Finally, we did not find conclusive evidence of differential impacts of naturally acquired antibodies by HIV status, and further investigation is warranted.

Availability of data and materials

All relevant data are included in the article or uploaded as additional information.

References

  1. Brianti P, De Flammineis E, Mercuri SR. Review of HPV-related diseases and cancers. New Microbiol. 2017;40(2):80–5.

    CAS  PubMed  Google Scholar 

  2. Cobos C, et al. The role of human papilloma virus (HPV) infection in non-anogenital cancer and the promise of immunotherapy: a review. Int Rev Immunol. 2014;33(5):383–401.

    Article  CAS  PubMed  Google Scholar 

  3. Black E, Richmond R. Prevention of cervical cancer in Sub-Saharan Africa: the advantages and challenges of HPV vaccination. Vaccines (Basel). 2018;6(3):61.

    Article  PubMed  Google Scholar 

  4. World Health O. WHO guideline for screening and treatment of cervical pre-cancer lesions for cervical cancer prevention. 2nd ed. Geneva: World Health Organization; 2021.

    Google Scholar 

  5. Crosbie EJ, et al. Human papillomavirus and cervical cancer. Lancet. 2013;382(9895):889–99.

    Article  PubMed  Google Scholar 

  6. Carter JR, Ding Z, Rose BR. HPV infection and cervical disease: a review. Aust N Z J Obstet Gynaecol. 2011;51(2):103–8.

    Article  PubMed  Google Scholar 

  7. Arbyn M, et al. Estimates of incidence and mortality of cervical cancer in 2018: a worldwide analysis. Lancet Glob Health. 2020;8(2):e191–203.

    Article  PubMed  Google Scholar 

  8. Rosalik K, Tarney C, Han J. Human papilloma virus vaccination. Viruses. 2021;13(6):1091.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Mboumba Bouassa RS, et al. Natural and vaccine-induced B cell-derived systemic and mucosal humoral immunity to human papillomavirus. Expert Rev Anti Infect Ther. 2020;18(6):579–607.

    Article  CAS  PubMed  Google Scholar 

  10. Cheng L, Wang Y, Du J. Human papillomavirus vaccines: an updated review. Vaccines. 2020;8(3):391.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. de Oliveira CM, Fregnani J, Villa LL. HPV vaccine: updates and highlights. Acta Cytol. 2019;63(2):159–68.

    Article  PubMed  Google Scholar 

  12. Garland SM, et al. Impact and effectiveness of the quadrivalent human papillomavirus vaccine: a systematic review of 10 years of real-world experience. Clin Infect Dis. 2016;63(4):519–27.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Malagón T, et al. Proportion of incident genital human papillomavirus detections not attributable to transmission and potentially attributable to latent infections: implications for cervical cancer screening. Clin Infect Dis. 2021;75:365–71.

    Article  PubMed Central  Google Scholar 

  14. Beachler DC, et al. An examination of HPV16 natural immunity in men who have sex with men (MSM) in the HPV in men (HIM) study. Cancer Epidemiol Biomark Prev. 2018;27(4):496–502.

    Article  CAS  Google Scholar 

  15. Mooij SH, et al. No evidence for a protective effect of naturally induced HPV antibodies on subsequent anogenital HPV infection in HIV-negative and HIV-infected MSM. J Infect. 2014;69(4):375–86.

    Article  PubMed  Google Scholar 

  16. Pamnani SJ, et al. Impact of serum antibodies to HPV serotypes 6, 11, 16, and 18 to risks of subsequent genital HPV infections in men: the HIM study. Cancer Res. 2016;76(20):6066–75.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Pierce Campbell CM, et al. Human papillomavirus (HPV) L1 serum antibodies and the risk of subsequent oral HPV acquisition in men: the HIM study. J Infect Dis. 2016;214(1):45–8.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Viscidi RP, et al. Serum immunoglobulin A response to human papillomavirus type 16 virus-like particles in human immunodeficiency virus (HIV)-positive and high-risk HIV-negative women. J Infect Dis. 2003;188(12):1834–44.

    Article  CAS  PubMed  Google Scholar 

  19. Viscidi RP, et al. Seroreactivity to human papillomavirus (HPV) types 16, 18, or 31 and risk of subsequent HPV infection: results from a population-based study in Costa Rica. Cancer Epidemiol Biomark Prev. 2004;13(2):324–7.

    Article  Google Scholar 

  20. Viscidi RP, et al. Human papillomavirus capsid antibody response to natural infection and risk of subsequent HPV infection in HIV-positive and HIV-negative women. Cancer Epidemiol Biomark Prev. 2005;14(1):283–8.

    Article  CAS  Google Scholar 

  21. Castellsague X, et al. Risk of newly detected infections and cervical abnormalities in women seropositive for naturally acquired human papillomavirus type 16/18 antibodies: analysis of the control arm of PATRICIA. J Infect Dis. 2014;210(4):517–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Lin SW, et al. HPV16 seropositivity and subsequent HPV16 infection risk in a naturally infected population: comparison of serological assays. PLoS ONE. 2013;8(1):e53067.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Rosillon D, et al. Risk of newly detected infections and cervical abnormalities in adult women seropositive or seronegative for naturally acquired HPV-16/18 antibodies. Cancer Med. 2019;8(10):4938–53.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Safaeian M, et al. Risk of HPV-16/18 infections and associated cervical abnormalities in women seropositive for naturally acquired antibodies: pooled analysis based on control arms of two large clinical trials. J Infect Dis. 2018;218(1):84–94.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Safaeian M, et al. Epidemiological study of Anti-HPV16/18 seropositivity and subsequent risk of HPV16 and -18 infections. J Natl Cancer Inst. 2010;102(21):1653–62.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Triglav T, et al. Effect of naturally acquired type-specific serum antibodies against human papillomavirus type 16 infection. J Clin Virol. 2017;90:64–9.

    Article  CAS  PubMed  Google Scholar 

  27. Wentzensen N, et al. A competitive serological assay shows naturally acquired immunity to human papillomavirus infections in the Guanacaste Natural History Study. J Infect Dis. 2011;204(1):94–102.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Wilson L, et al. Seroprevalence of 8 oncogenic human papillomavirus genotypes and acquired immunity against reinfection. J Infect Dis. 2014;210(3):448–55.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Beachler DC, et al. Natural acquired immunity against subsequent genital human papillomavirus infection: a systematic review and meta-analysis. J Infect Dis. 2016;213(9):1444–54.

    Article  CAS  PubMed  Google Scholar 

  30. Looker KJ, et al. Evidence of synergistic relationships between HIV and Human Papillomavirus (HPV): systematic reviews and meta-analyses of longitudinal studies of HPV acquisition and clearance by HIV status, and of HIV acquisition by HPV status. J Int AIDS Soc. 2018;21(6):e25110.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Denslow SA, et al. Incidence and progression of cervical lesions in women with HIV: a systematic global review. Int J STD AIDS. 2014;25(3):163–77.

    Article  PubMed  Google Scholar 

  32. Huchko MJ, et al. The time has come to make cervical cancer prevention an essential part of comprehensive sexual and reproductive health services for HIV-positive women in low-income countries. J Int AIDS Soc. 2015;18(Suppl 5):20282–20282.

    Article  PubMed  PubMed Central  Google Scholar 

  33. DerSimonian R, Kacker R. Random-effects model for meta-analysis of clinical trials: an update. Contemp Clin Trials. 2007;28(2):105–14.

    Article  PubMed  Google Scholar 

  34. Borenstein M, Higgins JPT. Meta-analysis and subgroups. Prev Sci. 2013;14(2):134–43.

    Article  PubMed  Google Scholar 

  35. Baseman JG, Koutsky LA. The epidemiology of human papillomavirus infections. J Clin Virol. 2005;32(Suppl 1):S16-24.

    Article  PubMed  Google Scholar 

  36. Gautam A, et al. Human papillomavirus genome based detection and typing: a holistic molecular approach. Curr Mol Med. 2019;19:237–46.

    Article  CAS  PubMed  Google Scholar 

  37. Wells, G. A., Shea, B., O'Connell, D., Peterson, J., Welch, V., Losos, M., Tugwell, P. The Newcastle-Ottawa Scale (NOS) for assessing the quality of nonrandomised studies in meta-analyses. http://www.ohri.ca/programs/clinical_epidemiology/oxford.asp.

  38. Egger M, et al. Bias in meta-analysis detected by a simple, graphical test. BMJ. 1997;315(7109):629–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. CRAN. meta: general package for meta-Analysis. The comprehensive R archive network. https://cran.r-project.org/web/packages/meta/meta.pdf.

  40. Kelly H, et al. Association of antiretroviral therapy with high-risk human papillomavirus, cervical intraepithelial neoplasia, and invasive cervical cancer in women living with HIV: a systematic review and meta-analysis. Lancet HIV. 2018;5(1):e45–58.

    Article  PubMed  Google Scholar 

  41. Beachler DC, et al. A longitudinal study of human papillomavirus 16 L1, e6, and e7 seropositivity and oral human papillomavirus 16 infection. Sex Transm Dis. 2015;42(2):93–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Kreimer AR, et al. Evidence for single-dose protection by the bivalent HPV vaccine—review of the Costa Rica HPV vaccine trial and future research studies. Vaccine. 2018;36(32, Part A):4774–82.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Lissouba P, Van de Perre P, Auvert B. Association of genital human papillomavirus infection with HIV acquisition: a systematic review and meta-analysis. Sex Transm Infect. 2013;89(5):350–6.

    Article  PubMed  Google Scholar 

  44. Houlihan CF, et al. Human papillomavirus infection and increased risk of HIV acquisition. A systematic review and meta-analysis. AIDS. 2012;26(17):2211–22.

    Article  PubMed  Google Scholar 

  45. Gravitt PE. The known unknowns of HPV natural history. J Clin Investig. 2011;121(12):4593–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Twisk DE, et al. Detection of incident anal high-risk human papillomavirus DNA in men who have sex with men: incidence or reactivation? J Infect Dis. 2018;218(7):1018–26.

    Article  PubMed  Google Scholar 

  47. Ranjeva SL, et al. Recurring infection with ecologically distinct HPV types can explain high prevalence and diversity. Proc Natl Acad Sci U S A. 2017;114(51):13573–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Ho GY, et al. Natural history of cervicovaginal papillomavirus infection in young women. N Engl J Med. 1998;338(7):423–8.

    Article  CAS  PubMed  Google Scholar 

  49. Woodman CB, et al. Natural history of cervical human papillomavirus infection in young women: a longitudinal cohort study. Lancet. 2001;357(9271):1831–6.

    Article  CAS  PubMed  Google Scholar 

  50. Eldridge RC, Pawlita M, Wilson L, Castle PE, Waterboer T, Gravitt PE, Schiffman M, Wentzensen N. Smoking and subsequent human papillomavirus infection: a mediation analysis. Annals Epidemiol. 2017;27(11):724–30.

    Article  Google Scholar 

  51. Robbins HA, Li Y, Porras C, Pawlita M, Ghosh A, Rodriguez AC, Schiffman M, Wacholder S, Kemp TJ, Gonzalez P, Schiller J, Lowy D, Esser M, Matys K, Quint W, van Doorn LJ, Herrero R, Pinto LA, Hildesheim A, Waterboer T, Safaeian M. Glutathione S-transferase L1 multiplex serology as a measure of cumulative infection with human papillomavirus. BMC Infect Dis. 2014;14:120.

    Article  PubMed  PubMed Central  Google Scholar 

  52. Lu B, Hagensee ME, Lee JH, Wu Y, Stockwell HG, Nielson CM, Abrahamsen M, Papenfuss M, Harris RB, Giuliano AR. Epidemiologic factors associated with seropositivity to human papillomavirus type 16 and 18 virus-like particles and risk of subsequent infection in men. Cancer Epidemiol Biomark Prev. 2010;19(2):511–6.

    Article  CAS  Google Scholar 

  53. Lu B, Viscidi RP, Wu Y, Lee JH, Nyitray AG, Villa LL, Lazcano-Ponce E, da Silva RJ, Baggio ML, Quiterio M, Salmeron J, Smith DC, Abrahamsen ME, Papenfuss MR, Stockwell HG, Giuliano AR. Prevalent serum antibody is not a marker of immune protection against acquisition of oncogenic HPV16 in men. Cancer Res. 2012;72(3):676–85.

    Article  CAS  PubMed  Google Scholar 

  54. Herrero R, Wacholder S, Rodríguez AC, Solomon D, González P, Kreimer AR, Porras C, Schussler J, Jiménez S, Sherman ME, Quint W, Schiller JT, Lowy DR, Schiffman M, Hildesheim A, Costa Rica Vaccine Trial Group. Prevention of persistent human papillomavirus infection by an HPV16/18 vaccine: a community-based randomized clinical trial in Guanacaste. Costa Rica Cancer Discov. 2011;1(5):408–19.

    Article  PubMed  Google Scholar 

  55. Konno R, Yoshikawa H, Okutani M, Quint W, Suryakiran P, Lin L, Struyf F. Efficacy of the human papillomavirus (HPV)-16/18 AS04-adjuvanted vaccine against cervical intraepithelial neoplasia and cervical infection in young Japanese women. Hum Vaccin Immunother. 2014;10(7):1781–94.

    Article  PubMed  PubMed Central  Google Scholar 

  56. Yao X, Chen W, Zhao C, Wei L, Hu Y, Li M, Lin Z, Lin B, Liu X, Hong Y, Li Q, Pan Q, Zhang X, Li M, Zhao Y, Zhang L, Xu H, Hu F, Zhao J, Huang Y, Sheng W, Zheng Y, Hu S, Su Y, Huang S, Pan H, Zhao F, Qiao Y, Wu T, Zhang J, Xia N. Naturally acquired HPV antibodies against subsequent homotypic infection: A large-scale prospective cohort study. Lancet Reg Health West Pac. 2021;13:100196.

    Article  PubMed  PubMed Central  Google Scholar 

  57. Beachler DC, Viscidi R, Sugar EA, Minkoff H, Strickler HD, Cranston RD, Wiley DJ, Jacobson LP, Weber KM, Margolick JB, Reddy S, Gillison ML, D’Souza G. A longitudinal study of human papillomavirus 16 L1, e6, and e7 seropositivity and oral human papillomavirus 16 infection. Sex Transm Dis. 2015;42(2):93–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. Szarewski A, Poppe WA, Skinner SR, Wheeler CM, Paavonen J, Naud P, Salmeron J, Chow SN, Apter D, Kitchener H, Castellsagué X, Teixeira JC, Hedrick J, Jaisamrarn U, Limson G, Garland S, Romanowski B, Aoki FY, Schwarz TF, Bosch FX, Harper DM, Hardt K, Zahaf T, Descamps D, Struyf F, Lehtinen M, Dubin G, HPV PATRICIA Study Group. Efficacy of the human papillomavirus (HPV)-16/18 AS04-adjuvanted vaccine in women aged 15–25 years with and without serological evidence of previous exposure to HPV-16/18. Int J Cancer. 2012;131(1):106–16.

    Article  CAS  PubMed  Google Scholar 

  59. Moscicki AB, Ma Y, Farhat S, Darragh TM, Pawlita M, Galloway DA, Shiboski S. Redetection of cervical human papillomavirus type 16 (HPV16) in women with a history of HPV16. J Infect Dis. 2013;208(3):403–12.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Windon MJ, Waterboer T, Hillel AT, Chien W, Best S, Stewart C, Akst L, Troy T, Bender N, Miles B, Ryan WR, Mandal R, Pitman K, Eisele DW, Fakhry C, D’Souza G. Sex differences in HPV immunity among adults without cancer. Hum Vaccin Immunother. 2019;15(7–8):1935–41.

    Article  PubMed  PubMed Central  Google Scholar 

  61. Strickler HD, Kirk GD, Figueroa JP, Ward E, Braithwaite AR, Escoffery C, Drummond J, Goebel B, Waters D, McClimens R, Manns A. HPV 16 antibody prevalence in Jamaica and the United States reflects differences in cervical cancer rates. Int J Cancer. 1999;80(3):339–44.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This study was funded by the Canadian Institutes of Health Research. MM-G research program is supported by a Tier II Canada Research Chair in Population Health Modeling. MC Boily acknowledges funding from the MRC Centre for Global Infectious Disease Analysis (reference MR/R015600/1), jointly funded by the UK Medical Research Council (MRC) and the UK Foreign, Commonwealth & Development Office (FCDO), under the MRC/FCDO Concordat agreement and is also part of the EDCTP2 programme supported by the European Union.

Author information

Authors and Affiliations

Authors

Contributions

KG, TM, MB, MCB, and MMG conceived the study. KG, TM, MCB, and MMG planned analyses. KY and KG developed the search strategy, and independently did all stages of screening and data extraction. KY performed all analyses. KY, KG, MCB and MMG interpreted the results. KY and KG drafted the manuscript. All authors read and approved the final version of the manuscript.

Corresponding author

Correspondence to Mathieu Maheu-Giroux.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

KG reports postdoctoral awards from the Fonds de recherche du Québec—Santé, the Canadian Institutes of Health Research (CIHR), and the CIHR Canadian HIV Trials Network during the conduct of the study.

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.

Association of naturally acquired type-specific HPV antibodies and subsequent HPV re-detection: Systematic review and meta-analysis.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Yokoji, K., Giguère, K., Malagón, T. et al. Association of naturally acquired type-specific HPV antibodies and subsequent HPV re-detection: systematic review and meta-analysis. Infect Agents Cancer 18, 70 (2023). https://doi.org/10.1186/s13027-023-00546-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13027-023-00546-3

Keywords