Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Patterns of Human Papillomavirus DNA and Antibody Positivity in Young Males and Females, Suggesting a Site-Specific Natural Course of Infection

  • Henrike J. Vriend ,

    rianne.vriend@rivm.nl

    Affiliations Centre for Infectious Disease Control, National Institute for Public Health and the Environment (RIVM), Bilthoven, The Netherlands, Department of Internal Medicine, Division of Infectious Diseases, Tropical Medicine and AIDS, Center for Infection and Immunity Amsterdam (CINIMA), Academic Medical Center (AMC), Amsterdam, The Netherlands

  • Johannes A. Bogaards,

    Affiliations Centre for Infectious Disease Control, National Institute for Public Health and the Environment (RIVM), Bilthoven, The Netherlands, Department of Epidemiology and Biostatistics, VU University Medical Center, Amsterdam, The Netherlands

  • Fiona R. M. van der Klis,

    Affiliation Centre for Infectious Disease Control, National Institute for Public Health and the Environment (RIVM), Bilthoven, The Netherlands

  • Mirte Scherpenisse,

    Affiliation Centre for Infectious Disease Control, National Institute for Public Health and the Environment (RIVM), Bilthoven, The Netherlands

  • Hein J. Boot †,

    † Deceased.

    Affiliation Centre for Infectious Disease Control, National Institute for Public Health and the Environment (RIVM), Bilthoven, The Netherlands

  • Audrey J. King,

    Affiliation Centre for Infectious Disease Control, National Institute for Public Health and the Environment (RIVM), Bilthoven, The Netherlands

  • Marianne A. B. van der Sande,

    Affiliations Centre for Infectious Disease Control, National Institute for Public Health and the Environment (RIVM), Bilthoven, The Netherlands, Julius Center, University Medical Center, Utrecht, The Netherlands

  • on behalf of the Medical Microbiological Laboratories, Municipal Health Services

Abstract

Background

To monitor the impact of human papillomavirus types 16 and 18 vaccine on HPV infection dynamics in the Netherlands, we started an ongoing study in sexually transmitted infection (STI) clinics in 2009. Here, we analyze baseline type-specific HPV DNA and HPV-specific antibody positivity rates.

Methods

We enrolled 3569 men and women, 16–24 years of age, from 14 STI clinics, and estimated genital and anal HPV DNA and antibody positivity rates of 7 main carcinogenic HPV types. Generalized estimating equations regression analyses were applied to determine risk factors for, and associations between, type-specific HPV DNA and antibody positivity.

Results

Genital HPV DNA positivity rates were higher in women than in men; anal HPV DNA was especially high in men who have sex with men (MSM). HPV antibody seropositivity rates were also highest in women and MSM. High-risk sexual behavior was predictive of both HPV DNA and antibody positivity. Despite a strong correlation in serological profiles for multiple HPV types, seropositivity was independently associated with homologous HPV DNA detection.

Conclusions

HPV DNA and antibody positivity rates are higher in women and MSM than in heterosexual men, but their association is similar across gender. This suggests a site-specific natural course of infection.

Introduction

Human papillomavirus (HPV) is a common sexually transmitted virus known for its causal relation to cervical cancer. There are more than 100 HPV genotypes, with more than 15 carcinogenic types [1], [2]. In many countries, HPV vaccination has been introduced in sexually naïve girls to prevent infections with HPV-16/-18, which are most commonly found in cervical cancers. It is not yet known what the impact of HPV vaccination will be on HPV dynamics in partially vaccinated populations. Monitoring of type-specific HPV prevalence in both vaccinated and nonvaccinated people is, therefore, of great importance.

HPV infection does not always induce an immune response that results in HPV-specific antibodies (Ab) [3], [4], [5]. Even if women are diagnosed with precancerous cervical lesions that test positive for HPV DNA, they might still be negative for serum HPV Ab [6]. Whether HPV infection will lead to seroconversion depends on several factors, such as specific HPV types, persistence of infection, HPV DNA viral load, and site of infection [3], [4], [7], [8], [9], [10], [11]. In contrast to natural infection, HPV vaccination induces an immune response with high concentrations of HPV Ab, by far exceeding the HPV Ab concentrations found in nonvaccinated populations [12]. In addition, studies showed that vaccination against HPV-16/-18 can result in cross-protection against phylogenetically related genotypes [13], [14]. Therefore, it is possible that vaccination may not only result in a decline in HPV-16/-18 prevalence, but also in a decline in phylogenetically related genotypes such as HPV-31, -33, and -45. Conversely, unwanted effects like type replacement, i.e., the potential for nonvaccine HPV types to occupy the vacated ecologic niches, can occur as a result of the elimination of HPV-16/-18 [15]. This hypothesis has neither been confirmed nor rejected by epidemiological studies. As a result of reduced exposure to HPV-16/-18, an effect can also be expected among nonvaccinated men and women [16], [17], [18], [19].

