Are the Two Human Papillomavirus Vaccines Really Similar? A Systematic Review of Available Evidence: Efficacy of the Two Vaccines against HPV.

BACKGROUND: When the bivalent and the quadrivalent HPV vaccines were marketed they were presented as having comparable efficacy against cervical cancer. Differences between the vaccines are HPV types included and formulation of the adjuvant.
METHOD: A systematic review was conducted to assess the efficacy of the two vaccines against cervical cancer. Outcomes considered were CIN2+, CIN3+, and AIS.
RESULTS: Nine reports (38,419 women) were included. At enrollment mean age of women was 20 years, 90% had negative cytology, and 80% were seronegative and/or DNA negative for HPV 16 or 18 (naïve women). In the TVC-naïve, VE against CIN2+ was 58% (95% CI: 35, 72); heterogeneity was detected, VE being 65% (95% CI: 54, 74) for the bivalent and 43% (95% CI: 23, 57) for the quadrivalent. VE against CIN3+ was 78% (95% CI: <0, 97); heterogeneity was substantial, VE being 93% (95% CI: 77, 98) for the bivalent and 43% (95% CI: 12, 63) for the quadrivalent. VE in the TVC was much lower. No sufficient data were available on AIS.
CONCLUSIONS: In naïve girls bivalent vaccine shows higher efficacy, even if the number of events detected is low. In women already infected the benefit of the vaccination seems negligible.

1. Introduction

In the seventies, Dr. Harald Zur Hausen firstly postulated the link between human papillomaviruses (HPVs) and cervical cancer: studies to develop an anticancer vaccine followed. Approximately 70% of cervical cancers worldwide are associated with two high-risk HPV types (16/18) [1, 2] and almost 90% of genital warts are associated with two low-risk HPV types (6/11). Each year around 500,000 women develop invasive cervical cancer worldwide, with 83% of new cases and 85% of deaths occurring in developing countries [3, 4]. Risk factors associated with HPV infection are younger age at first coitus, higher number of sexual partners, smoking cigarettes, and history of herpes simplex virus infection [5, 6]. Some 75% of sexually active women develop a HPV infection [7], more frequently soon after their sexual debut: the majority of these infections (between 70% and 90%) spontaneously clear [8-10]. A minority progress from acute infection to cervical cancer, a process taking decades and going through precancerous lesions named cervical intraepithelial neoplasia (CIN) of increasing severity, from CIN1 to CIN3; spontaneous regression of the lesions is possible at any point [4, 8]. Incidence of genital warts is less precisely known, due to lack of data on the general population, but it is estimated to be around 1% [8].

Based on the results of five randomized controlled trials (RCTs) involving 40,000 women [11-17] two HPV vaccines entered the market. The quadrivalent vaccine against HPV6/11/16/18 was approved by the FDA and the EMA in 2006 [18, 19], whereas the bivalent vaccine against HPV16/18 was first approved in Europe in 2007 [20] and then in the USA in 2009 [21]. Soon after, several western countries such as the USA, Australia, and five European states [22] started national immunization campaigns targeting adolescent girls. The number of countries adopting the vaccines has since increased: in April 2014, 23 out of 29 European countries were reported to have implemented it [23], budgetary constraints being one relevant obstacle for the remaining countries. The primary target of HPV vaccination is adolescent girls aged 11 to 13 years, with some minor differences in national recommendations: in the USA routine vaccination is recommended at age 11 or 12 years with quadrivalent or bivalent vaccine for females and with quadrivalent vaccine for males in a 3-dose schedule during a 6-month interval [24]; in the UK HPV vaccine is recommended for girls under 15 years of age and consists in two injections spaced at least six and not more than 24 months apart. The vaccine is also recommended for men who have sex with men [25]. No prior assessment with Pap testing or screening for existing HPV infection is required.

Both vaccines contain human papillomavirus L1 self-assembling virus-like particles and are not infectious. Differences between the two vaccines are the number of HPV types included and formulation of the adjuvant (Table 1) possibly leading to different vaccine efficacy (VE) [36-38]. Head-to-head comparisons between the two vaccines are still exclusively based on immunogenicity [38, 39], although an immune correlate of protection has not yet been established [40].

Table 1

Characteristics of bivalent and quadrivalent HPV vaccine.

	Quadrivalent vaccine	Bivalent vaccine
