Skip to main content
  • Research article
  • Open access
  • Published:

Monitoring vaccine and non-vaccine HPV type prevalence in the post-vaccination era in women living in the Basilicata region, Italy

Abstract

Background

A large free-of-charge quadrivalent HPV (qHPV) vaccination program, covering four cohorts annually (women 11, 14, 17 and 24 years), has been implemented in Basilicata since 2007. This study evaluated vaccine and non-vaccine HPV prevalence 5-7 years post-vaccination program implementation in vaccinated and unvaccinated women.

Methods

This population-based, cross-sectional study was conducted in the public screening centers of the Local Health Unit in Matera between 2012 and 2014. Cervical samples were obtained for Pap and HPV testing (HC2, LiPA Extra® assay) and participants completed a sociodemographic and behavioral questionnaire. Detailed HPV vaccination status was retrieved from the official HPV vaccine registry. HPV prevalence was described overall, by type and vaccination status. The association between HPV type-detection and risk/protective factors was studied. Direct vaccine protection (qHPV vaccine effectiveness [VE]), cross-protection, and type-replacement were evaluated in cohorts eligible for vaccination, by analyzing HPV prevalence of vaccine and non-vaccine types according to vaccination status.

Results

Overall, 2793 women (18-50 years) were included, 1314 of them having been in birth cohorts eligible for the HPV vaccination program (18- to 30-year-old women at enrolment). Among the latter, qHPV vaccine uptake was 59% (at least one dose), with 94% completing the schedule; standardized qHPV type prevalence was 0.6% in vaccinated versus 5.5% in unvaccinated women (P <0.001); adjusted VE against vaccine type infections was 90% (95% CI: 73%-96%) for all fully vaccinated women and 100% (95% CI not calculable) in women vaccinated before sexual debut. No statistically significant difference in overall high-risk HPV, high-risk non-vaccine HPV, or any single non-vaccine type prevalence was observed between vaccinated and unvaccinated women.

Conclusions

These results, conducted in a post-vaccine era, suggest a high qHPV VE and that a well-implemented catch-up vaccination program may be efficient in reducing vaccine-type infections in a real-world setting. No cross-protective effect or evidence of type-replacement was observed a few years after HPV vaccine introduction.

Peer Review reports

Background

Genital human papillomavirus (HPV) infection is the most common sexually transmitted viral infection worldwide, with approximately 75% of sexually active men and women exposed to HPV during their lives [1].

HPV types that can infect the genital tract can be categorized as high risk (HR) or low risk (LR) according to the degree of risk associated with the development of cervical cancer. Infection with HR-HPV is found in virtually all cases of cervical cancer and is considered a necessary cause of invasive cervical cancer [2]. Among the 12 HPV types considered carcinogenic, HPV16 and HPV18 are the most aggressively oncogenic [3], and are estimated to cause about 70% of invasive cervical cancer cases [4]. Other oncogenic HPV types, mainly HPV 31, 33, 45, 52, and 58, are estimated to be responsible for another 18% of cervical cancers [4, 5].

In Italy, cervical cancer is the 15th most frequent cancer overall, and the third most frequent among women aged 15-44 years [6]. Every year, about 3000 women are newly diagnosed with cervical cancer, of whom 1000 die from the disease. The overall HPV prevalence for Italy in women with normal cytology has been estimated at 9.7% (95% confidence interval [CI]: 9.4-10) [6].

Large randomized controlled trials (RCTs) have shown HPV vaccines to be well tolerated and highly efficacious against vaccine-type persistent HPV infection and precancerous cervical lesions (vaccine efficacy 93% to 100%) [7, 8]. HPV immunization programs are now implemented in many countries, although the methods of delivery, uptake, and monitoring systems vary [9], as does the choice of vaccine [10]. A national HPV vaccination program had been started in all 21 regions of Italy by 2008.

In Basilicata, a region in southern Italy, an organized vaccination program was launched in 2007 with the prophylactic quadrivalent HPV (qHPV) vaccine Gardasil® (Merck Sharp & Dohme Corp, Whitehouse Station, NJ). A unique multiple-cohort vaccine implementation strategy was used, inviting yearly four birth cohorts (girls and women during the 11th, 14th, 17th, and 24th year of age) to be vaccinated free of charge in the Local Health Unit (LHU) of Matera. The first cohorts were born in 1996, 1993, 1990, and 1983, respectively. In addition, HPV vaccinations given through this organized program are systematically reported into a computerized vaccination registry. Cervical screening targeting women aged between 25 and 64 years is also well organized in Basilicata, and the adherence rate was on average 57% for 2013 to 2014 (Ronco G, personal communication).

This cross-sectional population-based study aimed to assess vaccine and non-vaccine HPV prevalence 5-7 years post-vaccination program implementation in vaccinated and unvaccinated women.

Methods

Study participants and setting

All 15 public screening centers managed by the LHU of Matera took part in the study. These centers used computerized HPV vaccine and screening registries. Between 21 May 2012 and 14 February 2014, all women aged 25-50 years participating in the Basilicata cervical cancer screening program were asked to take part in the study. Women aged 18-24 years, who were not age-eligible for the screening program, were identified from municipality registries and all were also invited by individual letters to participate in the study. Those not yet sexually active, hysterectomized women, and women who had undergone colposcopy in the 24 months prior to the enrolment visit were not eligible for enrolment. Eligible women were enrolled after signing an informed consent form.

