Early Cost Effectiveness of Whole-Genome Sequencing as a Clinical Diagnostic Test for Patients with Inoperable Stage IIIB,C/IV Non-squamous Non-small-Cell Lung Cancer.

BACKGROUND: Advanced non-small-cell lung cancer (NSCLC) harbours many genetic aberrations that can be targeted with systemic treatments. Whole-genome sequencing (WGS) can simultaneously detect these (and possibly new) molecular targets. However, the exact added clinical value of WGS is unknown.
OBJECTIVE: The objective of this study was to determine the early cost effectiveness of using WGS in diagnostic strategies compared with currently used molecular diagnostics for patients with inoperable stage IIIB,C/IV non-squamous NSCLC from a Dutch healthcare perspective.
METHODS: A decision tree represented the diagnostic pathway, and a cohort state transition model represented disease progression. Three diagnostic strategies were modelled: standard of care (SoC) alone, WGS as a diagnostic test, and SoC followed by WGS. Treatment effectiveness was based on a systematic review. Probabilistic cost-effectiveness analyses were performed, and threshold analyses (using €80,000 per quality-adjusted life-year [QALY]) was used to explore the early cost effectiveness of WGS.
RESULTS: WGS as a diagnostic test resulted in more QALYs (0.002) and costs (€1534 [incremental net monetary benefit -€1349]), and SoC followed by WGS resulted in fewer QALYs (-0.002) and more costs (€1059 [-€1194]) compared with SoC alone. WGS as a diagnostic test was only cost effective if it was priced at €2000 per patient and identified 2.7% more actionable patients than SoC alone. Treating these additional identified patients with new treatments costing >€4069 per month decreased the probability of cost effectiveness.
CONCLUSIONS: Our analysis suggests that providing WGS as a diagnostic test is cost effective compared with SoC followed by WGS and SoC alone if costs for WGS decrease and additional patients with actionable targets are identified. This cost-effectiveness model can be used to incorporate new findings iteratively and to support ongoing decision making regarding the use of WGS in this rapidly evolving field.

Key Points for Decision Makers

Introduction

Advanced non-small-cell lung cancer (NSCLC) harbours many genetic aberrations [1-4] that can be targeted with systemic treatments [4, 5]. Whole-genome sequencing (WGS) is being investigated as a possible molecular diagnostic test to identify druggable targets for NSCLC, given its ability to simultaneously test for all known, and potentially new, molecular targets [6]. However, its added clinical value as a molecular diagnostic test is currently unknown.

Many molecular diagnostic tests exist that can detect a single actionable target, for example, immunohistochemistry (IHC) and fluorescence in situ hybridization (FISH) [2, 7]. These tests are often performed sequentially, though this is not ideal given the time required and limited available tissue samples [8, 9]. Next-generation sequencing (NGS), e.g. multi-gene panels, can overcome this problem by simultaneously testing for multiple actionable targets. However, these tests are often restricted to a specific subset of targets and performed in different combinations with other tests according to different hospitals [3, 10, 11]. WGS is a form of NGS that provides information about the whole genome and thus can detect any potential aberration present in the tumour DNA, including those that other techniques focussed only on specific parts of the genome would have missed [6, 12]. However, WGS is expensive, and interpretation of results is complex [11].

The clinical utility of a sequencing technology lies in its ability to accurately identify individuals with a given disease and to aid in personalised and targeted treatment selection. Targeted treatments have demonstrated significant improvement of progression-free survival (PFS) and overall survival (OS) for some patients with NSCLC [13]. Many actionable targets can be tested for (e.g. epidermal growth factor receptor [EGFR], anaplastic lymphoma kinase [ALK], ROS proto-oncogene 1 [ROS1], B-Raf proto-oncogene [BRAF]), and more are still being investigated [4]. Additionally, some aberrations validated in other tumour types (e.g. phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit alpha [PIK3CA], human epidermal growth factor receptor 2 [HER2], and BRCA) also exist in NSCLC but are not yet targetable [3]. With the rapid pace of genomic discoveries and treatment developments, there is a chance that currently used multi-gene panels aimed at selected genes could miss novel targets [14]. Furthermore, adding new targets to panels takes time, and they must be validated. These arguments imply that there is a potential use for WGS as a molecular diagnostic test. However, it is unclear how many actionable targets can be found and how much health gain using WGS would yield over traditional methods. At the same time, the costs for WGS are high but decreasing rapidly, from about €3000 per patient in the Netherlands in 2018 [15] to somewhere between €2500 and €2000 in 2021 (Cuppen E, personal communication). In essence, it is unclear how WGS should be implemented in the diagnostic pathway to provide the most health benefits at an accepted cost [16]. Early cost-effectiveness analysis is a useful tool to address this, as it determines the prerequisites for cost effectiveness and the potential for adoption of a new technology using currently available evidence [17]. To date, only limited evidence of the cost effectiveness of introducing multi-gene panels in molecular diagnostics exists [10, 18, 19]. Only one full economic evaluation of the use of WGS in oncology was identified, and this was specific to incidental findings in colorectal cancer [20]. Economic evaluations of WGS are urgently needed to support its translation and implementation into clinical practice in oncology [19]. Therefore, the objective of this study was to determine the early cost effectiveness of using WGS in diagnostic strategies compared with currently used molecular diagnostics for patients with inoperable stage IIIB,C/IV non-squamous NSCLC, from a Dutch healthcare perspective.

