Background: Gay and bisexual men (GBM) are at disproportionately high risk of anal cancer. The precursor lesions, high-grade squamous intraepithelial lesions (HSIL), are very common and it is evident that not all HSIL progresses to cancer. The serologic response to anal human papillomavirus (HPV) in GBM has not been well characterized.

Methods: The Study of the Prevention of Anal Cancer is an ongoing cohort study of GBM ages 35 years and older. At six visits over three years, anal samples are collected for cytology, HPV DNA testing, and histology. Baseline serum was tested for HPV L1, E6, and E7 antibodies for 10 HPV types. Seroprevalence and associated predictors were analyzed.

Results: A total of 588 of 617 participants were included in this analysis. A total of 436 (74.2%) were seropositive for at least one of the 10 HPV types. Almost half had L1 antibodies to HPV6 (48.5%), over a third to HPV11 (36.4%) and HPV16 (34.5%). HIV-positive men were more likely to be HPV L1 seropositive. HSIL detection was highest among participants who were HPV serology and DNA positive. There was a borderline significant association between presence of HPV16 E6 antibodies and prevalent HSIL (OR = 2.97; 95% confidence interval, 0.92–9.60; P = 0.068).

Conclusions: HPV L1 seropositivity was common in this cohort of older GBM. These results suggest that HPV L1 seropositivity, in conjunction with anal HPV DNA detection, predicts concurrent HSIL. The apparent association between HPV16 E6 antibodies and prevalent HSIL is a finding with potential clinical significance that needs further exploration.

Impact: HPV seropositivity with concurrent DNA detection predicted anal HSIL detection. Cancer Epidemiol Biomarkers Prev; 27(7); 768–75. ©2018 AACR.

More than two decades of study of the natural history of human papillomavirus (HPV) in the cervix has shown that persistent high-risk HPV (HRHPV) infection is the precursor to precancerous high-grade squamous intraepithelial lesions (HSIL). If left untreated, some HSIL will eventually progress to invasive cancer, at a rate of around 1% per year (1). It is generally accepted that a broadly similar natural history applies to HPV infection at all sites in the lower anogenital tract (2). However, in gay and bisexual men (GBM), the prevalence of anal HPV infection is substantially higher than that of HPV infection of the cervix in women. In fact, even though GBM experience very high rates of anal cancer, the prevalence of anal HRHPV and associated HSIL is so high in GBM that it is clear that many, if not most infections must clear or never progress to cancer (3). Thus, it is evident that persistence, regression, and progression of anal HPV-related disease are not fully analogous to cervical HPV–related disease.

A better understanding of the immune response to anal HPV is likely to improve our understanding of the natural history of HPV at this site. Type-specific serum antibodies against the late protein (L1) of the viral capsid and the early HPV oncoproteins (E6 and E7) may be produced after infection. Infection-induced antibodies against HPV L1 are considered markers of cumulative exposure to the virus, representing past and present infection. Antibodies targeting E6 and E7 proteins are potential markers for HPV-associated cancer (4, 5). Seropositivity for HPV16 E6 antibodies has been shown to be a highly sensitive and specific marker HPV-driven oropharyngeal and anal cancer but is less common in cervical cancer (6–8). As yet, there have been no studies published reporting the association of HPV16 E6 and E7 seropositivity with anal HSIL.

We aimed to characterize the relationship between the serologic response to HRHPV L1, E6, and E7 and prevalent anal HPV infections and anal HSIL in a community-based cohort of Australian GBM.

The Study of the Prevention of Anal Cancer (SPANC) is a prospective cohort study exploring the epidemiology of anal HPV infection and related cytologic and histologic anal abnormalities in a community-recruited cohort of GBM (9). It is an ongoing study based in Sydney, Australia, involving HIV-negative and HIV-positive GBM ages 35 years and older. Participants were predominantly men recruited from gay community social events and organizations, referrals from other participants, or from previous community-based studies conducted by the Kirby Institute. Participants attend six visits over a period of three years, including a baseline enrollment visit, four follow-up clinical visits, and a final visit to discuss all study results. At each visit, an anal swab is collected for cytology screening and for HPV DNA genotyping using the Roche Linear array (Roche Molecular Systems), to determine the presence of 37 HPV genotypes. Serum samples are collected annually. A high resolution anoscopy (HRA) is performed and any potential HPV-related abnormalities are biopsied. The baseline data collection commenced in September 2010 and the study was fully recruited by August 2015. Follow-up of participants will continue to early 2018.

This article reports the serologic results from the baseline serum samples from SPANC participants who provided informed consent to storage and testing of their sera. Serum samples were analyzed at the German Cancer Research Center (DFKZ), Heidelberg, Germany, for L1, E6, and E7 antibodies to 10 mucosal HPV types: HPV6, 11, 16, 18, 31, 33, 35, 45, 52, and 58. Samples were analyzed using multiplex serology, an antibody detection method based on glutathione S-transferase (GST) capture ELISA, in combination with fluorescent bead technology (10). Briefly, full-length viral proteins are expressed in Escherichia coli bacteria in fusion with an N-terminal GST domain. Glutathione cross-linked to casein is coupled to fluorescence-labeled polystyrene beads (SeroMap, Luminex Corp.) and GST fusion proteins are affinity purified on the beads. Each fusion protein is bound to a spectrally distinct bead set, and fusion protein–loaded bead sets are mixed. Sera are incubated with the mixed bead sets at a final dilution of 1:100, and bound antibodies are detected with biotinylated goat anti-human IgG secondary antibody and streptavidin-R-phycoerythrin. A Luminex analyser (Luminex Corp) is used to identify the internal color of the individual beads and to quantify their reporter fluorescence [expressed as median fluorescence intensity (MFI) of at least 100 beads per set per serum]. Glutathione-casein–coupled bead sets are loaded with their respective antigens in one batch. Autofluorescence of each bead set and background reactions resulting from binding of secondary reagents to the antigen-loaded beads are determined in one well per plate without human serum. Mean bead background values are subtracted from raw MFI values and antigen-specific reactivity is determined by subtraction of the MFI of GST alone from the MFI of the specific antigen. MFI values were dichotomized as antibody positive or negative, using previously established seropositivity cutoffs that were calculated independently for each HPV antigen (10, 11). In this report, HRHPV includes HPV16, 18, 31, 33, 35, 45, 52, and 58.

Type-specific HPV seroprevalence was calculated and the exact binomial method was used to calculate 95% confidence intervals (CIs). ORs with corresponding 95% CIs were calculated for seroprevalent HPV type–specific associations with age group and HIV status. Among HIV-positive participants, the association between reported current and nadir CD4 T-cell count and seropositivity was analyzed. Univariate and multivariate logistic regressions were performed to determine the association between type-specific anal HPV DNA detection, anal low-grade squamous intraepithelial lesions (LSIL), HSIL, and HPV L1 capsid seropositivity at baseline. Composite “normal,” LSIL, and HSIL endpoints were used, comprising cytologic and/or histologic LSIL and HSIL diagnoses (12). Composite normal included normal cytology, normal histology, and/or biopsies not taken due to normal HRA appearances. Composite LSIL included histologic LSIL and cytologic LSIL and atypical squamous cells of undetermined significance (ASC-US). Composite HSIL included histologic HSIL and cytologic HSIL and atypical squamous cells—cannot exclude HSIL (ASC-H; ref. 12). HSIL detection was analyzed by pooled HRHPV serostatus (for the 8 included HRHPV types) overall and then stratified by anal HRHPV DNA detection. The association between HPV E6 and E7 seropositivity and anal HSIL was analyzed by logistic regression. Data analyses were performed using STATA version 14 (Stata Corporation).