At the enrolment visit, two cervical samples (one each for Pap and HPV testing) were obtained, and participants completed self-administered sociodemographic and behavioral questionnaires. HPV vaccination history, HPV and cytological data, sociodemographic and behavioral data, including sexual habits and age at first sexual intercourse, were recorded. Detailed individual HPV vaccination status (including date of each vaccination and number of doses received) was also retrieved from the official computerized HPV vaccine registry of the LHU of Matera. This information was used for the statistical analysis. HPV analyses and genotyping were centralized in an accredited laboratory at the Cancer Prevention and Research Institute (ISPO). Quality control of HPV testing, cytology, genotyping, and data input procedures were implemented throughout the study, which was conducted in accordance with the Good Epidemiological Practice guidelines and the principles of the Declaration of Helsinki, and was approved by the Matera LHU ethics committee.

Laboratory methods

HC2 assay

HPV detection was performed by Hybrid Capture® 2 (HC2; QIAGEN, Hilden, Germany) using HR-HC2 (to detect HPV 16, 18, 31, 33, 35, 39, 45, 51, 52, 56, 58, 59, and 68) and LR-HC2 (to detect HPV 6, 11, 42, 43, and 44) probes with a semi-automated procedure (Rapid Capture® system, QIAGEN). As recommended by the manufacturer, samples with a relative light units (RLU) to cutoff value (CO) ratio ≥1.0 were classified as HPV positive.

HPV genotyping

Genotyping was performed for all HC2 positive samples (HR and/or LR positive; RLU/CO ratio ≥1), all “borderline” samples (RLU/CO ratio <1 and ≥0.5), and 10% of HC2 negative samples (RLU/CO ratio <0.5) chosen randomly. After DNA extraction (QIAamp DNA Mini Kit, QIAGEN), samples were amplified by INNO-LiPA® HPV Genotyping Extra Amp (Innogenetics, Ghent, Belgium) followed by detection by reverse hybridization with HPV type-specific probes for 28 different HPV types: HPV 16, 18, 31, 33, 35, 39, 45, 51, 52, 56, 58, 59 (classified by International Agency for Research on Cancer [IARC] as “carcinogenic to humans”); HPV 26, 53, 66, 68, 70, 73, 82 (classified as “possibly carcinogenic to humans”) and HPV 6, 11, 40, 43, 44, 54, 69, 71, 74 (not classifiable as to its carcinogenicity to humans). Results were interpreted using LiRAS® for LiPA HPV software (Fujirebio Europe, Ghent, Belgium) followed by manual verification. When the INNO-LiPA test did not exclude the presence of HPV 39, 52, 68, or 56 co-infection, an additional type-specific polymerase chain reaction (PCR) was performed. PCR with specific HPV42 primers was performed on LR-HC2 positive or borderline samples (HR-HPV HC2 ratio <1.0, ≥0.5) that were INNO-LiPA negative, as HC2 may detect HPV42 that cannot be detected by INNO-LiPA.

Statistical methods

Sociodemographic characteristics (birth place, marital status, educational level, occupational status) collected in the self-administered questionnaire were compared with data from the Italian National Institute of Statistics (ISTAT; www.demo.istat.it) to assess whether the study population was representative of the general female population of Basilicata [11].

Crude, age-stratified, and age-standardized prevalence (with related 95% CI) were estimated overall, as well as for HR types, HR non-vaccine types, HR vaccine types, LR types, LR non-vaccine types, LR vaccine types, all vaccine types, and all non-vaccine types, as well as by HPV type. The female Italian National Population (1 January 2012 census) was used as the standard population. For women born in cohorts eligible for the HPV vaccination program (1983 to 1996, corresponding to women aged 18-30 years within the study recruitment period), overall and type-specific data were stratified by vaccination status.

Information on vaccination status, derived from the official computerized Matera LHU vaccine registry, was defined at the date of informed consent. Women who received the first vaccine dose before the date of informed consent were considered as vaccinated (women who received the first vaccine dose at the time of informed consent were considered as unvaccinated). A sensitivity analysis was performed considering as vaccinated only women with a complete vaccination schedule with all three doses administered before the consent date. Age at first sexual intercourse was self-reported and used to identify women who had been vaccinated before sexual debut.

Direct vaccine protection (effect on qHPV vaccine types), cross-protection (effect on types closely related to vaccine types; i.e., HPV 31, 33, and 45), and potential type-replacement (effect on HR non-vaccine types; HR non-16,18), HR non-vaccine non-cross-protective types (i.e., HR non-16, 18, 31, and HR-non 16, 18, 31, 33, 45), and each single HR non-vaccine HPV types (provided N >10) were evaluated in cohorts eligible for vaccination, by analyzing HPV prevalence of vaccine and non-vaccine types according to vaccination status. Comparison of standardized HPV prevalence in vaccinated and unvaccinated participants was obtained using a Poisson regression model adjusted for 5-year age groups. The difference between HPV prevalence by 5-year age groups was tested using the chi-squared test. Logistic regression was used to evaluate the relationship between cervical HPV infection and potential risk or protective factors, including vaccination. A univariate analysis was used to identify associations between HPV infection and each potential risk/protective factor. Factors for which an association with HPV infection was found (P ≤0.1) were included in the multivariate analysis. In the final multivariate logistic regression model, vaccine effectiveness (VE) of a full vaccination schedule was evaluated for the prevention of qHPV vaccine type detection as the complement of the adjusted odds ratio (aOR) of the vaccination status (VE = (1 – aOR)*100), in cohorts eligible to the vaccination program (18- to 30-year-old women). Similar analyses compared HPV prevalence in women vaccinated before sexual debut and unvaccinated women. Due to the low number of women vaccinated before sexual debut in the 25- to 30-year group, these comparisons were restricted to those aged 18-24 years. Statistical analyses were performed using STATA® (version 11 or 12) software (StataCorp LP, College Station, Texas, US).