Methods

Overview

A probabilistic decision model was constructed in Microsoft Excel consisting of a decision tree and a state transition model (STM). A societal perspective from the Netherlands was adopted, and the model simulated a cohort with a starting age of 60 years, in 240 one-month cycles over a lifetime horizon.

Diagnostic Strategies

The decision tree is displayed in Fig. 1. Diagnostic strategy A: ‘standard of care’ (SoC) represents the current clinical practice of molecular testing, consisting of an NGS multi-gene panel, FISH, IHC, and Archer fusionPlex, to test for actionable targets. In cases where no actionable target was found, IHC was performed to determine programmed death-ligand 1 (PD-L1) expression levels on the tumour cell’s membrane, which are used in the selection of immunotherapy treatment. Using these tests, the most optimal and least costly testing strategy was created that detected all common actionable targets. This was based on guidelines and an expert’s opinion (PR) [21]. In diagnostic strategy B: ‘WGS as a diagnostic test’, WGS replaced all SoC tests and was used alongside IHC to test for PD-L1 (since WGS cannot detect PD-L1). In diagnostic strategy C: ‘SoC followed by WGS’, patients were first tested with SoC diagnostics and then with WGS if no targets were found. Biopsy and technical success rates were included in all diagnostic strategies, referring to biopsies containing enough tumour cells and the technical performance of the diagnostic test (Fig. 1, Table 1). Where no results were available, patients were considered the ‘target unknown’ subgroup. More details on diagnostics are presented in the Appendix in the electronic supplementary material (ESM).

Fig. 1

Diagnostic strategies in inoperable non-small-cell lung cancer. If a biopsy contained too few tumour cells for WGS as a diagnostic test (strategy B), patients could still receive SoC as in strategy A. If a biopsy also contained too few tumour cells for SoC, no test could be performed. If a diagnostic test technically failed twice, no results were available. In case of an unknown target, patients were treated as in the ‘PD-L1 unselected’ subgroup. The circles containing ‘A1’, ‘B1’, and ‘C1’ are truncations of the fully unfolded branch that is only presented the first time it occurs in the decision tree. Target X is a summary of the biomarkers that are currently missed in routine SoC (because of technical and/or analysis challenges) but are reliably detected by WGS. ALK anaplastic lymphoma kinase, BRAF B-Raf proto-oncogene, EGFR epidermal growth factor receptor, IHC immunohistochemistry, KRAS Kirsten rat sarcoma virus proto-oncogene, MET MET proto-oncogene, NSCLC non-small-cell lung cancer, NTRK neurotrophic tropomyosin receptor kinase, PDCT platinum-doublet chemotherapy, PD-L1 programmed death-ligand 1, RET rearranged during transfection proto-oncogene, ROS1 ROS proto-oncogene 1, SoC standard of care, WGS whole-genome sequencing

Table 1

Model input parameters

