Background: Evidence suggests that natural antibodies developed after HPV16 infection may protect some women but not men against subsequent HPV16 reacquisition. Less is known whether antibodies developed following HPV16 infection are protective among men who have sex with men (MSM).

Methods: Four hundred seventy-five MSM from the Human Papillomavirus Infection in Men (HIM) study were tested for serum antibodies to HPV16 L1 using enzyme-linked immunosorbent assays, and for anal and genital HPV16 DNA using PCR consensus primer system (PGMY 09/11). Adjusted Cox regression was used to evaluate whether baseline HPV16 seropositivity impacts subsequent genital or anal HPV16 DNA.

Results: The risk of subsequent genital HPV16 [aHR = 1.05, 95% confidence interval (CI) = 0.66–1.68] and anal HPV16 infections among MSM (aHR = 2.34, 95% CI = 0.92–5.98) was similar or nonsignificantly higher in HPV16-seropositive than HPV16-seronegative MSM. The risk of genital HPV16 was also similar between HPV16-seronegative and HPV16-seropositive MSM in the highest tertile of HPV16 antibody levels and when restricting to those with new sex partners during follow-up (P > 0.20). Among the 118 MSM who were HPV16 seropositive, 90% remained HPV16 seropositive up to 4 years later. When tested together, MSM with the highest antibody titers (top tertile) had similar levels to females (mean = 130.3 vs. 134.5 EU/mL, P = 0.84).

Conclusions: Despite years of HPV16 seropositivity persistence and antibody titers comparable with females, this study suggested no evidence of HPV16 natural antibodies protecting against subsequent genital or anal HPV16 infection in MSM.

Impact: This could help partially explain the high incidence of genital and anal HPV16 infection and related anal cancer seen in middle-aged and older MSM. Cancer Epidemiol Biomarkers Prev; 27(4); 496–502. ©2018 AACR.

Human papillomavirus (HPV) can infect and lead to cancer at several anatomic sites, including the cervix, anus, vulva, penis, and oropharynx (1). Although HPV is a common infection (2), most individuals are able to clear or control their infection (3–5). However, there is no evidence of lifelong immunity by HPV infection, as type-specific infections can often reappear in individuals previously infected with HPV (6).

While a majority of unvaccinated females who acquire an HPV infection develop a measurable type-specific serum antibody response against epitopes on the HPV L1 capsid protein (7, 8), a minority of HPV-infected men develop these responses (9–11). While a few studies have suggested a lack of natural immunity in females, a recent systematic review and meta-analysis found evidence that HPV16 IgG L1 seropositivity induced through natural HPV16 infection provides some protection against subsequent acquisition of genital HPV16 infection in females, but not in males (12). Another study suggested that HPV natural immunity may be restricted to females with the highest antibody levels (13), while no study has directly evaluated whether seropositive men reach these higher antibody levels potentially necessary for protection. Several studies in the infectious disease field have observed a more robust immune response among adult females compared with adult males (14). If “HPV natural immunity” only occurs in females, it would suggest that men are particularly prone to HPV reacquisition later in life.

Few studies have examined natural immunity among men who have sex with men (MSM). MSM may be more likely to acquire HPV natural immunity given that they are more likely to be exposed to HPV at an anatomic site with access to the mucosal immune system (the anus) and, as has been shown (10), are more likely to seroconvert after infection than men who have sex with women (MSW). Additionally, HPV16-seropositive MSM have higher antibody titers than HPV16-seropositive MSW (previously unpublished, Supplementary Fig. S1). However, prior studies have not suggested evidence of HPV natural immunity in MSM, but that may be because they were restricted to a limited number of MSM, often only included HIV-infected MSM, and utilized different methods of HPV16 antibody measurement compared with many of the female studies demonstrating natural immunity (15–18). In addition, the previous natural immunity studies in MSW and MSM have not examined the duration of antibody persistence in men, as a rapid seroreversion rate could also explain a potential lack of natural immunity.

Understanding the natural history and immunity of HPV in MSM is particularly important given the relatively low HPV vaccine uptake in men, MSM's lack of protection from herd immunity through female vaccination, and their high incidence of HPV16-associated anal cancer (19). Therefore, we conducted a study examining whether antibodies generated in response to natural HPV 16 infection provides protection against subsequent genital or anal HPV16 infection, and whether HPV16 seropositivity persisted over time among the MSM in the prospective Human Papillomavirus Infection in Men (HIM) study.

Study design and participants

The HIM study is a prospective multinational study of HPV infection in 4,074 men conducted in three locations: Tampa (Florida), Sao Paulo (Brazil), and Cuernavaca (Mexico; refs. 4, 20). Enrollment for the HIM cohort study occurred between June 2005 and September 2009, and participants contributed semiannual follow-up visits for a median of 4.2 years. Eligibility requirements for the HIM study included being a male between the ages 18 and 70 years, no prior diagnosis of penile or anal cancer, and no history of HIV or of HPV vaccination (15). For this analysis, we included all 494 men who ever reported having sex with a man except for 19 men who were anal and genital HPV16 DNA positive at HIM baseline, leaving 475 MSM available for this study. There were 29 MSM who were only genital HPV16 DNA positive at baseline that were excluded from the genital HPV16 analysis, and 127 MSM who were either anal HPV16 DNA positive at baseline or were not tested for anal HPV16 who were excluded from the anal HPV16 analysis.