The aim of our study was to describe HPV DNA and HPV-specific Ab detection rates of women, men who have sex with women only (MSW), and men who have sex with men (MSM), all of whom were without benefit of HPV vaccination. Furthermore, we explored associations between homologous and heterologous pairs of HPV DNA and HPV Ab types. This description will serve as a baseline measurement to which we can compare future monitoring rounds on HPV dynamics within the Netherlands.

Materials and Methods

Ethics Statement

The medical ethics committee of the University of Utrecht, the Netherlands, confirmed in writing that they waived the need for separate ethical approval and the need for written consent. This anonymous study used serum already collected for routine STI consultation, therefore no additional invasive procedures were needed. All eligible participants were informed about the purpose of the study prior to the regular STI consultation and full information was provided about the samples to be collected and the additional questionnaire to be administered. Only participants who consented verbally with all conditions were included in the study.

Study Population and Design

In 2009, the bivalent HPV vaccine was introduced in the Netherlands among 12- to 16-year-old girls. To monitor the effects of HPV-16/-18 vaccination on type-specific HPV dynamics in a young highly sexually active population, the PASSYON (PApillomavirus Surveillance among STI clinic YOungsters in the Netherlands) study was set up [20]. This biennial cross-sectional study includes 16- to 24-year-old male and female attendees of the sexually transmitted infection (STI) clinic. In 2009 and 2011, the first two rounds of this study took place in 14 STI clinics throughout the Netherlands; 10 STI clinics participated in both rounds.

A genital self-sample (vaginal or penile) of all participants was requested (sampling technique is described elsewhere [20]). In addition, MSM were requested for an anal sample as well. Procedures were similar in both rounds, however, in 2011, females and MSW also were invited to submit an anal sample. In addition, participation included: the use of serum already collected during routine consultation or collection of an extra blood sample if routine serum was not available; filling out a questionnaire; and permission to link the anonymous standardized data collected during the intake of the routine consultation with STI test results (e.g., Chlamydia trachomatis, Neisseria gonorrhoea, syphilis, HIV). The questionnaire covered sexual risk factors such as condom use, number of sex partners, and history of STI. Demographic data were also collected. Ethnicity was not recorded, but we asked the participants which population they belonged to, regardless of country of birth, and used their responses as proxies.

HPV DNA detection

All swabs were suspended in 1 ml universal transport medium buffer and stored at −20°C until processing. After thawing, swabs were vortexed and 200 µL of the sample was spiked with phocine herpes virus-1. DNA was subsequently extracted using the MagnaPure platform (Total Nucleic Acid Isolation Kit, Roche) and eluted in 100 µL elution buffer. HPV-DNA was amplified using the SPF10 primer set according to the manufacturer's instructions (DDL Diagnostic Laboratory, the Netherlands [21]). HPV-specific amplicons were detected using the DNA enzyme-linked immunoassay (HPV-DEIA, DDL Diagnostic Laboratory, the Netherlands). Amplicons of HPV-positive samples were subsequently analyzed in Line probe assay (HPV-LiPA, DDL Diagnostic Laboratory, the Netherlands). In this report, only the seven carcinogenic HPV types, which can also be detected in the Ab detection assay, are reported (i.e., HPV-16, -18, -31, -33, -45, -52, -58). Only persons positive for one of the seven HPV-types are reported as being HPV positive.

Serum antibody detection

Serum samples were stored at −80°C until analysis. HPV-specific serum antibodies against L1 virus-like particles (VLPs) for serotypes -16, -18, -31, -33, -45, -52 and -58 were assessed using a VLP-based multiplex immunoassay in which the VLPs were coupled to a set of seven distinct fluorescent microspheres. VLPs were kindly donated by GSK. Most samples were analysed with a 1/100 or 1/200 dilution, some in 1/2000 dilution. Serum samples were assumed to be antibody seropositive at cut-offs determined previously with this assay: 9, 13, 27, 11, 19, 14, 31 Luminex Units/millilitre (LU/ml) for HPV-16, -18, -31, -33, -45, -52, -58, respectively. The serological measurements and control steps have been described in more detail elsewhere [22].

Statistical Analyses

