Cost effectiveness analysis of HPV primary screening and dual stain cytology triage compared with cervical cytology.

OBJECTIVES: To assess the clinical and cost-effectiveness of human papillomavirus (HPV) primary screening triage with p16/Ki-67 dual stain cytology compared to cytology.
METHODS: We conducted an Excel®-based budget impact model to estimate the preinvasive and invasive cervical cancer cases identified, mortality rate, direct medical costs, quality-adjusted life years (QALYs) and the incremental cost-effectiveness analysis of two strategies from the healthcare payer perspective. The study population is a cohort of women 30-65 years of age presenting for cervical screening.
RESULTS: HPV primary screening triage with p16/Ki-67 dual stain showed higher sensitivity without losing specificity compared to conventional Pap smear. The improving the screening performance leads to decrease the prevalence of precancerous lesion, annual incidence and mortality of cervical cancer. The incidence of cervical cancer case detected by new algorithm compared with conventional method were 31,607 and 38,927, respectively. In addition, the new algorithm was more effective and more costly (average QALY 24.03, annual cost $13,262,693) than conventional cytology (average QALY 23.98, annual cost $7,713,251). The incremental cost-effective ratio (ICER) per QALY gained was $1,395. The sensitivity analysis showed if the cost of cytology and HPV test increased three times, the ICER would fall to $303/QALY gained and increased to $4,970/QALY gained, respectively.
CONCLUSION: Our model results suggest that screening by use of HPV genotyping test as a primary screening test combined with dual stain cytology as the triage of HPV positive women in Thai population 30-65 years old is expected to be more cost-effective than conventional Pap cytology.

INTRODUCTION

Cervical cancer remains one of the most common cancers among Thai women. Even the age-standardized incidence rate (ASR) continue decreasing, it is still high (17.8/100,000) [1]. Nowadays, a comprehensive cervical cancer screening in Thailand includes two types; a cytology-based screening, and co-testing or combination of cytology plus human papillomavirus (HPV) testing. Thai government aims to reduce the ASR through HPV vaccination in combination with more effective screening strategy with HPV DNA testing. The latest update data regarding the cervical cancer screening in Thailand reported that the most cost-effectiveness strategy is high-risk human papillomavirus (hrHPV) testing. Compare with screening by cytology strategy, screening with hrHPV testing strategy could not only decrease the cost, but also detect more cases of cervical intraepithelial neoplasia 2 or higher (CIN2+) [2]. The results of this study support the full-scale implementation of HPV testing as primary cervical cancer screening in Thailand. Regarding the HPV positive women in HPV genotype algorithm form that study, the women positive for oncogenic HPV16/18 direct to colposcopy, with reflex cytological triage for women with other oncogenic types and direct referral for those in this group with high-grade cytological findings. The result of the Primary ASC-US LSIL Marker Study (PALMS) revealed that the p16/Ki-67 dual stain cytology combines superior sensitivity and non-inferior specificity over Pap cytology for detecting CIN2+ [3]. Different triage strategies involving p16/Ki-67 dual staining were evaluated in several studies [4567]. Thomas et al assessed the performance of p16/Ki-67 dual stain cytology for triaging HPV-positive women undergoing primary HPV screening. The sensitivity, positive predictive value (PPV), and negative predictive value (NPV) were higher for dual staining than Pap cytology triaging HPV positive women with p16/Ki-67 dual stain cytology [4]. We assumed that using dual stain cytology instead of conventional pap smear for triaging other oncologic type HPV positive women would increase the detection rate of the precancerous lesion (CIN2/3). The objective of this study is to evaluate the clinical and cost-effectiveness of HPV primary screening triage with p16/Ki-67 dual stain cytology compared to primary conventional cytology algorithms from a payer perspective (Thai government).

MATERIALS AND METHODS

1. Modeling approach