Ethics approval for the study was granted by the St. Vincent's Hospital (Sydney, Australia), Australia Human Research Ethics Committee (file number 09/203). Written informed consent was obtained from all individuals before any study-specific procedures were performed. Separate, optional written informed consent was obtained for collection and storage of serum for future testing, including for HPV serology and data linkage studies.

A total of 617 men were enrolled in SPANC, with a median age of 49 years (range 35–79 years). Nearly all men (588, 95.3%) identified as gay or homosexual and 17 (2.9%) identified as bisexual. More than one-third (220, 35.7%) were HIV-positive. The great majority of these men (206, 93.6%) were currently receiving antiretroviral therapy, reported a viral load of less than 50 copies/mL (197, 89.5%) and had a CD4 T-cell count of more than 350 cells/μL (194, 88.0%). Of the 606 men with adequate anal samples, 525 (86.6%) had at least one anal HPV type detected by DNA genotyping. Among these, a single HPV type was detected in 101 participants (16.7%), while the rest had multiple HPV types detected. The prevalence of any anal HPV was 94.5% (207) in HIV-positive men and 82.2% (318) in HIV-negative men (P < 0.001). HPV16 was the most common HPV type (178, 29.4%), followed by HPV6 (109, 18.0%) and HPV45 (84, 13.9%; ref. 13).

Of the 617 men enrolled, 603 consented to storage and additional testing of collected specimens and had adequate serum samples for antibody analyses. Fifteen participants reported receiving prior HPV vaccination, four in the past year and 11 more than a year ago. Vaccinated participants had mean HPV16 L1 MFIs of 3,205 (vaccine in the past year) and 1,540 (more than a year ago), compared with a mean MFI of 569 for participants who reported never receiving HPV vaccination. These 15 participants were excluded from further analyses on L1 because their immune response was likely to be vaccine-induced rather than infection-induced.

Of the remaining 588, 436 (74.2%) were seropositive for at least one of the 10 HPV types analyzed. The seroprevalence for each HPV type is presented in Table 1. Almost half of the participants had L1 antibodies to HPV6 (48.5%) and over a third had L1 antibodies to HPV11 (36.4%) and HPV16 (34.5%). There was no significant association between increasing age and seropositivity. In fact, for every HPV type and for HPV overall, seroprevalence was numerically lower in men aged 65 and above than in younger men, but did not reach statistical significance for any HPV type. HIV-positive men were more likely to be L1 seropositive for each specific type, and this was significant for six of the 10 types (HPV16, P < 0.001; HPV18, P = 0.005; HPV31, P = 0.005; HPV35, P = 0.012; HPV45, P = 0.001 and HPV52, P = 0.017). HIV-positive men reporting a nadir CD4 T-cell count of less than 200 cells/μL were more likely to be HPV16 (P = 0.044) and HPV18 seropositive (P = 0.028). HIV-positive men with a higher current CD4 T-cell count were more likely to be seropositive for HPV45 (P = 0.042). This association was not seen with other HPV types.

Table 1.

Prevalence of LI HPV antibodies by HPV type (n = 588a), by age, HIV status, and CD4 count

PrevalenceHPV6 n (%)HPV11 n (%)HPV16 n (%)HPV18 n (%)HPV31 n (%)HPV33 n (%)HPV35 n (%)HPV45 n (%)HPV52 n (%)HPV58 n (%)Sero-positive to at least one type n (%)c
 285 (48.5) 214 (36.4) 203 (34.5) 159 (27.0) 97 (16.5) 51 (8.7) 125 (21.3) 162 (27.6) 67 (11.4) 68 (11.6) 152 (25.9) 
Age (yrs) 
 35–44 195 92 (48.7) 62 (32.8) 59 (31.2) 48 (25.4) 30 (15.9) 13 (6.9) 33 (17.5) 43 (22.8) 21 (11.1) 19 (10.1) 140 (75.1) 
 45–54 229 115 (52.3) 89 (40.5) 73 (33.2) 63 (28.6) 39 (17.7) 19 (8.6) 57 (25.9) 67 (30.5) 29 (13.2) 24 (10.9) 166 (75.5) 
 55–64 122 54 (44.3) 47 (38.5) 51 (41.8) 37 (30.3) 21 (17.2) 15 (12.3) 26 (21.3) 41 (33.6) 13 (10.7) 18 (14.8) 89 (73.0) 
 ≥65 57 24 (42.1) 16 (28.1) 20 (35.1) 11 (19.3) 7 (12.3) 4 (7.0) 9 (15.8) 11 (19.3) 4 (7.0) 7 (12.3) 41 (71.9) 
 Pb 0.269 0.992 0.152 0.868 0.743 0.364 0.915 0.495 0.467 0.307 0.758 
HIV 
 Neg 387 172 (45.7) 128 (34.0) 110 (29.3) 87 (23.1) 50 (13.3) 31 (8.2) 68 (18.1) 86 (22.9) 34 (9.1) 43 (11.4) 273 (72.6) 
 Pos 216 113 (53.3) 86 (40.6) 93 (43.9) 72 (34.0) 47 (22.2) 20 (9.4) 57 (26.9) 76 (35.9) 33 (15.6) 25 (11.8) 163 (76.9) 
 P 0.078 0.114 <0.001 0.005 0.005 0.623 0.012 0.001 0.017 0.897 0.255d 
PrevalenceHPV6 n (%)HPV11 n (%)HPV16 n (%)HPV18 n (%)HPV31 n (%)HPV33 n (%)HPV35 n (%)HPV45 n (%)HPV52 n (%)HPV58 n (%)Sero-positive to at least one type n (%)c
 285 (48.5) 214 (36.4) 203 (34.5) 159 (27.0) 97 (16.5) 51 (8.7) 125 (21.3) 162 (27.6) 67 (11.4) 68 (11.6) 152 (25.9) 
Age (yrs) 
 35–44 195 92 (48.7) 62 (32.8) 59 (31.2) 48 (25.4) 30 (15.9) 13 (6.9) 33 (17.5) 43 (22.8) 21 (11.1) 19 (10.1) 140 (75.1) 
 45–54 229 115 (52.3) 89 (40.5) 73 (33.2) 63 (28.6) 39 (17.7) 19 (8.6) 57 (25.9) 67 (30.5) 29 (13.2) 24 (10.9) 166 (75.5) 
 55–64 122 54 (44.3) 47 (38.5) 51 (41.8) 37 (30.3) 21 (17.2) 15 (12.3) 26 (21.3) 41 (33.6) 13 (10.7) 18 (14.8) 89 (73.0) 
 ≥65 57 24 (42.1) 16 (28.1) 20 (35.1) 11 (19.3) 7 (12.3) 4 (7.0) 9 (15.8) 11 (19.3) 4 (7.0) 7 (12.3) 41 (71.9) 
 Pb 0.269 0.992 0.152 0.868 0.743 0.364 0.915 0.495 0.467 0.307 0.758 
