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

A cost-utility analysis of cervical cancer vaccination in preadolescent Canadian females



Despite the fact that approximately 70% of Canadian women undergo cervical cancer screening at least once every 3 years, approximately 1,300 women were diagnosed with cervical cancer and approximately 380 died from it in 2008. This study estimates the effectiveness and cost-effectiveness of vaccinating 12-year old Canadian females with an AS04-adjuvanted cervical cancer vaccine. The indirect effect of vaccination, via herd immunity, is also estimated.


A 12-health-state 1-year-cycle Markov model was developed to estimate lifetime HPV related events for a cohort of 12-year old females. Annual transition probabilities between health-states were derived from published literature and Canadian population statistics. The model was calibrated using Canadian cancer statistics. From a healthcare perspective, the cost-effectiveness of introducing a vaccine with efficacy against HPV-16/18 and evidence of cross-protection against other oncogenic HPV types was evaluated in a population undergoing current screening practices. The base-case analysis included 70% screening coverage, 75% vaccination coverage, $135/dose for vaccine, and 3% discount rate on future costs and health effects. Conservative herd immunity effects were taken into account by estimated HPV incidence using a mathematical model parameterized by reported age-stratified sexual mixing data. Sensitivity analyses were performed to address parameter uncertainties.


Vaccinating 12-year old females (n = 100,000) was estimated to prevent between 390-633 undiscounted cervical cancer cases (reduction of 47%-77%) and 168-275 undiscounted deaths (48%-78%) over their lifetime, depending on whether or not herd immunity and cross-protection against other oncogenic HPV types were included. Vaccination was estimated to cost $18,672-$31,687 per QALY-gained, the lower range representing inclusion of cross-protective efficacy and herd immunity. The cost per QALY-gained was most sensitive to duration of vaccine protection, discount rate, and the correlation between probability of screening and probability of vaccination.


In the context of current screening patterns, vaccination of 12-year old Canadian females with an ASO4-ajuvanted cervical cancer vaccine is estimated to significantly reduce cervical cancer and mortality, and is a cost-effective option. However, the economic attractiveness of vaccination is impacted by the vaccine's duration of protection and the discount rate used in the analysis.

Peer Review reports


Despite a large decline (>60%) in cervical cancer incidence and mortality since the implementation of screening in Canada in the 1960s, an estimated 1,300 women were diagnosed with cervical cancer and approximately 380 died in 2008 [1]. The guidelines for screening for cancer of the cervix in Canada state that all women aged 18 and over be screened, initially with two smears one year apart. If these smears are satisfactory then rescreening every three years is advised until the age of 69 [2]. The best national data currently available show 1 year participation rates do not vary greatly among provinces, ranging from 37% in British Columbia and Ontario to 44% in Nova Scotia. The 3-year participation rate in Canada is approximately 70% [3]. Despite screening programs, cervical cancer mortality has declined only marginally since 1990 [4]. Strong evidence points to specific high-risk human papillomavirus (HPV) types as the causative agent of the vast majority of cancerous and precancerous cervical lesions [5]. Specifically, HPV 16 and HPV 18 are implicated in approximately 70% of cervical cancers [3], and other high-risk HPV types including 31, 45, 33, 52 have been found to contribute an additional 14.5% [6]. Canadian data show that 85% of all HPV infection costs are related to screening and follow-up procedures, amounting to $250-$300 million annually [7].

In February 2007, the National Advisory Committee on Immunization (NACI) recommended vaccination of women between 9 and 13 years of age[4]. According to the Canadian Immunization Committee(CIC)-NACI HPV Working Group, the goal of an HPV immunization program in Canada is to decrease the morbidity and mortality of cervical cancer, its precursors and other HPV-related cancers in women through a combination of primary and secondary prevention strategies [4]. Two HPV vaccines exist at present, a quadrivalent HPV 6/11/16/18 (Gardasil®, Merck Frosst) currently available in Canada, and an HPV 16/18 AS04-adjuvanted vaccine (Cervarix™, GlaxoSmithKline) under regulatory review by Health Canada. Both vaccines are highly efficacious and safe [813]. Immunity against HPV is thought to be primarily type-specific [1418], however there is some cross-protection between genetically related types. The HPV-16/18 AS04-adjuvanted vaccine is reported to provide high efficacy against CIN2+ associated with HPV-16/18 as well as non-vaccine oncogenic HPV individual types including HPV 45, HPV 31, HPV 33 [8, 9]. Broad protection against CIN2+ infection with combined oncogenic, non-vaccine HPV types, suggesting that cross-protection might extend beyond types 45 and 31 [8, 9]. In addition, substantial overall effect is reported in cohorts that are relevant to universal mass vaccination and catch-up programmes[8].

By reducing the amount of HPV transmission in the population, vaccination also reduces the risk of infection to unvaccinated individuals. This effect is termed herd immunity or herd protection[19, 20], and its inclusion in the assessment of HPV vaccination generally results in lower cost-effectiveness ratios under certain common conditions [2125]. Accounting for herd protection effects in a cost-utility analysis normally requires use of a dynamic model and thus abandonment of the spreadsheet-based cohort modeling framework often used for vaccination cost-effectiveness analysis[26].

Several studies have examined the cost-effectiveness of HPV vaccination in different countries. Though studies vary in model type, input data, and model assumptions, most of them have shown that HPV vaccination is likely to be cost-effective in pre-adolescent females[7, 2225, 2730, 30, 31, 3136]; however, the inclusion of male vaccination was generally not cost-effective[2225, 35] or in some cases less cost-effective and including catch-up programs was also less cost-effective[23, 25]. Sensitivity analyses demonstrated that cost-effectiveness estimates of HPV vaccine are largely influenced by factors such as the discount rate, vaccine cost, vaccine efficacy, coverage, and duration of vaccine protection[7, 22, 25, 2731, 33, 34, 36].

