Cost-Effectiveness Evaluation of Quadrivalent Human Papilloma Virus Vaccine for HPV-Related Disease in Iran

authors:

avatar Mohsen Khatibi a , avatar Hamid Reza Rasekh b , * , avatar Zohreh Shahverdi c , avatar Hamid Reza jamshidi d

Department of Pharmacoeconomics and Pharmaceutical Administration, School of Pharmacy, Shahid Beheshti University of Medical Sciences, Tehran, Iran.
Department of Pharmacoeconomics and Pharmaceutical Administration, School of Pharmacy, Shahid Beheshti University of Medical Sciences, Tehran, Iran.
Department of Obstetrics and Gynecology, Shahid Beheshti University of Medical Sciences, Tehran, Iran.
Department of Pharmacology , School of Medicine, Shahid Beheshti University of Medical Sciences, Tehran, Iran.

how to cite: Khatibi M, Rasekh H R, Shahverdi Z, jamshidi H R. Cost-Effectiveness Evaluation of Quadrivalent Human Papilloma Virus Vaccine for HPV-Related Disease in Iran. Iran J Pharm Res. 2014;13(Suppl):e125460. https://doi.org/10.22037/ijpr.2014.1478.

Abstract

Human Papilloma Virus (HPV) vaccine has been added recently to the Iran Drug List. So, decision makers need information beyond that available from RCTs to recommend funding for this vaccination program to add it to the National Immunization program in Iran. Modeling and economic studies have addressed some of those information needs in foreign countries. In order to determine the long term benefit of this vaccine and impact of vaccine program on the future rate of cervical cancer in Iran, we described a model, based on the available economic and health effects of human papilloma virus (HPV), to estimate the cost-effectiveness of HPV vaccination of 15-year-old girls in Iran. Our objective is to estimate the cost-effectiveness of HPV vaccination in Iran against cervical cancer based on available data; incremental cost-effectiveness ratio (ICER) calculations were based on a model comparing a cohort of 15-year-old girls with and without vaccination. We developed a static model based on available data in Iran on the epidemiology of HPV related health outcome. The model compared the cohort of all 15-year old girls alive in the year 2013 with and without vaccination. The cost per QALY, which was found based on our assumption for the vaccination of 15-years old girl to current situation was 439,000,000 Iranian Rial rate (IRR). By considering the key parameters in our sensitivity analysis, value varied from 251,000,000 IRR to 842,000,000 IRR. In conclusion, quadrivalent HPV vaccine (Gardasil) is not cost-effective in Iran based on the base-case parameters value.

Introduction

The human papilloma virus (HPV) is one of the most common sexually transmitted viruses. Chronic infection with certain subtypes of the HPV is the primary cause of cervical cancer and its precancerous lesions. At least 50% of the adult population is infected with this virus during their lifetime. Despite of screening programs for cervical cancer, it remains the second most common cause of cancer-related death among women worldwide (1, 2).

Vaccines are essential tools for preventing the diseases. They will protect the vaccinated individual and help to protect the community by reducing the spread of infectious agents (3).

Gardasil vaccine is part of the national immunization program in the United States, Canada and Australia. It is covered by insurance in Canada and Australia. Gardasil has been registered and used in 124 countries. Moreover, it is part of the national immunization program in 19 countries of the world which means it is administered to all girls aged 11-12 years (4).

In order to determine the long term benefit of this vaccine and impact of vaccine program on the future rate of cervical cancer, many pharmacoeconomists have used mathematical Models. Some models focused on cost-effectiveness of different strategies (5, 6, 7 and 8).

Human Papilloma Virus (HPV) can cause cervical cancer, genital warts, and some other anogenital cancers (9, 10, 11, 12 and 13). In Iran, HPV types 16 and 18 cause around 72%of cervical cancer, the attribution of HPV 16 and 18 are 56.4% and 15.6% respectively (14). HPV types 6 and 11 cause 90% of cases of genital warts (13, 15). The Estimated percentages of health outcomes attributable to various HPV types showed in Table 1.

Table1

Estimated percentages of health outcomes attributable to various HPV types

Health outcomeHPV 6,11HPV 16HPV 18Reference
Cervical Cancer (Iran)056.4%15.6% 13-
Genital Warts90%00( 20- 35- 36-)
Cervical Intraepithelial Neoplasia grade I (CIN I)6.30%19.40%9.20% 33-
Cervical Intraepithelial Neoplasia grade II (CIN II)045.80%10% 34-
Cervical Intraepithelial Neoplasia grade III (CIN III)045.80%10% 34-