HIV 
 Neg 387 172 (45.7) 128 (34.0) 110 (29.3) 87 (23.1) 50 (13.3) 31 (8.2) 68 (18.1) 86 (22.9) 34 (9.1) 43 (11.4) 273 (72.6) 
 Pos 216 113 (53.3) 86 (40.6) 93 (43.9) 72 (34.0) 47 (22.2) 20 (9.4) 57 (26.9) 76 (35.9) 33 (15.6) 25 (11.8) 163 (76.9) 
 P 0.078 0.114 <0.001 0.005 0.005 0.623 0.012 0.001 0.017 0.897 0.255d 

aFifteen participants who reported receiving HPV vaccination were excluded from the analysis.

bPtrend.

cThe OR for age for participants who were HPV seronegative was 0.96 (95% CI, 0.76–1.22).

dThe OR for HIV status for participants who were HPV seronegative was 1.26 (95% CI, 0.85–1.86).

The relationship between HPV L1 serology and DNA detection is shown in Table 2. Approximately 50% of participants were both L1 seronegative and DNA negative for each of HPV6 and HPV 16. For all other HPV types, more than 60% were L1 seronegative and anal HPV DNA negative (Table 2). In general, for each HPV type, more men who were L1 seropositive were also DNA positive. The exception was HPV45, for which 13% of the seropositive and 14% of the seronegative were DNA positive. For HPV16, anal HPV16 DNA was detected in 86 (22.8%) of the seronegative and 85 (42.7%) of the seropositive.

Table 2.

Prevalence of anal HPV DNA detection and HPV L1 serology in 577a participants with both serology and HPV DNA results available

Type-specific HPVL1 serostatus and type-specific anal HPV DNA status
L1 Seropositive DNA PositiveL1 Seropositive DNA NegativeL1 Seronegative DNA PositiveL1 Seronegative DNA Negative
HPV6 83 (14.4) 199 (34.5) 21 (3.6) 274 (47.5) 
HPV11 39 (6.8) 170 (29.5) 16 (2.8) 352 (61.0) 
HPV16 85 (14.7) 114 (19.8) 86 (14.9) 292 (50.6) 
HPV18 28 (4.9) 128 (22.2) 42 (7.3) 379 (65.7) 
HPV31 12 (2.1) 84 (14.6) 24 (4.2) 457 (79.2) 
HPV33 7 (1.2) 43 (7.5) 23 (4.0) 504 (87.4) 
HPV35 8 (1.4) 115 (19.9) 18 (3.1) 436 (75.6) 
HPV45 21 (3.6) 139 (24.1) 60 (10.4) 357 (61.9) 
HPV52 15 (2.6) 51 (8.8) 64 (11.1) 447 (77.5) 
HPV58 9 (1.6) 57 (9.9) 57 (9.9) 454 (78.7) 
Type-specific HPVL1 serostatus and type-specific anal HPV DNA status
L1 Seropositive DNA PositiveL1 Seropositive DNA NegativeL1 Seronegative DNA PositiveL1 Seronegative DNA Negative
HPV6 83 (14.4) 199 (34.5) 21 (3.6) 274 (47.5) 
HPV11 39 (6.8) 170 (29.5) 16 (2.8) 352 (61.0) 
HPV16 85 (14.7) 114 (19.8) 86 (14.9) 292 (50.6) 
HPV18 28 (4.9) 128 (22.2) 42 (7.3) 379 (65.7) 
HPV31 12 (2.1) 84 (14.6) 24 (4.2) 457 (79.2) 
HPV33 7 (1.2) 43 (7.5) 23 (4.0) 504 (87.4) 
HPV35 8 (1.4) 115 (19.9) 18 (3.1) 436 (75.6) 
HPV45 21 (3.6) 139 (24.1) 60 (10.4) 357 (61.9) 
HPV52 15 (2.6) 51 (8.8) 64 (11.1) 447 (77.5) 
HPV58 9 (1.6) 57 (9.9) 57 (9.9) 454 (78.7) 

aFifteen participants who reported receiving HPV vaccination were excluded from the analysis.

The relationships between HPV DNA and serologic positivity with anal HPV-associated disease are presented in Table 3. The four categories examined were: double negatives, serology positive/DNA negative, serology negative/DNA positive, and double positives. For most HPV types, HSIL detection increased across these four categories. Compared with men who were HPV16 seronegative/DNA negative, participants who were HPV16 seronegative/DNA positive were 4 fold more likely (OR = 4.26; 95% CI, 2.23–8.12) to have HSIL detected, whereas participants who were HPV16 seropositive/DNA positive were almost 13-fold more likely (OR = 12.93; 95% CI, 4.99–33.48) to have HSIL detected (Table 3).

Table 3.

Composite cytologic–histologic diagnoses by combinations of HPV L1 serology and anal HPV DNA status

Composite normala (%)Composite LSILb (%)Composite HSILc (%)OR (95% CI)dPe
HPV6     <0.001 
 DNA - SERO - 94 (34.4) 64 (24.2) 106 (40.2)  
 DNA - SERO + 57 (29.1) 42 (21.4) 97 (49.5) 1.51 (0.98–2.23)  
 DNA + SERO - 3 (14.3) 7 (33.3) 11 (52.4) 3.25 (0.88–12.01)  
 DNA + SERO + 6 (7.5) 29 (36.3) 45 (56.3) 6.65 (2.72–16.29)  
HPV11 
 DNA - SERO - 117 (34.3) 81(23.8) 142 (41.8) <0.001 
 DNA - SERO + 42 (25.3) 43(25.9) 81 (48.8) 1.59 (1.02–2.48)  
 DNA + SERO - 0 (0) 10(62.5) 6 (37.5) —  
 DNA + SERO + 1 (2.6) 8 (20.5) 30 (76.9) 24.72 (3.32–183.99)  
HPV16 
 DNA - SERO - 100 (35.6) 82 (29.2) 99 (35.2) <0.001 
 DNA - SERO + 41(36.9) 33 (29.7) 37 (33.3) 0.91 (0.54–1.54)  
 DNA + SERO - 14(16.7) 11(13.1) 59 (70.2) 4.26 (2.23–8.12)  
 DNA + SERO + 5 (5.9) 16 (18.8) 64 (75.3) 12.93 (4.99–33.48)  
HPV18 
 DNA - SERO - 116 (31.8) 94 (25.8) 155(42.5) <0.001 
 DNA - SERO + 40 (31.8) 34 (27.0) 52 (41.3) 0.97 (0.60–1.57)  
 DNA + SERO - 2 (4.8) 12 (28.6) 28 (66.7) 10.48 (2.45–44.87)  
 DNA + SERO + 2 (7.1) 2 (7.1) 24 (85.7) 8.98 (2.08–38.76)  
HPV31 
 DNA - SERO - 139 (31.3) 111 (25.0) 194 (43.7) <0.001 
 DNA - SERO + 19 (23.5) 26 (32.1) 36 (44.4) 1.36 (0.75–2.47)  
 DNA + SERO - 2 (8.3) 5 (20.8) 17 (70.8) 6.09 (1.38–26.79)  
 DNA + SERO + 0 (0) 0 (0) 12 (100) —  