Parameter	Base-case value	SE	Distribution	Sources
Model population
Sex, male	0.56	–	Fixed	[22]
Age, mean	60.00	–	Fixed	[22]
Proportions actionable subgroups
EGFR (exon 19/21/T790M)	0.070	0.008	Dirichlet	[23]
EGFR (non-classic/other)	0.010	0.001	Dirichlet	[23]
ALK	0.020	0.002	Dirichlet	[23]
ROS1	0.019	0.002	Dirichlet	[1]
BRAF (V600)	0.021	0.002	Dirichlet	[24]
MET (amp/exon 14 sk)	0.014	0.001	Dirichlet	[24]
RET	0.017	0.002	Dirichlet	[24]
NTRK (kinase 1, 2, 3)	0.005	0.001	Dirichlet	[2]
Target X (WGS)a	0.005	0.001	Dirichlet	EO
KRAS (exon 2,3,4)	0.280	0.029	Dirichlet	[23]
PD-L1 ≥ 50%b	0.250	0.038c	Beta	[25]
Success rates
Technical (SoC)d	0.943	0.014	Beta	[26]
Biopsy (SoC)	0.950	0.014	Beta	[26]
Technical (WGS)d	0.956	0.010	Beta	[26]
Biopsy (WGS)	0.850	0.027	Beta	[26]
Diagnostic test costs
IHC (PD-L1)	97	–	Fixed	[15]
SoC diagnosticse	850	–	Fixed	[15]
WGSf	2500	–	Fixed	Cuppeng
Unit price per model cycle (€)h
Treatments
Afatinib	2494	–	Fixed	[27, 28]
Osimertinib	6746	–	Fixed	[27, 28]
Crizotinib	5605	–	Fixed	[27, 28]
Alectinib	6518	–	Fixed	[27, 28]
Lorlatinib	6762	–	Fixed	[27, 28]
Dabrafenib	6921	–	Fixed	[27, 28]
Trametinib	6762	–	Fixed	[27, 28]
Larotrectinib	5916	–	Fixed	[27, 28]
‘Treatment X’i	10,000	1500c	Gamma	Assumption
Pembrolizumab	8292	–	Fixed	[27, 28]
Cisplatin	83	–	Fixed	[27, 28]
Pemetrexed	3546	–	Fixed	[27, 28]
Docetaxel	962	–	Fixed	[27, 28]
Best supportive care	1845	277c	Gamma	[29]
Treatments for SAEsj	Table 2 (ESM)	–	–	–
Direct medical costs
Drug administration	282	42c	Gamma	[29]
End of life	2282	342c	Gamma	[30]
Laboratory test	80	12c	Gamma	[31]
Tumour response assessment (CT/MRI)	421	63c	Gamma	[32]
Outpatient visits	87	13c	Gamma	[33]
Informal care (per hour)	15	2c	Gamma	[33]
Home care (per hour)	25	4c	Gamma	[33]
Travel + (parking)	6	1c	Gamma	[33]
Indirect costs
Productivityk	4228	634c	Gamma	[34]
Indirect medicall	Table 6 (ESM)	–	–	[35]
Utilitiesm
No progression	0.710	0.033	Beta	[36]
Progression first line	0.670	0.041	Beta	[36]
Progression second line	0.590	0.089	Beta	[36]
Disutilitiesn	Table 2 (ESM)	–	–	–

Costs and effects were discounted by a rate of 4.0 and 1.5%, respectively, in line with Dutch guidelines [33] All costs were based on the average consumer price index of 2020

ALK anaplastic lymphoma kinase, amp amplification, BRAF B-Raf proto-oncogene, CT computed tomography, EGFR epidermal growth factor receptor, EO expert opinion, ESM electronic supplementary material, IHC immunohistochemistry, KRAS Kirsten rat sarcoma virus proto-oncogene, MET MET proto-oncogene, MRI magnetic resonance imaging, NTRK neurotrophic tropomyosin receptor kinase, PD-L1 programmed death-ligand 1, RET rearranged during transfection proto-oncogene, ROS1 ROS proto-oncogene 1, SAEs serious adverse events, SE standard error, sk skipping, SoC standard of care, WGS whole-genome sequencing

aIt was assumed that WGS was capable of finding the same proportions of actionable targets as SoC diagnostic tests, with an additional 0.5% of rare molecular targets, i.e. ‘target X’ for on-label treatment (based on expert opinion)

bThe probability that PD-L1 is unselected is 1-p(PD-L1 ≥50%)

cWhen not available, the SE was assumed to be equal to 15% of the mean

dIt was assumed that the success rates for performing the diagnostic tests for the first and second time were equal

eCalculations for the average costs for SoC diagnostics can be found in Fig. 1 in the ESM

fThe cost of WGS was based on the NovaSeq 6000 illumine and concerns tumour-normal pair sequencing at 90× average base coverage for the tumour and 30× for the germline control and includes interpretation of the results by a molecular tumour board

