- Research article
- Open Access
- Open Peer Review
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 Infectious Diseases volume 16, Article number: 502 (2016)
Acquisition of genital human papillomavirus (HPV) infection is common among the young, sexually active population. Genital HPV infections do not always lead to seroconversion. We aimed to assess the association between cervico-vaginal high risk (hr) HPV DNA and type-specific antibodies in an ethnically diverse cohort of young women.
Women of Dutch, South-Asian Surinamese, African Surinamese, Ghanaian, Moroccan and Turkish origin participating in a large-scale multi-ethnic population-based cohort (the HELIUS study) provided vaginal self-samples and blood samples, and completed a questionnaire regarding demographics, lifestyle and sexual health. Vaginal swabs were tested for HPV using the highly sensitive SPF10-PCR DEIA/LiPA25 system (version1). Serum samples were tested for type-specific L1 antibodies against 7 hrHPV types (16,18,31,33,45,52,58) with multiplex serology. We assessed the association between vaginal HPV DNA and type-specific seropositivity with logistic and linear regression, using generalized estimating equations (GEE). We determined whether this association varies by ethnicity by adding an interaction term.
We selected 532 women who completed the questionnaire, provided a vaginal swab and a blood sample. Their median age was 27 years (interquartile range 24–31 years). Prevalence of DNA of any of the 7 hrHPV was 22 %; HPV-52 was most common. Prevalence of antibodies against one or more hrHPV types was 24 %; HPV-16 seropositivity was most common. In multivariable logistic regression analysis using GEE, adjusting for other determinants, vaginal HPV DNA detection was associated with type-specific HPV seropositivity (OR 1.53, 95 % CI 1.06-2.20). In multivariable linear regression analysis using GEE, the geometric mean of type-specific antibody reactivity was 1.15 (95 % CI 1.04-1.27) times higher in women positive for HPV DNA compared to HPV DNA-negative women. There was little evidence that ethnicity modified the association between HPV DNA, and type-specific seropositivity, or with antibody reactivities (p = 0.47 and p = 0.57, respectively).
In this multi-ethnic group of young women in Amsterdam, cervico-vaginal hrHPV DNA detection was an independent determinant of type-specific HPV seropositivity.
Due to increased globalization, European populations are becoming more ethnically diverse, and health risks and outcomes differ between ethnic groups . Acquisition of human papillomavirus (HPV) infection is common among the young, sexually active population. Oncogenic human papillomaviruses cause the majority of cervical, anal, and penile cancer, and in some countries also of oropharyngeal cancer . The incidence of cervical carcinoma differs between ethnic groups in the Netherlands, with higher incidence in those of the non-Dutch origin .
Most infections with oncogenic HPV are cleared spontaneously and pass without malignant sequelae. Estimates of proportions of women who seroconvert after infection range from 55–85 %, depending on study design and methods employed [4–6]. The median time to seroconversion after genital infection with HPV ranged from 8.3–31.3 months between studies [4–6]. In a large multi-national study 79 % of cases of cervical HPV DNA cleared within 24 months . Concordance of genital HPV detection and serum antibodies is not expected per se, because the former may constitute current infection, and the latter past infection or infection of longer duration. Some studies, examining the concordance between cervico-vaginal HPV DNA and type-specific seropositivity showed positive associations, while others showed no significant relationship [8–15].
We aimed to assess the association between prevalent cervico-vaginal HPV DNA of seven high risk HPV (hrHPV) types and type-specific antibodies in serum, and to evaluate whether this association varies by ethnicity.
The HEalthy LIfe in an Urban Setting (HELIUS) study is a multi-ethnic cohort study conducted in Amsterdam which has been described in detail elsewhere . Briefly, the study was initiated in 2011 and includes people aged 18–70 years from the 6 largest ethnic groups in Amsterdam, i.e. those of Dutch, South-Asian Surinamese, African Surinamese, Ghanaian, Moroccan and Turkish origin. Participants were randomly selected per ethnicity from the Amsterdam municipal registers. Data were collected by questionnaire and physical examination; collection of biological samples took place during physical examination. Information on demographic characteristics (age, education, marital status), health-related characteristics (parity, HPV vaccination, smoking), and sexual behavior (age at sexual debut, number of lifetime male sex partners, type of sexual partner in the preceding six months) were obtained via questionnaire.
The selection process for this study has been described in detail elsewhere (Alberts et al. 2015, manuscript submitted for review). Briefly, we selected all women aged 18–34 years from the HELIUS participants who were enrolled until December 2013, who completed the questionnaire on demographics, health behavioral characteristics, and sexual behavior, provided a self-collected vaginal swab, and of whom a blood sample was available (n = 959). We randomly selected a maximum of 7 women per life year per ethnic group, resulting in a selection of 610 women. Participants who were not tested for HPV antibodies (n = 60), who reported HPV vaccination (n = 10) and participants with an invalid HPV DNA test result (n = 8) were excluded, resulting in a study sample of 532 women. The study protocol was approved by The Academic Medical Center Ethics Review Board (Amsterdam, the Netherlands), and all participants provided written informed consent (reference number METC 10/100# 17.10.1729).
HPV DNA detection and genotyping
Self-collected vaginal swabs (COPAN Italia, Brescia, Italy) were stored at 4 °C for a maximum of 5–6 days at the research location and sent to the Public Health Service of Amsterdam where they were stored at −20 °C until further analysis. Swabs were eluted in 600 μL PBS and DNA was extracted from 200 μL of each specimen using the isopropanol method  with an end volume of 100 mL. The presence of HPV in the swabs was assessed by using the highly sensitive SPF10-PCR DEIA/LiPA25 system version 1 (DDL Diagnostics Laboratory, Rijswijk, the Netherlands) . Briefly, with this technique a 65 base-pair open reading frame of the L1 region was amplified with SPF10 primers resulting in biotinylated amplicons. The amplification products were tested for HPV DNA with a DNA enzyme immuno-assay by use of hybridization with a cocktail of probes that recognized at least 54 HPV types. Optical densities were measured and compared with cut-off values from internal controls. Positive samples were genotyped by a reverse hybridization line probe assay which allows simultaneous detection of HPV-6, 11, 16, 18, 31, 33–35, 39, 40, 42–45, 51–54, 56, 58, 59, 66, 68/73, 70, and 74. Prevalence of vaginal HPV DNA in this population and determinants of vaginal HPV DNA detection have been reported elsewhere (Alberts et al. 2015, manuscript submitted for review).
HPV antibody detection
Serum samples were stored at −80 °C at the Academic Medical Center of Amsterdam until sent for analysis to the German Cancer Research Center, Heidelberg, Germany. Antibodies to the L1 proteins of HPV-16, 18, 31, 33, 45, 52 and 58 were analyzed by multiplex serology assay based on glutathione S-transferase (GST) fusion protein capture on fluorescent beads (SeroMAP Microspheres; Luminex) . Antibody quantity was determined as the median R-phycoerythrin fluorescence intensity (MFI) from at least 100 beads of the same internal color. MFI values were dichotomized as antibody positive or negative based on predetermined cut-off values for each antibody by analyzing MFI values obtained from Korean HPV DNA-negative self-reported virgin females using an algorithm of the mean plus 5 SDs after exclusion of outliers . This resulted in cut-offs varying from 394–712 MFI.
Differences in baseline characteristics between ethnic groups were assessed using the Pearson’s chi-squared or the Kruskal-Wallis test as appropriate. The statistical analysis was restricted to HPV types 16, 18, 31, 33, 45, 52 and 58.
First, we assessed type-specific HPV seropositivity on a dichotomous scale among both HPV DNA positive women and HPV DNA negative women per individual HPV type by performing standard univariable logistic analyses. As a woman can be infected with multiple HPV types simultaneously we also analyzed data on all 7 HPV types with univariable and multivariable logistic regression models using generalized estimating equations (GEE). We used a binomial distribution and an exchangeable correlation structure with robust variance.
We then investigated the association between cervico-vaginal HPV DNA and type-specific antibody reactivities using linear regression. For this purpose, we transformed antibody reactivities to their natural logarithm. We first undertook standard univariable and multivariable linear regression for individual HPV types. Next, we proceeded with univariable and multivariable linear regression using GEE. Because the mean and standard deviation of antibody reactivities of the 7 hrHPV types differed, we standardized the natural logarithm of antibody reactivity levels by calculating Z-scores. We used a Gaussian distribution and an exchangeable correlation structure, with robust variance. Because we included a log-transformed outcome variable in the linear regression analysis, we back-transformed the regression coefficients in order to ease interpretation. As a result, the coefficients can be interpreted as a multiplicative factor of the geometric mean; in other words the geometric mean of an antibody reactivity in women who are positive for cervico-vaginal HPV DNA, is equal to that factor times the geometric mean of women who are HPV DNA negative.
All multivariable models were adjusted for an a priori defined set of variables, i.e. age, ethnicity, marital status, smoking in pack years, and number of lifetime male sexual partners. We included age as a continuous variable using restricted cubic splines with 4 knots. All variables other than age were included as categorical variables. No backward selection was done.
Only one sexual risk variable could be included in the model because of multicollinearity between these variables. Also, only smoking in pack years and not smoking status was included in the models because of multicollinearity. For all other variables multicollinearity did not appear to be a concern.
As a last step we investigated whether the association between HPV seropositivity and HPV DNA varied by ethnicity. To do this we tested for interaction by adding an interaction term to the multivariable model. Stata Statistical Software release 13.1 (StataCorp, College Station, TX, USA) was used for statistical analyses.
A total of 532 women were included in the analyses (Table 1). The average age of the participants did not differ between ethnic groups (median 27 years, interquartile range (IQR) 23–30 years). All other baseline characteristics differed across ethnic groups; we found significant differences in education, marital status, and smoking. Furthermore, lifetime male sexual partners, age of sexual debut and type of sexual contact in in the last six months varied significantly by ethnicity. Most women reported having one steady sexual partner in the preceding six months.
Vaginal HPV DNA of at least one of the 7 hrHPV types was detected in 115 (22 %) women of whom some had multiple infections with different HPV types; DNA of HPV-52 was detected most often (Table 2). Of the 532 women, 130 (24 %) were seropositive for at least one of the 7 hrHPV types; HPV-16 had the highest seroprevalence (14 %). Among women positive for DNA of an individual HPV type, seropositivity for the corresponding HPV type ranged from 13–33 %. Among HPV DNA negative women, seropositivity ranged from 6–13 %.
For all 7 HPV types we found that HPV DNA positive women were seropositive for the corresponding type more often than HPV DNA negative women (indicated by ORs > 1 for all), but this was only significant for HPV-33, 52 and 58 (Table 2). In univariable logistic regression analysis using GEE cervico-vaginal HPV DNA detection was significantly associated with type-specific seropositivity (OR 1.70, 95 % CI 1.22 – 2.36) (Table 3). After adjustment for age, ethnicity, marital status, smoking in pack years, and number of lifetime male sexual partners the association remained positive and statistically significant (OR 1.53, 95 % CI 1.06 – 2.20).
The geometric mean of the antibody reactivities in MFI was higher among women who were positive for cervico-vaginal HPV DNA. This was the case for all 7 hrHPV types, though only significant for HPV-31, 52 and 58 (Table 4). When we adjusted for age, ethnicity, marital status, smoking in pack years, and lifetime number of male sexual partners in multivariable linear regression models, DNA detection of HPV-31 remained significantly associated with a higher antibody reactivity for this type. For all other HPV types the back-transformed regression coefficient was >1, but the association was not significant.
Next, we investigated the association between cervico-vaginal HPV DNA and standardized antibody reactivities with linear regression analysis using GEE. In univariable linear analysis with GEE women with prevalent cervico-vaginal HPV DNA had higher type-specific antibody reactivities than women without cervico-vaginal HPV DNA (back-transformed regression coefficient 1.15 [95 % CI 1.04 – 1.27]) (Table 5). In the multivariable analysis the back-transformed regression coefficient remained unchanged (1.15 [95 % CI 1.04 – 1.27]). This indicates that women with cervico-vaginal HPV DNA have higher type-specific antibody reactivities than women without cervico-vaginal HPV DNA, independent of other determinants.
We examined whether the effect of HPV DNA on seropositivity differed by ethnicity by adding an interaction term to the models. Neither in logistic, nor in linear regression analyses using GEE did we detect significant effect modification (p = 0.47 and p = 0.57, respectively).
Our results show that in this multi-ethnic group of young women in Amsterdam, cervico-vaginal HPV DNA detection is an independent risk factor of type-specific HPV seropositivity and that detection of cervico-vaginal HPV DNA is associated with higher antibody reactivities of type-specific antibodies.
Variables that describe sexual behavior directly or that could be considered a proxy for sexual behavior were strongly associated with higher type-specific HPV antibody reactivities and HPV seropositivity. The latter has also been shown by several other studies [8, 20, 21].
Since antibody formation in response to natural infection is a slow process that could extend to over a year, we expected in this young study population that concordance between cervico-vaginal DNA and antibodies would be low. Our study population however had a significant positive association even after adjustment for other well-known risk factors. Our study population most likely represents women who contracted their HPV infection relatively recently. A small proportion of these infections will become persistent and the majority will be cleared [22, 23].
With these results we aim to show that the interpretation of HPV (sero-)prevalences may not be as straightforward as sometimes is depicted, i.e. even though there is a significant positive correlation between HPV DNA and HPV seropositive participants, being HPV DNA positive does not imply that a person inevitably will be or will become HPV seropositive, while being HPV seronegative does not necessarily mean an individual has never been exposed. Concordance between HPV DNA and seropositivity therefore has little clinical application as a diagnostic tool. However it may be valuable in describing the natural history of HPV infection.
Because seropositivity to HPV is a measure of lifetime exposure to HPV we expected that age would be associated with HPV seropositivity. In our sample we found a significant positive association with age 19 and seropositivity relative to the reference age of 24 years, whereas the odds for seropositivity seem to not vary much with 29 and 34 years. A similar trend has been reported by Scherpenisse et al. (2012); seropositivity percentages rise sharply in the early twenties, then stabilize and do not vary much until 59 years of age .
Both significant and non-significant associations between cervical HPV DNA and type-specific seropositivity have previously been reported [8–15]. We found a significant positive association in our sample between cervico-vaginal HPV DNA and type-specific seropositivity expressed on a dichotomous scale, and further strengthened this conclusion by confirming it in an analysis with the outcome measure on a continuous scale. Our finding of a positive association between cervico-vaginal HPV DNA and type-specific antibodies is in agreement with several previous studies [9, 12–15]. We complement and expand available knowledge on HPV by including seven hrHPV types in our analysis, where most previous studies included fewer types.
Some previous studies differ in methodology to ours in several respects, in that they were conducted amongst women at higher risk for HPV infections (more sexual risk behavior or patients with cervical dysplasia or recruitment in STI clinics) [9, 13, 14]. In other publications, type-specific positive associations have also been found in young women that were recruited from populations not at high risk of HPV infection [12, 15].
The fact that we could only include women who provided a cervico-vaginal swab might have introduced some selection bias into our study. Alberts et al. (2015) have shown that women of Turkish and Moroccan origin from the same study population who declined to provide a sample exhibited less sexual risk behavior (Alberts et al., 2015, submitted for review). However, no significant differences were observed between women who did and did not provide a self-swab for the other ethnicities suggesting limited bias for these groups. The ethnic diversity of our sample, and the concurrent analysis of 7 hrHPV types are strengths of this study. Because seropositivity is a measure of lifetime exposure to HPV and not necessarily the result of a current infection, an association between cervico-vaginal HPV DNA and seropositivity can be difficult to interpret.
In conclusion, we observed a positive association between cervico-vaginal DNA detection and type-specific seropositivity. Moreover, we demonstrated higher type-specific antibody reactivities in women who are HPV DNA positive compared to HPV DNA negative women in a multi-ethnic urban population of young women. These associations were independent of demographic and behavioral factors. There was little evidence that ethnicity constituted a significant effect modifier in our population which suggests that effect modifiers other than genetic variation are at play that result in differing HPV infection outcomes in ethnic minorities.
generalized estimating equations
HEalthy LIfe in an Urban Setting
high-risk human papillomavirus
median fluorescence intensity
Keppel KG, Pearcy JN, Heron MP. Is there progress toward eliminating racial/ethnic disparities in the leading causes of death? Public Health Rep. 2010;125(5):689–97.
de Martel C, Ferlay J, Franceschi S, Vignat J, Bray F, Forman D, Plummer M. Global burden of cancers attributable to infections in 2008: a review and synthetic analysis. Lancet Oncol. 2012;13(6):607–15.
Arnold M, Aarts MJ, van der Aa M, Visser O, Coebergh JW. Investigating cervical, oesophageal and colon cancer risk and survival among migrants in The Netherlands. Eur J Public Health. 2013;23(5):867–73.
Carter JJ, Koutsky LA, Hughes JP, Lee SK, Kuypers J, Kiviat N, Galloway DA. Comparison of human papillomavirus types 16, 18, and 6 capsid antibody responses following incident infection. J Infect Dis. 2000;181(6):1911–9.
Ho GY, Studentsov YY, Bierman R, Burk RD. Natural history of human papillomavirus type 16 virus-like particle antibodies in young women. Cancer Epidemiol Biomarkers Prev. 2004;13(1):110–6.
Steele J, Collins S, Wen K, Ryan G, Constandinou-Williams C, Woodman CB. 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. 2008;15(9):1387–90.
Jaisamrarn U, Castellsague X, Garland SM, Naud P, Palmroth J, Del Rosario-Raymundo MR, Wheeler CM, Salmeron J, Chow SN, Apter D, et al. Natural history of progression of HPV infection to cervical lesion or clearance: analysis of the control arm of the large, randomised PATRICIA study. PloS One. 2013;8(11), e79260.
Viscidi RP, Ahdieh-Grant L, Clayman B, Fox K, Massad LS, Cu-Uvin S, Shah KV, Anastos KM, Squires KE, Duerr A, et al. Serum immunoglobulin G response to human papillomavirus type 16 virus-like particles in human immunodeficiency virus (HIV)-positive and risk-matched HIV-negative women. J Infect Dis. 2003;187(2):194–205.
Vriend HJ, Bogaards JA, van der Klis FR, Scherpenisse M, Boot HJ, King AJ, van der Sande MA, Medical Microbiological Laboratories MHS. Patterns of human papillomavirus DNA and antibody positivity in young males and females, suggesting a site-specific natural course of infection. PloS One. 2013;8(4), e60696.
Paaso AE, Louvanto K, Syrjanen KJ, Waterboer T, Grenman SE, Pawlita M, Syrjanen SM. Lack of type-specific concordance between human papillomavirus (HPV) serology and HPV DNA detection in the uterine cervix and oral mucosa. J Gen Virol. 2011;92(Pt 9):2034–46.
Castro FA, Dominguez A, Puschel K, Van De Wyngard V, Snijders PJ, Franceschi S, Pawlita M, Ferreccio C. Serological prevalence and persistence of high-risk human papillomavirus infection among women in Santiago, Chile. BMC Infect Dis. 2014;14:361.
Clifford GM, Shin HR, Oh JK, Waterboer T, Ju YH, Vaccarella S, Quint W, Pawlita M, Franceschi S. Serologic response to oncogenic human papillomavirus types in male and female university students in Busan, South Korea. Cancer Epidemiol Biomarkers Prev. 2007;16(9):1874–9.
Bontkes HJ, de Gruijl TD, Walboomers JM, Schiller JT, Dillner J, Helmerhorst TJ, Verheijen RH, Scheper RJ, Meijer CJ. Immune responses against human papillomavirus (HPV) type 16 virus-like particles in a cohort study of women with cervical intraepithelial neoplasia. II. Systemic but not local IgA responses correlate with clearance of HPV-16. J Gen Virol. 1999;80(Pt 2):409–17.
Wu X, Zhang C, Feng S, Liu C, Li Y, Yang Y, Gao J, Li H, Meng S, Li L, et al. Detection of HPV types and neutralizing antibodies in Gansu province, China. J Med Virol. 2009;81(4):693–702.
Skjeldestad FE, Mehta V, Sings HL, Ovreness T, Turpin J, Su L, Boerckel P, Roberts C, Bryan J, Jansen KU, et al. Seroprevalence and genital DNA prevalence of HPV types 6, 11, 16 and 18 in a cohort of young Norwegian women: study design and cohort characteristics. Acta Obstet Gynecol Scand. 2008;87(1):81–8.
Stronks K, Snijder MB, Peters RJ, Prins M, Schene AH, Zwinderman AH. Unravelling the impact of ethnicity on health in Europe: the HELIUS study. BMC Public Health. 2013;13:402.
van Houdt R, Bruisten SM, Geskus RB, Bakker M, Wolthers KC, Prins M, Coutinho RA. Ongoing transmission of a single hepatitis B virus strain among men having sex with men in Amsterdam. J Viral Hepat. 2010;17(2):108–14.
Molijn A, Kleter B, Quint W, van Doorn LJ. Molecular diagnosis of human papillomavirus (HPV) infections. J Clin Virol. 2005;32 Suppl 1:S43–51.
Waterboer T, Sehr P, Michael KM, Franceschi S, Nieland JD, Joos TO, Templin MF, Pawlita M. Multiplex human papillomavirus serology based on in situ-purified glutathione s-transferase fusion proteins. Clin Chem. 2005;51(10):1845–53.
Markowitz LE, Sternberg M, Dunne EF, McQuillan G, Unger ER. 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. 2009;200(7):1059–67.
Nielsen A, Kjaer SK, Munk C, Iftner T. 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. 2008;35(3):276–82.
Hildesheim A, Schiffman MH, Gravitt PE, Glass AG, Greer CE, Zhang T, Scott DR, Rush BB, Lawler P, Sherman ME, et al. Persistence of type-specific human papillomavirus infection among cytologically normal women. J Infect Dis. 1994;169(2):235–40.
Castle PE, Schiffman M, Herrero R, Hildesheim A, Rodriguez AC, Bratti MC, Sherman ME, Wacholder S, Tarone R, Burk RD. A prospective study of age trends in cervical human papillomavirus acquisition and persistence in Guanacaste, Costa Rica. J Infect Dis. 2005;191(11):1808–16.
Scherpenisse M, Mollers M, Schepp RM, Boot HJ, Meijer CJ, Berbers GA, van der Klis FR, de Melker HE. Changes in antibody seroprevalence of seven high-risk HPV types between nationwide surveillance studies from 1995–96 and 2006–07 in The Netherlands. PloS One. 2012;7(11), e48807.
We gratefully acknowledge the AMC Biobank for their support in biobank management and high-quality storage of collected samples. We are most grateful to the participants of the HELIUS study and the management team, research nurses, interviewers, research assistants and other staff who have taken part in gathering the data of this study.
Members of the HPV HELIUS study group (R. A. Vos, H. Borgdorff, J. van Bergen, S. Bruisten, S.E. Geerlings, R. van Houdt, S. Morré, H. J.C. de Vries, and J. van de Wijgert) are gratefully acknowledged for their advice and suggestions on study design.
The HELIUS study is conducted by the Academic Medical Center Amsterdam and the Public Health Service of Amsterdam. Both organizations provided core support for HELIUS. The HELIUS study is also funded by the Dutch Heart Foundation (2010 T084), the Netherlands Organization for Health Research and Development (ZonMw) (200500003), and the European Union (FP-7) (278901). This particular study was funded by grant 7115 0001 and 204005002 (Academic Collaborative for Public Health) from ZonMw and an internal grant of the Public Health for Research & Development, the Netherlands.
Availability of data and materials
Data are available from the HELIUS study, a third party. Dr. Snijder and Dr. Prins are affiliated with the HELIUS study and are co-authors of this paper in accordance with the HELIUS requirements for collaboration. Dr. Snijder is the Scientific Coordinator of HELIUS and may be contacted with further questions (firstname.lastname@example.org). Additionally, researchers interested in further collaboration with HELIUS may see the following URL: http://www.heliusstudy.nl/nl/researchers/collaboration.
AK carried out statistical analyses and drafted the manuscript. LC and CJA contributed to data management and statistical analysis. MSvdL supervised statistical analyses. MBS and CJA ensured quality of the dataset. WV carried out analyses of HPV DNA in the swabs. TW and AM were responsible for the HPV serology. MBS, MP and MSvdL designed the HELIUS study. MSvdL CJA and AK designed the current study. All authors participated in review of the manuscript, and saw and approved the last version.
The institution of M. F. Schim van der Loeff received study funding from Sanofi Pasteur MSD; he is a co-investigator in a Merck-funded investigator-initiated study; he is an investigator on a Sanofi Pasteur MSD sponsored trial; he served on a vaccine advisory board of GSK; his institution received in-kind contribution for an HPV study from Stichting Pathologie Onderzoek en Ontwikkeling (SPOO); his institution receives research funding from Janssen Infectious Diseases and Vaccines. Other authors: no conflicts of interest.
Ethics approval and consent to participate
The study protocol was approved by The Academic Medical Center Ethics Review Board (Amsterdam, the Netherlands), and all participants provided written informed consent (reference number METC 10/100# 17.10.1729).