Numerous studies have been conducted to evaluate its cost-effectiveness by calculating the cost, which is necessary for one quality-adjusted life year (QALY). These include various models such as Markov model, decision model, dynamic model, transmission or a combination one (16).

Studies conducted so far have mostly addressed its impact on cervical cancer and cervical intraepithelial neoplasia (CIN) (17). Few studies have dealt with the cost-effectiveness of Gardasil regarding other malignancies caused by the HPV (18).

The national immunization program of girls prior to sexual activity has shown that it reduces the HPV-related mortality and morbidity and it will be cost-effective (19, 20, 21, 22 and 23). Since HPV infection is asymptomatic, it is growing silently (24).

Vaccination can cause immunity both directly and indirectly through herd immunity (25). Recently a HPV vaccine was added to the Iran Drug List, but it has not been added to the national immunization program till now. Little data exist on the epidemiology of cervical cancer or the prevalence of HPV infection in Iran due to unavailability of national screening program (14).

About 663 new cases of cervical cancer are diagnosed each year in Iran (26), there is no organized population-based cervical screening program, therefore most cervical cancers are diagnosed in late stages, Pap Smear will be done in Iran when they refer to the gynecologist (14). The incidence rate of cervical cancer is low in Iran and many other Muslim countries (27).

The mean cervical cancer ASMR (Age Specific Mortality Rate) for Iran was 1.04 per 100,000 (14). The prevalence of HPV infection was 76% in Iranian cervical cancer patients and 7% among healthy Iranian women (14). The HPV types isolated in Iranian cervical cancer patients includes: HPV 16 (54%), 18 (14%) and 31 (6%) (14).

In Iran, the HPV Prevalence and HPV type is similar to European countries but the rate of Cancer is very low (28).

The US Food and Drug Administration approved a quadrivalent (HPV 6, 11, 16, 18) vaccine (Gardasil, manufactured by Merck and Co.) in June 2006. It was approved for use in girls and women 9-26 years of age (13). The vaccine efficacy will be almost 100% if given to young women before sexual exposures (11, 13 and 29). Other indications have been mentioned in table 2.

Table 2

FDA Approval History for Gardasil

Dateapproved indications
Dec, 2010Gardasil Approved to Prevent Anal Cancer
Oct, 2009Gardasil Approved for Use in Boys and Young Men
Sep, 2008Gardasil Approved for Preventing Certain Vulvar and Vaginal Cancers
Jun, 2006Gardasil Approved for Prevention of Cervical Cancer other HPV-related Diseases

Several cost-effectiveness studies and HTAs (Health Technology Assessments) have been done all over the world for estimation of potential cost-effectiveness of this vaccine; these studies applied a Markov model, a decision model, a dynamic transmission model, or a combination thereof. There was another simplified model that estimated the cost-effectiveness of adding HPV vaccination of 12-year-old girls in the United States (30). This model is very useful for Iran because required data for estimating the cost-effectiveness of Quadrivalent HPV vaccine is available based on this model. Other models are complicated and need many data which are not available in Iran.

The objective of this study is to estimate the cost-effectiveness of HPV vaccination in Iran against cervical cancer based on available data; incremental cost-effectiveness ratio (ICER) calculations were based on a model comparing a cohort of 15-year-old girls with and without vaccination.

Experimental

We developed a static model based on available data in Iran on the epidemiology of HPV related health outcome. The model compared the cohort of all 15-year old girls alive in the year 2013 with and without vaccination.

Similar to the US STUDY MODEL approach, two strategies has been considered; the first one is reference (baseline) strategy including the current situation. As practised in Iran we didn’t consider any national cervical cancer screening. The second one is comparator strategy including the routine quadrivalent (16/18/6/11) HPV vaccination of girls at age 15 with 70% coverage at the year of vaccination.

In other word, the first strategy including the cohort of 15-year-old girl group started from 2013 based on Iranian population without vaccination and another strategy including the same group with vaccination strategy. In both strategies, we used an incidence-based model of the health and economic effects of HPV-related health outcomes including cervical cancer, CIN I. CIN II, CIN III and genital wart.