Results

Description of study population

All 8349 women aged 18-24 years living in the Matera area were invited to participate in the study. Among them, only 906 presented themselves to screening centers, and of those, 895 accepted and satisfied the eligibility criteria for enrolment (i.e., 10.7% of women initially contacted). The proportion of women aged 25-50 years attending for cervical cancer screening who agreed to participate in the study was >98% (1909/1941). Overall, 2793 women were included in the statistical analysis (11 enrolled women were excluded either because sample dates were inconsistent or because samples had not been sent for testing), 894 and 1899 in the 18-24 years and 25- to 50-year groups, respectively.

Study population characteristics are described in Additional file 1: Table 1a and b. The demographic profile of the study population was generally representative of the overall Italian population, and comparison with ISTAT data revealed a slightly higher than expected proportion of women born outside Italy among participants in both age groups (18-24 years and 25-50 years), especially in the younger group (study: 8.4%; ISTAT: 4.1%). The proportion of married women in the 25- to 50-year group (study: 68.1%) was similar to that expected (ISTAT: 67.6%), but the proportion was higher in the younger age group (study: 7.9%; ISTAT: 4.7%). In the 25- to 50-year group, the proportion of participants with a high school diploma or higher was 11% higher than expected (study: 76.3%; ISTAT: 65.5%), whereas in the younger age group the proportion was quite similar (study: 84.6%; ISTAT: 87.6%).

Table 1 HPV vaccination status by age group

Table 1 shows age distribution by vaccination status. Among women included in the analysis, 1314 were in birth cohorts eligible for the HPV vaccination program since its introduction (18- to 30-year-old women in study). The proportion of vaccinated women differed significantly by age class (P <0.001): 67% in the 18- to 24-year group, and 42% in the 25- to 30-year group (overall 59% in those 18-30 years), whereas it was near zero among women older than 30 years.

Among vaccinated women, 38.4% (301/783) had received their first dose before sexual debut. This value reached to ≈50% in the 18- to 24-year group (296/597) and then dropped sharply in the older age group (2.8% for those 25-30 years and 0% for those 31-55 years).

Almost all women (94%) who had received at least one dose had gone on to complete the full vaccination schedule of three doses. Schedule completion was higher in women who had been vaccinated before sexual debut, with 97% receiving all three doses.

Within the 18- to 30-year subgroup, the demographic profiles of vaccinated and unvaccinated women differed for some parameters: vaccinated women were on average younger than unvaccinated women, more frequently born in Italy (98% vs 84%), single (92% vs 71%), with a higher educational level (84% vs 72%), and more frequently unemployed (see Table 2). However, the two groups were comparable for the main risk factors for HPV infection: the number of sexual partners both in the last 6 months and over their lifetime (84% of both vaccinated and unvaccinated women had one sexual partner in the last 6 months and ≈61% of both vaccinated and unvaccinated women had more than one lifetime partner), and they had similar smoking habits (P =0.268) (see Table 3).

Table 2 Comparison between vaccinated and unvaccinated women (aged 18-30 years) by sociodemographic characteristics (N =1314)
Table 3 Comparison between vaccinated and unvaccinated women (aged 18-30 years) by behavioral characteristics (N =1314)

Overall HPV prevalence and type distribution

The overall standardized HPV prevalence in the study population (women aged 18-50 years) was 9.5% for all types, 7.6% for HR-HPV types, and 3.2% for LR-HPV types (Additional file 1: Table S2). Overall, HPV prevalence was higher in younger than older age groups, with a significant trend of decreasing prevalence with increasing age; this trend was similar for all HR types and all LR types.

The overall standardized prevalence of non-vaccine HPV types was 8.0% (4.9% for HR non-vaccine types) and 2.4% for the four qHPV vaccine types combined (HPV 6, 11, 16, and 18) (2.3% for HR-qHPV vaccine types) (Additional file 1: Table S3). The decreasing trend with increasing age was not observed for qHPV vaccine types. Instead, the prevalence of qHPV vaccine types was lower in women age-eligible for vaccination (2.5% for women aged 18-24 years, and 3.1% for those aged 25-30 years for all four HPV vaccine types) compared with women aged 31-35 years (5.1%) (P =0.027).

HPV prevalence by vaccination status – restricted to cohorts eligible for vaccination program (18- to 30-year group)

Overall HPV prevalence

Among the cohorts eligible for the vaccination program, the most frequent HPV type in unvaccinated women was HPV16, followed by types 42, 52, 31, and 51. In vaccinated women, HPV42 was the most frequent type, followed by types 51, 39, and 52, with the same prevalence as types 31 and 58 (Fig. 1). HPV 16 was the only HR type for which a statistically significant difference was observed between vaccinated and unvaccinated women (Fig. 1).

Fig. 1
figure 1

HPV type prevalence by vaccination status (women aged 18–30 y)

No statistically significant differences were observed between unvaccinated and vaccinated women in the overall standardized HPV prevalence (17.1% vs 15.8%, P =0.39), HR-HPV type prevalence (13.9% vs 12.1%, P =0.21), and LR-HPV type prevalence (6.0% vs 6.1%, P =0.96) (Table 4). HR-HPV prevalence was significantly higher in unvaccinated women compared with women vaccinated before sexual debut (16.6% vs 9.5%, P =0.01 in the 18- to 24-year old group) (Table 5). This difference was maintained after adjustment for age at enrolment, marital status, smoking status, number of sexual partners (within the past 6 months and lifetime), and sexually transmitted disease (past or present), although this was no longer statistically significant (aOR =0.62, P =0.09) (Table 5).