gCuppen E, personal communication

hOne model cycle has a length of 1 month

iThe cost of treatment X was assumed to be €10,000 per model cycle since this would most likely concern new and expensive targeted therapies

jIn presence of SAEs, treatment costs were counted for treating SAEs

kThe friction costs method was used

lIndirect medical costs were age dependent, ranging from €4760 to 32,070 based on an age of 60 to 90 years and were calculated using the PAID tool v3.0.

mIn the absence of SAEs, health utilities were assumed to be treatment independent

nDisutility scores were counted in presence of SAEs, that occurred as treatment side effects

Targets and Treatment Selection

The following targets were considered in the model: EGFR mutations (exon 19, 21, T790M, non-classic, and other); ALK rearrangements; ROS1 rearrangements; BRAF (V600) mutations; MET proto-oncogene (MET) (amplifications and exon 14 skipping); rearranged during transfection proto-oncogene (RET) rearrangement; neurotrophic tropomyosin receptor kinase (NTRK) (kinase 1, 2, and 3) gene fusions; KRAS mutation in exon 2, 3, and 4; and PD-L1 expression level. Molecular targets in the aforementioned genes with targeted therapies available in routine practice that are currently missed in SoC (because of technical and/or analysis challenges) but that can be reliably detected with WGS were labelled as ‘target X’ [26]. For the base-case analysis, WGS was assumed to be capable of identifying a proportion of 0.5% target X, based on expert opinion (PR). SoC tests were unable to detect target X because either the test design or the testing strategy was incomplete. All targets in the model were mutually exclusive.

The treatment strategies that were used in the model were a close representation of Dutch clinical practice and were determined through consultation with clinical oncology experts (JM, MP, JA) and the Dutch guidelines [21]. First-line treatments consisted of targeted therapies for the following actionable targets: EGFR, ALK, ROS1, BRAF, and NTRK (including target X). It was assumed that patients with target X were treated with a targeted therapy, summarised as ‘treatment X’. This treatment X was assumed to have the same effectiveness as first-line chemotherapy in patients with EGFR mutations and a cost of €10,000 per patient per model cycle. If no target or a MET, RET, or KRAS aberration was found, first-line treatment was based on PD-L1 and consisted of immunotherapy with or without platinum-doublet chemotherapy (PDCT) in accordance with clinical guidelines [21]. Patients in the unknown target subgroup were treated as PD-L1 unselected. It was assumed that patients switched to second-line treatment after disease progression in the first line. In the second line, patients could receive targeted therapies, immunotherapies with or without PDCT, or chemotherapy, according to the treatment they received in the first line. After disease progression in the second line, patients received best supportive care, which concerned symptom management. Lower detection rates of molecular actionable targets resulted in more patients treated with immunotherapies [21]. Immunotherapy has higher costs and is less effective than targeted therapy. The effectiveness and costs of treatment X will determine the cost effectiveness of WGS, as more patients detected with target X will lead to more patients receiving targeted therapy. Details of the treatment strategies can be found in the Appendix in the ESM.

Model Structure for Disease Progression

Patients’ disease progression was simulated in the STM with four health states: no progression, progression first line, progression second line, and death (Fig. 2).

Fig. 2

State transition model structure. The transition probabilities between the health states are determined by the treatment that was assigned in the decision tree. Note that p4 and p5 were assumed to be equal. p1 probability of progression in the first line (of treatment administration), p2 probability of dying in the first line, p3 probability of progression in the second line, p4 probability of dying in the second line, p5 probability of dying with best supportive care

Model Input

Identified Targets and Treatment Effects

Parameter values for the decision tree and STM are listed in Table 1. Results from a recent systematic review were used to model the OS and PFS of targeted therapies and immunotherapies in NSCLC [13, 37–51]. Treatment effects that were not available in the systematic review were replaced with the OS and PFS of chemotherapy of patients with EGFR mutations. Transition probabilities for first-line treatments were derived from parametric survival models using different distributions chosen according to goodness of fit. Transition probabilities for second-line treatments were derived using exponential distributions, with the assumption of no cycle-dependent transition probabilities in the model. Additionally, the background mortality derived from the Dutch population was applied. More details on model inputs are presented in the Appendix in the ESM.

Costs

Costs and resource use were based on the Dutch costing manual [33], the Dutch healthcare authority [32], and the literature [27]. In the presence of serious adverse events (SAEs) requiring treatment, treatment costs were accounted for. Costs and effects were discounted by a rate of 4.0 and 1.5%, respectively, in line with Dutch guidelines [33]. All costs were based on the average consumer price index of 2020.