A previous study examined HPV natural immunity against genital HPV16 among both MSW and 234 MSM from the HIM study (15). This follow-up study further investigates HPV16 natural immunity in MSM in the HIM study by increasing the power to examine the association (includes 475 MSM with longer follow-up), examines natural immunity against anal HPV16, using a standardized L1 VLP ELISA that has been extensively used to characterize antibody responses in women (13, 21), and examines the persistence of HPV16 seropositivity over time among those HPV16 seropositive at baseline.

Laboratory testing

Anal and genital samples were collected semiannually in the HIM study. Briefly, a clinician utilized three saline-wetted swabs to collect exfoliated cells from the penis (coronal sulcus, glans, and ventral and dorsal areas of the shaft) and scrotum (4). Then, utilizing a separate swab, 360 degrees of the anal epithelium was swabbed between the anal os and the anal canal dentate line (22). HPV DNA testing was performed utilizing a polymerase chain reaction consensus primer system (PGMY 09/11), as described elsewhere (23).

Baseline serum samples from 234 MSM were previously tested with HPV16 L1 VLP-ELISA, a standard measure of immunogenicity measuring polyclonal antibodies, by a laboratory at Johns Hopkins University (JHU; ref. 15). VLP-ELISAs have been widely used in natural immunity studies in females, and are considered a more sensitive measure of cumulative HPV16 exposure than other neutralizing assays such as the Competitive Luminex Immunoassay (cLIA; ref. 24). However, studies utilizing VLP-ELISA often have different procedures, such as virus-like particle production, that could theoretically impact results. This study utilized a validated VLP-ELISA from the National Cancer Institute (NCI)'s HPV Immunology lab, which uses the same laboratory procedures (such as virus-like particle production) as a natural HPV16 immunity study performed in females who demonstrated immunity (13). Baseline serum was tested for all 475 MSM, and subsequently serum from the MSM who were HPV16 seropositive at baseline was tested at the 1-, 2-, 3-, and 4-year follow-up visits.

We also evaluated whether HPV16 antibody titers differed between HIM MSM in the highest tertile of HPV16 seropositivity and a random sample of 42 women ages 18 to 25 from the control arm of the Costa Rica vaccine trial (CVT) who were in the highest tertile of HPV16 seropositivity (≥60 EU/mL). A previous study in the CVT suggested that having a high level of natural antibodies (>60 EU/mL) may protect against HPV infection (13). We reran the HIM and CVT samples together using the same assay and VLP batch, to avoid any laboratory related bias.

Statistical methods

Utilizing the JHU HPV16 VLP-ELISA results, we compared the HPV16 optical density of HPV16-seropositive MSM with HPV16-seropositive MSW using box plots, t test, and Mann–Whitney tests. The NCI and JHU HPV16 VLP-ELISA results were compared using percent agreement and kappa statistic. Characteristics between HPV16-seropositive and HPV16-seronegative individuals (by the NCI assay) were compared using χ2 and Fisher exact tests for categorical variables and t tests for continuous variables.

Cumulative incidence curves and Cox proportional hazards regression were utilized to compare incident genital HPV16 in HPV16-seropositive versus -seronegative individuals—adjusting for age, country, sexual orientation [MSM or men who have sex with men and women (MSMW)], circumcision status, and number of new (past six months) and lifetime sex partners (any type of sex and anal sex particularly). This analysis was repeated comparing incident anal HPV16 in HPV16-seropositive versus -seronegative individuals. Seropositivity was also described by antibody titer level—continuously and categorically. Antibody titers were categorized a priori into tertiles to match the approach of previous studies (13, 18).

We conducted stratified analyses by whether an individual had a new sexual partner during follow-up as previously described (25). Incident HPV infections detected among those without a new sexual partner may be more likely to represent reactivation of a latent HPV infection and may be less likely to be protected against through natural immunity (23). Considering that sexual behavior measures were self-reported and may be prone to misclassification, we examined the effect of adjusting for incident nononcogenic HPV infection in the models as separate marker of sexual behavior.

We conducted sensitivity analyses restricting to HPV16 infections that persisted for 6 or more months. For this analysis, we excluded the 19 individuals who had an incident HPV16 infection at their last visit or only contributed two visits of follow-up. We also described the time to clearance for HPV16 seropositivity utilizing Kaplan–Meier curves and by utilizing Cox proportional hazard regression. All statistical tests were two sided and were considered significant at the alpha = 0.05 level. Analyses were performed utilizing STATA v13.1 and SAS.

Characteristics of MSM in the HIM study

All 475 men who ever reported having sex with a man in the HIM study and who were not positive for anal and genital HPV16 DNA at baseline were included in this analysis. Among the 475 remaining MSM, the mean age was 33.1 years with over half of the participants being white, never smokers, and residing in Brazil (Table 1). Approximately a quarter of participants reported having two or more new sexual partners in the last six months and over 20 lifetime sexual partners. A majority (73%) of the MSM had some report of also having sex with a woman in their lifetime and will be termed MSMW henceforth.

Table 1.

Baseline characteristics of HIM MSM participants by HPV16 serologic status