We adopted a governmental perspective and included all direct medical costs and benefits regardless of who incurred the costs or received the benefits (31, 32). No other societal costs were included in the analysis. Our study question was “What is the cost per QALY gained by adding vaccination of 15-year-old girls to existing situation (no official screening) in Iran?”

According to a World Bank report published in 2012, the life expectancy of women in Iran is 74 years. We considered a hypothetical population of persons 15-74 years of age. The number of 15-year-old girls was based on 2011 sex-specific population estimates. The number of 16-year-old girls was calculated based on the number of 15-year-olds and the probability of survival from age 15 years to age 16 years. We continued this calculation based on the specific age related mortality for other groups till 74-years in an analogues manner.

Table 3

Parameters and Base Case Value

ParametersBase case value
Vaccine efficacy100%
Protection durationLifetime
Time horizon59-Ys
Vaccine coverage70%
Euro rate 33,500
Price of Gardasil (Euro) 85.5
Margins1.460
Vaccine cost 4,181,805
Vials of vaccine needed without considering the booster dose3
Vaccine cost per series12,545,415
Financial discount rate (Annual)3%
Health discount rate (Annual)3%
Life expectancy (Year)74

In this cohort we assumed that 3-doses of HPV vaccine will be administered to 15-year-old girls. This cohort will be started in year 15 and will be continued through year 74. Vaccination coverage assumption was 70%. Vaccination efficacy was assumed to be 100% for the HPV 16, 18, 11, and 6 related outcomes (11, 31 and 33). The duration of vaccine protection was assumed to be life-long, but in sensitivity analysis we assumed to administer another booster dose in the year of 10. Based on the proposed price of vaccine in Iran, the cost of full series of quadrivalent HPV vaccine considered 265.5 euro per series.

Table 4

Age Specific Incidence Rate (ASIR) per 100,000 in Iran

YearASIR(CC)ASIR (CIN I)ASIR(CIN II)ASIR (CIN III)ASIR(FGW)
10-14 Y0.0700043
15-19 Y01608030287
20-24 Y0.22510320130620
25-29 Y0.74140380410394
30-34 Y1.3240140180265
35-39 Y2.75240140180199
40-44 Y4.351205050139
45-49 Y6.561205050144
50-54 Y8.9170401092
55-59 Y7.4270401086
60-64 Y6.224010076
65-69 Y9.424010055
70-74 Y7.732001040
More than 7511.72001021
Reference (25-) (15-, 36-) (15-, 36-) (15-, 36-) (15-, 36-)

We examined the following HPV-related health outcomes: cervical cancer; CIN grades I, II, and III; and genital warts; the age-specific incidence rates of the HPV-related adverse health outcomes were used to estimate the potential outcomes that could be obtained through their life for both strategies.

Age-specific incidence rates of cervical cancer (ASIR CC) were extrapolated from 2008 population-based cancer registries in Iran (26). Age-specific incidence rates of CIN grades I, II, and III, and female genital warts (FGW) were based on estimates obtained from the literature (16, 37).

Table 5

Treatment cost of HPV adverse health outcomes

HPV adverse health outcomesTreatment Cost/case (IRR)
Cervical Cancer: 80,000,000.00
Genital Wart:2,500,000.00
CIN I: 5,000,000.00
CIN II: 10,000,000.00
CIN III: 30,000,000.00

Treatment cost of HPV adverse health outcomes

The cost of care of HPV adverse health outcomes calculated for patients in current clinical practice in Iran in 2013 for each stage of cervical cancer, CIN I, CIN II, CIN III and genital wart. To determine the different kind of cost averted by above mentioned adverse health outcomes, we evaluated the patients files archived in the hospitals and gynecologists private offices. Finally after finding the different kind of direct costs, we referred to the official list of tariffs (38) to calculate the cumulative cost of each HPV adverse health outcomes.

Costs averted and qalys saved by vaccination

After calculating the potential HPV adverse health outcomes, we defined the infected population in both strategies. In order to define the cost averted by vaccination, we calculated the cost of infected population via multiplying the infected people in both strategies to the estimated cost of each HPV related disease and then reduced the cost of vaccinated strategy to non-vaccinated strategy. In order to define the QALY averted by vaccination, we calculated the QALY loss of infected population via multiplying the infected people in both strategies to the estimated QALY loss of each HPV related disease and then reduced the QALY loss of vaccinated strategy to non-vaccinated strategy.