Utilities

Health state utilities for the states no progression and progression in the first and second line were obtained from the literature [36]. In the absence of SAEs, health utilities for the three states were assumed to be treatment independent. Disutilities were counted in the presence of treatment-related SAEs. It was assumed that SAEs occurred as soon as the treatment started and disappeared when treated. More details on (dis)utility weights, the studies they were derived from, and SAEs can be found in the Appendix in the ESM.

Model Validation

The model was checked for face validity by consulting clinical oncology and pathology experts (JM, MP, JA, RB, PR, KM), as were all model assumptions. Furthermore, the model was checked for validity using the AdViSHE (Assessment of the Validation Status of Health-Economic decision models) checklist [52] by one of the authors (ZM). The completed checklist is included in the Appendix in the ESM.

Analyses

Model Output

Model outcomes included costs in €, year 2020 values; life-years (LYs); and quality-adjusted LYs (QALYs) per actionable subgroup. Outcomes were averaged per diagnostic strategy. The incremental cost-effectiveness ratio (ICER) and net monetary benefit were calculated using a willingness-to-pay (WTP) threshold of €80,000 per QALY, as advised for this patient population in the Netherlands [53].

Sensitivity and Threshold Analyses

Deterministic sensitivity analysis (DSA) was performed to test the influence of different parameters on model outcomes. Parameter values were varied based on their standard error or 15% of the mean if not available. Specific parameter values were used for the upper and lower boundaries of the cost for WGS (€2000 [Cuppen E, personal communication] and €3040 [15]) and for the upper boundary of the cost for SoC: €1881 [54]. Probabilistic sensitivity analyses (PSA) were performed to evaluate the parameter uncertainty of all model inputs simultaneously. This was done by drawing random values from pre-specified distributions using Monte Carlo simulations with 2000 iterations. As only the list prices of systemic treatments were available, we performed a scenario analysis in which all treatment costs were halved to reflect the impact of price negotiations between government/hospitals and pharmaceutical companies. One- and three-way threshold analyses were performed to estimate the thresholds for costs for WGS; the proportion of patients with target X identified with WGS; the effectiveness of treatment X (applying the hazard ratio [HR] to the base case); the costs for treatment X; and the costs for SoC for WGS to become cost effective, considering the WTP threshold. A comparison of the strategies ‘WGS as diagnostic test’ and ‘SoC followed by WGS’ versus ‘SoC alone’ was performed in the threshold analyses.

Results

Base-Case Analysis

Table 2 lists the intermediate results of the base-case and scenario sensitivity analysis on the treatment costs and QALYs of the different treatment strategies for the actionable subgroups. Each diagnostic strategy resulted in different proportions of patients within the actionable subgroups. The actionable subgroup resulted in different LYs, QALYs, and costs, which were model averaged per diagnostic strategy. Table 3 lists the resulting overall results of the base-case analysis. SoC followed by WGS was dominated by SoC. WGS as a diagnostic test was not cost effective compared with SoC, considering the WTP threshold. Additional results are listed in Table 8 in the ESM.

Table 2

Treatment costs and quality-adjusted life-years of all treatment strategies and proportions of patients per actionable subgroup for each diagnostic strategy

Actionable subgroups	Treatment strategiesa	QALY (mean)	Total treatment cost		Proportion of patients per diagnostic strategy (%)
			Base case	SSAb	SoC	WGS as diagnostic test	SoC followed by WGS
ALK	Alectinib	2.52	224,007	126,552	1.89	1.90	1.89
EGFR	Osimertinib	2.19	170,197	103,011	6.64	6.64	6.64
EGFR	Afatinib	1.82	79,062	57,170	0.94	0.94	0.94
ROS1	Crizotinib	1.67	93,261	66,964	1.80	1.80	1.80
NTRK	Larotrectinib	1.67	95,262	67,966	0.47	0.47	0.47
Target Xc	Treatment X	1.67	121,461	81,064	0.00	0.43	0.23
BRAF	Dabrafenib + trametinib	1.67	184,536	113,509	1.99	1.99	1.99
PD-L1 ≥ 50%	Pembrolizumab	1.56	114,083	74,615	5.06	5.04	4.75
PD-L1 ≥ 50%	Pembrolizumab + PDCT	1.27	147,452	82,523	15.18	15.12	14.25
PD-L1 unselected	Pembrolizumab + PDCT	1.08	126,830	71,909	66.03	65.67	67.04