OverallHPV16 seronegativeHPV16 seropositive
Baseline characteristics of HIM's MSM participantsn = 475an = 347 (73%)bn = 128 (27%)bP
Age, mean (SD) 33.1 (10.1) 32.7 (10.5) 34.2 (9.3) 0.16 
 Younger than 30 191 (40%) 148 (78%) 43 (23%) 0.08 
 30—39 157 (33%) 105 (67%) 52 (33%)  
 40 or older 127 (26%) 94 (74%) 33 (26%)  
Country    0.03 
 United States 78 (16%) 60 (77%) 18 (23%)  
 Brazil 297 (63%) 205 (69%) 92 (31%)  
 Mexico 100 (21%) 82 (82%) 18 (18%)  
Race/ethnicity    0.19 
 White 263 (55%) 186 (71%) 77 (29%)  
 Black 91 (19%) 66 (73%) 25 (27%)  
 Mixed 93 (20%) 76 (82%) 17 (18%)  
 Other 21 (4.4%) 13 (62%) 8 (38%)  
 Refused 7 (1.5%) 6 (86%) 1 (14%)  
Cigarette smoker    0.86 
 Current 129 (27%) 96 (74%) 33 (26%)  
 Former 85 (18%) 63 (74%) 22 (26%)  
 Never 261 (55%) 188 (72%) 73 (28%)  
Number of new sex partners (last 6 months)    0.01 
 0 193 (41%) 153 (79%) 40 (21%)  
 1 106 (22%) 68 (64%) 38 (36%)  
 2 or more 109 (23%) 76 (70%) 33 (30%)  
 Not reported 67 (14%) 50 (75%) 17 (25%)  
Lifetime number of sex partners    <0.001 
 0 to 3 82 (17%) 69 (84%) 13 (16%)  
 4 to 19 193 (41%) 148 (77%) 45 (23%)  
 20 or more 143 (30%) 88 (62%) 55 (38%)  
 Not reported 57 (12%) 42 (74%) 15 (26%)  
Sexual orientation    <0.001 
 MSM 130 (27%) 78 (60%) 52 (40%)  
 MSMW 345 (73%) 269 (78%) 76 (22%)  
HIV infectionc    0.75 
 No 455 (96%) 331 (73%) 124 (27%)  
 Yes (acquired during study) 4 (0.8%) 3 (75%) 1 (25%)  
 Not reported 16 (3.4%) 13 (81%) 3 (19%)  
Circumcision    0.04 
 No 342 (72%) 241 (70%) 101 (30%)  
 Yes 133 (28%) 106 (80%) 27 (20%)  
OverallHPV16 seronegativeHPV16 seropositive
Baseline characteristics of HIM's MSM participantsn = 475an = 347 (73%)bn = 128 (27%)bP
Age, mean (SD) 33.1 (10.1) 32.7 (10.5) 34.2 (9.3) 0.16 
 Younger than 30 191 (40%) 148 (78%) 43 (23%) 0.08 
 30—39 157 (33%) 105 (67%) 52 (33%)  
 40 or older 127 (26%) 94 (74%) 33 (26%)  
Country    0.03 
 United States 78 (16%) 60 (77%) 18 (23%)  
 Brazil 297 (63%) 205 (69%) 92 (31%)  
 Mexico 100 (21%) 82 (82%) 18 (18%)  
Race/ethnicity    0.19 
 White 263 (55%) 186 (71%) 77 (29%)  
 Black 91 (19%) 66 (73%) 25 (27%)  
 Mixed 93 (20%) 76 (82%) 17 (18%)  
 Other 21 (4.4%) 13 (62%) 8 (38%)  
 Refused 7 (1.5%) 6 (86%) 1 (14%)  
Cigarette smoker    0.86 
 Current 129 (27%) 96 (74%) 33 (26%)  
 Former 85 (18%) 63 (74%) 22 (26%)  
 Never 261 (55%) 188 (72%) 73 (28%)  
Number of new sex partners (last 6 months)    0.01 
 0 193 (41%) 153 (79%) 40 (21%)  
 1 106 (22%) 68 (64%) 38 (36%)  
 2 or more 109 (23%) 76 (70%) 33 (30%)  
 Not reported 67 (14%) 50 (75%) 17 (25%)  
Lifetime number of sex partners    <0.001 
 0 to 3 82 (17%) 69 (84%) 13 (16%)  
 4 to 19 193 (41%) 148 (77%) 45 (23%)  
 20 or more 143 (30%) 88 (62%) 55 (38%)  
 Not reported 57 (12%) 42 (74%) 15 (26%)  
Sexual orientation    <0.001 
 MSM 130 (27%) 78 (60%) 52 (40%)  
 MSMW 345 (73%) 269 (78%) 76 (22%)  
HIV infectionc    0.75 
 No 455 (96%) 331 (73%) 124 (27%)  
 Yes (acquired during study) 4 (0.8%) 3 (75%) 1 (25%)  
 Not reported 16 (3.4%) 13 (81%) 3 (19%)  
Circumcision    0.04 
 No 342 (72%) 241 (70%) 101 (30%)  
 Yes 133 (28%) 106 (80%) 27 (20%)  

aOverall column represented by count (n), column percentage (%) unless otherwise specified.

bColumns represented by count and row percentage (%) unless otherwise specified.

cThe HIM study originally restricted to individuals without HIV at baseline; thus, those with HIV likely acquired it during the study.