To estimate the discounted QALYs lost per case of cervical cancer, CIN I, CIN II, CIN III, and genital warts were based on published estimates of the quality of life without adverse these health outcomes (39) and the estimated reduction in quality of life associated with these HPV-related health outcomes (40, 41, 42, 43), for cervical cancer was assumed to have of six outcomes at diagnosis:

1- Local lesion lead to survive: QALY loss assumed 0.27 for 4 months and 0.07 for long life till 74 years

2- Local lesion lead to death: QALY loss assumed 0.36 for 3 years and 1 for long life till 74 years

3- Regional metastatic lead to survive: QALY loss assumed 0.37 for 3 years and 0.1 for long life till 74 years

4- Regional metastatic lead to death: QALY loss assumed 0.41 for 3 years and 1 for long life till 74 years

5- Distant metastatic lead to survive: QALY loss assumed 0.45 for 3 years and 0.24 for long life till 74 years

6- Distant metastatic lead to death: QALY loss assumed 0.45 for 3 years and 1 for long life till 74 years.

Based on expert opinion, we assumed the distributions of cervical cancer stage at diagnosis as follows:

1- Local: 35% and the probability of survival 0.87

2- Regional: 40% and the probability of survival 0.5

3- Distant: 25% and the probability of survival 0.09

For CIN I; we assumed loss in quality of life of 0.03 for 18 months, and no quality loss after this 18 months (42).

For CIN II; we assumed loss in quality of life of 0.07 for 18 months, and no quality loss after these 18 months (42).

For CIN III; we assumed loss in quality of life of 0.07 for 2 years, and no quality loss after these 18 months (42). For genital warts in females, we assumed loss of quality of life and the duration of such loss were assumed to be one of the following four scenarios (42):

1- 0.05 loss for 3 months, with probability of 0.475

2- 0.1 loss for 6 months, with probability of 0.025

3- 0.15 loss for 3 months, with probability of 0.475

4- 0.15 loss for 6 months, with probability of 0.025

Incremental cost per QALY gained

Vaccination costs, averted treatment costs and the number of QALYs saved were calculated for each year over a 62-year period, discounted to present value by using an annual discount rate of 3% for both health outcomes and costs (31).

The incremental cost per QALY gained by using HPV vaccine to existing situation was calculated as (V-A)/Q, where V is the cost of vaccination, A is the averted treatment costs due to vaccination, and Q is the number of QALYs saved due to vaccination (32).

Sensitivity analysis