Commercial name	Gardasil/Silgard	Cervarix
Manufacturer	Sanofi Pasteur MSD SNC	GlaxoSmithKline Biologicals S.A.
HPV types	HPV 6 L1 protein 20 μg	HPV 16 L1 protein 20 μg
	HPV 11 L1 protein 40 μg	HPV18 L1 protein 20 μg
	HPV 16 L1 protein 40 μg
	HPV 18 L1 protein 20 μg
Common characteristics	L1 protein in the form of noninfectious virus-like particles produced by recombinant DNA technology
Differences in cellular culture	Yeast cells (Saccharomyces cerevisiae CANADE 3C-5 (strain 1895))	Hi-5 Rix4446 cells derived from Trichoplusia ni using a Baculovirus expression system
Differences in adjuvant	Amorphous aluminium hydroxyphosphate sulfate adjuvant, 225 μg	AS04 adjuvant system composed of aluminium hydroxide and 3-O-desacyl-4′-monophosphoryl lipid A, 50 μg
Therapeutic indications	Gardasil is a vaccine for use from the age of 9 years for the prevention of(i) premalignant genital lesions (cervical, vulvar, and vaginal) and cervical cancer causally related to certain oncogenic human papillomavirus (HPV) types; (ii) genital warts (condyloma acuminata) causally related to specific HPV types.	Cervarix is a vaccine for use from the age of 9 years for the prevention of  (i) premalignant cervical lesions and cervical cancer causally related to certain oncogenic human papillomavirus (HPV) types.
Efficacy data leading to registration	95.2% (87.2, 98.7)1	90.4% (53.4, 99.3)
Efficacy data in the latest publication	43% (13, 63)3	93% (79, 99)4
Mean follow-up of phase III trials	3.6 years1	4 years2

1Vaccine efficacy against CIN/AIS associated with vaccine related HPV in the TVC-naïve from EMA registration data [19].

2Vaccine efficacy against CIN2+ associated with vaccine related HPV in the TVC-naïve from EMA registration data [20].

3Vaccine efficacy against any CIN3 in the TVC-naïve from Muñoz et al. 2010 [26].

4Vaccine efficacy against any CIN3+ in the TVC-naïve from Lehtinen et al. 2012 [27].

Three meta-analyses, published when vaccines were marketed, showed comparably high efficacy of the two vaccines against precancerous lesions associated with HPV16/18 [30, 41, 42]. The meta-analyses subsequently published [43, 44] have confirmed the high efficacy of the vaccines against lesions associated with HPV16/18, but they did not provide information about VE against any cervical lesions irrespective of HPV type, nor did they discuss possible differences in terms of efficacy of the two vaccines. Thus, our systematic review that includes studies with longer follow-up aims to assess differences between the two vaccines from a public health perspective, considering all cervical lesions.

2. Methods

The systematic review was developed based on a prespecified protocol (protocol number FARM8N2ZFL) funded by the Italian Medicines Agency (AIFA) within a program of independent research on drugs [45]. The PRISMA (Preferred Reporting Items for Systematic Reviews and Meta-Analyses) statement guided the content and reporting of the review [46]. Published and unpublished RCTs comparing any of the two HPV vaccines versus placebo or any other control were considered for inclusion. We exclusively considered studies involving women, irrespective of age at enrolment.

As for protocol, primary outcome measures were cervical lesions (i.e., cervical cancer, CIN2, CIN3, and AIS) associated with any HPV type and cervical lesions exclusively associated with HPV16/18 occurring in three study populations: according to protocol population (ATP), the general population of vaccinated women (total vaccine cohort (TVC)) approximating all women regardless of status of HPV infection at vaccination, and a selected population of women seronegative for HPV16/18 and HPV DNA negative for 14 oncogenic HPV types, approximating the group of young adolescent girls targeted in the national immunization campaigns (TVC-naïve). In this paper only data related to cervical lesions associated with any HPV type occurring in the TVC and in the TVC-naïve population are reported, as the other data are not relevant from a public health perspective.

2.1. Literature Search

Trial identification: we searched the Cochrane Library (to Issue 3, 2014), MEDLINE (to March 2014), and EMBASE (to March 2014) using keywords and MeSH terms as reported in Annex 1 online, in combination with a highly sensitive filter for identifying RCT [47]. There were no language or time restrictions. Reference lists of relevant papers were also examined to identify additional studies. For unpublished RCTs, we searched the Internet for prepublication study presentations at conferences or meetings. Moreover experts in the field and vaccine manufacturers were contacted for further information (unpublished studies and single patient data). Clinical trial registers were searched for ongoing studies. Two review authors independently screened abstracts of potential studies and retrieved full articles for those deemed eligible.