The baseline HPV16 seroprevalence in this population was 27% (n = 128/475). Among the 128 HPV16-seropositive MSM, the median and mean antibody titer were 29.0 EU/mL and 63.3 EU/mL (SD = 83.0), respectively. Baseline HPV16 seropositivity was considerably more common among those with more recent sex partners, more lifetime sex partners, and in MSM compared with MSMW (all P ≤ 0.01). Additionally, HPV16 seropositivity was modestly more common in men who were uncircumcised (30% vs. 20%; P = 0.04) and those from Brazil (31% vs. 20%; P = 0.03). HPV16 seropositivity was not associated with age, race, smoking status, or HIV status (Table 1).

Comparisons between serologic assays

Results from the VLP-ELISAs were compared across the Hopkins and NCI laboratories and between duplicates of the NCI lab. All samples were tested twice by the NCI laboratory with strong agreement (Kappa = 0.95, 95% confidence interval (CI) = 0.92–0.98; Supplementary Fig. S2). The intra-assay CV between duplicates from the NCI laboratory with detectable antibody was 6.2% and the interassay CV between batches was 4.1%. Results among the 234 MSM tested by the JHU and NCI laboratories demonstrated substantial but lower agreement (Kappa = 0.68, 95% CI = 0.60–0.75; Supplementary Fig. S2). The percent agreement between two HPV16 VLP assays was 87.6%, as 95 of the 120 (79%) individuals classified as HPV16 seropositive by the Hopkins lab were classified as HPV16 seropositive by the NCI labs and 95 of the 129 (74%) individuals classified as HPV16 by the NCI lab were classified as HPV16 seropositive by the Hopkins lab. All subsequently described serologic analyses were tested via the NCI lab.

Evaluation of HPV16 natural immunity against genital HPV16 infection

Among the 446 eligible MSM, the incidence of genital HPV16 infection was similar in HPV16-seropositive MSM and HPV16-seronegative MSM (7.18 vs. 5.92 per 100 person-years; HR = 1.24, 95% CI = 0.80–1.93; Table 2 and Fig. 1). After adjustment for age, country, sexual orientation (MSM vs. MSMW), lifetime number of sex partners, number of recent sex partners, and circumcision status, the risk of HPV16 infection remained similar between HPV16-seropositive and -seronegative MSM (aHR = 1.05, 95% CI = 0.66–1.68; Table 2).

Table 2.

Natural HPV16 immunity against genital HPV16 in 446 MSM in the HIM studya

Genital HPV16 incidence
MeasureanGenital HPV16 infectionsPerson-time (person-years)Incidence rate (infs/100 p-ys)Unadjusted HRAdjusted HRb
HPV 16 L1 seropositive 
 No 321 63 1064 5.92 REF REF 
 Yes 125 29 404 7.18 1.24 (0.80–1.93) 1.05 (0.66–1.68) 
  T1: (8.0–16.2 EU/mL) 40 10 143 6.99 1.20 (0.61–2.35) 1.12 (0.56–2.25) 
  T2: (16.3–54.5 EU/mL) 41 129 6.98 1.21 (0.60–2.44) 0.92 (0.45–1.90) 
  T3: (≥54.5 EU/mL) 44 10 132 7.58 1.31 (0.67–2.57) 1.12 (0.56–2.23) 
  Continuous     0.67 0.97 
Restricting to persistent (6+ mo) HPV16 infectionsc 
 HPV16 seronegative 308 19 1028 1.85 REF REF 
 HPV16 seropositive 119 389 1.54 0.82 (0.33–2.06) 0.73 (0.28–1.89) 
Among those with a new sex partner during follow-up 
 HPV16 seronegative 254 55 857 6.42 REF REF 
 HPV16 seropositive 108 27 355 7.61 1.18 (0.74–1.87) 1.13 (0.69–1.85) 
Among those with no new sex partners during follow-up 
 HPV16 seronegative 67 207 3.86 REF REF 
 HPV16 seropositive 17 49 4.08 1.42 (0.29–7.05) 1.01 (0.09–11.43) 
Genital HPV16 incidence
MeasureanGenital HPV16 infectionsPerson-time (person-years)Incidence rate (infs/100 p-ys)Unadjusted HRAdjusted HRb
HPV 16 L1 seropositive 
 No 321 63 1064 5.92 REF REF 
 Yes 125 29 404 7.18 1.24 (0.80–1.93) 1.05 (0.66–1.68) 
  T1: (8.0–16.2 EU/mL) 40 10 143 6.99 1.20 (0.61–2.35) 1.12 (0.56–2.25) 
  T2: (16.3–54.5 EU/mL) 41 129 6.98 1.21 (0.60–2.44) 0.92 (0.45–1.90) 
  T3: (≥54.5 EU/mL) 44 10 132 7.58 1.31 (0.67–2.57) 1.12 (0.56–2.23) 
  Continuous     0.67 0.97 
Restricting to persistent (6+ mo) HPV16 infectionsc 
 HPV16 seronegative 308 19 1028 1.85 REF REF 
 HPV16 seropositive 119 389 1.54 0.82 (0.33–2.06) 0.73 (0.28–1.89) 
Among those with a new sex partner during follow-up 
 HPV16 seronegative 254 55 857 6.42 REF REF 
 HPV16 seropositive 108 27 355 7.61 1.18 (0.74–1.87) 1.13 (0.69–1.85) 