Treatment strategies are ranked by effectiveness, i.e. mean total QALYs produced. The PD-L1 subgroups include patients with MET, RET, KRAS and no target. Costs are presented as €

ALK anaplastic lymphoma kinase, BRAF B-Raf proto-oncogene, EGFR epidermal growth factor receptor, KRAS Kirsten rat sarcoma virus proto-oncogene, MET MET proto-oncogene, NTRK neurotrophic tropomyosin receptor kinase, PDCT platinum-doublet chemotherapy, PD-L1 programmed death-ligand 1, QALY quality-adjusted life-year, RET rearranged during transfection proto-oncogene, ROS1 ROS proto-oncogene 1, SoC standard of care, SSA scenario sensitivity analysis, WGS whole-genome sequencing

aThe different treatment strategies are described by the treatment administered in the first line, but costs and QALYs are produced by both first and second line

bScenario sensitivity analysis in which systemic treatments cost were halved as result of price negotiations

c‘Target X’ is a summary of the biomarkers that are currently missed in routine SoC (because of technical and/or analysis challenges) but are reliably detected by WGS

Table 3

Results of the base-case cost-effectiveness analyses

Diagnostic strategy	LY (95% CI)	QALYs (95% CI)	Costs (95% CI)	Strategy comparison	IncrementalLY (95% CI)	Incremental QALYs (95% CI)	Incremental costs (95% CI)	ICERa	iNMBb
SoC	1.878 (1.755 to 2.011)	1.234 (1.077 to 1.400)	149,703 (141,726 to 158,225)	–	–	–	–	–	–
SoC followed by WGS	1.878 (1.753 to 2.009)	1.232 (1.072 to 1.393)	150,777 (142,730 to 159,324)	vs. SoC	–0.002 (–0.008 to 0.002)	–0.002 (–0.005 to 0.001)	1059 (761 to 1284)	Dominated	−1194
WGS as a diagnostic test	1.882 (1.758 to 2.014)	1.236 (1.079 to 1.405)	151,237 (143,456 to 160,196)	vs. SoC	0.004 (–0.033 to 0.042)	0.002 (–0.022 to 0.027)	1534 (111 to 2929)	657,572	−1349

Diagnostic strategies are sorted by costs. Costs are presented in €

CI confidence interval, ICER incremental cost-effectiveness ratio, iNMB incremental net monetary benefit, LY life-year, QALY quality-adjusted life-year, SoC standard of care, WGS whole-genome sequencing

aA diagnostic strategy was dominated by another diagnostic strategy if it resulted in fewer QALYs and more costs

bA diagnostic strategy was considered cost effective compared with SoC if the iNMB was equal to or above 0, based on a willingness-to-pay threshold of 80,000 per QALY

Sensitivity Analyses

Results of the DSA showed that the model results were robust (Fig. 2 in the ESM). Results of the PSA are displayed in Fig. 3 and in Fig. 3 in the ESM. The 95% confidence intervals of both WGS strategies indicated higher costs and uncertainty regarding the QALY gain compared with SoC.

Fig. 3

Incremental cost-effectiveness plane of the compared diagnostic strategies. The ellipses represent the 95% confidence interval of 2000 iterations of the diagnostic strategy comparisons. The point estimates represent the mean results of the probabilistic sensitivity analysis. The size of the 95% confidence interval ellipses is explained by the differences between diagnostic strategies caused by parameter uncertainty located in the decision tree of the model. WGS as a diagnostic test had the most differences in diagnostic tests and SoC followed by WGS the least compared with SoC alone. The dotted line represents the willingness-to-pay threshold of €80,000 per QALY. QALY quality-adjusted life-years, SoC standard of care, WGS whole-genome sequencing

WGS as a diagnostic test had a 3.0% probability of being cost effective compared with SoC, and SoC followed by WGS had a 0% probability, given a WTP threshold of €80,000 per QALY.

Threshold Analyses