Excel-based budget impact model was constructed based on a Markov decision-analytic model in order to estimate the number of accumulated cases of CIN2, CIN3, cervical cancer and budget impact of each screening program. The Markov model was created based on the natural history of disease, including persistent high-risk HPV infection, CIN (low grade, high grade), invasive cervical cancer (ICC) and death. All women enter to screening model at the same time and exit from the model when they are diagnosed with CIN2, CIN3 or ICC. For women with undetected disease (either CIN or ICC) in the first round, they will be sent to routine screening or follow up testing. Some women may be lost to follow up and return unscheduled. The disease may regress, persist or progress spontaneously over time. Natural history and epidemiological data of HPV infection and cancer are based on published literature referenced in the model (Table 1). A model of the patient flow through a cervical cancer screening was showed in Fig. 1. The model considers 2 hypothetical cohorts, comparing costs and outcomes of a cohort screened with a new proposed algorithm (HPV primary screening triage with p16/Ki-67 dual stain cytology) with a control cohort of a conventional pap smear. Screening with the proposed algorithm is assumed to increase the detection rate of the precancerous lesion (CIN2/3) which leads to decrease the prevalence of cervical cancer, improve life expectancy and quality of life (QoL). The primary endpoint for the analysis was the incremental cost-effective ratio (ICER) per quality-adjusted life year (QALY) gained. ICER was calculated by additional annual cost divided by average QALY gained between two strategies.

Table 1

Data of HPV infection and cancer based on published reference

Clinical parameters				Input value
The performance of screening test [2425]
	Cytology (threshold = ASCUS)
		Sensitivity of cytology for CIN2		53.20%
		Sensitivity of cytology for CIN3		57.70%
		Sensitivity of cytology for ICC		57.70%
		Specificity of cytology		73.40%
	HPV testing
		Sensitivity of pooled hrHPV testing for CIN2		86.40%
		Sensitivity of pooled hrHPV testing for CIN3		89.90%
		Sensitivity of pooled hrHPV testing for ICC		89.90%
		Specificity of pooled hrHPV testing		62.70%
		Sensitivity of genotyping 16/18 for CIN2		43.60%
		Sensitivity of genotyping 16/18 for CIN3		53.40%
		Sensitivity of genotyping 16/18 for ICC		59.20%
		Specificity of genotyping 16/18		91.90%
	Dual staining (pooled HPV triage)
		Sensitivity for CIN2		86.80%
		Sensitivity for CIN3		89.80%
		Specificity for CIN2+		71.40%
		Sensitivity for ICC		93.80%
Epidemiology data [1124262728293031]
	Prevalence of hrHPV			5.6%
	Prevalence of HPV16 and 18			1.7%
	Prevalence of CIN1			0.6%
	Prevalence of CIN2			0.3%
	Prevalence of CIN3			0.8%
	Prevalence of invasive cervical cancer			0.075%
	% of HSIL+ population that is HPV+			88.4%
	% of LSIL population that is HPV+			61.5%
	% of ASCUS population that is HPV+			21.4%
	% CIN1 that are hrHPV 16/18			13.6%
	% CIN2 that are hrHPV 16/18			23.1%
	% CIN3 that are hrHPV 16/18			50.3%
	% of ICC that are hrHPV 16/18			75.0%
	General population annual death rate			0.800%
Natural history parameters
	Progression [17263233343536373839]
		Well to hrHPV infection		3.20%
		Transformation from hrHPV (12 types)
			to CIN1	9.10%
			to CIN2	0.10%
			to CIN3	0.10%
		Transformation from hrHPV 16/18
			to CIN1	7.30%
			to CIN2	2.20%
			to CIN3	2.00%
		Progression from CIN1
			to CIN2	3.10%
			to CIN3	0.90%
		Progression from CIN2 (base case assumes CIN2 does not progress directly to ICC)
			to CIN3	4.20%
			to ICC	0.00%
		CIN3 to ICC		4.50%
		Annual mortality rate for cervical cancer		8.30%
	Regression [32353940]
		Regression from hrHPV (12 types) to..
			with NORMAL smear to well	58.60%
			with BORDERLINE/MILD smear to well	45.60%
		Regression from hrHPV 16/18 to..
			with NORMAL smear to well	43.80%
			with BORDERLINE/MILD smear to well	21.80%
		Regression from CIN1
			to well	21.20%
			to hrHPV	2.40%
		Regression from CIN2
			to well	9.40%
			to CIN1	9.40%
		Regression from CIN3
			to well	3.80%
			to CIN1	1.60%