Since no significant differences were observed in HPV DNA or Ab positivity rates between the first two rounds (2009 and 2011), we chose to combine these rounds in this study to gain more power. We included only HIV-negative persons having a urogenital DNA sample as well as a serum sample. Women who reported HPV vaccination were excluded in these baseline analyses. We did not exclude men who reported to be HPV vaccinated (n = 17), since Ab concentrations of HPV-16/-18 were not significantly different from those of persons who reported not to be HPV vaccinated, and since men are not (yet) vaccinated in the Netherlands for HPV, it was unlikely that they actually had been vaccinated. Type-specific positivity rates for the seven detectable HPV DNA and HPV Ab types were estimated, stratified by women, MSW, and MSM. Positivity rates were estimated by dividing the number of type-specific HPV DNA positives by the total number tested, with 95% confidence intervals (CI) based on Wilson's score method. The Student's t-test was used to compare positivity rates. When the number testing positive was below five the Fisher's exact test was used.

Generalized estimating equations (GEE) analyses were performed to assess independent risk factors associated with the presence of HPV DNA, and separate analyses were performed to asses independent risk factors associated with HPV Ab seropositivity. GEE models are marginal and take into account multiple measurements, making them highly suitable for estimating population-averaged risks associated with repeated measures [23]. GEE models estimate conditional log odds ratios (βx) between an independent risk factor (x) and a multivariate outcome variable, while assuming a particular correlation structure between the repeated measures, in this case multiple HPV types [24]. In assessing independent risk factors for the presence of HPV DNA and for HPV Ab seropositivity, we used an ‘exchangeable correlation’ structure with a common odds ratio alpha, which is appropriate when the association between pairs of HPV types can be assumed to be constant. Analyses were stratified by gender (i.e., women, MSW, MSM), since univariable analyses showed that all three outcome variables (genital HPV DNA, anal HPV DNA, and HPV Ab) differed significantly between women, MSW, and MSM. Multivariable analyses were performed by including the following variables into a GEE model: age, ethnicity, years sexually active, number of lifetime sex partners, number of sex partners in the last six months, condom use in the last six months with a steady partner, condom use in the last six months with a casual partner, a chlamydia or gonorrhea infection in the past, a current urogenital chlamydia infection, a current anogenital chlamydia or gonorrhea infection, passive anal sex in the past six months, and the seven different HPV genotypes. Subsequently, all nonsignificant (p≥0.05) variables were removed by backward selection.

GEE models were also used to estimate possible type-specific associations between the presence of HPV DNA and HPV Ab. These analyses were also stratified by gender, but could not be performed for MSM because of the small group size. For these analyses, we used a saturated model for the independent risk factors and started with a fully parameterized model for the pair-wise log odds ratios between all 14 genital HPV DNA and HPV Ab outcome variables. The starting log odds ratio model was composed of 14×13/2 = 91 linear parameter combinations, one for each unique binary pair of type-specific DNA and/or Ab of the seven HPV types. All nonsignificant parameters were subsequently deleted from the model using Bonferroni-corrected p values based on Wald tests, and this procedure was repeated until all remaining parameters were significant. Final estimates were obtained after orthonormal transformation of the log odds ratio estimation matrix. In modelling the association between pairs of responses, we made use of the alternating logistic regression algorithm of the SAS GENMOD procedure [25].

Results

Characteristics study population

In 2009 and 2011, 1637 and 1932 participants, respectively, were tested for a genital HPV infection. Of these 3569 participants, 3439 (96%) were also tested for HPV Ab. After exclusion of 6 HIV-positive persons and 105 self-reported HPV-vaccinated women, a study population of 3328 participants remained, of whom two-thirds were females (n = 2233). A subgroup of 118 (5%) women, 56 (6%) MSW, and 124 (72%) MSM also had a test result for anal HPV DNA infection.

Characteristics stratified for women, MSW, and MSM are given in Table 1. The median age of the 3328 participants was 21 years (IQR: 20–23 years). Of all men, 16% (n = 173) reported having sex with men. Sexual behavioral factors showed a median age of first intercourse of 16 years (IQR: 15–17 years) and a median number of lifetime partners of 7 (IQR: 4–12) (6, 9, and 12 lifetime partners for women, MSW, and MSM, respectively). The percentage of persons with a current or past anogenital chlamydia or gonorrhoea infection is given in Table 2. Among women and MSW, 14% had a current urogenital chlamydia and 2% had a current gonorrhea infection. For MSM this was 5% and 4%, respectively. A current anal chlamydia and/or gonorrhea infection was reported for 1% of women and 14% of MSM. No anal infections were seen in MSW.

thumbnail
Table 1. Demographics and characteristics of the study population: women, MSW, and MSM.

https://doi.org/10.1371/journal.pone.0060696.t001