2.2. Study Selection

Two reviewers carried out independent assessment of citations retrieved. When more than one publication reported the same trial, the one with a longer follow-up was selected. Reasons for exclusion were recorded. Quality of trials was assessed using the criteria outlined in the Cochrane Handbook [47] and included the assessment of (i) generation of the randomization sequence, (ii) quality of the allocation concealment, (iii) completeness of follow-up, and (iv) blinding of the outcome assessment. Based on quality assessment the risk of bias of the included trials was defined as low, high, or unclear. Differences in opinion were resolved through discussion involving a third author if needed.

2.3. Statistical Analysis

Two review authors extracted the data independently using a data extraction form. Statistical analyses were carried out using the STATA software version 11. For time to event data the hazard ratio (HR) was used as a measure of association. The results were summarized by using the inverse of variance method and the random effects model. Point estimates as well as their 95% confidence intervals (95% CI) were calculated and represented by the forest plot. Vaccine efficacy (VE) was calculated as VE(%) = (RU − RV)/RU × 100, where RU is the rate of disease in the unvaccinated and RV is the rate in the vaccinated [48]. The equation can be rewritten to use the HR in the following way: VE(%) = (1 − HRv/u) × 100, where HRv/u is the HR of the vaccinated versus the unvaccinated. Heterogeneity among studies was assessed using the I 2 statistics [47].

2.4. Subgroup Analyses

The following prespecified subgroup analyses were planned:

Geographical areas (Europe, Africa, Asia, North America, and South and Central America): the literature suggests that prevalence and circulation of HPV high risk types varies according to geographical areas [49-51].

Vaccine type (bivalent, quadrivalent): data on VE can be influenced by type of vaccine used as differences between the two formulations could be not negligible.

3. Results

3.1. Selection of Studies

Study identification and selection process is outlined in Figure 1. Of the 726 records initially identified, 3 were duplicates and 670 were excluded based on title and abstract assessment. The most common reasons for exclusion were the following: reports were not RCTs, did not assess the outcome of interest, were not related to oncogenic HPV or to vaccine administration, or were studies reporting exclusively laboratory or immunogenicity data. We assessed the full text of 53 articles and excluded 44 (reasons reported in Figure 1). Thus nine reports [26–29, 31–35], corresponding to five registered protocols [52-56], were included in the systematic review: in three trials (20,797 women) the bivalent vaccine was used (PATRICIA trial being the larger study) [27-29] and in two trials (17,622 women) the quadrivalent vaccine was assessed (FUTURE I and FUTURE II trials) [26]. Five publications reporting subset of data of the above-mentioned trials relating to specific geographical areas were also included [31-35].

Figure 1

Study selection flowchart.

The data of the phase III trials included in our meta-analysis have a mean follow-up of 4 years for the bivalent trial [27] and 3.6 years for the quadrivalent trial [26], longer than that of the systematic reviews with meta-analysis published up to now.

3.2. Quality of Included Trials and Characteristics of Trial Participants

Characteristics of the trials are summarized in Table 2.

Table 2

Characteristics of the five trials selected.

Study [reference]	Protocol	Number of study sites	Countries included	Year of study enrolment	Funding source	Inclusion criteria		Exclusion criteria	Vaccine type∗	Control	Women enrolled (TVC)	Women assessed (TVC)	Length of follow-up
						Age years	Sexual partners	Any of the following
GSK [11, 12, 28]	NCT00120848  Phase II trial	27	Brazil, Canada, and USA	2001	GSK	15–25	≤6	Abnormal cervical cytology, HPV16/18 seropositivity,  DNA positivity for 14 oncogenic HPV types,  history of abnormal Pap test.	B	Aluminium containing placebo	560 V  553 P	505 V  497 P	5.9 years (average)  6.4 years (maximum)
PATRICIA [13, 27]	NCT00122681  Phase III trial	135	14 (Asia-Pacific, Europe, Latin America, and North America)	2004-2005	GSK	15–25	≤6	History of abnormal Pap test, pregnancy or breastfeeding, chronic disease, autoimmune disease, immunodeficiency.	B	HAV vaccine	9319 V 9325 P	8694 V 8708 P	4 years (average)
Konno et al. [29]	NCT00316693  Phase II trial	13	Japan	2006	GSK	20–25	Not specified	History of abnormal Pap test or genital warts, pregnancy,  previous vaccination with HPV or HAV vaccine, MPL administration, hepatitis A infection.	B	HAV vaccine	519 V 521 P	501 V 501 P	13.6 months in Konno 2010
FUTURE I [16, 26, 30]	NCT00092521  Phase III trial	62	16 (Asia-Pacific, Europe, and America)	2001–2003	Merck	16–24	≤4	History of abnormal Pap test or genital warts, pregnancy, being not healthy.	Q	Aluminium containing placebo	2723 V 2732 P	7980 V 7236 P Follow-up FUTURE I + FUTURE II trial	3.6 years (average) 4.9 years (maximum)
FUTURE II [17, 26, 30]	NCT00092534  Phase III trial	90	13 (as above)	2002-2003	Merck	15–26	≤4	History of abnormal Pap test or genital warts, pregnancy, being not healthy.	Q	Aluminium containing placebo	6087 V 6080 P	As above	As above