Table 4 HPV standardized prevalence (%) stratified by vaccination status and adjusted odds ratio evaluating association between HPV prevalence and vaccination status (women 18-30 years)
Table 5 HPV standardized prevalence (%) stratified by vaccination status and adjusted odds ratio evaluating association between HPV prevalence and vaccination status (women 18-24 years)

qHPV vaccine type prevalence and vaccine effectiveness

The standardized prevalence of qHPV vaccine types (16, 18, 6, 11) was significantly higher in unvaccinated than in vaccinated women (5.5% vs 0.6%, P <0.001) (Table 4). This was true for both HR-qHPV vaccine types (5.2% vs 0.6%, P <0.001) and LR-qHPV vaccine types (0.7% vs 0.0%, P =0.02). The corresponding adjusted VE against vaccine type infections was 90% (95% CI: 73%-96%). No qHPV vaccine type infections were detected among women vaccinated before sexual debut (Table 5), corresponding to a VE of 100% against vaccine type infections in this population compared with unvaccinated women in the same age group (18-24 years).

Non-vaccine type prevalence: cross-protection and type-replacement

Even though we observed a higher prevalence point estimate for HPV 31, and HPV 31, 33, and 45 combined in unvaccinated versus vaccinated women, these differences were not statistically significant (2.1% vs 1.3%, P =0.58, and 2.7% vs 1.5%, P =0.30, respectively; Table 4). No statistically significant differences in HPV prevalence between unvaccinated and vaccinated women were observed in these cohorts for HR non-vaccine HPV types; HR non-16, 18, 31 HPV types; HR non-16, 18, 31, 33, 45 HPV types, or for any single non-vaccine HPV type (Table 4).

Similar results were observed when women vaccinated before sexual debut were compared with unvaccinated women of the same age-group and after adjustment for confounding (Tables 4 and 5). All analyses were repeated considering as vaccinated only women who completed all three doses before the consent date; similar results were found (data not shown).

Discussion

Our study is one of the largest in regard to the number of women enrolled in a post-vaccine era (1314 women age-eligible for the vaccination program were included). Due to this large sample size, the wide age range of women targeted by the vaccination program in Basilicata (11- to 24-year-old women each year), as well as the very good overall adherence to the vaccination program (about 59% vaccine uptake among all age ranges; 94% adherence to a full vaccination schedule), this region provides a unique setting to evaluate the real-life effect of a well-implemented vaccination program on viral HPV circulation, particularly on HPV types under the selective pressure of qHPV vaccine introduction.

This is also the first study analyzing overall HPV and type-specific prevalence in a post-vaccine era in Italy. By performing a comprehensive genotyping assay for surveillance, rather than focusing on HPV 16/18 only, we aimed to gain insights into cross-protection and type-replacement, in addition to vaccine effectiveness, as well as to provide a good baseline to monitor the impact of the possible future introduction of a nonavalent HPV vaccine.

A strength of this study was the high quality of the information gathered on individual vaccination status, which was issued from the official computerized vaccine registry of the LHU of Matera ensuring differentiation between vaccinated and unvaccinated individuals and allowing for a robust vaccine effectiveness evaluation. The collection of sociodemographic and behavioral data, using an anonymous questionnaire, is another strength of this study, as well as facilitated adjusted comparisons between vaccinated and unvaccinated women with several confounding factors for HPV infection in multivariate analysis. In particular, the self-declared age at first sexual intercourse allowed differentiation between women vaccinated before sexual debut, who had the potential to fully benefit from the vaccination, and those who were likely to have been exposed to HPV before vaccination.

One limitation of the study is the low participation rate (10.8%) of women not age-eligible for the organized screening program (aged 18-24 years). This participation rate, lower than that expected based on previous studies in Italy [12, 13], was mainly due to the large proportion of women in this age group who were attending university outside of Matera province for most of the year. This may have limited the number of the youngest study participants, leading to a lower than planned number of participants in cohorts age-eligible for the vaccination program. Thus, the results concerning the overall and type-specific HPV prevalence stratified by vaccination status, particularly the evaluation of type-replacement, must be interpreted with caution, as the study was not sufficiently powered to analyze the impact of vaccination on single HPV types. With a sample size equal to 1314 women (41% unvaccinated and 59% vaccinated) and a power of 80%, we were able to show a statistically significant reduction in prevalence of 2.8% for HPV 16/18 and 4.6% for HR-HPV types between unvaccinated and vaccinated women.

Another limitation is that our study was conducted in a post-vaccine era only. In absence of HPV prevalence data collected in a similar manner in a pre-vaccine era (same centers, same region, same analysis method, etc.), we cannot fully assess the impact of the vaccine program introduction. Nevertheless, our data seem overall consistent with HPV prevalence data observed in pre-vaccine eras elsewhere in Italy. For example, the HPV prevalences observed in 2007 in a screened population aged 25-65 years living in central and southern Italy [14] were very close to those observed in our population (18-50 years): 11.9% versus 11.3% for overall HPV; 9.2% versus 8.9% for HR-HPV, and 4.2% versus 4.1% for LR-HPV prevalence, respectively. The HPV prevalence observed among unvaccinated women aged 18-24 years in our study is consistent with figures from an HPV study conducted in southern Italy in 2006 to 2007 [13] and data from studio PREGIO conducted in 2007 to 2009 [12]. Lastly, it is important to underline that in the Basilicata region, the HPV vaccination program was limited only to female birth cohorts (11, 14, 17, and 24 years), and not for the whole population. Thus, also considering that the HPV prevalence data for unvaccinated women in our study were highly comparable to other Italian pre-vaccine studies, a herd immunity effect is unlikely.