thumbnail
Table 2. Past and current (ano)genital chlamydia and gonorrhea infections: women, MSW, and MSM.

https://doi.org/10.1371/journal.pone.0060696.t002

HPV DNA positivity and HPV-specific antibody seropositivity

Women had the highest genital HPV DNA positivity rate (45%) with the highest type-specific rate for HPV-16 (17% of all women). Among men, MSM had lower positivity rates for genital (penile) DNA than MSW, 16% versus 26%, respectively. Although low genital HPV DNA positivity rates were seen in MSM, anal HPV DNA positivity rates were high. Within the subgroups that also provided anal samples (124 MSM, 118 women, and 56 MSW), similar anal HPV DNA positivity rates were seen in MSM (33%) and women (32%), whereas anal HPV DNA positivity was rare in MSW (4%) (Figure 1). HPV seropositivity was highest in women (49%), followed by MSM (34%) and MSW (19%).

thumbnail
Figure 1. Positivity rate for genital and anal DNA detection in women, MSW and MSM.

Black gives the positivity rate for persons with a HPV DNA detection only at the genitals, dark grey gives the positivity rate for persons with an HPV DNA infection genital as well as anal, and light grey are the positivity rates when persons are only detected with HPV DNA anally. The figures includes women (N = 118), MSW (N = 56) and MSM (N = 124) of whom a genital and an anal sample were available. MSW =  men who have sex with women only; MSM = men who have sex with men.

https://doi.org/10.1371/journal.pone.0060696.g001

Comparing HPV DNA positivity and HPV Ab seropositivity rates by gender, the main differences were the lower positivity rates for HPV Ab -16, -18, -31, and -52 compared to homologous HPV DNA positivity rates in MSW (Figure 2a–c). Women and MSM, however, had generally higher positivity rates for HPV Ab in comparison to homologous HPV DNA positivity rates. When comparing anal and genital infections in MSM, the anal HPV DNA positivity rates were higher than the genital HPV DNA positivity rates, although not significant (p = 0.25) (Figure 2c–d). Subsequently, the percentage of simultaneous detection of HPV DNA and homologous HPV Ab increased. The difference between anal HPV DNA positivity and HPV Ab seropositivity in MSM was similar to the genital rates observed in women. When observing women and MSM with both a genital as well as an anal sample, differences were seen for HPV Ab seropositivity among those infected with HPV DNA at one anatomical site only. Although numbers were small, women tended to be more often HPV Ab seropositive when only a genital HPV DNA infection was present compared to MSM (45% versus 18%, respectively (p = 0.16)). Whereas, HPV Ab seropositivity rates when only an anal HPV DNA infection was present were quite similar between women and MSM (29% versus 21%, respectively (p = 0.71)).

thumbnail
Figure 2. Genital (and anal) HPV DNA and HPV-specific antibody seropositivity in women, MSW, and MSM.

Figure 2a–c gives the genital HPV DNA (light grey), HPV-specific antibody seropositivity (dark grey), as well as the percentage of persons detected with genital HPV DNA as well as HPV antibodies (black) with 95% confidence intervals in women (N = 2233), MSW (N = 922), and MSM (N = 173) per HPV genotype. Figure 2d gives the anal HPV DNA (light grey), HPV-specific antibody seropositivity (dark grey) as well as the percentage of persons detected with anal HPV DNA as well as HPV antibodies (black) with 95% confidence intervals in MSM (N = 124). Ab = antibody; MSM = men who have sex with men; MSW = men who have sex with women only.

https://doi.org/10.1371/journal.pone.0060696.g002

Simultaneous detection of type-specific HPV DNA and homologous HPV Ab genotype was most frequent for genotype 16 among all genders and sites. From Figure 2a–b, the odds ratio for HPV-16 Ab seropositivity together with HPV-16 DNA positivity is directly estimated at 2.9 for women and at 3.3 for MSW. However, given the presence of HPV-16 DNA, the probability of testing positive for HPV-16 Ab is much higher for women than for MSW, 41% versus 11%, respectively (p<0.01).

Determinants of HPV DNA positivity and HPV-specific antibody seropositivity

Multivariable analyses resulted in different risk factors for genital HPV DNA positivity, anal HPV DNA positivity, and HPV Ab seropositivity. Most noticeable was the significant association of ethnicity with HPV Ab seropositivity for women as well as MSM, whereas ethnicity was not significantly associated with HPV DNA positivity. Persons who reported they were non-Dutch were more likely to be HPV Ab seropositive. Furthermore, in women the site of the current chlamydia infection corresponded with the site of the associated HPV DNA positivity: a genital chlamydia infection was associated with genital HPV DNA infection and an anal chlamydia infection with an anal HPV DNA infection. In MSM, reporting passive anal sex in the past six months was associated with HPV Ab seropositivity. All significant risk factors resulting from the multivariable GEE analysis are shown in Table 3, 4 and 5.