Among those with no new sex partners during follow-up 
 HPV16 seronegative 67 207 3.86 REF REF 
 HPV16 seropositive 17 49 4.08 1.42 (0.29–7.05) 1.01 (0.09–11.43) 

aThere were 29 MSM who were genital HPV16 DNA positive at baseline that were excluded from this analysis. These individuals were anal HPV16 DNA positive at baseline and thus were included on that analysis.

bAdjusted for age, country, sexual orientation (MSM vs. MSMW), lifetime number of sex partners, number of recent sex partners, and circumcision status.

cThis analysis excluded 19 individuals with an incident HPV16 infection at their last visit or only contributed two visits of follow-up.

Figure 1.

Cumulative incidence of genital HPV16 infection comparing HPV16-seronegative and HPV16-seropositive MSM in the HIM study. This figure presents the proportion of HPV16-seropositive and HPV16-seronegative MSM that acquired a genital HPV16 infection during the follow-up of the HIM study.

Figure 1.

Cumulative incidence of genital HPV16 infection comparing HPV16-seronegative and HPV16-seropositive MSM in the HIM study. This figure presents the proportion of HPV16-seropositive and HPV16-seronegative MSM that acquired a genital HPV16 infection during the follow-up of the HIM study.

Close modal

Risk of genital HPV16 infection was also similar among HPV16-seropositive men in the highest tertile of HPV16-seropositive with levels above 54.4 EU/mL (aHR = 1.12, 95% CI = 0.56–2.23) compared with HPV16-seronegative MSM; Table 2). Similarly, HPV16 seropositivity was not associated with risk of genital HPV16 infection when restricting to those with new sex partners, when further adjusting for incident nononcogenic HPV infection, or when excluding individuals who were anal HPV16 positive at baseline (each P > 0.20; Table 2). Results were similar by age groups, and there was no significant evidence of HPV16 natural immunity protection against subsequent infection among MSM in their 20s (aHR = 1.09, 95% CI = 0.55–2.20), 30s (aHR = 0.98, 95% CI = 0.43–2.27), or among MSM 40 years or older (aHR = 0.62, 95% CI = 0.17–2.19).

There were only 25 incident genital HPV16 infections that persisted ≥6 months, with 19 detected among HPV16-seronegative MSM (1.85 per 100 person-years) and 6 detected among HPV16-seropositive MSM (1.54 per 100 person-years). No differences in risk of subsequent incident persistent infection was present by sero-status (aHR = 0.73, 95% CI = 0.28–1.89) or by antibody titer (P = 0.36); however, none of the six persistent genital HPV16 infections detected among the seropositive MSM were among those with the highest HPV16 antibody levels (≥54.4 EU/mL).

Evaluation of HPV16 natural immunity against anal HPV16 infection

Among the 348 eligible MSM, the incidence of anal HPV16 was modestly lower than genital HPV16 (4.38 vs. 6.27 per 100 person-years; P < 0.01). There were 22 incident anal HPV16 infections detected. The incidence of anal HPV16 was nonsignificantly higher in MSM who were HPV16 seropositive at baseline (6.45 per 100 person-years) than HPV16 seronegative (6.45 vs. 3.46 per 100 person-years; HR = 1.81, 95% CI = 0.78–4.21; Table 3). Similarly, after adjustment the risk of anal HPV16 infection was nonsignificantly higher in HPV16-seropositive MSM (aHR = 2.34, 95% CI = 0.92–5.98; Table 3). Results remained nonsignificantly elevated when restricting to those with the highest tertile of seropositivity (>54.5 EU/mL), when also adjusting for incident nononcogenic HPV type during follow-up, and when excluding those who had a genital HPV16 infection at baseline (P > 0.20).

Table 3.

Natural HPV16 immunity against anal HPV16 in 348 MSM in the HIM study

Anal HPV16 incidence
MeasurenAnal HPV16 infectionsPerson-time (person-years)Incidence rate (infs/100 p-ys)Unadjusted HRAdjusted HRa
HPV 16 L1 seropositive 
 No 244 12 347 3.46 REF REF 
 Yes 104 10 155 6.45 1.81 (0.78–4.21) 2.34 (0.92–5.98) 
  T1: (8.0–16.2 EU/mL) 38 56 3.57 1.02 (0.23–4.55) 1.40 (0.81–10.23) 
  T2: (16.3–54.5 EU/mL) 30 42 9.52 2.76 (0.89–8.57) 2.89 (0.81–10.22) 
  T3: (≥54.5 EU/mL) 36 57 7.02 1.91 (0.61–5.97) 2.74 (0.81–9.35) 
  Continuous     0.86 0.52 
Anal HPV16 incidence
MeasurenAnal HPV16 infectionsPerson-time (person-years)Incidence rate (infs/100 p-ys)Unadjusted HRAdjusted HRa
HPV 16 L1 seropositive 
 No 244 12 347 3.46 REF REF 
 Yes 104 10 155 6.45 1.81 (0.78–4.21) 2.34 (0.92–5.98) 
  T1: (8.0–16.2 EU/mL) 38 56 3.57 1.02 (0.23–4.55) 1.40 (0.81–10.23) 
  T2: (16.3–54.5 EU/mL) 30 42 9.52 2.76 (0.89–8.57) 2.89 (0.81–10.22) 
  T3: (≥54.5 EU/mL) 36 57 7.02 1.91 (0.61–5.97) 2.74 (0.81–9.35) 
  Continuous     0.86 0.52 