ASCUS, atypical squamous cells of undetermined significance; CIN, cervical intraepithelial neoplasia; HPV, human papillomavirus; hrHPV, high-risk human papillomavirus; ICC, invasive cervical cancer.

Fig. 1

A model of the patient flow through a cervical cancer screening.

CIN, cervical intraepithelial neoplasia; FU, follow up; HPV, human papillomavirus.

2. Study population and time horizon

The study population is a cohort of women 30–65 years of age presenting for cervical screening. According to our estimations, the total population in Thailand (in 2015) is 67,959,360 [8]. Approximately 26.6% of those are women aged 30–65 years and only 50% attend to screening program [910]. After excluding ineligible individuals, 10.2% pregnant women and 12% women with HIV or hysterectomized women, the total number of the population eligible for screening is 7,953,963 [11]. The age of screening population and interval to rescreening were based on the current screening program of National Health Security Office and Ministry of Public Health of Thailand and also aligned with our previous studies [212].

3. Screening scenarios

The current practice in Thailand, conventional cytology method, provided the baseline reference for comparison with a new comparative practice, HPV DNA test with genotyping 16&18 triage with p16/Ki-67.

1. Cytology: Screening by a cytology-based program using Papanicolaou standard cytology followed by colposcopy of the result is atypical squamous cells of undetermined significance (ASCUS) or worse. The women with negative cytology return to screening in 5 years (Fig. 2A).

Fig. 2

Screening model. (A) Cytology: screening with conventional cytology. (B) HPV/dual stain: HPV DNA test with genotyping 16&18 plus triage with p16/Ki-67 dual stain.

ASCUS, atypical squamous cells of undetermined significance; CIN, cervical intraepithelial neoplasia; HPV, human papillomavirus; ICC, invasive cervical cancer.

*Women with negative colposcopy return to routine screening, women with CIN or ICC were referred to treatment.

2. HPV/dual stain: HPV 16/18 genotyping is used as primary screening then refer to colposcopy if the result is positive for HPV 16 or 18. The dual staining was performed in cases of other 12 high risks HPV positive and those with positive result undergo colposcopy. Women with negative result of HPV genotyping return to screening in 5 years (Fig. 2B).

Compliance, treatment, and follow up protocol

The proportions of attendance at re-testing for management of both algorithms was 61.9%. We used the same compliance rate for management of both algorithms, thus it did not affect the result. The model assumed that all CIN1 cases had followed up visits without treatment, all CIN3 cases and 50% of CIN2 cases received treatment. Women with CIN1 or posttreatment CIN2+ returned to follow up every 6 months. A number of negative follow up visits before returning to routine screening population was 2, 4 times for CIN1 and CIN2-ICC woman respectively [13]. The recurrent rate after treatment for CIN2, CIN3 was assumed to be 4% and for ICC was assumed to be 20% [14]. For women treated with no recurrence within one routine screening cycle, they would re-enter to next screening interval or may loss follow up and return unscheduled. For women treated with recurrence, they would stay in the current status or progress to more severe states in the model.

4. Model parameters and assumptions

Clinical input