thumbnail
Table 3. Adjusted odds ratio and 95% confidence interval for the detection of genital and anal HPV DNA and HPV antibodies using multilevel analyses (generalized estimating equations [GEE)]) in women.

https://doi.org/10.1371/journal.pone.0060696.t003

thumbnail
Table 4. Adjusted odds ratio and 95% confidence interval for the detection of genital and anal HPV DNA and HPV antibodies using multilevel analyses (generalized estimating equations [GEE)]) in men who have sex with women only (MSW).

https://doi.org/10.1371/journal.pone.0060696.t004

thumbnail
Table 5. Adjusted odds ratio and 95% confidence interval for the detection of genital and anal HPV DNA and HPV antibodies using multilevel analyses (generalized estimating equations [GEE)]) in men who have sex with men (MSM).

https://doi.org/10.1371/journal.pone.0060696.t005

Association between genital HPV DNA and HPV-specific antibody genotypes

The log odds ratio model for genital HPV DNA positivity and HPV Ab seropositivity could be reduced from a total of 91 parameters to as few as 7 independent parameters in women. Table 6 shows the resulting estimates of the common odds ratios for any pair of the following: (1) homologous HPV DNA and HPV Ab genotypes; (2) HPV DNA genotypes; (3) HPV Ab genotypes; and (4) heterologous DNA and Ab genotypes. In addition to these associations, we found significantly increased odds ratios for DNA and Ab of HPV-58; DNA of HPV-31 and HPV-33; and Ab detection for a subset of genotypes (i.e., HPV-33, -45, -52, -58). For example, seropositivity for HPV-33 Ab in combination with HPV-45 Ab is a frequently occurring test result, as well as seropositivity for HPV-33 Ab with HPV-52 Ab. For MSW, the number of independent parameters that could be estimated was smaller than for women, which is in line with the smaller sample size and lower positivity rate in MSW as compared to women. Associations were similar between the sexes (Table 6), except that a much stronger association was observed between pairs of HPV Ab genotypes in MSW compared with women.

thumbnail
Table 6. Association (odds ratio and 95% confidence interval) between genital HPV DNA and HPV-specific antibody genotypes, for women and MSW.

https://doi.org/10.1371/journal.pone.0060696.t006

The results of the log odds ratio model provided a justification for the assumption of an ‘exchangeable correlation’ structure in assessing independent risk factors for the presence of HPV DNA and for HPV Ab seropositivity. The associations between HPV DNA genotypes or between HPV Ab genotypes estimated by the log odds ratio model corresponded well with the alphas estimated in the multivariable GEE analyses (Table 3, 4 and 5). Note that the common odds ratios estimated by assuming an exchangeable correlation structure are slightly higher (with exception of the odds ratio for HPV Ab in MSW), because these incorporate the small elevations related to certain HPV types in the log odds ratio model. Having verified the assumptions underlying an exchangeable correlation structure, we then compared the associations between HPV DNA genotypes or HPV Ab genotypes in MSM and those in women and in MSW. We found that MSM have similar associations between anal HPV DNA genotypes as women (Table 3 and 5). Associations between HPV Ab genotypes of MSM overlap with those of women and MSW when taking the 95% confidence intervals into account.

Discussion

Positivity rates were high in this young and sexually active population, without the benefit of HPV vaccination. HPV DNA positivity and HPV Ab seropositivity were higher in women than in MSW, with MSM in between, but the association between detection of type-specific DNA and serum Ab was similar across gender. The PASSYON study sheds new light on differences in the natural history of HPV infection between men and women, indirectly providing evidence for a site-specific natural course of infection.