aAdjusted for age, country, sexual orientation (MSM vs. MSMW), lifetime number of sex partners, number of recent sex partners, and circumcision status.

Longitudinal pattern of HPV16 seropositivity in MSM from the HIM and women from the Costa Rica Vaccine Trial

One-hundred and eighteen MSM who were HPV16 seropositive at baseline were tested for HPV16 seropositivity at HIM follow-up visits. At the 1-year follow-up visit, 94% (111/118) remained HPV16 seropositive, and at the 4-year follow-up visit, 90% remained HPV16 seropositive. Among the eight MSM who “HPV16 seroreverted” during the study and had a serology test at a visit after the one that defined seroreversion, there were five (62.5%) who tested seropositive at a subsequent annual visit.

We evaluated whether HPV16 antibody titers differed between MSM in highest tertile of HPV16 seropositivity (>54 EU/mL) compared with a random sample of 42 women from the control arm of the Costa Rica vaccine trial (CVT) who were in the highest tertile of HPV16 seropositivity (≥60 EU/mL; ref. 20). We retested the high HPV16 antibody CVT and HIM samples in the same run with the NCI's HPV16 VLP-ELISA, and found that the mean HPV16 antibody titer level between women from the CVT and MSM from the HIM to be similar (CVT: mean = 134.5 EU/mL, HIM: mean = 130.3 EU/mL, P = 0.84). The median HPV16 antibody titer level was 98.9 EU/mL in women from the CVT and 81.1 EU/mL in MSM from the HIM (P = 0.23).

This study, one of the largest evaluations of natural immunity in MSM, observed that HPV16 seropositivity is common among MSM and appears to persist for multiple years. However, there is no evidence that HPV16 natural antibodies protect against subsequent genital or anal HPV16 infection in MSM despite the fact that antibody levels observed were comparable with those seen in women for whom protection from reinfection have previously been documented (13). If confirmed, this may partially explain the high incidence of genital and anal HPV16 infection and related anal cancer seen in MSM throughout the lifespan (19).

The lack of natural HPV16 immunity seen in MSM in this study provides further support to previous findings that found natural immunity did not provide protection among MSW or MSM (15–18). Similarly, a meta-analysis found no evidence that natural immunity provided protection in men (MSW and MSM; pooled RR = 1.22, 95% CI = 0.67–1.77) against genital HPV16, despite evidence for moderate natural immunity in women (pooled RR = 0.65, 95% CI = 0.50–0.80, P heterogeneity = 0.05) against cervical HPV (12). We originally questioned whether the lack of natural immunity in MSW would apply to MSM for several reasons. First, MSM are more likely to acquire HPV natural immunity given that they are more likely to be exposed to HPV at an anatomic site with further access to the mucosal immune system (the anus), and are more likely to seroconvert after infection than MSW (10). Second, it appears, based on the preliminary results of this study, that antibodies often persist for years in HPV16-seropositive MSM and that the antibody levels are higher than those in HPV16-seropositive MSW. Third, previous studies of natural immunity in MSM had a limited number of MSM, often only included HIV-infected MSM, and utilized different methods of HPV16 antibody measurement compared with many of the female studies demonstrating natural immunity.

No protection against subsequent genital or anal HPV16 was observed in the current study despite the fact that the study included a larger number of HIV-uninfected MSM and utilized the same methods of antibody methods as a previous natural immunity study that observed protection in females (13). Although we did observe a fair bit of disagreement between the two labs with different methods of VLP production, suggesting that differences in the serologic assay methods could have an impact on the results. Results also remained similar in a number of sensitivity analyses including no evident protection among individuals with the highest level of antibody titers (≥55 EU/mL) at baseline, a level that has been suggested to be most protective in women (13).

There are a number of potential hypotheses as to why natural immunity against genital or anal HPV16 was not observed in this study. First, it is possible that MSM (and other men) do not obtain natural immunity against HPV infection, while a subset of women may obtain natural immunity. Our initial examination of those with higher antibody titers suggested a relatively similar antibody titer level among HIM MSM and women from the CVT (although there are several cohort differences between the two cohorts such age, sexual behavior, and geographic location that preclude a direct comparison by sex in this study). However, the mechanism of natural immunity is not well understood and antibody measures may only be a marker of other immune functions that may correlate with antibody titers in females but not males. Second, it is possible that the HPV DNA endpoints that were detected at follow-up visits may represent deposition rather true infections, particularly for genital HPV infection, which has a very low likelihood of having clinical significance in men. While results remained similar when restricting those with new sex partners and when adjusting for nononcogenic HPV infection (another measure of sexual activity), it is still possible that results are biased to the null in men due to the higher proportion of HPV deposition. Indeed, while there were a limited number of persistent (6 month+) genital and anal HPV16 infections identified, future studies should consider utilizing this endpoint given its lower likelihood to represent deposition and our nonsignificant point estimate (HR = 0.73), suggesting natural immunity may be evident in larger studies with persistent endpoints.