HPV33 
 DNA - SERO - 149 (30.5) 125 (25.6) 215 (44.0) 0.004 
 DNA - SERO + 7 (16.7) 16 (38.1) 19 (45.2) 1.88 (0.77–4.59)  
 DNA + SERO - 3 (13.0) 1 (4.4) 19 (82.6) 4.39 (1.28–15.10)  
 DNA + SERO + 1 (14.3) 0 (0) 6 (85.7) 4.16 (0.50–34.90)  
HPV35 
 DNA - SERO - 131 (31.0) 103 (24.4) 189 (44.7) 0.027 
 DNA - SERO + 25 (22.3) 35 (31.3) 52 (46.4) 1.44 (0.85–2.44)  
 DNA + SERO - 2 (11.1) 3 (16.7) 13 (72.2) 4.51 (0.10–20.30)  
 DNA + SERO + 2 (25.0) 1 (12.5) 5 (62.5) 1.73 (0.33–9.07)  
HPV45 
 DNA - SERO - 108 (31.2) 93 (26.8) 146 (42.1) 0.001 
 DNA - SERO + 40 (29.4) 33 (24.3) 63 (46.3) 1.17 (0.73–1.86)  
 DNA + SERO - 11 (19.3) 9 (15.8) 37 (64.9) 2.49 (1.21–5.10)  
 DNA + SERO + 1 (4.8) 7 (33.3) 13 (61.9) 9.62 (1.24–74.63)  
HPV52 
 DNA - SERO - 136 (31.3) 110 (25.4) 188 (43.2) 0.004 
 DNA - SERO + 10 (20.4) 16 (32.7) 23 (46.9) 1.66 (0.77–3.61)  
 DNA + SERO - 11 (17.5) 14 (22.2) 38 (60.3) 2.50 (1.23–5.06)  
 DNA + SERO + 3 (20.0) 2 (13.3) 10 (66.7) 2.41 (0.65–8.93)  
HPV58 
 DNA - SERO - 134 (30.5) 109 (24.8) 196 (44.7) 0.004 
 DNA - SERO + 18 (32.1) 20 (35.7) 18 (32.1) 0.68 (0.34–1.36)  
 DNA + SERO - 7 (12.3) 10 (17.5) 40 (70.2) 3.91 (1.70–8.98)  
 DNA + SERO + 1 (11.1) 3 (33.3) 5 (55.6) 3.42 (0.39–29.59)  
Composite normala (%)Composite LSILb (%)Composite HSILc (%)OR (95% CI)dPe
HPV6     <0.001 
 DNA - SERO - 94 (34.4) 64 (24.2) 106 (40.2)  
 DNA - SERO + 57 (29.1) 42 (21.4) 97 (49.5) 1.51 (0.98–2.23)  
 DNA + SERO - 3 (14.3) 7 (33.3) 11 (52.4) 3.25 (0.88–12.01)  
 DNA + SERO + 6 (7.5) 29 (36.3) 45 (56.3) 6.65 (2.72–16.29)  
HPV11 
 DNA - SERO - 117 (34.3) 81(23.8) 142 (41.8) <0.001 
 DNA - SERO + 42 (25.3) 43(25.9) 81 (48.8) 1.59 (1.02–2.48)  
 DNA + SERO - 0 (0) 10(62.5) 6 (37.5) —  
 DNA + SERO + 1 (2.6) 8 (20.5) 30 (76.9) 24.72 (3.32–183.99)  
HPV16 
 DNA - SERO - 100 (35.6) 82 (29.2) 99 (35.2) <0.001 
 DNA - SERO + 41(36.9) 33 (29.7) 37 (33.3) 0.91 (0.54–1.54)  
 DNA + SERO - 14(16.7) 11(13.1) 59 (70.2) 4.26 (2.23–8.12)  
 DNA + SERO + 5 (5.9) 16 (18.8) 64 (75.3) 12.93 (4.99–33.48)  
HPV18 
 DNA - SERO - 116 (31.8) 94 (25.8) 155(42.5) <0.001 
 DNA - SERO + 40 (31.8) 34 (27.0) 52 (41.3) 0.97 (0.60–1.57)  
 DNA + SERO - 2 (4.8) 12 (28.6) 28 (66.7) 10.48 (2.45–44.87)  
 DNA + SERO + 2 (7.1) 2 (7.1) 24 (85.7) 8.98 (2.08–38.76)  
HPV31 
 DNA - SERO - 139 (31.3) 111 (25.0) 194 (43.7) <0.001 
 DNA - SERO + 19 (23.5) 26 (32.1) 36 (44.4) 1.36 (0.75–2.47)  
 DNA + SERO - 2 (8.3) 5 (20.8) 17 (70.8) 6.09 (1.38–26.79)  
 DNA + SERO + 0 (0) 0 (0) 12 (100) —  
HPV33 
 DNA - SERO - 149 (30.5) 125 (25.6) 215 (44.0) 0.004 
 DNA - SERO + 7 (16.7) 16 (38.1) 19 (45.2) 1.88 (0.77–4.59)  
 DNA + SERO - 3 (13.0) 1 (4.4) 19 (82.6) 4.39 (1.28–15.10)  
 DNA + SERO + 1 (14.3) 0 (0) 6 (85.7) 4.16 (0.50–34.90)  
HPV35 
 DNA - SERO - 131 (31.0) 103 (24.4) 189 (44.7) 0.027 
 DNA - SERO + 25 (22.3) 35 (31.3) 52 (46.4) 1.44 (0.85–2.44)  
 DNA + SERO - 2 (11.1) 3 (16.7) 13 (72.2) 4.51 (0.10–20.30)  
 DNA + SERO + 2 (25.0) 1 (12.5) 5 (62.5) 1.73 (0.33–9.07)  
HPV45 
 DNA - SERO - 108 (31.2) 93 (26.8) 146 (42.1) 0.001 
 DNA - SERO + 40 (29.4) 33 (24.3) 63 (46.3) 1.17 (0.73–1.86)  
 DNA + SERO - 11 (19.3) 9 (15.8) 37 (64.9) 2.49 (1.21–5.10)  
 DNA + SERO + 1 (4.8) 7 (33.3) 13 (61.9) 9.62 (1.24–74.63)  
HPV52 
 DNA - SERO - 136 (31.3) 110 (25.4) 188 (43.2) 0.004 
 DNA - SERO + 10 (20.4) 16 (32.7) 23 (46.9) 1.66 (0.77–3.61)  
 DNA + SERO - 11 (17.5) 14 (22.2) 38 (60.3) 2.50 (1.23–5.06)  
 DNA + SERO + 3 (20.0) 2 (13.3) 10 (66.7) 2.41 (0.65–8.93)  
HPV58 
 DNA - SERO - 134 (30.5) 109 (24.8) 196 (44.7) 0.004 
 DNA - SERO + 18 (32.1) 20 (35.7) 18 (32.1) 0.68 (0.34–1.36)  
 DNA + SERO - 7 (12.3) 10 (17.5) 40 (70.2) 3.91 (1.70–8.98)  
 DNA + SERO + 1 (11.1) 3 (33.3) 5 (55.6) 3.42 (0.39–29.59)  