In this study, we estimate the effectiveness and cost-effectiveness of an AS04-adjuvanted cervical cancer vaccine in a hypothetical cohort of 12-year old females, in the context of current cervical cancer screening programs. We use a recently-developed methodology based on equations derived from dynamic models to conservatively account for herd immunity effects in the context of an existing spreadsheet-based Markov cohort model for HPV vaccination in Canadian women.


Model structure

A simple and transparent Markov model was developed in Microsoft Excel that allows for rapid and flexible use with a primary focus on cervical cancer. This model was validated against a more complex model and adapted for its application to the Canadian setting[37, 38]. Figure 1 displays the model schematics.

Figure 1
figure 1

Cohort model structure. 12-disease state Markov process model with 12 month cycle, simulates natural history of HPV infection and cervical cancer. Transition rates vary by age, HPV type, screening, treatment, vaccination. Simple screening module. Detected: subjects with disease detected through screening have same pathways but different probabilities. CIN: Cervical Intraepithelial Neoplasia.

The Markov model replicates the natural history of cervical disease over the lifetime of the cohort in one-year cycles and in the following health states: Normal, HPV infection (high-risk, low-risk), Cervical Intraepithelial Neoplasia (CIN), Cervical Cancer, Cervical Cancer Death, and Other Death. Transition probabilities were age and HPV-type (high-risk, low-risk) specific and derived through the literature. Women in the model can progress and regress from states up to persistent CIN2,3 and thereafter may progress to cervical cancer. Women in the cancer state can either remain in cancer or move to a 'cured' cancer state and remain there over their lifetime. A simplified screening module is incorporated, with no follow-up testing after a negative cytology result, and colposcopy and biopsy with treatment after a positive cytology result. Regular (i.e., 70% women are screened every 3 years), irregular (i.e., 18% various screening patterns), and never screened (i.e. 12%) are taken into account[39]. Sensitivity for screening is incorporated to estimate the proportion of patients with their disease detected and treated. Natural immunity after HPV infection was not accounted for in the model. The number of sexual partners, age at first pregnancy, race and socioeconomic status are surrogate markers for infection with HPV[40] but were not included in the model.


The model was calibrated by manually manipulating input parameters to match the model output to observed epidemiological data while keeping transition probabilities within pre-determined ranges of plausibility. Age-specific cervical cancer incidence and mortality rates[41], as well as the overall number of cases and deaths estimated by the model were matched to observed data during this process. To validate the matching between the observed epidemiological data and model-predicted data, we conducted a chi-square test for significance at the level of p < 0.05.

Vaccine efficacy

Vaccine efficacy is modelled by assuming it causes a reduction in the probability of acquiring HPV infection. The model uses an overall average vaccine effectiveness against all oncogenic types, computed as the average efficacy against HPV 16/18, and other non-vaccine oncogenic HPV types weighted by prevalence of cervical cancer attributed to each HPV-type[6]. Based on recently published data, we used 98% efficacy against HPV 16 and 18 infections and 37% for other non-vaccine oncogenic types[8, 9]. Long-term data on vaccine efficacy is not currently available. Recently data shows strong and sustained antibody response for HPV 16 and HPV 18 up to 7.3 years post-vaccination [10, 42] and mathematical models have predicted the long term vaccine duration to be at least 20 years[43]. Lifelong vaccine protection was assumed for the base-case analysis, and shorter duration times were tested in the sensitivity analysis. To prevent over-estimation of the vaccine benefit on the incidence of CIN events, CIN1 and CIN2,3 were thereafter added to take into account the difference between the proportion of HPV in CIN and cancer.

Herd immunity

Herd protection effects were included within the cohort model through a pseudo-dynamic approximation based on the adaptation of a standard compartmental dynamic model, the Susceptible-Infectious-Recovered (SIR) model[44]. The model was calibrated using published estimates of HPV transmissibility[45, 46]. The model equations were solved to determine the risk of infection to unvaccinated individuals in a population with a given level of vaccine coverage. This was combined with age-stratified data on sexual mixing patterns to determine the age-specific change in infection risk over time for those members of an aging cohort of women not successfully vaccinated under the vaccination program[24, 47]. This change in infection risk was represented by a multiplicative factor rherd, representing the proportionate reduction in infection risk due to herd immunity effects. For instance, r herd = 0.5 means that susceptible individuals have their infection risk reduced by 50% due to vaccination in other members of the population). r herd can be set to 1 when herd immunity is not taken into account through this pseudo-dynamic approximation, i.e., when infection risk is assumed constant as under conventional Markov cohort models of vaccination. Details of the methodology appear in [47].

Costs and Quality Adjusted Life Years (QALYs)

Total costs, total QALYs, number of cervical cancer cases, and number cervical cancer deaths were derived over the lifetime of one birth cohort of girls after assigning specific costs and QALY to each health state and cumulating those over the lifetime. Utilities for pre-cancer and cervical cancer were taken from the literature. The disutility of pre-cancer and cancer health states were subtracted from age-specific utility in the general population.

Data on costs were obtained from Canadian literature, though some of the data cited in the publications were derived from the US [7, 48]. When necessary, costs were inflated to 2006 using the health and personal care component of the Consumer Price Index[49]. The cost of the vaccine was assumed to be $135 per dose. Cost of vaccine administration was not included in the base case analysis, though this may be represented in the upper bound of the vaccine cost in the sensitivity analysis. Cost for adverse events was not included, though HPV vaccine has been shown to have an excellent safety profile. Future costs and health effects were discounted at 3% and varied in sensitivity analyses.