This study has several strengths. It was one of the largest evaluations of natural immunity in MSM, a group at high risk of HPV16-associated cancer. It also included robust serologic assays to measure natural HPV16 antibodies utilizing processes that have demonstrated HPV16 natural immunity in women (13). However, there are also limitations. While one of the larger studies in MSM, this study is still underpowered to observe endpoints such as persistent genital or anal HPV infection or HPV-associated disease. The primary endpoint used in this study, one time detection of genital and anal HPV16 is prone to deposition (a potential measurement bias). While there were fewer anal HPV16 infections detected during follow-up, the lower 95% CI of 0.92 excludes substantial protection against anal HPV16. In fact, other natural immunity studies in men have also suggested a potential higher risk of anal or oral HPV infection in HPV-seropositive men, potentially due to residual confounding by sexual risk behaviors (17, 26). Additionally, this study only examined natural immunity against HPV16, and it is possible that the level of antibody needed to confer protection may differ by HPV type (27).

Understanding the natural history and immunity of HPV in MSM is particularly important given the relatively low HPV vaccine uptake in men, MSM's lack of protection from herd immunity through female vaccination, and their high incidence of HPV16-associated anal cancer (19). Further research into vaccination and screening for anal cancer screening in middle-aged MSM is warranted.

A.G. Nyitray reports receiving a commercial research grant from Merck & Co. and has received renumeration from EUROGIN. No potential conflicts of interest were disclosed by the other authors.

Conception and design: D.C. Beachler, L.A. Pinto, A. Hildesheim, R. Viscidi, A.R. Kreimer, A.R. Giuliano

Development of methodology: D.C. Beachler, R. Viscidi, A.R. Kreimer

Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): L.A. Pinto, T.J. Kemp, A.G. Nyitray, R. Viscidi, A.R. Kreimer, A.R. Giuliano

Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): D.C. Beachler, L.A. Pinto, A. Hildesheim, J. Schussler, A.R. Kreimer, A.R. Giuliano

Writing, review, and/or revision of the manuscript: D.C. Beachler, L.A. Pinto, A.G. Nyitray, A. Hildesheim, R. Viscidi, J. Schussler, A.R. Kreimer, A.R. Giuliano

Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): D.C. Beachler, L.A. Pinto, T.J. Kemp, J. Schussler

Study supervision: A. Hildesheim, A.R. Kreimer, A.R. Giuliano

We would like to thank the participants of the HIM and CVT studies. Genotyping of the anal canal samples in the HIM Study—PI: A.G. Nyitray: NIAID-NIH (1R21AI101417-01), NCI-NIH (1R03 CA134204-01), Merck & Co., Inc. (IISP 33707). HIM study infrastructure was supported through an NCI grant to ARG (CA R01CA098803). D.C. Beachler was supported by NCI's cancer prevention fellowship. The Costa Rica HPV Vaccine Trial is a long-standing collaboration between investigators in Costa Rica and the NCI. The trial is sponsored and funded by the NCI (contract N01-CP-11005).

The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