B: bivalent vaccine; Q: quadrivalent vaccine.

Studies were all double blind RCTs, with an adequate sample size. They were all manufacturer-sponsored. Eligible participants were healthy, not pregnant women aged between 15 and 26 years, with 6 or less lifetime sexual partners and no history of abnormal Pap smear at enrolment. Almost 90% of women had normal cytology at study entry. Risk of bias was low in 4 out of 5 trials: generation of the randomization sequence, quality of allocation concealment, completeness of follow-up, and blinding of outcome assessment were adequate. Risk of bias was unclear (generation of the randomization sequence and allocation concealment were not described) in one smaller trial conducted in Japan [29]; as the Japanese trial provided only data on CIN2+ lesions exclusively associated with HPV16/18, its results are not reported here. Characteristics of the women enrolled in the trials used in this meta-analysis are reported in Table 3.

Table 3

Baseline characteristics of women enrolled in the four studies presented in the meta-analysis.

Study [reference]	Age mean (SD)		HPV 16 positivity at enrolment (%)				HPV 18 positivity at enrolment (%)				Lifetime number of sexual partners (median)			Smoking status (%)		Chlamydia trachomatis (%)		Hormonal contraceptive use (%)		Cytological abnormality at entry (%)
			DNA		Serology		DNA		Serology
	V	C	V	C	V	C	V	C	V	C		V	C	V	C	V	C	V	C	V	C
GSK [11, 12, 28]	20 ± 3	21 ± 3	NR	NR	NR	NR	NR	NR	NR	NR		NR	NR	NR	NR	NR	NR	NR	NR	NR	NR
PATRICIA [13, 27]	20 ± 3	20 ± 3	14	14	2	1	11	10	1	1	0 ∗	4%	4%	30	30	6	5	59	61	10	9
											1 ∗	74%	74%
											2 ∗	14%	15%
											≥3 ∗	8%	8%
FUTURE I [16, 26, 30]	20 ± 2	20 ± 2	9	8	12	12	3	3	3	3		2	2	26	26	4	5	58	57	11	12
FUTURE II [17, 26, 30]	20 ± 2	20 ± 2	9	9	11	11	4	4	4	4		2	2	NR	NR	4	4	59	60	12	11

V: vaccine group; C: control group; NR: not reported.

Results are presented as percentage of women stratified by number of sexual partners.

3.3. Pooling of the Data

Pooling of the data was possible for CIN2+ [26-28], CIN3+ [26, 27], and AIS [26, 28]. Data for each outcome are presented below for the TVC and the TVC-naïve cohorts.

CIN2+. The pooled HRs for CIN2+ lesions associated with any type of HPV in the TVC and in the TVC-naïve are reported in Figure 2. Corresponding values of VE were 26% (95% CI: 11, 39) in the TVC [26, 27] and 58% (95% CI: 35, 72) in the TVC-naïve [26-28]. Results suggested substantial heterogeneity among bivalent and quadrivalent vaccines: I 2 test was 68.7% and 66.4% in the TVC and in the TVC-naïve, respectively.

Figure 2

Vaccine efficacy against CIN2+ lesions, in total vaccine (a) and total vaccine naïve cohort (b), any HPV type.

CIN3+. The pooled HRs for CIN3+ lesions associated with any type of HPV in the TVC and in the TVC-naïve are reported in Figure 3. Corresponding values of VE were 32% (95% CI: <0, 56) in the TVC [26, 27] and 78% (95% CI: <0, 97) in the TVC-naïve [26, 27]. Results suggested substantial heterogeneity among bivalent and quadrivalent vaccines: I 2 test was 86.3% and 90.7% in the TVC and in the TVC-naïve, respectively.

Figure 3