aComposite normal included normal cytology, normal histology, and/or biopsies not taken due to normal macroscopic appearance.

bComposite LSIL included histologic LSIL and cytologic LSIL and ASC-US.

cComposite HSIL included histologic HSIL and cytologic HSIL and ASC-H.

dOR comparing composite HSIL with composite normal diagnosis.

ePheterogeneity

Men who were L1 seropositive for at least one HRHPV type were significantly more likely to have HSIL detected (P = 0.032) than men who were seronegative to all. Among men with detectable HRHPV DNA, there was a nonsignificant trend toward more HSIL being detected among seropositive than seronegative participants (P = 0.088). There was no difference in HSIL detection by serostatus in participants with no detectable HRHPV DNA (P = 0.349; Table 4).

Table 4.

Anal HRHPV DNA, L1 serology status and association with composite cytologic–histologic HSIL diagnosis

No cHSIL (%)cHSIL (%)P
Any HRHPV DNA positive 
 HRHPV sero - 77 (44.0) 98 (56.0) 0.088 
 HRHPV sero + 70 (35.4) 128 (64.7)  
HR HPV DNA negative 
 HRHPV sero - 90 (52.3) 82 (47.7) 0.349 
 HRHPV sero + 15 (62.5) 9 (37.5)  
No cHSIL (%)cHSIL (%)P
Any HRHPV DNA positive 
 HRHPV sero - 77 (44.0) 98 (56.0) 0.088 
 HRHPV sero + 70 (35.4) 128 (64.7)  
HR HPV DNA negative 
 HRHPV sero - 90 (52.3) 82 (47.7) 0.349 
 HRHPV sero + 15 (62.5) 9 (37.5)  

The association of E6 and E7 antibodies with HSIL is presented in Table 5. The association of HPV16E6 seropositivity with composite HSIL approached statistical significance (OR = 2.97; 95% CI, 0.92–9.60; P = 0.068). Of the 14 men who had composite HSIL and tested HPV16 E6 positive, three participants were also positive for E7 (14). One participant was HPV33 E6 and E7 double seropositive. He had possible LSIL on cytology and normal histology. No participant was HPV18, 31, 35, 45, 52, or 58 E6 and E7 double seropositive.

Table 5.

Anal HPV E6 and E7 serology and composite cytologic–histologic HSIL diagnosis

HPV E6 and E7 serostatusNo composite HSIL (%) (n = 305)Composite HSILa (%) (n = 263)OR (95% CI)P
HPV6 E6 
 No 305 (53.7) 263 (46.3) N/A N/A 
 Yes   
HPV6 E7 
 No 302 (53.5) 263 (46.6) 0.87 (0.74–1.03) 0.101 
 Yes 3 (100) 0 (0)   
HPV11 E6 
 No 299 (53.7) 258 (46.3) 0.97 (0.29–3.20) 0.955 
 Yes 6 (54.6) 5 (45.5)   
HPV11 E7 
 No 304 (53.7) 262 (46.3) 1.16 (0.07–18.64) 0.916 
 Yes 1 (50.0) 1 (50.0)   
HPV16 E6 
 No 301 (54.3) 253 (45.7) 2.97 (0.92–9.60) 0.068 
 Yes 4 (28.6) 10 (71.4)   
HPV16 E7 
 No 290 (53.1) 256 (46.9) 0.53 (0.21–1.32) 0.171 
 Yes 15 (68.2) 7 (31.8)   
HPV18 E6 
 No 304 (53.6) 267 (46.4) 3.51 (0.36–33.93) 0.2785 
 Yes 1 (25.0) 3 (75.0)   
HPV18 E7 
 No 304 (53.8) 261 (46.2) 2.33 (0.21–25.84) 0.491 
 Yes 1 (33.3) 2 (66.7)   
HPV31 E6 
 No 301 (53.9) 261 (46.2) 1.46 (0.39–5.49) 0.577 
 Yes 4 (44.4) 5 (55.6)   
HPV31 E7 
 No 303 (53.5) 263 (46.5) N/A N/A 
 Yes 2 (100) 0 (0)   
HPV33 E6 
 No 300 (53.6) 260 (46.4) 0.69 (0.16–2.92) 0.617 
 Yes 5 (62.5) 3 (37.5)   
HPV33 E7 
 No 285 (54.1) 242 (45.9) 1.24 (0.65–2.34) 0.513 
 Yes 20 (48.8) 21 (51.2)   
HPV35 E6 
 No 303 (53.5) 263 (46.5) N/A N/A 
 Yes 2 (100) 0 (0)   
HPV35 E7 
 No 300 (53.9) 259 (46.1) 0.93 (0.25–3.49) 0.910 
 Yes 5 (55.6) 4 (44.4)   
HPV45 E6 
 No 296 (53.6) 253 (46.4) 1.30 (0.52–3.25) 0.575 
 Yes 9 (47.4) 10 (52.6)   
HPV45 E7 
 No 301 (53.6) 261 (46.4) 0.58 (0.10–3.17) 0.527 
 Yes 4 (66.7) 2 (33.3)   
HPV52 E6 
 No 301 (53.6) 261 (46.4) 0.58 (0.10–3.17) 0.527 
 Yes 4 (66.7) 2 (33.3)   
HPV52 E7 
 No 296 (53.5) 257 (46.5) 0.77 (0.27–2.19) 0.621 
 Yes 9 (60.0) 6 (40.0)   
HPV58 E6 
 No 299 (53.6) 259 (46.4) 0.77 (0.21–2.76) 0.688 
 Yes 6 (60.0) 4 (40.0)   
HPV58 E7 
 No 299 (53.8) 257 (46.2) 1.16 (0.37–3.65) 0.795 
 Yes 6 (50.0) 6 (50.0)   
HPV E6 and E7 serostatusNo composite HSIL (%) (n = 305)Composite HSILa (%) (n = 263)OR (95% CI)P
HPV6 E6 
 No 305 (53.7) 263 (46.3) N/A N/A 
 Yes   
HPV6 E7 
 No 302 (53.5) 263 (46.6) 0.87 (0.74–1.03) 0.101 
 Yes 3 (100) 0 (0)   
HPV11 E6 
 No 299 (53.7) 258 (46.3) 0.97 (0.29–3.20) 0.955 
 Yes 6 (54.6) 5 (45.5)   
HPV11 E7 
 No 304 (53.7) 262 (46.3) 1.16 (0.07–18.64) 0.916 
 Yes 1 (50.0) 1 (50.0)   
HPV16 E6 
 No 301 (54.3) 253 (45.7) 2.97 (0.92–9.60) 0.068 
 Yes 4 (28.6) 10 (71.4)   
HPV16 E7 
 No 290 (53.1) 256 (46.9) 0.53 (0.21–1.32) 0.171 
 Yes 15 (68.2) 7 (31.8)   
HPV18 E6 
 No 304 (53.6) 267 (46.4) 3.51 (0.36–33.93) 0.2785 
 Yes 1 (25.0) 3 (75.0)   
HPV18 E7 
 No 304 (53.8) 261 (46.2) 2.33 (0.21–25.84) 0.491 
 Yes 1 (33.3) 2 (66.7)   