We applied 1-way sensitivity analysis in which we varied one of parameter values while holding other parameters at their base-case values. The parameters we varied included the cost of the vaccine series, the cost per case of all HPV-related health outcomes (±25% of their base-case values); the discount rate (0%, 7.2%); the incidence rates of health outcomes (±25% of their base-case values for CIN 1, CIN 2, CIN 3, and genital warts, and duration of vaccine protection (injection of one booster dose in year 10).

Results

To calculate the number of woman at each stage, we used the age specific survival rate. Table 6 shows the woman at each stage calculated based on the age specific survival rate.

Life expectancy of women in Iran is 74 years. Therefore, the life expectancy of 15- Years old girls was 59 years.

Table 6

Number of woman at each stage calculated based on the age specific survival rate

Year1390Year1390
10-14 Y2,783,04745-49 Y2,003,143
15-19 Y3,259,60750-54 Y1,762,295
20-24 Y4,212,92255-59 Y1,353,485
25-29 Y4,318,02060-64 Y981,945
30-34 Y3,456,09665-69 Y700,389
35-39 Y2,720,78570-74 Y558,821
40-44 Y2,420,370More than 75895,799

The economic output of interest from the model included total discounted costs and the incremental cost per QALY gained ratio. Both costs and QALYs were discounted at a 3% annual rate (31). We measured the cost-per-QALY ratio as the incremental cost difference between the two strategies divided by the incremental QALY difference between the two strategies. Under base-case parameter values, the estimated cost per QALY gained by adding vaccination of 15-year-old girls to existing cervical cancer screening was around 440,000,000 IRR (Iranian Rial Rate).

There is much controversy about the choice of discount rates, and whether costs and benefits should be discounted at the same or different rates. In recognition of this, Drummond and Jefferson (44) suggest that sensitivity analysis be done using alternative discount rates, including zero. One reason that zero discount rates have been suggested in both the health and the environmental area is that the benefits are likely to be felt in the future, whereas many costs are incurred now. Considering the zero percent discount rate of cost, the ICER is 414,000,000 IRR.

For the HPV vaccine, the exact total duration of protection is not known yet, because the current maximum length of clinical trials is around 6 years. Consequently, it could be argued that base-case analysis on the cost-effectiveness of the HPV vaccine should not use durations of protection beyond 6 years, let alone lifelong protection. By considering 10 years protection, the ICER is around 599,000,000 IRR.

When the cost of treatment for HPV-related complications (Cervical cancer, CIN I, CIN II, CINII, and Genital wart) varied from -25% to +25%, the cost per QALY gained ranged from 447,000,000 IRR to 428,000,000 IRR. When the time horizon varied from 59 to 85 years, the cost per QALY gained ranged from 440,000,000 IRR to 250,000,000 IRR. Changes in the other parameter values such as QALYs that is associated with HPV-related health outcomes also affected the results: when health related outcome varied from -25% to +25%, the cost per QALY gained from 489,000,000 IRR to 387,000,000 IRR. When we excluded the other HPV related complications such as CIN I, CIN II, CIN III and genital wart, the ICER changed to 842,000,000 IRR. In other words, if we consider the prevention of cervical cancer, only the cost per QALY is 842,000,000 IRR.

Discussion

National decision makers and reimbursement committees prefer to utilize cost utility Analysis. In some countries, decision makers have established ICER thresholds to determine whether a healthcare technology is cost-effective. For example, in the United States, USD50, 000 per QALY gained has been considered as a threshold of cost-effectiveness analysis (45, 46). Similarly in England, the National Institute for Health and Clinical Excellence (NICE) has used ICER of 30,000£ per QALY gained as a threshold of cost effective for the National Health Service (NHS) (47, 48).

For developing countries like Iran, World Health Organization (WHO) has recommended a cost-effectiveness threshold indicating that a healthcare technology is cost effective if the ICER is less than three times the GDP (Gross Domestic Production) per capita (49) WHO’s recommendation about threshold of developing countries considers ICER less than triplet of GDP of Iran for 2012 is 5,810 $. Based on WHO recommendation, ICER less than 17,430 USD per QALYs could be considered cost-effective. Since the official exchange rate of USD to IRR is 24,500 in November 2013, the cost-effective threshold in Iran would be around 427,000,000 IRR. In Iran, the availability of data and information for this kind of disease is very limited. Therefore, we used a very simple model in order to estimate the cost-effectiveness of HPV vaccine versus current situation. There is no formal and national screening program for cervical cancer diagnosis in Iran. The cost per QALY, which was found based on our assumption for the vaccination of 15-years old girl to current situation was 439,000,000 IRR. By considering the key parameters in our sensitivity analysis, value varied from 251,000,000 IRR to 842,000,000 IRR. Table 7 shows the summaries of our assumptions and related cost per QALY.

Table 7

Summaries of our assumptions and related Cost per QALY

Assumption Cost of Vaccinated (IRR)QALY loss of VaccinatedCost of Non VaccinatedQALY loss of Non VaccinatedCost effectiveness Threshold IranICERIs it Cost Effective?
Base Case5,527,698,780,568 12,295 1,060,714,746,135 22,468427,035,000 439,092,468No
7.2% Financial Discount Rate 5,272,538,157,362 12,295 636,625,601,020 22,468427,035,000 455,697,686No
0% Financial Discount Rate 5,902,143,602,544 12,295 1,688,572,121,292 22,468427,035,000 414,182,699Yes
Booster Dose 7,157,041,160,945 12,295 1,060,714,746,135 22,468427,035,000 599,252,425No
Cost+25% 5,399,764,452,681 2,295 848,571,796,908 22,468427,035,000 447,369,948No
Cost-25% 5,687,616,690,428 12,295 1,325,893,432,669 22,468427,035,000 428,745,618No
Cervical Cancer Only 5,931,983,028,333 16,687 1,060,714,746,135 22,468427,035,000 842,587,744No
Health Related HPV-25 5,403,066,718,428 10,975 855,225,834,442 20,270427,035,000 489,286,661No
Health Related HPV+25 5,683,488,858,243 13,944 1,317,575,885,751 25,216427,035,000 387,351,396Yes
100 Y Life Expectancy5,527,698,780,56819,775 1,060,714,746,13537,541427,035,000251,434,341Yes

Comparing the cost-effectiveness threshold of Iran to base-case scenario shows that this vaccine is not cost-effective. According to the sensitivity analysis, the most important parameter, which had more and substantial influence, was limiting the indications of Quadrivalent HPV vaccine to cervical cancer only. The time horizon and effectiveness of vaccination during the long life time also had substantial effect on the results. The variation of the cost of treatment had the minimum effect.

Only for the below three scenarios, quadrivalent HPV vaccine is cost-effective in Iran:

1- If the life expectancy of Iranian women would be 100 years instead of 74 years

2- If the HPV related adverse effect considered 25% more than what we assumed

3- If the financial discount rate assumed zero percent

This study is the first cost-effectiveness study of HPV vaccine and vaccination in Iran and the issue of declining immunity was addressed by assigning a higher cost per vaccination (one dose only) as in the sensitivity analysis to show the cost of a booster. Another advantage of this study is fewer assumptions were needed than more complex Markov and hybrid models did. Since the base-case is established by the incidence rate of different diseases, there is no need to have the following domestic information for modeling: the possible progression from HPV infection to CIN I, CIN II, CIN III, cervical cancer, the probability of HPV infection, the probability of HPV transmission, cervical cancer screening and sexually transmitted disease prevention activities, but the potential benefits of preventing vaginal, anal, vulvar, and oropharyngeal cancers in our estimated cost-effectiveness of HPV vaccine were not considered and the QALY loss related to the initial distress of receiving an abnormal Pap smear result (50) was not considered in this study. Thus, the benefits of HPV vaccination may be underestimated because of reducing the number of abnormal Pap smear results (51, 52). The potential benefits of herd immunity, boys and men vaccination were not considered in this study. As the domestic data for the incidence of CIN I, CIN II, CIN III, and genital warts was not available to us, the international incidence rate has been assumed for mentioned diseases.

In conclusion, quadrivalent HPV vaccine is not cost-effective in Iran based on the base-case parameters values.

Acknowledgements

References

  • 1.

    Smith JS, Lindsay L, Hoots B, Keys J, Franceschi S, Winer R, Clifford GM. Human papillomavirus type distribution in invasive cervical cancer and high-grade cervical lesions: a meta-analysis update. Int. J. Cancer. 2007;121:621-632. [PubMed ID: 17405118].

  • 2.

    Parkin DM, Bray F. Chapter 2: the burden of HPV-related cancers. Vaccine. 2006;24:11-25.

  • 3.

    Tavajohi S, Rastegar H, Ostad SN, Rezayat SM, Ghahramani MH. Evaluation of Potency of Measles Vaccine used in Iran: Comparison of WHO and NIBSC Method in Cell Culture. Iran. J. Pharm. Res. 2005;3:155-160.

  • 4.

    Koulova A, Tsui J, Irwin K. Country recommendations on the inclusion of HPV vaccines in national immunization programmes among high –income countries. Vaccine. 2008;26:6529-6541. [PubMed ID: 18805453].

  • 5.

    Garnett GP, Kim JJ, French K. Chapter 21: modeling the impact of HPV vaccines on cervical cancer and screening programmes. Vaccine. 2006;24:178-186.

  • 6.

    Dasbach EJ, Elbasha EH, Insinga RP. Mathematical models for predicting the epidemiologic and economic impact of vaccination against human papillomavirus infection and disease. Epidemiol. Rev. 2006;28:88-100. [PubMed ID: 16740585].

  • 7.

    Newall AT, Beutels P, Wood JG. Cost-effectiveness analyses of human papillomavirus vaccination. Lancet Infect. Dis. 2007;7:289-296. [PubMed ID: 17376386].

  • 8.

    Goldie S, Goldhaber-Fiebert JD, Garnett GP. Chapter 18: public health policy for cervical cancer prevention. The role of decision science, economic evaluation, and mathematical modeling. Vaccine. 2006;24:155-163.

  • 9.

    Institute of Medicine. Vaccines for the 21st century: a tool for decision making. Washington: National Academy of Sciences; 2000.

  • 10.

    Parkin DM, Bray F. the burden of HPV-related cancers. Vaccine. 2006;24:11-25.

  • 11.

    Centers for Disease Control, Prevention. Quadrivalent human papillomavirus vaccine: recommendations of the Advisory Committee on Immunization Practices (ACIP). MMWR Recomm. Rep. 2007;56:1-24.

  • 12.

    Clifford GM, Smith JS, Aguado T, Franceschi S. Comparison of HPV type distribution in high-grade cervical lesions and cervical cancer: a meta-analysis. Br. J. Cancer. 2003;89:101-105. [PubMed ID: 12838308].

  • 13.

    Dunne EF, Markowitz LE. Genital human papillomavirus infection. Clin. Infect. Dis. 2006;43:624-629. [PubMed ID: 16886157].

  • 14.

    Khorasanizadeh F, Hassanloo J, Khaksar N, Mohammad Taheri S, Marzaban M, H Rashidi B, Akbari Sari A, Zendehdel K. Epidemiology of cervical cancer and human papilloma virus infection among Iranian women - analyses of national data and systematic review of the literature. Gynecol. Oncol. 2013;128:277-281. [PubMed ID: 23200918].

  • 15.

    Greer CE, Wheeler CM, Ladner MB, BeutnerK, Coyne MY, Liang H. Human papillomavirus (HPV) type distribution and serological response to HPV type 6 virus-like particles in patients with genital warts. J. Clin. Microbiol. 1995;33:2058-2063. [PubMed ID: 7559948].

  • 16.

    Insinga RP, Dasbach EJ, Elbasha EH. Assessing the annual economic burden of preventing and treating anogenital human papillomavirus-related disease in the US: analytic framework and review of the literature. Pharmacoeconomics. 2005;23:1107-1122. [PubMed ID: 16277547].

  • 17.

    Fleurence RL, Dixon JM, Milanova TF. Review of the economic and quality-of-life burden of cervical human papillomavirus disease. Am. J. Obstet. Gynecol. 2007;196:206-212. [PubMed ID: 17346523].

  • 18.

    Hu D, Goldie S. The economic burden of noncervical human papillomavirus disease in the United States. Am. J. Obstet. Gynecol. 2008;198:1-7. [PubMed ID: 18166295].

  • 19.

    Sanders GD, Taira AV. Cost-effectiveness of a potential vaccine for human papillomavirus. Emerg. Infect. Dis. 2003;9:37-48. [PubMed ID: 12533280].

  • 20.

    Kulasingam SL, Myers ER. Potential health and economic impact of adding a human papillomavirus vaccine to screening programs. JAMA. 2003;290:781-789. [PubMed ID: 12915431].

  • 21.

    Goldie SJ, Kohli M, Grima D. Projected clinical benefits and cost-effectiveness of a human papillomavirus 16/18 vaccine. J. Natl. Cancer Inst. 2004;96:604-615. [PubMed ID: 15100338].

  • 22.

    Brisson M, Van de Velde N, De Wals P. The potential cost-effectiveness of prophylactic human papillomavirus vaccines in Canada. Vaccine. 2007;25:5399-5408. [PubMed ID: 17561316].

  • 23.

    Kulasingam S, Connelly L, Conway E. A costeffectiveness analysis of adding a human papillomavirus vaccine to the Australian national cervical cancer screening program. Sex Health. 2007;4:165-175. [PubMed ID: 17931529].

  • 24.

    Giuliano AR. Human papillomavirus vaccination in males. Gynecol. Oncol. 2007;107:24-26.

  • 25.

    Garnett GP. Role of herd immunity in determining the effect of vaccines against sexually transmitted disease. J. Infect. Dis. 2005;191:97-106.

  • 26.

    Ministry of Health and medical education of Iran. Iranian Annual National Cancer Registration Report. 2008-2009.

  • 27.

    Arbyn M, Castellsague X, de Sanjose S, Bruni L, Saraiya M, Bray F. Worldwide burden of cervical cancer in 2008. Ann. Oncol. 2011;22:2675-2686. [PubMed ID: 21471563].

  • 28.

    Hamkar R, Mokhtari_Azad T. Prevalence of variouse types of HPV among cervical cancer and normal biopsy in north of Iran. Iran. J. Infec. Dis. Tropical Medicine. 2004;8:22-29.

  • 29.

    Elbasha EH, Dasbach EJ, Insinga RP. Model for assessing human papillomavirus vaccination strategies Emerg. Infect. Dis. 2007;13:28-41.

  • 30.

    Chesson HW, Ekwueme DU, Saraiya M, Dunne EF, Markowitz LE. The cost-effectiveness of male HPV vaccination in the United States. Vaccine. 2011;29:8443-8450. [PubMed ID: 21816193].

  • 31.

    Gold MR, Siegel JE, Russell LB, Weinstein MC. Cost-effectiveness in health and medicine. New York: Oxford University Press; 1996.

  • 32.

    Haddix AC, Teutsch SM, CorsoPS. . Prevention effectiveness: a guide to decision analysis and economic evaluation. New York: Oxford University Press; 2003.

  • 33.

    Harper DM, Franco EL, Wheeler CM, Moscicki AB, Romanowski B, Roteli-Martins CM. 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. [PubMed ID: 16631880].

  • 34.

    US Cancer Statistics Working Group. United States cancer statistics: 2003 incidence and mortality. Atlanta: Centers for Disease Control and Prevention and National Cancer Institute; 2006.

  • 35.

    Hankey BF, Ries LA, Edwards BK. The surveillance, epidemiology, and end results program: a national resource. Cancer Epidemiol. Biomarkers Prev. 1999;8:1117-1121. [PubMed ID: 10613347].

  • 36.

    Watson M, Saraiya M, Weir H, Ahmed F. Leading partnerships to assess the burden of HPV-related cancers. 2007 CDC Cancer. 2. Atlanta: Centers for Disease Control and Prevention; 2007. p. 1-5.

  • 37.

    Insinga RP, Glass AG, Rush BB. Diagnoses and outcomes in cervical cancer screening: a population-based study. Am. J. Obstet. Gynecol. 2004;191:105-113. [PubMed ID: 15295350].

  • 38.

    Iran ministry of health and medical education. Tariff of healthcare services in public and private sectors in Iran (2012). Tehran: Iran ministry of health and medical education press; 2012.

  • 39.

    Gold MR, Franks P, McCoy KI, Fryback DG. Toward consistency in cost-utility analyses: using national measures to create condition-specific values. Med. Care. 1998;36:778-792. [PubMed ID: 9630120].

  • 40.

    Kulasingam SL, Myers ER. Potential health and economic impact of adding a human papillomavirus vaccine to screening programs. JAMA. 2003;290:781-789. [PubMed ID: 12915431].

  • 41.

    Goldie SJ, Kohli M, Grima D, Weinstein MC, Wright TC, Bosch FX. Projected clinical benefits and cost-effectiveness of a human papillomavirus 16/18 vaccine. J. Natl. Cancer Inst. 2004;96:604-615. [PubMed ID: 15100338].

  • 42.

    Institute of Medicine. Vaccines for the 21st century: A tool for decision making. Washington, DC: National Academy of Sciences; 2000.

  • 43.

    Kim JJ, Wright TC, Goldie SJ. Cost-effectiveness of alternative triage strategies for atypical squamous cells of undetermined significance. JAMA. 2002;287:2382-2390. [PubMed ID: 11988059].

  • 44.

    Drummond MF, Jefferson TO. Guidelines for authors and peer reviewers of economic submissions to the BMJ. British Med. J. 1996;313:275-283.

  • 45.

    Eichler HG, Kong SX, Gerth WC, Mavros P, Jonsson B. Use of cost-effectiveness analysis in health-care resource allocation decision-making: How is cost-effectiveness thresholds expected to emerge? Value in Health. 2004;7:518-528. [PubMed ID: 15367247].

  • 46.

    Grosse SD. Assessing cost-effectiveness in healthcare: History of the $50,000 per QALY threshold. Expert Rev. Pharmacoeconomics Res. 2008;8:165-178.

  • 47.

  • 48.

    Devlin N, Parkin D. Does NICE have a cost-effectiveness threshold and what other factors influence its decisions? A binary choice analysis. Health Econ. 2004;13:437-452. [PubMed ID: 15127424].

  • 49.

  • 50.

    Kulasingam SL, Myers ER, Lawson HW, McConnell KJ, Kerlikowske K, Meinikow J. Cost-effectiveness of extending cervical cancer screening intervals among women with prior normal pap tests. Obstet. Gynecol. 2006;107:321-328. [PubMed ID: 16449119].

  • 51.

    Schiffman M. Integration of human papillomavirus vaccination, cytology, and human papillomavirus testing. Cancer. 2007;111:145-153. [PubMed ID: 17487850].

  • 52.

    Haupt RM. GARDASIL update. CDC Cancer Conference Program of Events. Atlanta: Centers for Disease Control and Prevention; 2007.