In women, HPV DNA positivity rates were lower for each HPV genotype than HPV Ab seropositivity rates, whereas in MSW, type-specific HPV DNA positivity rates were higher than HPV Ab seropositivity rates with the exceptions of HPV Ab-33, -45, and -58. These specific Ab types had high positivity rates in all genders, which was confirmed by the results of the GEE analyses. Significantly higher association between Ab-33 and -45 and between Ab-33 and -58 was found, meaning that when an individual is positive for one of these Ab types, he or she is also frequently found positive for the other type. The gender differences we observed in HPV DNA positivity and HPV Ab seropositivity may be related to the type of tissue being infected. The dry keratinized tissue of the penis is much harder for the virus to infect than the soft mucosal tissue of the vagina or the anus, therefore men are less likely than women to develop a humoral response [8], [9], [10], [11]. Results in MSM support the hypothesis that tissue could play an important role in HPV Ab seropositivity, since high HPV Ab seropositivity rates in combination with high anal HPV DNA positivity rates were found in MSM, whereas genital HPV DNA positivity was low. The ratios between HPV Ab and anal HPV DNA positivity rates in MSM were similar to the ratios of HPV Ab and vaginal/anal HPV DNA positivity rates in women. To examine whether gender-specific factors other than anatomical site could influence seroconversion, we also compared the HPV Ab seropositivity rates within women and MSM who provided both genital and anal samples. Women with only a genital HPV DNA infection were more often seropositive for the homologous HPV Ab serotype as compared to MSM with only a penile HPV DNA infection, whereas HPV Ab seropositivity rates were comparable for women and MSM with only an anal HPV DNA infection. These findings support the hypothesis that penile HPV infections will result less often in seroconversion than vaginal or anal HPV infections. However, our study was not adequately powered to detect differences in seropositivity rate between women and MSM who were infected at one anatomical site only, and we cannot rule out the possibility that additional factors (e.g. genetic, behavioural) might contribute to gender-specific differences in seroconversion. We can only conjecture that anatomical site is a determinant in the natural course of HPV infection. Obviously, a cross-sectional design is inherently limited in drawing conclusions about dynamic associations and it would be worthwhile to further investigate the relation between site-specific HPV infection and the humoral immune response in a longitudinal study.

Multivariable GEE analyses helped clarify risk factors for HPV positivity. One risk factor for a HPV DNA infection in women was a current chlamydia infection. The site-specific relation between a chlamydia infection and an HPV infection was noticeable: a current genital chlamydia infection was associated with genital HPV DNA positivity and a current anal chlamydia infection with anal HPV DNA positivity. The underlying mechanism could be less efficient clearance because of this concurrent infection, resulting in a persistent infection, or easier invasion into the body. There were also behavioral factors, such as reported number of lifetime sex partners, associated with HPV DNA positivity and HPV Ab seropositivity. These results are in line with previous studies [11], [26], [27], [28], [29], [30].

The added value of GEE analyses is the possibility to include all type-specific HPV DNA and HPV Ab outcomes in one model, thereby making optimal use of the available data. In addition, including both HPV DNA and HPV Ab test results in one model gave us the option to estimate associations between homologous and heterologous pairs of HPV DNA and HPV Ab genotypes. Despite the difference in HPV DNA positivity and HPV Ab seropositivity between women and MSW, no significant gender difference was seen in the association between a homologous or heterologous pair of HPV DNA and HPV Ab types. Equal odds ratios for the associations between HPV DNA positivity and HPV Ab seropositivity in women and MSW do not imply that the natural course of the infection would be similar between the sexes. As shown, women have a higher probability to be seropositive when HPV DNA positive then MSW. A higher degree of persistence could be an explanation for the higher seropositivity in women. Furthermore, a higher association between any two HPV Ab genotypes in MSW was observed as compared to women. This higher association in MSW could be due to the less frequent occurrence of seroconversion. It is conceivable that characteristics related to seroconversion are restricted to specific persons such that, on a population level, it will lead to a stronger clustering of seropositivity in MSW than in women, where seroconversion is more common.

We also found type-specific HPV DNA and HPV Ab associations and found some associations between homologous HPV DNA and HPV Ab types, and between two heterologous HPV DNA types.

These associations could have occurred because of the high level of multiple HPV infections. However, the HPV types tested in the VLP-MIA are phylogenetically related to each other, and belong to the alpha-7 species (HPV-18, -45) and alpha-9 species (HPV-16, -31, -33, -52, -58) [31]. Since a high percentage of homology is found between HPV types of the same species, cross-reactive HPV-specific antibodies could account for the type-specific associations.

To summarize, the higher HPV DNA and Ab positivity rates in women and MSM as compared to heterosexual men, and the similar HPV DNA and Ab associations across gender suggest a site-specific natural course of infection. These findings in combination with the determinants for HPV DNA positivity and HPV Ab seropositivity can help us interpret the results of future rounds with more accuracy and will enable us to detect changes in HPV infection dynamics, both of vaccine and nonvaccine types.

Acknowledgments

The authors thank the following for their valuable contributions to the design or execution of the study: Mariet Feltkamp, Titia Heijman, Elske van Logchem, Adam Meijer, Chris Meijer, Hester de Melker, Martijn van Rooijen, Rutger Schepp, Peter Snijders, Jan Sonsma, and Hans van Vliet. Furthermore, the STI clinics, including all nurses and physicians, within the municipal health centers and the hospitals are acknowledged for their permission to collect data from their patients and their effort. The authors acknowledge the medical microbiologic laboratories and the analysts for storage and testing of the samples.