All model parameters were shown in Table 1. Demographic data were taken from national data for Thailand [15]. The screening performance inputs were taken from the Addressing THE Need for Advanced HPV Diagnostics HPV: high-risk Human Papillomavirus (ATHENA) trial and PALMS trial [316]. This model did not take into account of HPV vaccination. Colposcopy was considered to be a gold standard for diagnosis of CIN or ICC with 100% sensitivity and specificity. Disease progression and the proportion of fast-growing cancers were assumed not to be age-specific. The risk of progression and regression was assumed to be constant over time. Average probabilities of progression and regression of CIN were stratified by HPV status and showed in Table 1. The QoL or burden of disease was presented by the quality-adjusted life year or QALY. The QALY method could estimate the number of years lived and the QoL during those years that can be attributed to an intervention. Each year in perfect health was assigned the value of 1.0 down to a value of 0.0 for being dead. The QoL weight at the stage of ICC was the weighted average of cervical cancer at different stages (localized cervical cancer was 0.76; regional cervical cancer was 0.67, and distant cervical cancer is 0.48) (Supplementary Table 1) [1718].

Cost input

The total direct medical costs were calculated from the cost of screening, diagnosis and treatment cost in provider perspective (Table 2). Screening costs of cytology including lab fee, profession fee, diagnostic costs and treatment cost for CIN and ICC were derived from our previous studies [219] which based on the cost of Center of Health Assurance at King Chulalongkorn Memorial Hospital. Cost of HPV DNA test and p16/Ki-67 were obtained from Roche Diagnostic (Thailand). As the treatment options for ICC depend on the International Federation of Gynecology and Obstetrics (FIGO) staging system, the cost of treatment was the weighted average of cervical cancer at different stages. Primary surgery is the treatment option for FIGO stage I to IIA cervical cancer. Concurrent chemotherapy and radiotherapy is the treatment of stage IIB to IVA cervical cancer. All cost and clinical input were discounted at an annual rate 3.5% [20].

Table 2

Details of total direct medical costs

Cost parameters			Input value (USD*)
Screening costs [41]
	Office visit (routine/repeat screening)		2.00
	Cytology test (lab fee)		5.30
	Cytology test (professional fee)		3.00
	HPV DNA test		17.00
	P16/Ki-67 Dual staining		35.00
Diagnosis costs [41]
	Office visit (diagnostic follow-up)		12.86
	Colposcopy plus biopsy		21.42
	CINtec® p16 Histology		25.37
Treatment costs [414243]
	Treatment for CIN2/CIN3		1,292.00
	Treatment for ICC†		7,403.00
		• Stage IA1	1,206.29
		• Stage IA2–IIA	2,904.94
		• Stage IIB–IVA	163,334.63
		• Stage IVB	9,168.74
	End of life cancer treatment cost		10,019.00
Discounting rate [20]
	Discount rate for cost		0.035
	Discount rate for health outcomes		0.035

CIN, cervical intraepithelial neoplasia; HPV, human papillomavirus; ICC, invasive cervical cancer.

*The currency used was US dollar (US Dollar exchange rate on May 3,2018; 1 USD = 35 THB); †The treatment for invasive cervical cancer cost was the weighted average of cervical cancer at different stages (stage I, 0.37; stage II, 0.19; stage III,0.33; and stage IV, 0.11) [43].

5. Sensitivity analysis

We performed a one-way sensitivity analysis to estimate the impact of uncertainty in different parameters of cost inputs. Because of the skewness in cost distribution, the ranges of costs were varied 10% below and 3 times above the base-case estimation.

RESULTS

1. Base case analysis

The screening performance, costs, average QALY and cost-effectiveness results based on a time horizon of 10 years were showed in Table 3. A proposed algorithm with HPV test as primary screening combined with p16/Ki-67 dual stain cytology showed better performance compared to Pap cytology by increasing the sensitivity by 48% (82.65% vs. 55.85%, respectively) without losing specificity (96.47 vs. 95.48%, respectively). This strategy was also lower the false positive by 22% results in less number of patients referred to colposcopy. Due to the gain in accuracy, the higher detection rate of precancerous cases leads to higher number of treated cases and reduce the number of cancer cases detected by 18.8%. Considering only the cost of treatment for invasive cervical cancer, the avoided treatment cost of the benefit of programmed as a consequence of reduced use of healthcare services is 54,189,960 USD (7320 cases × 7,403 USD). Moreover, the long-term outcome of the model showed that using HPV testing with dual stain cytology could decrease the prevalence of precancerous lesion (CIN2, CIN3), annual incidence and mortality of cervical cancer (Fig. 3). Under base case assumptions, the total and annual cost of the HPV testing with dual stain cytology model costed up to 71% more than cytology method. (total cost: $1,326,269,261 vs. $771,325,070, annual cost: $13,262,693 vs $7,713,251). The most influential direct cost component was the cost of screening (Supplementary Fig. 1). Using HPV testing with dual stain cytology for primary screening cervical cancer provided a more effective option at an ICER of $1,395 per QALY.