Our population had a high prevalence of HPV42. It is important to note that most prevalence studies in Italy did not genotype LR types or the genotype systems used did not detect HPV 42. In the present study, all genotyping analyses were performed by INNO-LiPA Genotyping Extra, which is not able to detect HPV 42, one of the LR-HPV HC2 target types. Consequently, we evaluated HPV 42 prevalence by type-specific PCR, for all LR-HC2 positive samples that were negative for specific HPV genotypes as detected by INNO-LiPA. This may explain why the prevalence of HPV 42 is so high in our study.

The high prevalence of HPV infection in the youngest age groups and decreasing HPV prevalence with increasing age was expected from HPV natural history [15,16,17,18]. Still, it was interesting to note that this trend was not observed for qHPV vaccine types. Indeed, the prevalence of qHPV vaccine types was statistically significantly lower in women age-eligible for vaccination (2.5% for women aged 18-24 years; 3.1% for those aged 25-30 years) compared with the prevalence observed in the 31- to 35-year group (5.1%), indicating a probable impact of this vaccination program at the population level.

In our study, VE against qHPV vaccine types was 90% considering all vaccinated women (18-30 years), regardless of sexual exposure at the time of vaccination. These figures are close to those observed by Tabrizi [19], with a VE of 73% against qHPV vaccine-types among women aged 18-24 years (89% after adjusting for age and use of hormonal contraceptives), and more recently by Markowitz et al. (VE of 89% against qHPV vaccine types in women aged 14-24 years) [20]. HPV vaccination is not therapeutic [21], and women vaccinated after sexual debut may have already been exposed to vaccine HPV types [22]. In the present study, a substantial proportion of women (60%) were sexually active at vaccination. However, the overall VE was high, suggesting that a well-implemented catch-up vaccination program may be successful. No qHPV vaccine type infections were detected in women vaccinated before sexual debut, among which almost all (97%) had completed the full vaccination schedule. The high VE of 100% in this population confirms that the full benefit of vaccination is obtained when women are vaccinated with a complete vaccination schedule before sexual debut. These data also confirm the results observed in qHPV vaccine RCTs [7, 23] in a real-life setting. Due to the low number of women who had not completed a full vaccination schedule in our study (19 women received only one dose, 29 women received only two doses), we did not directly assess the VE of an incomplete vaccination schedule.

Although not statistically significant, lower HPV prevalence point estimates of HPV 31, 33, and 45 were observed in vaccinated versus unvaccinated women, suggesting that the qHPV vaccine might provide some degree of cross-protection against non-vaccine HPV types. This seems in line with evidence from other clinical [24,25,26] and observational studies [26], even though the limited power of the current study precludes firm conclusions on this point.

Vaccination is expected to reduce the prevalence of HPV vaccine types, but there is a theoretical concern that vaccine introduction may affect the distribution of other oncogenic types and induce type-replacement [27]. Similar to other observational studies [20, 28,29,30], our data provide no evidence of type-replacement a few years after HPV vaccine introduction. In their meta-analysis. Drolet et al. [29] found no statistically significant difference in the prevalence of HR-non-vaccine types between pre- and post-vaccination periods in any of the age groups studied (13-19 years and 20-24 years). In women aged 20-24 years, a small, but not statistically significant, increase in non-vaccine HR-HPV types (relative risk: 1.09, 95% CI: 0.98-1.22) was found to be negatively associated with increasing vaccination coverage (P =0.03) [29]. A recent meta-analysis by Mesher et al. [31] shows a significant increase in some HR non-vaccine types (HPV 39, HPV 52) in the post-vaccine era. These results, as well as the HPV types found to be significantly increased, are inconsistent between age groups, vaccine type (bivalent vs quadrivalent), and according to the methodological quality of studies, suggests random fluctuation. When only studies with a low potential for bias were considered, no statistically significant difference was observed. Having considered the potential limitations and uncertainties of this approach to evaluate type-replacement, Mesher et al. concluded that their meta-analysis provided no clear evidence for type-replacement [31]. This finding needs further exploration.

“Unmasking,” a potential diagnostic artifact of broad-spectrum assays, may give rise to an apparent increase in non-vaccine types following vaccination [27, 32], and may partially explain small observed increases. In the presence of high concentration types, broad-spectrum assays can miss types that are present in much lower concentrations. As HPV 16/18 are usually the most prevalent types, other types are less likely to be detected when present in the same samples. In vaccinated subjects, non-vaccine types that may not have been detected in the absence of vaccination due to the presence of HPV 16/18, may be “unmasked” by vaccination and thus be detected more broadly in vaccinated versus unvaccinated subjects.

“Unmasking” should be differentiated from type-replacement. It is equally important to evaluate the test used and ensure that it performs adequately. The World Health Organization HPV LabNet provides blinded “proficiency panels” designed to evaluate whether assays can detect a monoinfection equally well in the presence of other HPV types. Comparison of results from more than 100 laboratories worldwide that have used a variety of HPV assays has shown that underestimation of some HPV types, when other types are present in the same sample, is a problem for some assays [33, 34]. In this regard, continuing monitoring to ensure the adequate performance of assays used for surveillance is important.