Medical Microbiological Laboratories: D.S. Luijt (Laboratory for Infectious Diseases Groningen), J.W.A. Rossen (St. Elisabeth Hospital Tilburg), M. Schutten (Erasmus Medical Center Rotterdam), R. Schuurman (University Medical Center Utrecht), A.G.C.L. Speksnijder (Public Health Laboratory Amsterdam), P. Wolffs (University Hospital Maastricht);

Municipal Health Services: G. Aalfs (MHS Drenthe), H. van Buel-Bruins (MHS IJsselland), P.J. Cornelissen (MHS Rivierenland), T.E. Doorn (University Medical Center Utrecht), H.M. Götz (MHS Rotterdam-Rijnmond), F. de Groot (MHS Groningen), C.J.P.A. Hoebe (MHS Zuid Limburg), P. Hut-van Vliet (MHS Fryslân), C.J.G. Kampman (MHS Twente), H. van Kruchten (MHS Hart voor Brabant), A.P. van Leeuwen (MHS Amsterdam), M. Pelgrim (MHS Gelderland-Midden).

Author Contributions

Data collection: HV. Contributed to the serum antibody detection methods: FVDK MS. Contributed to the DNA detection methods: HB AK. Conceived and designed the experiments: MVDS HV HB FVDK. Analyzed the data: HV JB. Wrote the paper: HV JB FVDK MS HB AK MVDS.