HPV31 E6 
 No 301 (53.9) 261 (46.2) 1.46 (0.39–5.49) 0.577 
 Yes 4 (44.4) 5 (55.6)   
HPV31 E7 
 No 303 (53.5) 263 (46.5) N/A N/A 
 Yes 2 (100) 0 (0)   
HPV33 E6 
 No 300 (53.6) 260 (46.4) 0.69 (0.16–2.92) 0.617 
 Yes 5 (62.5) 3 (37.5)   
HPV33 E7 
 No 285 (54.1) 242 (45.9) 1.24 (0.65–2.34) 0.513 
 Yes 20 (48.8) 21 (51.2)   
HPV35 E6 
 No 303 (53.5) 263 (46.5) N/A N/A 
 Yes 2 (100) 0 (0)   
HPV35 E7 
 No 300 (53.9) 259 (46.1) 0.93 (0.25–3.49) 0.910 
 Yes 5 (55.6) 4 (44.4)   
HPV45 E6 
 No 296 (53.6) 253 (46.4) 1.30 (0.52–3.25) 0.575 
 Yes 9 (47.4) 10 (52.6)   
HPV45 E7 
 No 301 (53.6) 261 (46.4) 0.58 (0.10–3.17) 0.527 
 Yes 4 (66.7) 2 (33.3)   
HPV52 E6 
 No 301 (53.6) 261 (46.4) 0.58 (0.10–3.17) 0.527 
 Yes 4 (66.7) 2 (33.3)   
HPV52 E7 
 No 296 (53.5) 257 (46.5) 0.77 (0.27–2.19) 0.621 
 Yes 9 (60.0) 6 (40.0)   
HPV58 E6 
 No 299 (53.6) 259 (46.4) 0.77 (0.21–2.76) 0.688 
 Yes 6 (60.0) 4 (40.0)   
HPV58 E7 
 No 299 (53.8) 257 (46.2) 1.16 (0.37–3.65) 0.795 
 Yes 6 (50.0) 6 (50.0)   

aComposite HSIL included histologic HSIL and cytologic HSIL and ASC-H.

In this cohort of GBM, three quarters of men had L1 antibodies to at least one HPV type, and HPV6, HPV11, and HPV16 L1 antibodies were most commonly detected. HIV-positive men were more likely to be HPV L1 seropositive, and a nadir CD4 T-cell count of less than 200 cells/μL was predictive of HPV16 and HPV18 L1 seropositivity. Type-specific anal DNA was more prevalent in men who were L1 seropositive for the same HPV type. For most HRHPV types, HSIL detection was highest among participants who were both HPV serology and DNA positive. The presence or absence of HPV L1 antibodies was not significantly associated with HSIL among participants who were anal HPV DNA negative. There was a borderline significant increased risk of HSIL in those who were HPV16 E6 seropositive.

The high rates of HPV L1 seroprevalence we have found are comparable with findings from other serologic studies involving GBM. HPV L1 seroprevalence ranged from 31% to 86% of participants in these studies (15–21). In a Dutch population-based cross-sectional survey, men who reported sexual attraction to men had an HPV16 L1 seroprevalence of 31.4% compared with 16.8% among heterosexual men (22). Similarly, in U.S. cross-sectional national surveys, HPV16 L1 seroprevalence was 37.7% among men who reported ever having had sex with a man (23), and was 30% among 169 Austrian HIV-positive men, 59 of whom were homosexual (24).

Inconsistent relationships of HPV L1 antibody status with age have been described. Increasing HPV seroprevalence up to age 35 years was described in the multinational HIM study (18, 25). In the Amsterdam-based H2M study, seroprevalence increased across the three age groups of less than 35, 35–44, and 45 and above. The HYPER study, which enrolled 200 GBM aged between 16 and 20 years in Melbourne, Australia found very low seroprevalence (3.0% for HPV18 to 12.5% for HPV6) for the quadrivalent vaccine HPV types (26). In a San Francisco study of HIV-positive MSM, no association with age was found. Other studies group men into age group inclusive of 45–70 years or older than 45 years (15, 25) or include separate groups of less than 10 men older than 60 (17) and thus cannot identify whether a decline occurs at ages older than 60 years. In SPANC, HPV seroprevalence was nonsignificantly lower in men ages 65 years and older.

HPV L1 seroprevalence in most studies of GBM is more common among HIV-positive men. In a GBM population recruited from a sexual health clinic in the Amsterdam-based H2M study, HIV-positive men were more likely than HIV-negative men to be HPV16 seropositive (65% vs. 31%; OR = 4.05; 95% CI, 2.84–5.77; ref. 27). This is a larger difference than seen in SPANC (44% vs. 30%) despite both studies reporting higher numbers of lifetime and recent sexual partners among HIV-positive participants (9, 27).

HPV16 L1 seroprevalence was associated with concurrent detection of HPV16 DNA in the anal canal. Although L1 seropositivity is believed to reflect lifetime history of HPV infection, in our study, type-specific seropositivity was associated with current infection and with current disease (HSIL). This suggests that HPV seropositivity may predict persistent infection, but this requires confirmation in longitudinal data. These results correlate with findings on HPV infection in the vagina in women. The HELIUS study from the Netherland reported that in adjusted analyses, vaginal HPV DNA detection was significantly associated with type-specific HPV seropositivity (28).

The SPANC study is the first study to report on the relationship between HPV seropositivity and HPV-associated noninvasive anal disease. We found that HPV16 L1 seroprevalence was a strong predictor of concurrent HSIL. Compared with men who were HPV16 seronegative/DNA negative, participants who were HPV16 seronegative/DNA positive were 4-fold more likely to have HSIL detected, whereas participants who were HPV16 seropositive/DNA positive were 13-fold more likely to have HSIL detected. Thus, HPV16 seropositivity added incrementally to the likelihood of HSIL detection in participants who were anal HPV16 DNA positive. In addition, men who were seropositive for antibodies against at least one HRHPV type were significantly more likely to have HSIL detected (P = 0.018) than HRHPV seronegative men.

Previous studies have clearly demonstrated that HPV16 E6 seropositivity is associated with anal cancer. A nested case–control study within the European Prospective Investigation into Cancer and Nutrition (EPIC) found that 7 of 24 (29%) people who subsequently developed anal cancer and 0.6% of controls (4 of 718) were HPV16 E6 seropositive, translating to a 75-fold increase in cancer risk (7). In a case–control study nested in the Swiss HIV Cohort antibodies against HPV16 L1 (OR = 4.5) and E6 (OR = ∞; 95% CI, 4.6–∞) were strongly associated with anal cancer. In that study, antibodies against HPV16 E6 were 100% specific, in that the prevalence was zero in people without anal cancer, but were a relatively insensitive marker for anal cancer (29). HPV16 E6 seropositivity has been shown to be rare in cancer-free populations (30). A pooled analysis of almost 5,000 controls from cancer studies found that only 32 people (0.7%) were HPV16 E6 seropositive (31). The prevalence of HPV16 E6 antibodies in SPANC participants, presumed to be cancer-free on enrollment, was 2.3%, slightly higher than in the other cancer-free populations described above. However, there are few data on HPV16 E6 seropositivity in populations of homosexual men (32). The association between HPV16 E6 seropositivity and HSIL that we found was of borderline significance, with a 3-fold increase in the likelihood of anal HSIL (OR, 2.97; 95% CI, 0.92–9.60). There have been no previous studies of the association between HPV E6 antibodies and anal HSIL. Although the seroprevalence of HPV16 E6 is low (4% in participants with HSIL), the association should be explored further, to determine whether HPV16 E6 serology has a role as a biomarker in an anal cancer screening program.