Table 1 provides a summary of the input parameters and data sources. Base case estimates and ranges were determined using a combination of literature or best assumption. Data was Canadian-specific whenever possible.

Table 1 Summary of input parameters into health economic model

Base case

Vaccination is evaluated in a hypothetical cohort of 100,000 12-year old Canadian girls. This age was chosen as the target group for vaccination because it represents women prior being at risk for HPV infection (i.e. before sexual debut). HPV vaccination rates differ across Canada by province. According to a Canadian report, HPV vaccination rates among school-based cohorts was 80% in Atlantic provinces, 85% in Newfoundland, 85-87% in Quebec, and 53% in Ontario [50]. Based on this variability we assumed a base-case rate of 75%. It was assumed that all adolescents would receive three doses of the vaccine and be fully immunized after 1 year. Vaccination is compared with the Canadian status quo for cervical cancer prevention, cytological (Pap) screening. For this analysis, we assume that current screening patterns will not change with implementation of vaccination. Four scenarios are included in base case analyses based on the inclusion or exclusion of herd protection and cross-protective efficacy. Sensitivity analysis was conducted on all four scenarios.

We conducted a cost-utility analysis from the perspective of the healthcare system (i.e., only direct costs to the healthcare system were taken into account), evaluating the incremental cost per QALY gained.

The incremental cost-utility ratio is calculated as:


where program A is current cervical cancer screening, and program B is cervical cancer vaccination plus current screening. As a threshold to decide whether the estimated incremental cost-effectiveness ratio (ICER) could be deemed cost-effective, we use the gross domestic product (GDP) per capita for Canada, $37,000 (2006 est.)[51].

Sensitivity analysis

To deal with uncertainty in parameter values, sensitivity analyses were conducted to determine the impact of varying parameter values on the cost per QALY-gained. Several input parameters were varied within a pre-specified range of values based on published literature or best/worst case assumptions These included vaccine parameters (efficacy, duration of protection, coverage rate, age at vaccination), screening parameters (frequency), economic parameters (vaccine cost, screening costs, disease costs, discount rate), and quality of life parameters (utilities). We also tested a non-zero correlation between probability of screening and probability of vaccination, based on a Canadian survey of 500 adult women that found that contrary to their initial hypothesis, a higher proportion of females who had less regular screening tests were willing to recommend an HPV vaccine to their daughters/nieces [52]. For vaccine waning scenarios, we assume that vaccine efficacy does not wane until a specific time period after initial vaccination (e.g. 10 years, 30 years) after which it begins to wane linearly over a period of 5 years. Upper bound screening were varied by 50% from base case, as a recent review suggests the base case costs we used are conservative [48].



The model predicted and observed rates of cervical cancer incidence and deaths are compared in Figure 2. The modeled number of cumulative cervical cases over the lifetime of a 12 year old cohort is well matched to observed Canadian data (822 and 824, respectively). The model reflects, with a slight overestimation, the average Canadian cervical cancer death rates with a predicted cumulative 350 deaths compared with 319 observed over the lifetime of the cohort. The total number of observed vs. modeled cervical cancer cases and deaths were tested for differences using a Chi-square test for significance and evaluated at the p = 0.05 level. No statistically significant differences were found. Predictions are well matched to Canadian estimates when looking across five-year age groups as shown in figure 2.

Figure 2
figure 2

Model Calibration: Observed vs. model predicted cases and deaths. Note: The total number of observed vs. modeled cervical cancer cases and deaths were tested for differences using a Chisquare test for significance and evaluated at the p = 0.05 level. No statistically significant differences were found.

In addition, the model predicts a cumulative 105,960 cases of CIN1 and 26,100 CIN2,3 compared with the 105,922 CIN1 and 25,294 CIN2,3 expected using the same epidemiologic data [41].

Effectiveness of Vaccination

With a 75% vaccination rate, undiscounted base case scenario results estimate 390-457 cervical cancer cases and 168-196 cervical-related deaths would be prevented over the life-time a cohort of 100,000 girls vaccinated at 12 years of age compared with an unvaccinated cohort (Table 2). Discounted cases prevented ranged from 118-138 and 38-44 for cervical cancer and deaths, respectively. This represents at least a 47% reduction in cases and deaths. Taking into account herd protection effects, the estimated reduction in cervical cancer cases and deaths increases up to 77%-78%

Table 2 Undiscounted lifetime risk-reduction in burden of disease for bases case analyses, by duration of vaccine protection

Under a scenario where protection from the vaccine wanes after 10 years, the estimated reduction in undiscounted lifetime cancer cases and deaths is at least 9% and 4%, respectively (Table 2). With a vaccine with 30 years duration of protection the estimated reduction in undiscounted lifetime cancer cases and deaths is at least 22% and 13%, respectively. When herd protection effects are included, the estimated reductions in cervical cancer cases and deaths increase quite significantly to at least of 21% and 34% for 10 and 30 year duration of vaccine protection, respectively.

Cost-utility of vaccination

Table 3 shows predicted cost-utility results in a vaccinated cohort under the base case scenarios compared with an unvaccinated cohort undergoing the same screening and treatment scenario. Compared with no vaccination, vaccinating 100,000 girls aged 12 years resulted in incremental QALYs gained of 0.0083-0.0095 per person and an incremental cost of $256 - $263 over the lifetime of the cohort. The resulting discounted cost per QALY gained for vaccination was in the range of $26,947 - $31,687 depending on whether or not cross-protection was included. The inclusion of herd immunity reduced the cost per QALY gained by 12% - 31%.