Vaccine efficacy against CIN3+ lesions, in total vaccine (a) and total vaccine naïve cohort (b), any HPV type.

AIS. The pooled HR for lesions associated with any type of HPV was 0.31 (95% CI: 0.14–0.70) in the TVC [26, 27], corresponding values of VE being 69% (95% CI: 30, 86).

AIS cases in the TVC-naïve cohort were zero in the vaccine group and ten in the placebo; thus only an approximate estimate of the efficacy was possible. The pooled HR for lesions associated with any type of HPV was 0.01 (95% CI: 0.01–0.22) resulting in a VE of 99% (95% CI: 78, 99) [26, 27]. For all the comparisons I 2 suggested low heterogeneity.

3.4. Analysis by Geographical Area

Although formally required by the scientific advisory unit of ECDC in Stockholm, the two manufacturers did not provide single patient data (the authors received only partial data unsuitable for the analysis from GSK and no answer at all from Sanofi Pasteur MSD).

Five published papers [31-35] reported data according to geographical areas. Pooling of the data was not appropriate since geographical areas definition differed (Table 4).

Table 4

Vaccine efficacy against CIN2+ in different geographical regions (Latin America, Asia-Pacific, and Europe).

Study [reference]	Vaccine type∗	Geographical region	Protocols included	Women enrolled	Outcomes assessed	Lesion n/N		VE% (95% CI)
						Vaccine	Control
Perez et al. [31]	Q	Latin America	V501-007; V501-0013; V501-0015; V501-0016; V501-0018	6400	CIN2+ HPV 6/11/16/18 in
					ATP cohort	1/2415	21/2377	95.3 (71.0, 99.9)
					TVC-naïve	3/2671	26/2681	88.5 (62.5, 97.8)
					TVC	45/2718	67/2725	33.1 (1.0, 55.2)
Tay et al. [32]	Q	Asia-Pacific	V501-0013; V501-0015; V501-0016	814	CIN2+ HPV 6/11/16/18 in
					ATP cohort	0/302	5/312	100.0 (−12.4, 100.0)
Majewski et al. [33]	Q	Europe	V501-007; V501-0013; V501-0015; V501-0016	9265	CIN2+ HPV 6/11/16/18 in
					ATP cohort	0/4043	38/4043	100.0 (89.9, 100.0)
					CIN2+ any type in
					TVC-naïve	23/2470	54/2527	56.6 (28.0, 74.6)
Barr et al. [34]	Q	North America	V501-005; V501-007; V501-0013; V501-0015; V501-0016	5996	CIN2+ HPV 16/18 in
					TVC-naïve	0/2100	35/2116	100.0 (89.0, 100.0)
					TVC	19/2313	57/2356	66.4 (42.7, 81.1)
					CIN2+ any type in
					TVC	72/2313	108/2356	33.0 (8.9, 51.0)
de Carvalho et al. [35]	B	Brazil	NCT00689741 HPV001; NCT00120848 HPV007; NCT00518336 HPV023	431	CIN2+ HPV 16/18 in
					TVC	0/219	3/212	100.0 (−129.8, 100.0)
					CIN2+ any type in
					TVC	5/219	8/212	40.6 (−106.0, 84.7)

CIN: cervical intraepithelial neoplasia; ATP: according to protocol; TVC: total vaccine cohort.

Vaccine type Q = quadrivalent; B = bivalent.

4. Discussion

Since their introduction into the market, the effectiveness of the two vaccines against cervical cancer, based on first published data [11-17], has been subject of debate. Enthusiastic positions assumed that if vaccines are immunogenic and prevent infections associated with HPV16/18 they also prevent cervical cancer and therefore should be widely used [57-59]. Uncertainty was related to the following issues: correlation between immune response and clinical outcomes, need for a booster dose, replacement with other oncogenic strains, and possible reduction of Pap-test screening among the vaccinated [60-65].