The SPANC study is a large prospective study and one of the few studies globally to perform anal cytology, HRA screening, and HPV genotyping at the same visit on all participants. Serologic results were available for the vast majority of participants. Recruitment for the study was mostly from community-based settings, with broad inclusion criteria and thus we believe that results are likely to be generalizable to the GBM population in Australia.

There were several limitations in this research. Differences in serologic methods and variations between cut-off definitions for positivity between laboratories limit direct comparisons with other studies. The results presented are cross sectional and thus we cannot determine whether HPV anal infections or anal HSIL are acute or persistent or whether the serologic response is related to past or current infection with a specific HPV type. This will also be examined in subsequent analyses of SPANC longitudinal data. In addition, multiple HPV types were often present and it was not possible to conclusively attribute a HPV type based on a cytology specimen to an HSIL lesion.

In summary, HPV L1 seropositivity was common in this cohort of GBM, particularly among HIV-positive men. Concurrent detection of type-specific anal DNA was twice as common in men who were L1 seropositive for the same type and HSIL detection was highest among participants who were both HPV serology and DNA positive. These results suggest that HPV L1 seropositivity, in conjunction with anal HPV DNA, predicts concurrent HSIL and could be evaluated as part of an algorithm for current HSIL risk. The apparent association between HPV16 E6 antibodies and prevalent HSIL is a finding with potential clinical significance that needs further exploration.

A. Cornall has received speakers bureau honoraria from Seqirus (formerly bioCSL). S.M. Garland reports receiving commercial research grants from Merck, Glaxo Smith Kline, CSL, and the Commonwealth Department of Health; has received speakers bureau honoraria from MSD and SPMSD; and is a consultant/advisory board member for Merck. A.E. Grulich reports receiving a commercial research grant from Gilead as well as other commercial research support from Viiv and Seqirus; has received speakers bureau honoraria from Merck; and is a consultant/advisory board member for Viiv. No potential conflicts of interest were disclosed by the other authors.

The views expressed in this publication do not necessarily represent the position of the Australian Government.

Conception and design: F. Jin, R.J. Hillman, C.K. Fairley, A.E. Grulich

Development of methodology: T. Waterboer, F. Jin, R.J. Hillman, S. Tabrizi, A.E. Grulich

Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): I.M. Poynten, T. Waterboer, D.J. Templeton, R.J. Hillman, C. Law, A. Cornall, S. Tabrizi, J.M. Roberts

Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): I.M. Poynten, T. Waterboer, F. Jin, D.J. Templeton, S. Tabrizi, J.M. Roberts, S.M. Garland, A.E. Grulich

Writing, review, and/or revision of the manuscript: I.M. Poynten, T. Waterboer, F. Jin, D.J. Templeton, R.J. Hillman, A. Cornall, S. Tabrizi, J.M. Roberts, S.M. Garland, C.K. Fairley, A.E. Grulich

Study supervision: A.E. Grulich

This research is supported by a National Health and Medical Research Council Program New Investigator Grant (to I.M. Poynten). The SPANC study is funded by a National Health and Medical Research Council Program Grant (Sexually Transmitted Infections: Causes, Consequences and Interventions Grant #568971; to A.E. Grulich and C.K. Fairley) and a Cancer Council New South Wales Strategic Research Partnership Program Grant (Preventing Morbidity and Mortality from Anal Cancer Grant #13-11; I.M. Poynten, F. Jin, R.J. Hillman, D.J. Templeton, S. Tabrizi, S.M. Garland, AF, C.K. Fairley). Cytologic testing materials were provided by Hologic (Australia) Pty Ltd. The Kirby Institute is affiliated with the Faculty of Medicine, University of New South Wales and funded by the Australian Government of Health and Ageing.

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.
McCredie
M
,
Sharples
K
,
Paul
C
,
Baranyai
J
,
Medley
G
,
Jones
R
, et al
Natural history of cervical neoplasia and risk of invasive cancer in women with cervical intraepithelial neoplasia 3: a retrospective cohort study
.
Lancet Oncol
2008
;
9
:
425
34
.
2.
Darragh
T
,
Colgan
T
,
Cox
J
,
Heller
D
,
Henry
M
,
Luff
R
, et al
The lower anogenital squamous terminology standardization project for HPV-associated lesions: background and consensus recommendations from the College of American Pathologists and the American Society for Colposcopy and Cervical Pathology
.
J Lower Genit Tract Dis
2012
;
16
:
205
42
.
3.
Machalek
D
,
Poynten
I
,
Jin
F
,
Fairley
C
,
Farnsworth
A
,
Garland
S
, et al
Anal human papillomavirus infection and associated neoplastic lesions in men who have sex with men: a systematic review and meta-analysis
.
Lancet Oncol
2012
;
13
:
487
500
.
4.
Gravitt
P
. 
The known unknowns of HPV natural history
.
J Clin Invest
2011
;
121
:
4593
9
.
5.
Combes
J
,
Pawlita
M
,
Waterboer
T
,
Hammouda
D
,
Rajkumar
T
,
Vanhems
P
, et al
Antibodies against high-risk human papillomavirus proteins as markers for invasive cervical cancer
.
Int J Cancer
2014
;
135
:
2453
61
.
6.
Lang Kuhs
K
,
Pawlita
M
,
Gibson
S
,
Schmitt
N
,
Trivedi
S
,
Argiris
A
, et al
Characterization of human papillomavirus antibodies in individuals with head and neck cancer
.
Cancer Epidemiol
2016
:
46
52
.
7.
Kreimer
A
,
Brennan
P
,
Lang Kuhs
K
,
Waterboer
T
,
Clifford
G
,
Franceschi
S
, et al
Human papillomavirus antibodies and future risk of anogenital cancer: a nested case-control study in the European prospective investigation into cancer and nutrition study
.
J Clin Oncol
2015
;
33
:
877
84
.
8.
Reuschenbach
M
,
Waterboer
T
,
Wallin
K
,
Einenkel
J
,
Dillner
J
,
Hamsikova
E
, et al
Characterization of humoral immune responses against p16, p53, HPV16 E6 and HPV16 E7 in patients with HPV-associated cancers
.
Int J Cancer
2008
;
123
:
2626
31
.
9.
Machalek
D
,
Grulich
A
,
Hillman
R
,
Jin
F
,
Templeton
D
,
Tabrizi
S
, et al
The Study of the Prevention of Anal Cancer (SPANC): design and methods of a three-year prospective cohort study
.
BMC Public Health
2013
;
13
.
10.
Waterboer
T
,
Sehr
P
,
Michael
K
,
Franceschi
S
,
Nieland
J
,
Joos
T
, et al
Multiplex human papillomavirus serology based on in situ-purified glutathione s-transferase fusion proteins
.
Clin Chem
2005
;
51
:
1845
53
.
11.
Clifford
G
,
Shin
H
,
Oh
J
,
Waterboer
T
,
Ju
Y
,
Vaccarella
S
, et al
Serologic response to oncogenic human papillomavirus types in male and female university students in Busan, South Korea
.
Cancer Epidemiol Biomark Prev
2007
;
16
:
1874
9
.
12.
Machalek
D
,
Poynten
I
,
Jin
F
,
Hillman
R
,
Templeton
D
,
Law
C
, et al
A Composite cytology-histology endpoint allows a more accurate estimate of anal high-grade squamous intraepithelial lesion prevalence
.
Cancer Epidemiol Biomarkers Prev
2016
;
25
:
1134
43
.
13.
Poynten
I
,
Tabrizi
S
,
Jin
F
,
Templeton
D
,
Machalek
D
,
Cornall
A
, et al
Vaccine-preventable anal human papillomavirus in Australian gay and bisexual men
.
Papillomavirus Res
2017
;
3
:
80
84
.
14.
Holzinger
D
,
Wichmann
G
,
Baboci
L
,
Michel
A
,
Höfler
D
,
Wiesenfarth
M
, et al
Sensitivity and specificity of antibodies against HPV16 E6 and other early proteins for the detection of HPV16-driven oropharyngeal squamous cell carcinoma
.
Int J Cancer
2017
;
140
:
2748
57
.
15.
Mooij
S
,
van der Klis
F
,
van der Sande
M
,
Schepp
R
,
Speksnijder
A
,
Bogaards
J
, et al
Seroepidemiology of high-risk HPV in HIV-negative and HIV-infected MSM: the H2M study
.
Cancer Epidemiol Biomarkers Prev
2013
;
22
:
1698
708
.
16.
van Rijn
V
,
Mooij
S
,
Mollers
M
,
Snijders
P
,
Speksnijder
A
,
King
A
, et al
Anal, penile, and oral high-risk HPV infections and HPV seropositivity in HIV-positive and HIV-negative men who have sex with men
.
PLoS One
2014
;
9
:
e92208
.
17.
Sharma
R
,
Efird
J
,
Chein
A
,
Holly
E
,
Krajden
M
,
Berry
J
, et al
Prevalence and risk factors for neutralizing antibodies to human papillomavirus types 16 and 18 in HIV-positive men who have sex with men
.
J Acq Immune Def Syndr
2013
;
64
:
479
87
.
18.
Lu
B
,
Viscidi
R
,
Lee
J
,
Wu
y
,
Villa
L
,
Lazcano-Ponce
E
, et al
Human papillomavirus (HPV) 6, 11, 16, and 18 seroprevalence is associated with sexual practice and age: results from the multinational HPV Infection in Men Study (HIM Study)
.
Cancer Epidemiol Biomarkers Prev
2011
;
20
:
990
1002
.
19.
Lu
B
,
Viscidi
R
,
Wu
Y
,
Lee
J
,
Nyitray
A
,
Villa
L
, 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
.
20.
Poynten
I
,
Jin
F
,
Templeton
D
,
Prestage
G
,
Donovan
B
,
Pawlita
M
, et al
Prevalence, incidence, and risk factors for human papillomavirus 16 seropositivity in Australian homosexual men
.
Sex Transm Dis
2012
;
39
:
726
32
.
21.
Poynten
I
,
Waterboer
T
,
Jin
F
,
Templeton
D
,
Prestage
G
,
Donovan
B
, et al
Human papillomavirus types 6 and 11 seropositivity: risk factors and association with ano-genital warts among homosexual men
.
J Infect
2013
;
66
:
503
11
.
22.
Heiligenberg
M
,
Michael
K
,
Kramer
M
,
Pawlita
M
,
Prins
M
,
Coutinho
R
, et al
Seroprevalence and determinants of eight high-risk human papillomavirus types in homosexual men, heterosexual men, and women: a population-based study in Amsterdam
.
Sex Transm Dis
2010
;
37
:
672
80
.
23.
Stone
K
,
Karem
K
,
Sternberg
M
,
McQuillan
G
,
Poon
A
,
Unger
E
, et al
Seroprevalence of human papillomavirus type 16 infection in the United States
.
J Infect Dis
2002
;
186
:
1396
402
.
24.
Höpfl
R
,
Petter
A
,
Thaler
P
,
Sarcletti
M
,
Widschwendter
A
,
Zangerle
R
. 
High prevalence of high risk human papillomavirus-capsid antibodies in human immunodeficiency virus-seropositive men: a serological study
.
BMC Infect Dis
2003
;
3
.
25.
Rahman
S
,
Pierce Campbell
C
,
Rollison
D
,
Wang
W
,
Waterboer
T
,
Michel
A
, et al
Seroprevalence and associated factors of 9-valent human papillomavirus (HPV) Types among men in the multinational HIM study
.
PLoS One
2016
;
11
:
e0167173
.
26.
Zou
H
,
Tabrizi
S
,
Grulich
A
,
Hocking
J
,
Garland
S
,
Bradshaw
C
, et al
Antibody responses following incident anal and penile infection with human papillomavirus in teenage men who have sex with men
.
Int J Cancer
2016
;
139
:
639
46
.
27.
Alberts
C
,
Michel
A
,
Bruisten
S
,
Snijder
M
,
Prins
M
,
Waterboer
T
, et al
High-risk human papillomavirus seroprevalence in men and women of six different ethnicities in Amsterdam, the Netherlands: The HELIUS study
.
Papillomavirus Res
2017
.
28.
Kovaleva
A
,
CJ
A
,
Waterboer
T
,
Michel
A
,
Snijder
M
,
Vermeulen
W
, et al
A cross-sectional study on the concordance between vaginal HPV DNA detection and type-specific antibodies in a multi-ethnic cohort of women from Amsterdam, the Netherlands - the HELIUS study
.
BMC Infect Dis
2016
;
16
:
502
.
29.
Bertisch
B
,
Franceschi
S
,
Lise
M
,
Vernazza
P
,
Keiser
O
,
Schöni-Affolter
F
, et al
Risk factors for anal cancer in persons infectedWith HIV: a nested case-control study in the Swiss HIV Cohort Study
.
Am J Epidemiol
2013
;
178
:
877
84
.
30.
Beachler
D
,
Waterboer
T
,
Campbell
C
,
Ingles
D
,
Kuhs
K
,
Nyitray
A
, et al
HPV16 E6 seropositivity among cancer-free men with oral, anal or genital HPV16 infection
.
Papillomavirus Res
2016
;
2
:
141
4
.
31.
Lang Kuhs
K
,
Anantharaman
D
,
Waterboer
T
,
Johansson
M
,
Brennan
P
,
Michel
A
, et al
Human Papillomavirus 16 E6 Antibodies in Individuals without Diagnosed Cancer: A Pooled Analysis
.
Cancer Epidemiol Biomarkers Prev
2015
;
24
:
683
9
.
32.
Combes
J
,
Clifford
G
,
Egger
M
,
Cavassini
M
,
Hirsch
H
,
Hauser
C
, et al
Human papillomavirus antibody response following HAART initiation among MSM
.
AIDS
2017
;
31
:
561
9
.