The results of our study will be an important basis for the evaluation of future changes in the epidemiology of HPV infection in Italy, particularly in view of the upcoming 9-valent HPV vaccine that protects against five additional oncogenic HPV types [35]. Ongoing monitoring of the incidence of cervical pre-cancerous lesions and cancers, as well as other HPV-related lesions, such as anal and oropharyngeal cancers, will also be essential to understand the overall cost-effectiveness and population level benefit of HPV vaccines on cancer prevention. This study provides a preliminary evaluation of the potential effect of the qHPV vaccination program in Basilicata, showing that it significantly reduces the prevalence of qHPV vaccine types in vaccinated populations in a real-life setting compared with the prevalence in unvaccinated women.

Conclusion

The results of this study, conducted shortly after the implementation of the qHPV vaccination program in Basilicata, provide insight into the real-life population effect of HPV vaccination programs in those 18-50 years. They show a high VE of the qHPV vaccine on vaccine type prevalence, and found no evidence of type-replacement. This effectiveness is expected to be even more important in young girls who are HPV naive at vaccination, and will be examined in the routine cohort of 12-/13-year-olds, in whom vaccine uptake is higher and who will attend their first cervical screening visit in 2021.

Abbreviations

aOR:

Adjusted odds ratio

CI:

Confidence interval

CO:

Cutoff value

HC2:

Hybrid Capture® 2

HPV:

Human papillomavirus

HR:

High risk

LHU:

Local health unit

LR:

Low risk

NE:

Not evaluable

PCR:

Polymerase chain reaction

qHPV:

Quadrivalent HPV

RCTs:

Randomized controlled trials

RLU:

Relative light units

VE:

Vaccine effectiveness