One-way threshold analyses showed that the cost for WGS would need to decrease from €2500 to €998 for WGS as a diagnostic test and to €61 for SoC followed by WGS to become cost effective (vs. SoC alone). The proportion of patients with target X would need to increase from 0.5 to 3.8% for WGS as a diagnostic test and to 5.9% for SoC followed by WGS to become cost effective. When assuming that the effect of treatment X on OS and PFS was equal to that of alectinib (applying an HR of 0.21), the WGS strategies were still not cost effective. To become cost effective in each treatment strategy, costs for treatment X would need to decrease to below zero for WGS as a diagnostic test and for SoC followed by WGS. Furthermore, costs for SoC alone would need to increase from €850 to €2334 for WGS as a diagnostic test and to below zero for SoC followed by WGS to become cost effective.

Figure 4 displays the results of the three-way threshold analysis. Figure 4a shows that, if costs for WGS remained at €2500 but treatment X had a better treatment effect on OS and PFS (applying an HR of 0.50), the proportion of patients with target X would need to increase from 0.5 to 7.4% for WGS as a diagnostic test and to 11.4% for SoC followed by WGS to become cost effective. Figure 4b shows that, if costs for treatment X dropped below €4069, whereas the proportion of patients with target X increased from 0.5 to 2.1% for WGS as a diagnostic test and to 3.3% for SoC followed by WGS, a better treatment effect of treatment X resulted in a higher probability of WGS as a diagnostic test and of SoC followed by WGS becoming cost effective. Figure 4c shows that, if costs for WGS remained at €2500, whereas costs for SoC decreased to €619, the proportion of patients with target X would need to increase from 0.5 to 4.3% for WGS as a diagnostic test and to 5.9% for SoC followed by WGS to become cost effective. If costs for SoC increased to €1881, the proportion of patients with target X would need to increase from 0.5 to 1.5% for WGS as a diagnostic test and to 5.9% for SoC followed by WGS to become cost effective. Figure 4d shows that, if costs for treatment X remained at €10,000, whereas costs for WGS decreased to €2000, the proportion of patients with target X would need to increase from 0.5 to 2.7% for WGS as a diagnostic test and to 4.8% for SoC followed by WGS to become cost effective. If costs for WGS increased to €3040, the proportion of patients with target X would need to increase from 0.5 to 5.0% for WGS as a diagnostic test and to 7.1% for SoC followed by WGS to become cost effective.

Fig. 4

Three-way threshold analysis of the cost-effectiveness of WGS as a diagnostic test and SoC followed by WGS vs. SoC alone. The lines represent the incremental cost-effectiveness ratio of €80,000 per quality-adjusted life-year for the parameter value combinations of the x- and y-axis and a threshold value for a third parameter that is represented by the different lines. Parameter value combinations that are on the other side of the line in contrast to the area that is indicated with ‘SoC preferred’, result in the cost effectiveness of the other strategies, depending on which line is crossed. The following variables were varied simultaneously: a, cost for WGS, proportion of patients with target X found by WGS and two threshold values for treatment effect of treatment X; b, cost for treatment X, proportion of patients with target X found by WGS and two threshold values for treatment effect of treatment X; c, cost for WGS, proportion of patients with target X found by WGS and two threshold values for cost for SoC; d, cost for treatment X, proportion of patients with target X found by WGS and two threshold values for cost for WGS. HR hazard ratio, SoC standard of care, WGS whole-genome sequencing

Discussion

The early cost effectiveness of using WGS in diagnostic strategies compared with SoC to identify targeted therapies for patients with inoperable NSCLC was determined. This model used a Dutch healthcare perspective and was informed by inputs derived from the literature. Our findings show that both diagnostic strategies including WGS were not cost effective compared with most optimal SoC. The most efficient use of WGS appears to be as a diagnostic test for all patients, supported by SoC diagnostics in cases where the biopsy contains too few tumour cells. Using WGS as a diagnostic test resulted in a modest increase in QALYs (0.002) and higher costs (€1534) per patient. This was mainly because of the few additional targets in the model that could be treated with cost-effective targeted therapies. For example, we found that WGS as a diagnostic test could become cost effective if costs decreased to €2000 and 2.7% additional patients with actionable targets for routine targeted therapy were identified (assuming treatment X costs €10,000 per month). However, this is dependent on the cost effectiveness of treatment(s) targeted at the additional targets identified with WGS. If costs for WGS remained at €2500, identifying 2.7% more patients eligible for targeted therapy could result in cost effectiveness if monthly treatment costs to treat these additional patients decreased substantially (31.1%).