Table 3 Discounted cost, QALYs, and cost per QALY results for base case analyses, 3% discount, 75% vaccination coverage.

Sensitivity Analysis

Figure 3 shows results of the sensitivity analysis for vaccination with HPV 16/18 with cross-protective efficacy, with and without herd immunity. Results for HPV 16/18 only were similar, and therefore not shown here. HPV vaccination remained cost-effective for the most part (the cost per QALY gained remained below $37,000). The cost per QALY gained was less sensitive to changes in parameter values when herd immunity was included in the model.

Figure 3
figure 3

Sensitivity analysis results for vaccination with HPV 16/18 + cross-protection, with and without herd immunity.

The cost per QALY gained was generally not sensitive to changes in disease and screening costs, coverage, and age of vaccinated cohort. The cost per QALY gained was somewhat sensitive to changes in vaccine efficacy, and vaccine cost, though values still remained under the cost-effectiveness threshold.

Overall, the parameters that had the most impact on the cost per QALY-gained included the discount rate, duration of vaccine protection, and the correlation between probability of screening and probability of vaccination. Varying the discounting rates from 0% to 5% showed the greatest impact on the cost-utility ratio. The cost per QALY-gained varied from a low of $2,000 - $4,000 when a 0% discount rate was applied, to $50,000 - $69,000 when the discount rate was set to 5%. Given the lack of data on long-term duration of vaccine protection, we tested various waning immunity and booster scenarios. A scenario of linear waning of vaccine-derived immunity starting after 10 years of constant protection was the least effective strategy and resulted in cost per QALY-gained $68,000 and $100,000 with and without a booster, respectively. When herd protection effects are included, the cost per QALY-gained remains under the cost-effectiveness threshold. Results show that when protection is assumed to last 30 years, the cost per QALY does not change significantly when a booster is added, particularly when herd protection effects are included.

Our results indicate that if screening frequency and vaccination rates have a positive correlation (i.e. women who are infrequently or never screened are less likely to get vaccinated and vice versa), the cost per QALY-gained increases to $42,000. When herd immunity is taking into account, the cost per QALY-gained remains under the cost-effective threshold. If a negative correlation is assumed (i.e. women who are infrequently or never screened are more likely to get vaccinated and vice versa), the cost per QALY-gained decreases from our base case analysis, but does not go under $15,000 with or without inclusion of herd immunity.

Though no single parameter tested has resulted in a cost-savings situation, a multivariate sensitivity analysis that included varying only those parameters that had a positive impact on the cost per QALY-gained (i.e. lifelong vaccine protection, 0% discounting, lower vaccine cost), resulted in a very low cost per QALY gained of $1,000 - $2,000. Varying only the parameters that had a negative impact on the cost per QALY-gained (i.e. 10 year vaccine protection, 5% discounting, low vaccination rates among under-screened or never screened women) resulted in a cost per QALY-gained of $84,000 when herd protection was included to >$100,000 when herd immunity was not included.


The results from this cost-utility study suggest that within the context of current cervical cancer screening programs in Canada, cervical cancer vaccination of pre-adolescent/adolescent females is likely to be cost-effective from a healthcare perspective. These results were robust at the threshold of $37,000 per QALY-gained and under a wide range of assumptions and parameter uncertainties. Including herd immunity effects was estimated to increase base case vaccine effectiveness by approximately 8% - 23% and cost-effectiveness by 12% - 31%. Increases in vaccine effectiveness and cost-effectiveness are consistent with other economic analyses of HPV vaccination that have accounted for herd immunity effects [2225, 29]. Our base case cost-utility ratios are consistent with other Canadian analyses evaluating quadrivalent HPV female vaccination in a cohort model [7] and HPV 16/18 female vaccination in a dynamic model [53], showing cost per QALY-gained ranging from $21,000 - $33,115 per QALY-gained.

Our results indicate increased effectiveness and cost-effectiveness with a bivalent HPV 16/18 vaccine that confers cross-protection against other non-vaccine oncogenic HPV types (e.g. HPV 31, 33, 45) [8, 9], though longer term clinical trial data is needed to understand vaccine efficacy for these individual HPV types. The inclusion of cross-protective efficacy led to greater disease and mortality reductions. Cross-protective efficacy alone without accounting for herd immunity resulted in an additional 67 cancer cases prevented and 28 deaths prevented.

Consistent with other cost-effectiveness studies of HPV vaccination [7, 25, 29, 33, 34, 53], our results were sensitive to the discount rate and duration of vaccine protection. Because of the long latency between HPV vaccination and cervical cancer averted due to vaccination, we found that varying the discounting rates showed the greatest impact on the cost-utility ratio. According to much of the health economic literature, health effects should be discounted at the same rate as costs, however; some argue that health effects should be discounted at a lower rate than costs [54]. Others have suggested that discounting for vaccines should start from the time of risk reduction rather than the time of intervention [55, 56].

Our results also suggested that duration of vaccine protection is important in the effectiveness and cost-effectiveness of vaccination. A vaccine where protection starts to wane before 30 years would result in less than a 26% reduction in cumulative cancer cases and 15% in deaths. A booster shot for older females may be warranted in this case, however vaccination compliance is often lower for adult programs than for school-based programs. For instance, compliance for adult Hepatitis A vaccination programs in Canada has been estimated to be 25% [26]. Under herd immunity, our results suggest that individuals who lose their vaccine-derived immunity later in life may still derive some indirect protection through herd immunity effects. However given the uncertainty and sensitivity of this parameter, surveillance of populations in coming years to determine the long-term duration of protection from HPV vaccines will be needed.