References

  1. 1. Bouvard V, Baan R, Straif K, Grosse Y, Secretan B, et al. (2009) A review of human carcinogens–Part B: biological agents. Lancet Oncol 10: 321–322.
  2. 2. Munoz N, Castellsague X, de Gonzalez AB, Gissmann L (2006) Chapter 1: HPV in the etiology of human cancer. Vaccine 24 Suppl 3S3/1–10.
  3. 3. Carter JJ, Koutsky LA, Hughes JP, Lee SK, Kuypers J, et al. (2000) Comparison of Human Papillomavirus Types 16, 18, and 6 Capsid Antibody Responses Following Incident Infection. The Journal of Infectious Diseases 181: 1911–1919.
  4. 4. Ho GY, Studentsov YY, Bierman R, Burk RD (2004) Natural history of human papillomavirus type 16 virus-like particle antibodies in young women. Cancer Epidemiol Biomarkers Prev 13: 110–116.
  5. 5. Steele J, Collins S, Wen K, Ryan G, Constandinou-Williams C, et al. (2008) Measurement of the humoral immune response following an incident human papillomavirus type 16 or 18 infection in young women by a pseudovirion-based neutralizing antibody assay. Clin Vaccine Immunol 15: 1387–1390.
  6. 6. Vonka V, Hamsikova E, Kanka J, Ludvikova V, Sapp M, et al. (1999) Prospective study on cervical neoplasia IV. Presence of HPV antibodies. Int J Cancer 80: 365–368.
  7. 7. Carter JJ, Madeleine MM, Shera K, Schwartz SM, Cushing-Haugen KL, et al. (2001) Human Papillomavirus 16 and 18 L1 Serology Compared across Anogenital Cancer Sites. Cancer Res 61: 1934–1940.
  8. 8. Edelstein ZR, Carter JJ, Garg R, Winer RL, Feng Q, et al. (2011) Serum Antibody Response Following Genital α9 Human Papillomavirus Infection in Young Men. Journal of Infectious Diseases 204: 209–216.
  9. 9. Lu B, Viscidi RP, Lee JH, Wu Y, Villa LL, et al. (2011) 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 20: 990–1002.
  10. 10. Newall AT, Brotherton JML, Quinn HE, McIntyre PB, Backhouse J, et al. (2008) Population Seroprevalence of Human Papillomavirus Types 6, 11, 16, and 18 in Men, Women, and Children in Australia. Clinical Infectious Diseases 46: 1647–1655.
  11. 11. Stone KM, Karem KL, Sternberg MR, McQuillan GM, Poon AD, et al. (2002) Seroprevalence of Human Papillomavirus Type 16 Infection in the United States. Journal of Infectious Diseases 186: 1396–1402.
  12. 12. De Carvalho N, Teixeira J, Roteli-Martins CM, Naud P, De Borba P, et al. (2010) Sustained efficacy and immunogenicity of the HPV-16/18 AS04-adjuvanted vaccine up to 7.3 years in young adult women. Vaccine 28: 6247–6255.
  13. 13. Brown DR, Kjaer SK, Sigurdsson K, Iversen O-E, Hernandez-Avila M, et al. (2009) The Impact of Quadrivalent Human Papillomavirus (HPV; Types 6, 11, 16, and 18) L1 Virus-Like Particle Vaccine on Infection and Disease Due to Oncogenic Nonvaccine HPV Types in Generally HPV-Naive Women Aged 16–26 Years. The Journal of Infectious Diseases 199: 926–935.
  14. 14. Paavonen J, Naud P, Salmerón J, Wheeler CM, Chow SN, et al. (2009) Efficacy of human papillomavirus (HPV)-16/18 AS04-adjuvanted vaccine against cervical infection and precancer caused by oncogenic HPV types (PATRICIA): final analysis of a double-blind, randomised study in young women. The Lancet 374: 301–314.
  15. 15. Durham DP, Poolman EM, Ibuka Y, Townsend JP, Galvani AP (2012) Reevaluation of epidemiological data demonstrates that it is consistent with cross-immunity among human papillomavirus types. J Infect Dis 206: 1291–1298.
  16. 16. Bogaards JA, Coupe VM, Xiridou M, Meijer CJ, Wallinga J, et al. (2011) Long-term impact of human papillomavirus vaccination on infection rates, cervical abnormalities, and cancer incidence. Epidemiology 22: 505–515.
  17. 17. Brisson M, Van de Velde N, Boily MC (2011) Different population-level vaccination effectiveness for HPV types 16, 18, 6 and 11. Sex Transm Infect 87: 41–43.
  18. 18. Donovan B, Grulich AE (2011) The quadrivalent HPV vaccine is effective prophylaxis against HPV-related external genital lesions in young men. Evid Based Med 16: 157–158.
  19. 19. Smith MA, Lew JB, Walker RJ, Brotherton JM, Nickson C, et al. (2011) The predicted impact of HPV vaccination on male infections and male HPV-related cancers in Australia. Vaccine 29: 9112–9122.
  20. 20. Vriend HJ, Boot HJ, van der Sande MA (2012) Type-Specific Human Papillomavirus Infections Among Young Heterosexual Male and Female STI Clinic Attendees. Sex Transm Dis 39: 72–78.
  21. 21. Molijn A, Kleter B, Quint W, van Doorn LJ (2005) Molecular diagnosis of human papillomavirus (HPV) infections. J Clin Virol 32 Suppl 1S43–51.
  22. 22. Scherpenisse M, Mollers M, Schepp RM, Boot HJ, de Melker HE, et al.. (2012) Seroprevalence of seven high-risk HPV types in The Netherlands. Vaccine jn press.
  23. 23. Fitzmaurice GM, Laird NM, Ware JH (2004) Applied Longitudinal Analysis. New York: Wiley-Interscience.
  24. 24. Xue X, Gange SJ, Zhong Y, Burk RD, Minkoff H, et al. (2010) Marginal and mixed-effects models in the analysis of human papillomavirus natural history data. Cancer Epidemiol Biomarkers Prev 19: 159–169.
  25. 25. Carey V, Zeger SL, Diggle PJ (1993) Modelling multivariate binary data with alternating logistic regressions. Biometrika 80: 517–526.
  26. 26. Markowitz LE, Sternberg M, Dunne EF, McQuillan G, Unger ER (2009) Seroprevalence of human papillomavirus types 6, 11, 16, and 18 in the United States: National Health and Nutrition Examination Survey 2003–2004. J Infect Dis 200: 1059–1067.
  27. 27. Nielsen A, Kjaer SK, Munk C, Iftner T (2008) Type-specific HPV infection and multiple HPV types: prevalence and risk factor profile in nearly 12,000 younger and older Danish women. Sex Transm Dis 35: 276–282.
  28. 28. Nielson CM, Harris RB, Dunne EF, Abrahamsen M, Papenfuss MR, et al. (2007) Risk factors for anogenital human papillomavirus infection in men. J Infect Dis 196: 1137–1145.
  29. 29. Nyitray AG, da Silva RJ, Baggio ML, Lu B, Smith D, et al. (2011) The prevalence of genital HPV and factors associated with oncogenic HPV among men having sex with men and men having sex with women and men: the HIM study. Sex Transm Dis 38: 932–940.
  30. 30. Shikary T, Bernstein DI, Jin Y, Zimet GD, Rosenthal SL, et al. (2009) Epidemiology and risk factors for human papillomavirus infection in a diverse sample of low-income young women. J Clin Virol 46: 107–111.
  31. 31. de Villiers EM, Fauquet C, Broker TR, Bernard HU, zur Hausen H (2004) Classification of papillomaviruses. Virology 324: 17–27.