References

  1. Trottier H, Franco EL. Human papillomavirus and cervical cancer: burden of illness and basis for prevention. Am J Manag Care. 2006;12:S462–72.

    PubMed  Google Scholar 

  2. Baseman JG, Koutsky LA. The epidemiology of human papillomavirus infections. J Clin Virol. 2005;32:S16–24.

    Article  PubMed  Google Scholar 

  3. Khan MJ, Castle PE, Lorincz AT, Wacholder S, Sherman M, Scott DR, et al. The elevated 10-year risk of cervical precancer and cancer in women with human papillomavirus (HPV) type 16 or 18 and the possible utility of type-specific HPV testing in clinical practice. J Natl Cancer Inst. 2005;97(14):1072–9.

    Article  PubMed  Google Scholar 

  4. de Sanjose S, Quint WG, Alemany L, Geraets DT, Klaustermeier JE, Lloveras B, et al. Retrospective International Survey and HPV Time Trends Study Group. Human papillomavirus genotype attribution in invasive cervical cancer: a retrospective cross-sectional worldwide study. Lancet Oncol. 2010;11(11):1048–56.

    Article  CAS  PubMed  Google Scholar 

  5. Serrano B, Alemany L, Tous S, Bruni L, Clifford GM, Weiss T, et al. Potential impact of a nine-valent vaccine in human papillomavirus related cervical disease. Infect Agent Cancer. 2012;7(1):38.

    Article  PubMed  PubMed Central  Google Scholar 

  6. ICO Information Centre on HPV and Cancer. www.hpvcentre.net (data updated in June 2016). Available: http://www.hpvcentre.net/statistics/reports/ITA_FS.pdf and http://www.hpvcentre.net/statistics/reports/ITA.pdf.

  7. Garland SM, Hernandez-Avila M, Wheeler CM, Perez G, Harper DM, Leodolter S, et al. Females United to Unilaterally Reduce Endo/Ectocervical Disease (FUTURE) I Investigators. Quadrivalent vaccine against human papillomavirus to prevent anogenital diseases. N Engl J Med. 2007;356(19):1928–43.

    Article  CAS  PubMed  Google Scholar 

  8. Paavonen J, Naud P, Salmerón J, Wheeler CM, Chow SN, Apter D, et al. HPV PATRICIA Study Group. 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. Lancet. 2009;374(9686):301–14.

    Article  CAS  PubMed  Google Scholar 

  9. Markowitz LE, Tsu V, Deeks SL, Cubie H, Wang SA, Vicari AS, et al. Human papillomavirus vaccine introduction–the first five years. Vaccine. 2012;30(Suppl 5):F139–48.

    Article  PubMed  Google Scholar 

  10. Romanowski B. Long term protection against cervical infection with the humanpapillomavirus: review of currently available vaccines. Hum Vaccin. 2011;7(2):161–9.

    Article  PubMed  Google Scholar 

  11. Italian National Institute of Statistics. Demography in figures. Available at: www.demo.istat.it. Accessed 11 August, 2017.

  12. Giambi C, Donati S, Carozzi F, Salmaso S, Declich S, Atti ML, et al. A cross-sectional study to estimate high-risk human papillomavirus prevalence and type distribution in Italian women aged 18–26 years. BMC Infect Dis. 2013;13:74.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Ammatuna P, Giovannelli L, Matranga D, Ciriminna S, Perino A. Prevalence of genital human papilloma virus infection and genotypes among young women in Sicily, South Italy. Cancer Epidemiol Biomarkers Prev. 2008;17(8):2002–6.

    Article  CAS  PubMed  Google Scholar 

  14. Giorgi Rossi P, Bisanzi S, Paganini I, Di Iasi A, Angeloni C, Scalisi A, et al. HPV Prevalence Italian Working Group. Prevalence of HPV high and low risk types in cervical samples from the Italian general population: a population based study. BMC Infect Dis. 2010;10:214.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Bosch FX, Burchell AN, Schiffman M, Giuliano AR, de Sanjose S, Bruni L, et al. Epidemiology and natural history of human papillomavirus infections and type-specific implications in cervical neoplasia. Vaccine. 2008;19(26 Suppl 10):K1–16.

    Article  Google Scholar 

  16. Ripabelli G, Grasso GM, Del Riccio I, Tamburro M, Sammarco ML. Prevalence and genotype identification of human papillomavirus in women undergoing voluntary cervical cancer screening in Molise, central Italy. Cancer Epidemiol. 2010;34(2):162–7.

    Article  PubMed  Google Scholar 

  17. Masia G, Mazzoleni AP, Contu G, Laconi S, Minerba L, Montixi S, et al. Epidemiology and genotype distribution of human papillomavirus (HPV) in women of Sardinia (Italy). Vaccine. 2009;27(Suppl 1):A11–6.

    Article  CAS  PubMed  Google Scholar 

  18. Agarossi A, Ferrazzi E, Parazzini F, Perno CF, Ghisoni L. Prevalence and type distribution of high-risk human papillomavirus infection in women undergoing voluntary cervical cancer screening in Italy. J Med Virol. 2009;81(3):529–35.

    Article  CAS  PubMed  Google Scholar 

  19. Tabrizi SN, Brotherton JM, Kaldor JM, Skinner SR, Cummins E, Liu B, et al. Fall in human papillomavirus prevalence following a national vaccination program. J Infect Dis. 2012;206(11):1645–51.

    Article  PubMed  Google Scholar 

  20. Markowitz LE, Liu G, Hariri S, Steinau M, Dunne EF, Unger ER. Prevalence of HPV after introduction of the vaccination program in the United States. Pediatrics. 2016;137(3):e20151968.

    Article  PubMed  Google Scholar 

  21. Haupt RM, Wheeler CM, Brown DR, Garland SM, Ferris DG, Paavonen JA, et al. FUTURE I and II Investigators. Impact of an HPV6/11/16/18 L1 virus-like particle vaccine on progression to cervical intraepithelial neoplasia in seropositive women with HPV16/18 infection. Int J Cancer. 2011;129(11):2632–42.

    Article  CAS  PubMed  Google Scholar 

  22. Muñoz N, Kjaer SK, Sigurdsson K, Iversen OE, Hernandez-Avila M, Wheeler CM, et al. Impact of human papillomavirus (HPV)-6/11/16/18 vaccine on all HPV-associated genital diseases in young women. J Natl Cancer Inst. 2010;102(5):325–39.

    Article  PubMed  Google Scholar 

  23. Castellsagué X, Muñoz N, Pitisuttithum P, Ferris D, Monsonego J, Ault K, et al. End-of-study safety, immunogenicity, and efficacy of quadrivalent HPV (types 6, 11, 16, 18) recombinant vaccine in adult women 24–45 years of age. Br J Cancer. 2011;105(1):28–37.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Wheeler CM, Kjaer SK, Sigurdsson K, Iversen OE, Hernandez-Avila M, Perez G, et al. 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 sexually active women aged 16–26 years. J Infect Dis. 2009;199(7):936–44.

    Article  PubMed  Google Scholar 

  25. Wheeler CM, Castellsagué X, Garland SM, Szarewski A, Paavonen J, Naud P, et al. HPV PATRICIA Study Group. Cross-protective efficacy of HPV-16/18 AS04-adjuvanted vaccine against cervical infection and precancer caused by non-vaccine oncogenic HPV types: 4-year end-of-study analysis of the randomised, double-blind PATRICIA trial. Lancet Oncol. 2012;13(1):100–10.

    Article  CAS  PubMed  Google Scholar 

  26. Malagón T, Drolet M, Boily MC, Franco EL, Jit M, Brisson J, et al. Cross-protective efficacy of two human papillomavirus vaccines: a systematic review and meta-analysis. Lancet Infect Dis. 2012;12(10):781–9.

    Article  PubMed  Google Scholar 

  27. Tota JE, Ramanakumar AV, Jiang M, Dillner J, Walter SD, Kaufman JS, et al. Epidemiologic approaches to evaluating the potential for human papillomavirus type-replacement postvaccination. Am J Epidemiol. 2013;178(4):625–34.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Tabrizi SN, Brotherton JM, Kaldor JM, Skinner SR, Liu B, Bateson D, et al. Assessment of herd immunity and cross-protection after a human papillomavirus vaccination program in Australia: a repeat cross-sectional study. Lancet Infect Dis. 2014t;14(10):958–66.

    Article  PubMed  Google Scholar 

  29. Drolet M, Bénard É, Boily MC, Ali H, Baandrup L, Bauer H, et al. Population-level impact and herd effects following human papillomavirus vaccination programs: a systematic review and meta-analysis. Lancet Infect Dis. 2015;15(5):565–80.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Kavanagh K, Pollock KG, Potts A, Love J, Cuschieri K, Cubie H, et al. Introduction and sustained high coverage of the HPV bivalent vaccine leads to a reduction in prevalence of HPV 16/18 and closely related HPV types. Br J Cancer. 2014;110(11):2804–11.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Mesher D, Soldan K, Lehtinen M, Beddows S, Brisson M, Brotherton JML, et al. Population-level effects of human papillomavirus vaccination programs on infections with nonvaccine genotypes. Emerg Infect Dis. 2016;22(10):1732–40.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Choi YH, Chapman R, Gay N, Jit M. Potential overestimation of HPV vaccine impact due to unmasking of non-vaccine types: quantification using a multi-type mathematical model. Vaccine. 2012;30(23):3383–8.

    Article  PubMed  Google Scholar 

  33. Eklund C, Zhou T. Dillner J; WHO Human Papillomavirus Laboratory Network. Global proficiency study of human papillomavirus genotyping. J Clin Microbiol. 2010;48(11):4147–55.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Eklund C, Forslund O, Wallin KL, Zhou T, Dillner J. WHO Human Papillomavirus Laboratory Network. The 2010 global proficiency study of human papilloma virus genotyping in vaccinology. J Clin Microbiol. 2012;50(7):2289–98.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Joura EA, Giuliano AR, Iversen OE, Bouchard C, Mao C, Mehlsen J, et al. Broad Spectrum HPV Vaccine Study. A 9-valent HPV vaccine against infection and intraepithelial neoplasia in women. N Engl J Med. 2015;372(8):711–23.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

The authors take full responsibility for the content of this contribution and thank Graham Joint, Synercom Ltd, Macclesfield, UK (supported by Sanofi-Pasteur MSD) and Meredith Rogers, The Lockwood Group, Stamford, CT, USA (supported by Merck & Co.) for editorial support.

The authors would like to thank the participants, as well as the investigators and their study-site personnel, for their contribution to the study. Also, they specifically would like to thank:

• Sophie Carles, François Simondon, and Alfred Mahr (former SPMSD employees) for their scientific contribution in study design

• Samantha Athrux-Tallau (former SPMSD employee) for her help in study implementation and initial study coordination

• Xavier Cornen (SPMSD) and Perrine Bertet (former SPMSD employee) for data management

• Stéphane Thomas (SPMSD) and Marie-Laure Kurzinger (former SPMSD employee) for their contribution to statistical analyses

• Alessandra Salati (ISPO) for her contribution to data collection, HPV sample management and testing; Gianfranco Manneschi (ISPO) for his contribution to the Study database creation, validation, and data management

• The Theorem team: V. Bass, G. Friso, I. Penne, M. Penne, C Ulbricht

• The LHU team: C. Annona, R. Basile, M.A. Bruno, L. Carbone, E. Cetani, A. De Filippis, F. D’Eufemia, A.M. Di Tursi, C. Fiorino, B. Frascella, R. Giannace, L. Gravina, C. Guerrieri, M.M. La Piccirella, L. Libardi, I. Lo Franco, R. Montemurro, T. Moramarco, C. Maremonti, M. Panetta, L. Pietromatera, C. Quintano, A.V. Scarciolla, A. Sorge, L. Taccardi, A. Venezia, M. Vignola for participants enrolment, samples management, pap-test reading, and data collection

Funding

Sanofi-Pasteur MSD, the study sponsor, provided financial support for the conduct of the research and preparation of the article, and was involved in the study design, in the collection, analysis, and interpretation of data and in the writing of the study trial report.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Author information

Authors and Affiliations

Authors

Contributions

FC, as scientific coordinator of the study, contributed the critical review and validation of all study documents (including study design, study protocol, laboratory manual and statistical analysis), laboratory analysis, data analysis and interpretation of study results, and critically reviewed all versions of the manuscript and approved the final version. DP contributed to study design, protocol development and revision, data analysis, critical interpretation of the study results, wrote the first manuscript draft, and critically reviewed all versions of the manuscript and approved the final version. CO contributed to study documents, data analysis, critical interpretation of the study results, wrote first manuscript draft, and critically reviewed all versions of the manuscript and approved the final version. EB contributed to study documents, data collection and interpretation, participated in writing the first manuscript draft, and critically reviewed manuscript drafts and approved the final version. MC, PM, and MZ contributed to study design, protocol development and revision, contributed to data analysis and interpretation, and critically reviewed all versions of the manuscript and approved the final version. PSA contributed to protocol development/revision, participants’ enrolment, data collection and data interpretation, and critically reviewed manuscript drafts and approved the final version. EAM contributed to protocol development/revision, participants’ enrolment, data collection and data interpretation, and critically reviewed manuscript drafts and approved the final version. GDF contributed to data analysis and result interpretation, wrote the manuscript outline, and critically reviewed all versions of the manuscript and approved the final version. EP contributed to study coordination, critical review of all study documents (including study design, study protocol, laboratory manual and statistical analysis), and critically reviewed the manuscript drafts and approved the final version. LS contributed to study design, protocol development and revision, critically reviewed all versions of the manuscript and approved the final version.

Corresponding author

Correspondence to Géraldine Dominiak-Felden.

Ethics declarations

Ethics approval

All the patients gave a proper informed consent to participate in the study. The study was approved by the Matera LHU ethics committee.

Consent for publication

Not applicable.

Competing interests

The authors declare the following interests:

FC, DP, CO, EB, MC, PM, and MZ report funding provided by Sanofi-Pasteur MSD, the study sponsor, to their employing institution (ISPO), to perform study-related activities.

EAM and PSA report funding provided by Sanofi-Pasteur MSD, the study sponsor, to their employing institution, to perform study-related activities.

GDF, LS, and EP were employees of Sanofi-Pasteur MSD at the time the study was conducted.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional file

Additional file 1: Table S1

A. Description of the sociodemographic characteristics of the study population (N =2793). B. Description of the behavioral characteristics of the study population (N =2793). Table S2. Age-stratified, crude, and standardized HPV prevalence (%). Table S3. Age-stratified, crude, and standardized HPV prevalence (%) by vaccine and non-vaccine types. (DOCX 23 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Carozzi, F., Puliti, D., Ocello, C. et al. Monitoring vaccine and non-vaccine HPV type prevalence in the post-vaccination era in women living in the Basilicata region, Italy. BMC Infect Dis 18, 38 (2018). https://doi.org/10.1186/s12879-018-2945-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12879-018-2945-8

Keywords