It has been suggested that there are more targets available or under investigation in incurable cancers that are not detected by SoC and that can be targeted with systemic treatments [3, 6, 11, 26, 55, 56]. Additionally, WGS can detect non-hotspot targets within known genes, e.g. EGFR, that are missed by SoC and can be treated. However, these data are immature and not yet validated so were not incorporated in current analyses. If the number of actionable targets eventually increases, simultaneous testing using NGS techniques will likely be the most optimal strategy in terms of costs [3, 10, 18], probably because more (cost-)effective targeted therapies become available over time [10]. In many cases, newly found targets have led to off-label treatments or use as experimental drugs in clinical trial settings. This is expensive and not always accessible to all potentially eligible patients [6, 55, 56]. Therefore, finding more actionable targets with WGS does not necessarily lead to a more efficient use of healthcare resources [10], although benefits for patients might increase [6]. Additionally, the ‘added value’ of WGS goes beyond the direct health benefits we captured in our analysis. For example, the scope of using WGS data for treatment decision support could be extended to an earlier disease stage, allowing for previous WGS results to be used for selecting novel treatments that emerge for future patients [57]. The value of WGS could encapsulate wider public benefits, workability, macro-economic value for diagnostic laboratories, and other social factors [57]. However, these values cannot always be easily assessed in health economic evaluations.

A strength of this research is that the diagnostic pathways were defined based on guidelines and extensive clinical expert consultation. We analysed two different positions of WGS in the diagnostic pathway and performed threshold analyses on both strategies. This enabled us to explore many different scenarios and analyse their impact on the early cost effectiveness of WGS. This can be further extended to new scenarios that can help us better understand what to focus future research on to reduce uncertainty. We created a detailed model that can be easily adjusted when new data become available. In addition, we believe that, with adjustments, the model is also capable of simulating other forms of advanced cancers beyond NSCLC.

A limitation of this research is the inability of a cohort STM to properly capture the complexity of personalised medicine. In this context, it is often better to use a patient-level model. However, such a model may have only limited advantages regarding implementing multiple time scales in an efficient way. In this early stage, understanding how to use WGS as diagnostic test should be investigated first. This entails capturing the boundaries of essential parameter values, which can also be done with a cohort model. Because of the limited data, we had to make some model assumptions. The most influential assumption was regarding the treatment effectiveness of some therapies, most notably of treatment X. We also assumed a most optimal SoC strategy in which all the standard routine biomarkers are being tested with a comprehensive multi-gene panel at the DNA and RNA level. If the SoC strategy was not as optimal (which is most likely when using real-world data), diagnostic costs would likely increase, and actionable targets could be missed. However, we did explore the impact of most assumptions with extensive sensitivity analyses and by consulting clinical experts. Costs for keeping SoC tests up to date and delays caused by required technical adaptations and validations were not included in our model.

To our knowledge, this research provides the first set of evidence of the most optimal use of WGS as a diagnostic test from a health economic perspective. Our findings provide direction for future research into and design of the diagnostic pathway for WGS. Optimising (adapting and updating) SoC diagnostics to reduce costs and increase efficiency could be considered but requires time and resources. In contrast, optimising WGS is less often necessary and requires fewer resources, as it provides information about the whole genome and concerns only one test. Additionally, we provided threshold values of several factors that need to be achieved to make future research worthwhile. Future research should focus on the clinical value of WGS, specifically understanding how to use genomic information for treatment decision support and biomarker development. Furthermore, novel predictive biomarkers for clinical response to immunotherapy are needed and could lead to more health gains and cost savings [5]. Currently, only a subset of patients experiences clinical benefit from immunotherapy, resulting in over treatment with expensive therapies.

Conclusions

Our analysis suggests that providing WGS as a diagnostic test would be cost effective compared with SoC followed by WGS and SoC alone if costs for WGS decreased and additional patients with actionable targets were identified. This cost-effectiveness model can be used to incorporate new findings iteratively and to support ongoing decision making regarding the use of WGS in this rapidly evolving field.

Below is the link to the electronic supplementary material.

Supplementary file1 (PDF 831 kb)

This analysis suggests that providing whole-genome sequencing (WGS) as a diagnostic test results in more quality-adjusted life-years (QALYs) and costs, that standard of care (SoC) followed by WGS results in fewer QALYs and more costs, and that both strategies are not cost effective compared with SoC alone.

If costs for sequencing decrease to €2000, WGS as a diagnostic test and SoC followed by WGS are cost effective compared with SoC alone when 2.7 and 4.8% additional patients with actionable targets are detected, respectively.