Table 3

Screening performance, cost, average QALY and ICER per QALY gained

Variables			Cytology	HPV/dual stain
Screening performance based on colposcopy population
	Number of screening cycles		20	20
	Total colposcopy population (per screening cycle)		270,487	220,535	−18.5%
	False positive (per screening cycle)		249,800	194,562	−22.1%
	False negative (per screening cycle)		16,350	5,452	−66.7%
	Sensitivity (≥CIN2)		55.85%	82.65%	48.0%
	Specificity (≥CIN2)		95.48%	96.47%	1.0%
Screening performance and total number of cancer/precancer cases detected
	Screening performance (%)
		Cervical cancer detected	57.7%	88.9%	54.1%
		CIN3 detected	57.7%	85.2%	47.7%
		CIN2 detected	53.2%	79.2%	49.0%
	Total number of cancer/precancer cases detected
		Cervical cancer detected	38,927	31,607	−18.81%
		CIN3 Detected	213,218	257,188	20.62%
		CIN2 Detected	161,582	230,669	42.76%
Total annual cost, average QALY and ICER per QALY gained
	Total cost		$771,325,070	$1,326,269,261
	Annual cost		$7,713,251	$13,262,693
	Per person per year (over total screening population)		$1	$2
	Per member per year (over total population)		$0	$0
	Average QALY		23.98	24.03
	ICER per QALY gained		-	$1,395

CIN, cervical intraepithelial neoplasia; HPV, human papillomavirus; ICER, incremental cost-effectiveness ratio; QALY, quality-adjusted life year.

Fig. 3

The prevalence of preinvasive cervical cancer, cervical cancer, and mortality rate from cancer.

CIN, cervical intraepithelial neoplasia; HPV, human papillomavirus.

2. Sensitivity analysis

One-way sensitivity analysis for ICER was performed by changing the parameter of cost inputs. The values for sensitivity analysis and results were presented as tornado diagram (Supplementary Table 2, Supplementary Fig. 2). The analysis showed that the costs of HPV and cytology had the greatest impact on ICER. The ICER fell to $303/QALY gained at a cytology cost of $15.9 and increased to $4,970/QALY gained at a HPV testing cost of $51. No other parameters changed the ICER by more than 25%. All of the variables had no impact on model conclusions.

DISCUSSION

The superior accuracy of HPV primary testing over Pap screening is well established [21]. Our earlier cost-effectiveness study also supported full scale implementation of HPV testing as a primary cervical cancer screening in Thailand [2]. In this study, we found that using dual-stained cytology as triage part of HPV algorithm could improve screening performance. The gaining in sensitivity without decreasing specificity is consistent with the results of Uijterwaal confirming that dual stain cytology is suitable for triage women with HPV positive with normal cytology to colposcopy [5]. Wentzensen et al. [6] compared the detection rate of CIN2+ cases in HPV positive women with different triage strategies; between PAP smear and dual stain cytology. The result revealed that the dual stain had similar sensitivity and statistically higher specificity compared with cytology. As a result of improvement in diagnostic accuracy for CIN2+ of HPV screening with dual stain triage algorithm, the number of women with precancerous detected and treated in the earlier round was higher. The prevalence and mortality rate of cancer case is lower. The number of women who return at next screening interval in this strategy might be more than the conventional cytology strategy. Subsequently, the cost of screening, treatment and follow up might be increased. The diagnostic cost calculated from the cost of colposcopy plus biopsy was slightly fewer than cytology method because of less patients referred to colposcopy. In addition, the new algorithm consisted of multiple steps and used more expensive tests. Therefore, the screening cost was much higher than other costs and had the most impact on the total cost.