The cost-utility ratios are similar for various possible ages of vaccination between the ages 12-25. However, one should expect that vaccination at a younger age should generally result in lower cost-utility ratios, since more infection is thereby averted. In our results, cost-utility ratios may be similar across this age range because (a) HPV prevalence in Canada peaks at age 20-25, therefore vaccination before age 20 will prevent only a small number of additional events relative to vaccination at age 12, and (b) future cancer cases averted are discounted more heavily for vaccination at a younger age, which reduces estimated QALYs gains from a given expenditure on vaccination. This is confirmed by analyzing undiscounted results, in which the most favorable cost-effectiveness ratio is shown in the younger age cohort.

Our results also indicated that cost-effectiveness improves when women who are under-screened or never screened have higher vaccination rates than those women who are more frequent screeners. This occurs because the women who are vaccinated are those most at-risk for cervical cancer (i.e. under-screened or never screened), but fewer women in the total cohort are being vaccinated and therefore the total vaccine costs for cohort are reduced. However, the most effective strategy in terms of disease reduction occurs when a high coverage rate is to the entire cohort, irrespective of screening frequency, and this strategy remains below the cost-effectiveness threshold.

There are several limitations that should be noted. The model slightly overestimated the average Canadian cervical cancer death rate. It is possible that this could bias the study results in favour of the vaccine program. However, as noted earlier, our results are consistent with other economic analyses of cervical cancer vaccination in Canada and elsewhere. Varying levels of HPV infection in a population will have an impact on effectiveness and cost-effectiveness of vaccination programs and it is important to note that this was not taken into account here. Also, our base case value of vaccine cost was based on an assumption and did not include vaccine administration costs. Vaccine administration costs will vary depending on the vaccination strategy (e.g., new school-based program versus adding HPV vaccine to a pre-existing school-based program), therefore, it is important to look at the range of the cost-utility ratios in the sensitivity analysis under vaccination costs. Even in our upper range estimate of vaccine cost, the cost per QALY-gained was within the range of what is considered cost-effective.

We used a simple model to estimate the effectiveness and cost-effectiveness of cervical cancer vaccination in 12-year old females in Canada. This spreadsheet model requires fewer data and programming calculations which allows for rapid and flexible use. The value of such a simple model is that data requirements are more modest, which may suit its application to policy development in jurisdictions where data availability and/or quantity are limited. At the same time, this model has been validated against a more complex model and predicts similar cost-effectiveness [37]. We conservatively incorporated herd immunity effects in this cohort model of cervical cancer vaccination. The fullest and most accurate way to account for herd immunity effects is through a dynamic model, however, the approach described here allows us to approximate herd immunity on effects in the context of an existing spreadsheet-based cohort model under certain conditions. This may be useful as a means of preliminary incorporation of herd immunity effects preceding development of a full dynamic model.

This type of model can help inform policy until more complex models (i.e. sexual network transmission models) are available to answer and understand more in-depth policy questions. Some of the policy questions that will need to be explored are the long-term effectiveness and cost-effectiveness of current and possible modified cervical cancer screening policies (e.g. longer screening intervals and new technologies such as liquid-based cytology, HPV DNA testing) in the context of HPV vaccination in Canada.


Using a recently-developed methodology to conservatively account for herd immunity effects in the context of an existing cohort model, our results suggest that immunization of pre-adolescent Canadian females with an ASO4-adjuvanted cervical cancer vaccine in 12 year old females is both effective at reducing cervical disease and deaths and is cost-effective when compared with cervical screening alone. However, the cost-effectiveness is highly influenced by the duration of protection of the vaccine and the discount rate used in the analysis.



cervical cancer


cervical intraepithelial neoplasia


gross domestic product


human papillomavirus


incremental cost-effectiveness ratio