1.
IARC Working Group on the Evaluation of Carcinogenic Risks to Humans
. 
Human papillomaviruses
.
IARC Monogr Eval Carcinog Risks Hum
2007
;
90
:
1
636
.
2.
Baseman
JG
,
Koutsky
LA
. 
The epidemiology of human papillomavirus infections
.
J Clin Virol
2005
;
32
Suppl 1
:
S16
24
.
3.
Ho
GYF
. 
Natural history of cervicovaginal papillomavirus infection in young women
.
N Engl J Med
1998
;
338
:
423
8
.
4.
Giuliano
AR
,
Lee
JH
,
Fulp
W
,
Villa
LL
,
Lazcano
E
,
Papenfuss
MR
, et al
Incidence and clearance of genital Human Papillomavirus infection in Men (HIM): a cohort study
.
Lancet
2011
;
377
:
932
40
.
5.
Giuliano
AR
,
Nyitray
AG
,
Kreimer
AR
,
Pierce Campbell
CM
,
Goodman
MT
,
Sudenga
SL
, et al
EUROGIN 2014 roadmap: differences in human papillomavirus infection natural history, transmission and human papillomavirus-related cancer incidence by gender and anatomic site of infection
.
Int J Cancer
2015
;
136
:
2752
60
.
6.
Rodriguez
AC
,
Schiffman
M
,
Herrero
R
,
Hildesheim
A
,
Bratti
C
,
Sherman
ME
, et al
Low risk of type-specific carcinogenic HPV re-appearance with subsequent cervical intraepithelial neoplasia grade 2/3
.
Int J Cancer
2012
;
131
:
1874
81
.
7.
Carter
JJ
,
Koutsky
LA
,
Hughes
JP
,
Lee
SK
,
Kuypers
J
,
Kiviat
N
, et al
Comparison of human papillomavirus types 16, 18, and 6 capsid antibody responses following incident infection
.
J Infect Dis
2000
;
181
:
1911
9
.
8.
Tong
Y
,
Ermel
A
,
Tu
W
,
Shew
M
,
Brown
D
. 
Association of HPV types 6, 11, 16, and 18 DNA detection and serological response in unvaccinated adolescent women
.
J Med Virol
2013
;
85
:
1786
93
.
9.
Stone
KM
,
Karem
KL
,
Sternberg
MR
,
McQuillan
GM
,
Poon
AD
,
Unger
ER
, et al
Seroprevalence of human papillomavirus type 16 infection in the United States
.
J Infect Dis
2002
;
186
:
1396
402
.
10.
Lu
B
,
Viscidi
RP
,
Wu
Y
,
Nyitray
AG
,
Villa
LL
,
Lazcano-Ponce
E
, et al
Seroprevalence of human papillomavirus (HPV) type 6 and 16 vary by anatomic site of HPV infection in men
.
Cancer Epidemiol Biomarkers Prev
2012
;
21
:
1542
6
.
11.
Giuliano
AR
,
Viscidi
RP
,
Torres
BN
,
Ingles
DJ
,
Sudenga
SL
,
Villa
LL
, et al
Seroconversion following anal and genital HPV infection in men: the HIM study
.
Papillomavirus Res
2015
;
1
:
109
15
.
12.
Beachler
DC
,
Jenkins
G
,
Safaeian
M
,
Kreimer
AR
,
Wentzensen
N
. 
Natural acquired immunity against subsequent genital human papillomavirus infection: a systematic review and meta-analysis
.
J Infect Dis
2016
;
213
:
1444
54
.
13.
Safaeian
M
,
Porras
C
,
Schiffman
M
,
Rodriguez
AC
,
Wacholder
S
,
Gonzalez
P
, et al
Epidemiological study of anti-HPV16/18 seropositivity and subsequent risk of HPV16 and -18 infections
.
J Natl Cancer Inst
2010
;
102
:
1653
62
.
14.
Klein
SL
,
Flanaga
KL
. 
Sex differences in immune responses
.
Nat Rev Immunol
2016
;
16
:
626
38
.
15.
Lu
B
,
Viscidi
RP
,
Wu
Y
,
Lee
JH
,
Nyitray
AG
,
Villa
LL
, et al
Prevalent serum antibody is not a marker of immune protection against acquisition of oncogenic HPV16 in men
.
Cancer Res
2012
;
72
:
676
85
.
16.
Lu
B
,
Hagensee
ME
,
Lee
JH
,
Wu
Y
,
Stockwell
HG
,
Nielson
CM
, et al
Epidemiologic factors associated with seropositivity to human papillomavirus type 16 and 18 virus-like particles and risk of subsequent infection in men
.
Cancer Epidemiol Biomarkers Prev
2010
;
19
:
OF1
6
.
17.
Mooij
SH
,
Landen
O
,
van der Klis
FR
,
van der Sande
MA
,
de Melker
HE
,
Coutinho
RA
, 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
:
375
86
.
18.
Beachler
DC
,
Viscidi
R
,
Sugar
EA
,
Minkoff
H
,
Strickler
HD
,
Cranston
RD
, 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
:
93
7
.
19.
Palefsky
JM
. 
HPV infection in men
.
Dis Markers
2007
;
23
:
261
72
.
20.
Giuliano
AR
,
Lazcano
E
,
Villa
LL
,
Flores
R
,
Salmeron
J
,
Lee
JH
, et al
Circumcision and sexual behavior: factors independently associated with human papillomavirus detection among men in the HIM study
.
Int J Cancer
2009
;
124
:
1251
7
.
21.
Pinto
LA
,
Kemp
TJ
,
Torres
BN
,
Isaacs-Soriano
K
,
Ingles
D
,
Abrahamsen
M
, et al
Quadrivalent human papillomavirus (HPV) vaccine induces HPV-specific antibodies in the oral cavity: results from the mid-adult male vaccine trial
.
J Infect Dis
2016
;
214
:
1276
83
.
22.
Nyitray
AG
,
Carvalho da Silva
RJ
,
Baggio
ML
,
Lu
B
,
Smith
D
,
Abrahamsen
M
, et al
Age-specific prevalence of and risk factors for anal human papillomavirus (HPV) among men who have sex with women and men who have sex with men: the HPV in men (HIM) study
.
J Infect Dis
2011
;
203
:
49
57
.
23.
Giuliano
AR
,
Lazcano-Ponce
E
,
Villa
LL
,
Flores
R
,
Salmeron
J
,
Lee
JH
, et al
The human papillomavirus infection in men study: human papillomavirus prevalence and type distrubtion among men residing in Brazil, Mexico, and the United States
.
Cancer Epidemiol Biomark Prev
2008
;
17
:
2036
43
.
24.
Robbins
HA
,
Li
Y
,
Porras
C
,
Pawlita
M
,
Ghosh
A
,
Rodriguez
AC
, et al
Glutathione S-transferase L1 multiplex serology as a measure of cumulative infection with human papillomavirus
.
BMC Infect Dis
2014
;
14
:
120
.
25.
Wilson
L
,
Pawlita
M
,
Castle
PE
,
Waterboer
T
,
Sahasrabuddhe
V
,
Gravitt
PE
, et al
Seroprevalence of eight oncogenic human papillomaviruses and acquired immunity against re-infection
.
J Infect Dis
2014
;
210
:
448
55
.
26.
Pierce Campbell
CM
,
Viscidi
RP
,
Torres
BN
,
Lin
HY
,
Fulp
W
,
Abrahamsen
M
, 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
:
45
8
.
27.
Pamnani
SJ
,
Sudenga
SL
,
Viscidi
R
,
Rollison
DE
,
Torres
BN
,
Ingles
DJ
, 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
:
6066
75
.