Our systematic review highlights that, for precancerous lesions (CIN), the only available proxy of cervical cancer and heterogeneity among pooled studies is substantial (Figures 2 and 3): the bivalent vaccine shows higher efficacy against precancerous lesions. We focus our comments on the TVC-naïve cohort, as VE in the TVC is confirmed to be much lower and HPV vaccination is not universally offered to women already sexually active and data on AIS are too sparse to make sensible comments. For CIN2+ lesions estimates of efficacy of the two vaccines in the TVC-naïve cohort differ but the wide limits of the confidence intervals partially overlap (VE 65%; 95% CI: 54, 74 for the bivalent and VE 43%; 95% CI: 23, 57 for the quadrivalent), whereas for CIN3+ lesions estimates of efficacy largely differ and the limits of the confidence intervals do not overlap (VE 93%; 95% CI: 77, 98 for the bivalent and VE 43%; 95% CI: 12, 63 for the quadrivalent). However, due to limited number of patients with lesions detected, leading to wide confidence intervals of our estimates of effect, our conclusions should be interpreted with caution. The heterogeneity observed might be due to higher efficacy of the bivalent vaccine against cervical cancer possibly related to the specific adjuvant used (ASO4 adjuvant system), as suggested in two recent immune response head-to-head studies that consistently showed a higher neutralizing antibody production [38, 39] and a higher CD4+ T cell response [38] in bivalent than in quadrivalent vaccine recipients. Heterogeneity can also be due to baseline differences between the populations enrolled in the two trials, although such differences were not reported in the two trials (Table 3). However, misclassification of naïve women cannot be ruled out since the two manufacturers used different laboratory tests to measure immune response and to identify naïve girls. In fact cLIA test was used in the FUTURE trials whereas ELISA test that has a higher sensitivity than cLIA [66] was used in the PATRICIA trial. Moreover, data are often differently and poorly reported in the published trials [26–29, 31–35]; thus our ability to make meaningful comparisons and further analysis (e.g., assess the possible effect modification by smoking status or age) is hindered. We asked the manufacturers to provide individual patient data, but we did not receive a positive answer.

Another relevant point is that the length of the follow-up in the trials assessed seems insufficient to detect information relevant to the public and to policy-makers: as time interval from HPV infection to cervical cancer development is approximately 20 years, all information gathered in a much shorter period of time is not conclusive. Nevertheless the PATRICIA trial has a planned follow-up of 4 years [27] and longer follow-up data on bivalent vaccine are only available for 436 Brazilian women enrolled in a previous phase II trial [67]. We will have more information in 2020, when the results of the Finnish study that extended the follow-up for Finnish girls enrolled in the PATRICIA study will be published [68] and when the extension studies for the FUTURE II trial assessing the quadrivalent vaccine will also be available. In the meantime an open debate in this respect is urgently needed: national health agencies should set up a surveillance system to provide data on actual vaccines efficacy in the field and to allow international comparisons. At the moment, this comparison is not possible and we still do not know how to choose between the two vaccines [69]. Contrary to what is previously reported in other meta-analyses [30, 41–44], our systematic review suggests that the quadrivalent and the bivalent vaccines differ in terms of efficacy. This could be attributable to the different adjuvants contained in the two vaccines. Such difference in efficacy has not been widely recognized by national health agencies. For example, in Italy HPV vaccines are chosen and purchased through tendering schemes organized by regional health authorities that are based on the lowest price [70].

Apart from cervical cancer prevention, quadrivalent vaccine is known to effectively prevent genital warts [71, 72], whereas the bivalent vaccine can only marginally impact on these benign but distressing lesions [73]. The UK and a few Italian regions have recently substituted the bivalent with the quadrivalent vaccine, assuming comparable efficacy of the two vaccines against cervical cancer and giving an additional value to the activity against genital warts of the quadrivalent vaccine [10, 74]. Unfortunately, it is not possible to anticipate the consequences of the UK's and Italian choice: if different vaccine efficacy of the bivalent and the quadrivalent vaccine is going to be confirmed, lack of equivalence in terms of cervical cancer prevention could become an issue. The availability in the coming years of new broader spectrum HPV prophylactic vaccines could provide more insight into the current debate [75, 76].

5. Conclusions

In conclusion, we acknowledge that this systematic review has some limitations due to a low number of women with events and a high heterogeneity among trials that suggest caution in the interpretation of results. However, our conclusions are consistent with those from recent immunogenicity head-to-head studies [38, 39] and that provides strength to our interpretation. Our systematic review suggests that after nine years since HPV vaccines were introduced, their estimates of efficacy seem to diverge over time. The decision to consider the two vaccines similar in terms of cervical cancer prevention seems challenged by our longer term follow-up analyses. This might have implications for policy and pragmatic choices and deserves an open and comprehensive discussion. Moreover, international comparisons of the actual effectiveness of the two vaccines used in the field can add valuable information.

Finally, regulatory agencies should encourage the pharmaceutical companies to provide data across trials to assess all relevant outcomes in a comparable way, to reduce uncertainty, and to support health policy-makers that have to choose between alternative options [77, 78].

The supplementary material shows the search strategy that was adapted for each database searched.