quality adjusted life year




  1. Canadian cancer statistics 2008. []

  2. Health Canada: Programmatic Guidelines for Screening for Cancer of the Cervix in Canada. 1998

    Google Scholar 

  3. Cervical Cancer Screening in Canada: 1998 Surveillance Report. []

  4. Dobson S, Deeks SMD: An Advisory Committee Statement (ACS), National Advisory Committee on Immuzation (NACI): Statement on human papillomavirus vaccine. Canada Communicable Disease Report - Relevé des maladies transmissibles au Canada (CCDR-RMTC). 2007, 33: 1-32. []

    Google Scholar 

  5. Bosch FX, Manos MM, Munoz N, Sherman M, Jansen AM, Peto J, Schiffman MH, Moreno V, Kurman R, Shah KV: Prevalence of human papillomavirus in cervical cancer: a worldwide perspective. International biological study on cervical cancer (IBSCC) Study Group. J Natl Cancer Inst. 1995, 87: 796-802. 10.1002/ijc.20244.

    Article  CAS  PubMed  Google Scholar 

  6. Munoz N, Bosch FX, Castellsague X, Diaz M, de SS, Hammouda D, Shah KV, Meijer CJ: Against which human papillomavirus types shall we vaccinate and screen? The international perspective. Int J Cancer. 2004, 111: 278-285. 10.1016/j.vaccine.2007.04.086.

    Article  CAS  PubMed  Google Scholar 

  7. Brisson M, Van dV, De WP, Boily MC: The potential cost-effectiveness of prophylactic human papillomavirus vaccines in Canada. Vaccine. 2007, 25: 5399-5408. 10.1016/S0140-6736(09)61248-4.

    Article  PubMed  Google Scholar 

  8. Paavonen J, Naud P, Salmeron J, Wheeler C, Chow SN, Apter DL, Kitchener HC: 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: 301-314. 10.1016/S0140-6736(07)60946-5.

    Article  CAS  PubMed  Google Scholar 

  9. Paavonen J, Jenkins D, Bosch FX, Naud P, Salmeron J, Wheeler CM, Chow SN, Apter DL, Kitchener HC, Castellsague X, et al: Efficacy of a prophylactic adjuvanted bivalent L1 virus-like-particle vaccine against infection with human papillomavirus types 16 and 18 in young women: an interim analysis of a phase III double-blind, randomised controlled trial. Lancet. 2007, 369: 2161-2170. 10.1016/S0140-6736(06)68439-0.

    Article  CAS  PubMed  Google Scholar 

  10. Harper DM, Franco EL, Wheeler CM, Moscicki AB, Romanowski B, Roteli-Martins CM, Jenkins D, Schuind A, Costa Clemens SA, Dubin G: Sustained efficacy up to 4.5 years of a bivalent L1 virus-like particle vaccine against human papillomavirus types 16 and 18: follow-up from a randomised control trial. Lancet. 2006, 367: 1247-1255. 10.1016/S0140-6736(04)17398-4.

    Article  CAS  PubMed  Google Scholar 

  11. Harper DM, Franco EL, Wheeler C, Ferris DG, Jenkins D, Schuind A, Zahaf T, Innis B, Naud P, De Carvalho NS, et al: Efficacy of a bivalent L1 virus-like particle vaccine in prevention of infection with human papillomavirus types 16 and 18 in young women: a randomised controlled trial. Lancet. 2004, 364: 1757-1765. 10.1016/S1470-2045(05)70101-7.

    Article  CAS  PubMed  Google Scholar 

  12. Villa LL, Costa RL, Petta CA, Andrade RP, Ault KA, Giuliano AR, Wheeler CM, Koutsky LA, Malm C, Lehtinen M, et al: Prophylactic quadrivalent human papillomavirus (types 6, 11, 16, and 18) L1 virus-like particle vaccine in young women: a randomised double-blind placebo-controlled multicentre phase II efficacy trial. Lancet Oncol. 2005, 6: 271-278. 10.1056/NEJMoa020586.

    Article  PubMed  Google Scholar 

  13. Koutsky LA, Ault KA, Wheeler CM, Brown DR, Barr E, Alvarez FB, Chiacchierini LM, Jansen KU: A controlled trial of a human papillomavirus type 16 vaccine. N Engl J Med. 2002, 347: 1645-1651. 10.1056/NEJMoa020586.

    Article  CAS  PubMed  Google Scholar 

  14. Carter JJ, Koutsky LA, Wipf GC, Christensen ND, Lee SK, Kuypers J, Kiviat N, Galloway DA: The natural history of human papillomavirus type 16 capsid antibodies among a cohort of university women. J Infect Dis. 1996, 174: 927-936. 10.1016/j.vaccine.2005.09.002.

    Article  CAS  PubMed  Google Scholar 

  15. Stanley M: Immune responses to human papillomavirus. Vaccine. 2006, 24: S16-S22. 10.1016/j.vaccine.2005.09.002.

    Article  PubMed  Google Scholar 

  16. Unckell F, Streeck RE, Sapp M: Generation and neutralization of pseudovirions of human papillomavirus type 33. J Virol. 1997, 71: 2934-2939.

    CAS  PubMed  PubMed Central  Google Scholar 

  17. Wang Z, Christensen N, Schiller JT, Dillner J: A monoclonal antibody against intact human papillomavirus type 16 capsids blocks the serological reactivity of most human sera. J Gen Virol. 1997, 78: 2209-2215. 10.1006/viro.1994.1649.

    Article  CAS  PubMed  Google Scholar 

  18. Christensen ND, Kirnbauer R, Schiller JT, Ghim SJ, Schlegel R, Jenson AB, Kreider JW: Human papillomavirus types 6 and 11 have antigenically distinct strongly immunogenic conformationally dependent neutralizing epitopes. Virology. 1994, 205: 329-335. 10.1177/0272989X02239651.

    Article  CAS  PubMed  Google Scholar 

  19. Brisson M, Edmunds WJ: Economic evaluation of vaccination programs: the impact of herd-immunity. Med Decis Making. 2003, 23: 76-82. 10.1086/425271.

    Article  CAS  PubMed  Google Scholar 

  20. Garnett GP: Role of herd immunity in determining the effect of vaccines against sexually transmitted disease. J Infect Dis. 2005, 191: S97-106. 10.3201/eid1301.060438.

    Article  PubMed  Google Scholar 

  21. The cost-effectiveness of HPV vaccination in the United States. []

  22. Elbasha EH, Dasbach EJ, Insinga RP: Model for assessing human papillomavirus vaccination strategies. Emerg Infect Dis. 2007, 13: 28-41. 10.1016/j.vaccine.2007.10.056.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Insinga RP, Dasbach EJ, Elbasha EH, Puig A, Reynales-Shigematsu LM: Cost-effectiveness of quadrivalent human papillomavirus (HPV) vaccination in Mexico: a transmission dynamic model-based evaluation. Vaccine. 2007, 26: 128-139. 10.1016/j.vaccine.2007.10.056.

    Article  PubMed  Google Scholar 

  24. Taira AV, Neukermans CP, Sanders GD: Evaluating human papillomavirus vaccination programs. Emerg Infect Dis. 2004, 10: 1915-1923. 10.1071/SH07043.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Kulasingam S, Connelly L, Conway E, Hocking JS, Myers E, Regan DG, Roder D, Ross J, Wain G: A cost-effectiveness analysis of adding a human papillomavirus vaccine to the Australian National Cervical Cancer Screening Program. Sex Health. 2007, 4: 165-175. 10.1016/j.vaccine.2007.10.001.

    Article  PubMed  Google Scholar 

  26. Bauch CT, Anonychuk AM, Pham BZ, Gilca V, Duval B, Krahn MD: Cost-utility of universal hepatitis A vaccination in Canada. Vaccine. 2007, 25: 8536-8548. 10.1017/S0266462307080026.

    Article  CAS  PubMed  Google Scholar 

  27. Bergeron C, Largeron N, McAllister R, Mathevet P, Remy V: Cost-effectiveness analysis of the introduction of a quadrivalent human papillomavirus vaccine in France. Int J Technol Assess Health Care. 2008, 24: 10-19. 10.1016/j.vaccine.2007.05.061.

    Article  PubMed  Google Scholar 

  28. Boot HJ, Wallenburg I, de Melker HE, Mangen MJ, Gerritsen AA, Maas van der NA, Berkhof J, Meijer CJ, Kimman TG: Assessing the introduction of universal human papillomavirus vaccination for preadolescent girls in The Netherlands. Vaccine. 2007, 25: 6245-6256. 10.3201/eid1402.070499.

    Article  PubMed  Google Scholar 

  29. Chesson HW, Ekwueme DU, Saraiya M, Markowitz LE: Cost-effectiveness of human papillomavirus vaccination in the United States. Emerg Infect Dis. 2008, 14: 244-251. 10.1016/j.vaccine.2006.05.116.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Garnett GP, Kim JJ, French K, Goldie SJ: Chapter 21: Modelling the impact of HPV vaccines on cervical cancer and screening programmes. Vaccine. 2006, 24 Suppl 3: S3/178-S3/186. 10.1016/j.vaccine.2007.07.018.

    Google Scholar 

  31. Ginsberg GM, Fisher M, Ben-Shahar I, Bornstein J: Cost-utility analysis of vaccination against HPV in Israel. Vaccine. 2007, 25: 6677-6691. 10.1016/j.vaccine.2007.07.018.

    Article  PubMed  Google Scholar 

  32. Goldie SJ, Kohli M, Grima D, Weinstein MC, Wright TC, Bosch FX, Franco E: Projected clinical benefits and cost-effectiveness of a human papillomavirus 16/18 vaccine. J Natl Cancer Inst. 2004, 96: 604-615. 10.1016/j.vaccine.2007.05.058.

    Article  PubMed  Google Scholar 

  33. Goldie SJ, Kim JJ, Kobus K, Goldhaber-Fiebert JD, Salomon J, O'shea MK, Xavier BF, de SS, Franco EL: Cost-effectiveness of HPV 16, 18 vaccination in Brazil. Vaccine. 2007, 25: 6257-6270. 10.1016/j.vaccine.2007.05.058.

    Article  PubMed  Google Scholar 

  34. Sanders GD, Taira AV: Cost-effectiveness of a potential vaccine for human papillomavirus. Emerg Infect Dis. 2003, 9: 37-48. 10.1038/sj.bjc.6604023.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Kim JJ, ndres-Beck B, Goldie SJ: The value of including boys in an HPV vaccination programme: a cost-effectiveness analysis in a low-resource setting. Br J Cancer. 2007, 97: 1322-1328. 10.1186/1478-7547-6-4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Kulasingam SL, Benard S, Barnabas RV, Largeron N, Myers ER: Adding a quadrivalent human papillomavirus vaccine to the UK cervical cancer screening programme: A cost-effectiveness analysis. Cost Eff Resour Alloc. 2008, 6: 4-10.1016/j.vaccine.2008.02.040.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Debicki D, Ferko N, Demarteau N, Gallivan S, Bauch C, Anonychuk A, Mantovani L, Capri S, Chou CY, Standaert B, et al: Comparison of detailed and succinct cohort modelling approaches in a multi-regional evaluation of cervical cancer vaccination. Vaccine. 2008, 26: F16-F28. 10.1038/sj.bjc.6603501.

    Article  PubMed  Google Scholar 

  38. Kohli M, Ferko N, Martin A, Franco EL, Jenkins D, Gallivan S, Sherlaw-Johnson C, Drummond M: Estimating the long-term impact of a prophylactic human papillomavirus 16/18 vaccine on the burden of cervical cancer in the UK. Br J Cancer. 2007, 96: 143-150. 10.1038/sj.bjc.6603501.

    Article  CAS  PubMed  Google Scholar 

  39. Statistics Canada: Pap smear, by age group, women aged 18 to 69 years, Canada, provinces, territories, health regions (January 2000 boundaries) and peer groups, every 2 years. Canadian Community Health Survey: . 2000

    Google Scholar 

  40. Mitchell MF, Tortolero-Luna G, Wright T, Sarkar A, Richards-Kortum R, Hong WK, Schottenfeld D: Cervical human papillomavirus infection and intraepithelial neoplasia: a review. J Natl Cancer Inst Monogr. 1996, 17-25.

    Google Scholar 

  41. Canadian Cancer Registry, Public Health Agency of Canada, 1990-2001. 2009

  42. De Carvalho NS, Roteli-Martins CM, Teixeira J, Naud P, de Borba P, Zahaf T, Sanchez N, Schuind A: Immunogenicity and safety of HPV-16/18 AS04-Adjuvanted vaccine up to 7.3Y. 25th International papillomavirus conference; 8. 2009, May

    Google Scholar 

  43. David MP, Van HK, Hardt K, Tibaldi F, Dubin G, Descamps D, Van DP: Long-term persistence of anti-HPV-16 and -18 antibodies induced by vaccination with the AS04-adjuvanted cervical cancer vaccine: Modeling of sustained antibody responses. Gynecol Oncol.

  44. Hethcote HW, Lewis MA, van den Driessche P: An epidemiological model with a delay and a nonlinear incidence rate. J Math Biol. 1989, 27: 49-64.

    Article  CAS  PubMed  Google Scholar 

  45. Barnabas RV, Laukkanen P, Koskela P, Kontula O, Lehtinen M, Garnett GP: Epidemiology of HPV 16 and cervical cancer in Finland and the potential impact of vaccination: mathematical modelling analyses. PLoS Med. 2006, 3: e138-10.1371/journal.pmed.0030138.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Anderson RM, May RM, Anderson B: Infectious diseases of humans: dynamics and control. 1992, New York, USA: Oxford University Press

    Google Scholar 

  47. Bauch CT, Anonychuk AM, Effelterre TV, Pham BZ, Merid MF: Incorporating Herd Immunity Effects into Cohort Models of Vaccine Cost-Effectiveness. Med Decis Making. 2009

    Google Scholar 

  48. Krahn M, McLachlin M, Pham B, Rosen B, Sander B, Grootendorst P, Tomlinson G, John-Baptiste A, Frikemerid M, Hong Chen M, et al: Liquid-based techniques for cervical cancer screening: a systematic review and cost-effectiveness analysis. Ottawa. 2008, []

    Google Scholar 

  49. Prices and prices indexes: Consumer price indexes. []

  50. Colucci R, Hryniuk W, Savage C: HPV vaccination programs in Canada: Are we hitting the mark?. Report Card on Cancer in Canada. 2008, 7-10. 10.2165/00019053-200422180-00002.

    Google Scholar 

  51. GDP per capita. []

  52. Sauvageau C, Duval B, Gilca V, Lavoie F, Ouakki M: Human papilloma virus vaccine and cervical cancer screening acceptability among adults in Quebec, Canada. BMC Public Health. 2007, 7: 304-10.1016/j.vaccine.2007.09.003.

    Article  PubMed  PubMed Central  Google Scholar 

  53. Is it cost-effective to vaccinate girls and boys with the HPV vaccine?. []

  54. Beutels P, Edmunds WJ, Antonanzas F, De Wit GA, Evans D, Feilden R, Fendrick AM, Ginsberg GM, Glick HA, Mast E, et al: Economic evaluation of vaccination programmes: a consensus statement focusing on viral hepatitis. Pharmacoeconomics. 2002, 20: 1-7. 10.1056/NEJMoa071430.

    Article  PubMed  Google Scholar 

  55. Bos JM, Beutels P, Annemans L, Postma MJ: Valuing prevention through economic evaluation: some considerations regarding the choice of discount model for health effects with focus on infectious diseases. Pharmacoeconomics. 2004, 22: 1171-1179. 10.2165/00019053-199915040-00004.

    Article  PubMed  Google Scholar 

  56. Crott R: Economic analysis of HPV-vaccines: not so simple?. Vaccine. 2007, 25: 7717-10.1016/j.vaccine.2007.09.003.

    Article  PubMed  Google Scholar 

  57. Mayrand MH, Duarte-Franco E, Rodrigues I, Walter SD, Hanley J, Ferenczy A, Ratnam S, Coutlee F, Franco EL: Human papillomavirus DNA versus Papanicolaou screening tests for cervical cancer. N Engl J Med. 2007, 357: 1579-1588. 10.1056/NEJMoa071430.

    Article  CAS  PubMed  Google Scholar 

  58. Mittmann N, Trakas K, Risebrough N, Liu BA: Utility scores for chronic conditions in a community-dwelling population. Pharmacoeconomics. 1999, 15: 369-376. 10.2165/00019053-199915040-00004.

    Article  CAS  PubMed  Google Scholar 

  59. Health Outcomes Data Repository. []

  60. Insinga R: Health State transitions following an abnormal pap smear: implications for health utility assessment in cost-effectiveness analyses. 22nd International Papillomavirus Conference and Clinical Workshop. 2005

    Google Scholar 

  61. Appendix 11: Human Papillomavirus. Vaccines for the 21st century: a tool for decisionmaking. Edited by: Stratton KR, Durch JS, Lawrence S. 2000, Washington D.C.: National Academy Press, 213-221.

Pre-publication history

Download references


The authors would like to thank Baudouin Standaert for his guidance on this project. The authors also thank Reid Robson, Rachel Emerson, Laure Delbecque, Parisa Airia and Mariam Mir for their assistance with this paper.

Cervarix™ is a trade mark of the GlaxoSmithKline group of companies.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Andrea M Anonychuk.

Additional information

Competing interests

Funding for this project was provided by GlaxoSmithKline Canada. AA, ND, GVK are employed by GlaxoSmithKline. CTB and MFM have held various research contracts with GlaxoSmithKline Canada.

Authors' contributions

AMA contributed to the development of the research question and methodology, analysis, interpretation of results, writing of the report, and project management. CTB contributed to the development of the research question, developed the methodology, and contributed to the interpretation of results and writing of the report. MFM and ND contributed to the methodology, analysis, interpretation of results. GVK contributed to the interpretation and writing of the report. All authors reviewed this paper.

Andrea M Anonychuk, Chris T Bauch contributed equally to this work.

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Anonychuk, A.M., Bauch, C.T., Merid, M.F. et al. A cost-utility analysis of cervical cancer vaccination in preadolescent Canadian females. BMC Public Health 9, 401 (2009).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: