WOLVES is a long-term, population-based, prospective, cohort study to measure changes in the prevalence of HPV infections and associated diseases after the introduction of HPV vaccines among the general female population from 2009 to 2020. Women with a first residency (registered address) in the Wolfsburg area (about 125,000 citizens) were invited by letter to attend cervical cancer screening in one of 18 private gynecologic practices every year. Voluntary participation in the WOLVES study was open to women irrespective of existing HPV status and all participants gave written informed consent. The study was approved by the local ethics committee (Bo/07/2009).
Participants were recruited in cohorts, defined according to their birth date (1983/84, 1988/89, and 1993/94) and the year of analysis (2009/10 or 2014/15). Four cohorts are reported: cohort 1 = participants born in 1983/84 and analyzed in 2009/10; cohort 2 = participants born in 1988/89 and analyzed in 2009/10; cohort 3 = participants born in 1988/89 and analyzed in 2014/15; and cohort 4 = participants born in 1993/94 and analyzed in 2014/15 (Figure 1). At the time of sample collection, participants were aged 19−22 years in cohorts 2 and 4, and aged 24−27 years in cohorts 1 and 3. Women born in cohort 1 had a single cross-sectional medical examination in 2009/10, whereas those in cohorts 2 and 3 had annual medical examinations from 2009/10 to 2014/15. Participants in cohort 4 were analyzed in 2014/15 and will have annual visits until the end of 2020.
Participants of WOLVES do not differ from the general population of the same age, living in Wolfsburg area city, in terms of education, migration background, or parity. Town registry does not include any medical information, which is the main limitation of this study.
All participants were asked to complete a standardized medical questionnaire that included questions on education, birth country, marital status, pregnancies, parity, contraception, smoking, number of sex partners, age at first intercourse, history of abnormal Pap smears, vaccination status, sexually transmitted infections, and genital warts. During each medical examination, women had a pelvic examination with visualization of the uterine cervix. Pap smears were taken using a spatula and endocervical brush. A second sample was then obtained with a Qiagen Cervical Sampler (Medscan, Uppsala, Sweden) and suspended in 1 ml of specimen transport medium (STM/ Qiagen Inc., Hilden, Germany) for HPV DNA testing. From 2011, samples were taken for liquid cytology using a Cervex brush and ThinPrep (Hologic) transport medium. The gynecologist collected data on HPV vaccination status by checking the certificate of vaccination, which includes all vaccinations. For HPV vaccinations, the certificate lists batch numbers and the date of vaccination. A fully vaccinated status, based on number of doses administered according to the recommended schedule, was defined as three doses given at months 0, 2, and 6 (5−13 months).
Women were referred to colposcopy if they had genital warts or an abnormal high-grade Pap smear, or they had an abnormal Pap smear classified as borderline or low-grade and tested positive for high-risk HPV infection. All colposcopy examinations were done at the Klinikum Wolfsburg, the single certified colposcopy unit in the region.
Women were not included if they were not living in the Wolfsburg area, had a diagnosis of cervical or genital cancer, another malignancy, an organ transplant, or were undergoing an immunosuppressive therapy.
Diagnosis of genital warts was based on colposcopy and histological biopsies. Genital warts were classified as: (i) typical condylomata acuminata for lesions with typical acuminate morphology, which persist on mucosal epithelium or cornifying (cutaneous) skin and have highly vascularized dermal cores producing typical punctuated or cauliflower-like patterns; typical condylomata acuminata tend to be pigmentless and are mostly seen on pigmented skin (11); (ii) flat condylomata for genital papillomas with a more hyperkeratotic and pigmented surface and flat condylomata with a smooth surface and non-pigmented papules; or (iii) seborrheic wart-like lesions of the cutaneous skin of the external anogenital area. The focus of this paper is on low-risk HPV types and associated disease therefore we do not report on vulvar intraepithelial neoplasia, Bowenoid papulosis, Naevi or Mollusca contagiosa (2). These patients were excluded from the analysis of genital warts.
HPV DNA testing
HPV testing was done using the Hybrid Capture 2 assay (HC2; Qiagen Inc., Hilden Germany). Cervical Pap smear samples were analyzed for the presence of the 13 high-risk HPV types (16, 18, 31, 33, 35, 39, 45, 51, 52, 56, 58, 59, and 68) and the five low-risk HPV types (6, 11, 42, 43, and 44). Samples were diagnosed as positive if they attained or exceeded the FDA-approved threshold of 1.0 pg HPV DNA/ml, which corresponds to 1.0 relative light units (2).
HPV genotyping using SPF-10-PCR, followed by Reverse Line Probe Assay LiPA Extra (SPF-10-PCR) (12) was done on all HC2-positive samples and 10% of HC2-negative samples. Total DNA was isolated from the cervical samples using a QiaSymphony device (Qiagen, Hilden, Germany) and analyzed by the INNO-LiPA Extra HPV prototype assay (FujiReBioGent, Belgium) according to the manufacturer´s instructions. All PCR reactions were done with 10 µl input DNA in a final of 50 µl using reagents provided by Innogenetics according to the following sequence: 10 min 37°C, 9 min 94°C, and 40 cycles of 30 sec of denaturation at 94°C, followed by 45 sec at 52°C unchanging temperature and 45 sec of extension at 72°C run on a MJ Thermocycler PCT 200. The PCR product was then denaturated and a 10 µl aliquot hybridized to one strip at 49°C for 60 min, followed by multiple washing steps. The INNO-LiPA Extra test allowed identification of 13 established high-risk HPV types (16, 18, 31, 33, 35, 39, 45, 51, 52, 56, 58, 59, and 68), five known or potential high-risk types (26, 53, 66, 73, and 82), seven low-risk HPV types (6, 11, 40, 43, 44, 54, and 70), additional non-differentiated HPV types, and types with undefined risk (69, 71, and 74). The strips were analyzed on a flatbed scanner using LiRAS prototype software (Innogenetics, Inc), which displays the patterns and relative intensity of positive bands as arbitrary grey-tone values between 0.1 and 1.0 and allows direct data transfer to Excel spreadsheets. All HC2-low-risk positive samples, which were negative in the LiPA Extra, were retested with the CP4/5 PCR (13). The positive PCR products were purified with the Qiagen Gel Extraction Test QIAquick and then directly sequenced using the CP4 primer.
Study endpoints
The predefined endpoints were the rates of full HPV vaccination coverage (three doses), the prevalence of HPV 6 and 11 infection, and of genital warts (condylomata acuminata). We reported changes in these prevalence rates according to vaccination status and lifestyle factors (sexual history and smoking history) among participants in the four cohorts.
Statistical analysis
Vaccination coverage and prevalence rates of HPV 6/11 infection and genital warts are reported for participants aged 19−22 years (cohorts 2 and 4) and 24−27 years (cohorts 1 and 3) at the time of sample collection. Statistical analyses were done to compare similarly aged participants (cohort 1 versus 3 and cohort 2 versus 4) using 2 x 2 contingency tables for the prevalence of HPV 6/11 infection (positive versus negative) and genital warts (no versus yes). The Röhmel-Mansmann unconditional exact test was used to test for difference, with a two-side alpha of 0.05. Statistical analyses were performed using Testimate V.6.5.14.