The HPV and dual stain model was more effective and more costly than conventional cytology method with ICERs was 1,395 per QALY. All of the results including clinical outcome, QoL and budget impact are in line with cost-effectiveness analysis study in Belgium, conducted by Tjalma et al. [22]. The author compared the 2 strategies, primary cytology (the standard of care in Belgium) and the HPV primary screening and triage with dual stain cytology, same as our study. The result showed the latter strategy could decrease the prevalence of CIN and cervical cancer significantly and reduce the screening budget by 21%, resulting in saving 5.3 million euro a year in Belgium. The ICERs per QALY gained for the groups age 25–65 years and 30–65 years were -260 and -302 euro, respectively. However, the population cohort, time horizon, demographic data, cost data, and gross domestic product (GDP), are different. The population in model of Tjalma's study [22] has three different age groups for cervical cancer screening: women between ages 20–24, women between ages 25–65 and women between ages 30–65. The screening intervals were 3 years for primary cytology strategy and 5 years for HPV primary with reflex dual stain cytology.

Due to uncertain cost data, the sensitivity analysis was performed and found that the cost of cytology and HPV testing had an impact to ICER result. If HPV test cost was lower than base case 10%, it would decrease ICER by 12.7%. On the other hand, the ICER might rise to $4,970 per QALY gained if the cost of HPV test increased at least three times. Nevertheless, it was still more cost-effective than cytology method. The effect of cost variation in the HPV testing was similar to the previous study compared between cytology and HPV testing with cytology triage in Thailand [2]. In this study, the cost of HPV primary screening and dual stain does not play as the major parameter, even the cost of HPV test increase 3 times. HPV screening and dual stain triage still showed benefit over the cytology. We assumed that an ICER of less than three times the per capita GDP would be considered cost effective (current per capita GDP of Thailand $5,901) [23].

As far as we know, this is the first study using p16/Ki-67 dual stain cytology as triage part of HPV primary screening in Asia. We reported both intermediate and long-term clinical outcome, budget impact and cost-effectiveness data. Most of the data input derived from studies conducted in our country or in the setting that were similar to our practice. The results of study may differ by nations. However, the sensitivity analysis was conducted to display the effect of result due to variation in the cost of the screening tests. There are some limitations that need to be discussed further. This model did not take into account the impact of HPV vaccination which could affect the accuracy of HPV testing and prevalence of patients with precancerous lesions. According to the appropriate cervical cancer screening among vaccinated women, HPV test revealed more accurate than cytology [19]. Then the ICER would probably be decrease. We assumed that the rate of cervical lesion regression or progression were the same for all age groups and the treatment costs of cervical cancer in this model were the average costs across all stage. In the real practice, patients with advanced stage cervical cancer received concurrent chemoradiation which its cost is more expensive than early stage cancer. It is reasonable to assume that using difference cost would not change the result or may reduce ICER due to early detection cancer case and decreasing advance stage in the arm of new strategy (HPV/dual stain triage).

To find the best strategy, the future studies include new marker or triage with different tests are needed and the impact of HPV vaccination should be considered in the model. Moreover, our study was designed to focus on only direct medical cost this time. The effect of non-direct medical cost of social viewpoint is an important and should be included for future research.

This result is in line with our earlier study on health economic model with Thai context, primary screening with HPV-DNA (with partial genotyping 16&18 and 12 other high-risk HPV) demonstrated cost-effective compare to cytology screening. The p16/Ki-67 dual stain cytology could improve screening performance as a part of HPV algorithm. The new proposed HPV primary strategy allows early detection of cervical precancerous and cancer cases, reduction in mortality and lower treatment costs. We presented new information that should be considered to future national screening policy.