Cost-effectiveness of sleeping sickness elimination campaigns in five settings of the Democratic Republic of Congo

Gambiense human African trypanosomiasis (gHAT) is marked for elimination of transmission by 2030, but the disease persists in several low-income countries. We couple transmission and health outcomes models to examine the cost-effectiveness of four gHAT elimination strategies in five settings – spanning low- to high-risk – of the Democratic Republic of Congo. Alongside passive screening in fixed health facilities, the strategies include active screening at average or intensified coverage levels, alone or with vector control with a scale-back algorithm when no cases are reported for three consecutive years. In high or moderate-risk settings, costs of gHAT strategies are primarily driven by active screening and, if used, vector control. Due to the cessation of active screening and vector control, most investments (75-80%) are made by 2030 and vector control might be cost-saving while ensuring elimination of transmission. In low-risk settings, costs are driven by passive screening, and minimum-cost strategies consisting of active screening and passive screening lead to elimination of transmission by 2030 with high probability.


Introduction
Gambiense human African trypanosomiasis (gHAT) has caused three decades-long epidemics in Western, Central, and Eastern-Sub-Saharan Africa during the twentieth century [1]. Caused by the parasite Trypanosoma brucei gambiense and transmitted by tsetse, the disease is nearly always fatal if untreated [2]. No vaccine exists against gHAT, and current treatments are incompatible with mass drug administration strategies, but there is a relatively diverse toolbox available to control the infection [3]: a range of diagnostics for screening in at-risk villages by mobile teams, termed as active screening (AS), as well as diagnostics to test symptomatic individuals in fixed health facilities, termed as passive screening (PS). A new oral therapeutic -fexinidazole -became available in 2020, expanding the treatments available for outpatient care [4]. In addition, vector control (VC) can now be deployed across moderate scales to reduce tsetse density [5][6][7][8][9]. This toolbox has enabled a vast reduction in burden from 37,000 reported cases at the peak of the epidemic in 1998 to 864 reported cases in 2019 [10][11][12].
In 2012, under 15 years after the peak of the last epidemic, gHAT was targeted for elimination of transmission (EOT) by 2030, and by 2014 case reduction outpaced the intermediate goals set by the World Health Organization (WHO) [3,13]. Consequently, the pursuit of EOT brings new questions to the fore: because gHAT activities take place in settings facing both resource constraints and competing health needs, what are the resource implications of further pursuing gHAT EOT by 2030? Which are the most efficient strategies toward gHAT suppression?
The key economic arguments for and against disease elimination programmes may appear at odds: (i) the long-term benefits of elimination or eradication and the subsequent savings borne out of scaled-back disease control measures are substantial, yet (ii) the endgame is comprised of increasingly more expensive activities on a per-case basis without the guarantee of success. Smallpox eradication is heralded as an example of (i): eradication activities saved $1.35 billion (US$) but had a cost of $100 million [14]. Whereas polio -for which there were 42 cases reported in 2016 -exemplifies (ii), costing $3.3 billion each year, or $1 billion more than a suitable control strategy without eradication goal. [15].
To date, there has been one cost-effectiveness analysis of control and elimination strategies for gHAT, showing that new technologies to facilitate diagnosis, treatment, and curtail vector transmission could make elimination feasible at a moderate cost-effectiveness ($400-1500 per disability-adjusted life year (DALY) averted) for high-and medium-transmission settings, but at low cost-effectiveness for low-transmission settings [16]. However, it is worth reconsidering these questions while taking into account specific, local transmission dynamics, which is now possible thanks to recent developments in model calibration [17], and new cost estimates [18,19]. Moreover, including realistic levels of screening based on regional data allows us to consider expenses and cost-effectiveness in these real-world settings.
In the current study, we undertake an economic evaluation of four gHAT control and elimination strategies in distinct transmission settings in five health zones of the Democratic Republic of Congo (DRC), the country that reported 74% of global cases in 2018 [11] and 71% of global cases in 2019 [20]. We adopt a modelling framework in order to examine the comparative efficiency of gHAT elimination strategies while taking into account the short-and long-term interplay of epidemiological and economic factors of the disease.

Results
We selected five health zones in DRC, described in Table 1, spanning the spectrum of the WHO's incidence categories [3]. For each health zone, we simulated four strategies, depicted in Figure 1, using a variant of the "Warwick gHAT model", SEIRS-type deterministic dynamic model that explicitly simulates the transmission between humans via tsetse (see supplement S1.2.1). The model was fitted to case data from the health zones and projections were simulated for 2020-2050 [17,21]. The same underlying model framework has also been used to examine the epidemiology of gHAT in DRC and Chad [6,[22][23][24][25].
Strategies are configurations of the following interventions as delineated by the national policy and in line with WHO recommendations (see Figure 1 and section S1.2.2 for more details) [3,26]. Alongside passive screening (PS) in fixed health facilities, the strategies included active screening (AS) at (i) average or (ii) high-coverage levels, alone (i-ii) and in tandem with vector control (VC; iii-iv). Status quo strategies, or the current practice, are considered to be equivalent to Strategy 1, except in Yasa Bonga, where VC has been in place since mid-2015, and therefore consider the status quo equivalent to Strategy 3 and we omit strategies without VC in the analysis of this health zone. Additionally, the model has a transition from the suppression to post-elimination phase, devised to simulate cessation of AS and VC when no cases are reported for three consecutive years, allowing as well for re-activation of AS should a case be detected in a fixed health facility (i.e. through PS). Therefore, by relying on a realistic index of disease suppression -detected cases -our model shows not only whether strategies lead to earlier or later cessation of activities, but also whether strategies lead to mistakenly early cessation if screening activities are insufficient.

Feasibility of gHAT elimination and health impact
The feasibility of EOT and cessation of AS and VC are shown in Table 2. The probability of elimination was calculated as the proportion of the 10,000 model iterations where there were no new infections (detected or not) by 2030, and the probability of reactive screening (RS) was calculated as the proportion of the 10,000 model iterations where there was at least one case detected by PS after AS had been ended. The year that AS ends is the last year of AS activities before cessation; RS may occur after. While the risk category of each health zone (Table 1) influences the year when EOT is expected -places with higher incidence likely meet EOT later than places with lower incidence -the implementation of VC is predicted to substantially expedite EOT across all moderate-and high-risk settings and lower the probability of RS operations.
In low-risk Budjala, EOT appears imminent with any strategy, although employing Max AS over Mean AS would marginally reduce the uncertainty. In the moderate-risk settings of Mosango and Boma Bungu, EOT may occur by 2030 even with only Mean AS (79% and >99% probability, respectively) and EOT is almost certain by 2040. With additional AS coverage, the predicted probability of EOT by 2030 increases from 79% to 92% in Mosango, but VC would both ensure EOT by 2030 and bring forward the mean time to elimination by five years. In Yasa Bonga, where VC activities began in 2015 and scaled up completely by 2017, the model predicts that EOT may have been achieved by 2017 if VC was present in all areas of ongoing transmission, although predicted case detections are still expected in the early 2020s followed by AS cessation in 2024 (95% CI: 2021-2028). In high-risk Kwamouth, EOT is predicted to be almost impossible by 2030 (<0.01) in the absence of VC, and unlikely by 2040 (11-13% depending on AS coverage), but adding VC is predicted to bring forward EOT by more than two decades. Zero detections are more informative as an EOT proxy when VC is in situ; if AS is stopped after three years of zero detections, there is up to a 62% probability that RS would be necessary (in Kwamouth) in the absence of VC, but at most 19% of VC simulations result in RS (in Budjala).

3
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review)  In strategies 3 and 4, vector control (VC) is simulated assuming an 80% tsetse density reduction in 1 year, except in Yasa Bonga, where the reduction of tsetse density was estimated at 90% in the literature [9]. PS is in place under all strategies.

Health outcomes
The outputs of the transmission model were inputs in a probability tree model of disease outcomes (see Methods and Figure 4), from which cases, deaths and disability-adjusted life-years were calculated by standard conventions(for further elaboration on DALYs, see supplement S1.3) [27]. The health impact and net costs of each strategy between 2020 and 2040 are shown in Table 3.
Yasa Bonga, Boma Bungu, and Budjala are each predicted to have an average of ≤ 5 cases and gHAT-related deaths over the next 20 years. Mosango is predicted to have more cases (≤ 23) and deaths (≤ 13) in the absence of VC. Kwamouth has the most predicted cases and deaths, although the burden may be cut by three-quarters with the deployment of VC. In terms of DALYs, Kwamouth sustains the worst burden even under the most ambitious strategy (1,718 DALYs, 95% PI: 562-3,656) compared to the least ambitious strategy in moderate-risk Mosango (426 DALYs, 95% PI: 32-1,418), and any strategy in low-risk Boma Bungu (11)(12)(13)(14)(15)(16)(17). Supplemental outcomes of treatment are found in section S1.2.3; 4 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review)  97% of treated stage-1 patients and 93% of treated stage-2 patients will be cured, an additional 1% and 3%, respectively, will be cured after sustaining severe adverse events, and between 2% and 4% of patients will need rescue treatment.
Fewer than 1% of all treated cases will die. Further details of the outcomes under different time horizons are found in supplementary tables S18-S19 and on the projectwebsite:https://hatmepp.warwick.ac.uk/5HZCEA/v2/.

Costs
We devised cost functions to calculate costs from the perspective of the healthcare system as a whole (including all levels of government and donors) in 2018 US$. The total mean costs in each location range between $490,000 (Boma Bungu) and $5.43 million (Kwamouth), and on a yearly per-capita basis, costs range from $0.20 (Budjala) to $1.97 (Kwamouth) ( Table 3). The expected costs of strategy components are shown in Figure 2; these calculations are derived using the mean costs and the probability that activities have not ceased by a particular year. For all locations except Kwamouth, additional costs in the latter half of the 2030s arise from PS, as most simulations indicate that AS and VC would have ceased by this period. In Yasa Bonga, Mosango, and Kwamouth, about three-quarters of the cumulative costs are expected to be spent within the first five years of the 2020s and then stabilise when AS and VC cease ( Figure S4). Costs with the least ambitious (default) strategy in Kwamouth and Mosango are expected to overtake the costs of a strategy with VC by 2040, indicating that such a strategy would be economically neutral or even cost-saving in 20 years. In Boma Bungu and Budjala, where scant screening is undertaken and VC would cease quickly, costs are slow to accumulate throughout the period of study. The cost breakdowns across different investment horizons are shown in Figure S3 and show similar patterns to those in Figure 2. Interested readers may explore the costs and their breakdowns on the projectwebsite:https://hatmepp.warwick.ac.uk/5HZCEA/v2/.
5 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted April 8, 2021.  Table 3: Summary of effects and costs 2020-2040. Two differences should be noted between these estimates and those used for decision analysis shown in table 4. First, these estimates are not discounted. Second due to asymmetric distributions, a naive difference in mean costs would not equal the mean differences in costs across simulations -the metric we used in decision analysis. Undetected cases are reflected in deaths, as very few deaths (<1 percent) originate from treated cases. Estimates shown are means and 95% prediction intervals (PI) of the cases, deaths, DALYs, and costs across iterations of the dynamic transmission model.

Cost-effectiveness
The cost-effectiveness of each strategy in each health zone is displayed in Table 4 and select features are illustrated in Figure 3. The standard measure of cost-effectiveness (or comparative efficiency) is the incremental cost-effectiveness ratio (ICER), defined as the additional cost to avert an additional DALY of disease compared to the next-best strategy. We considered the health impacts and costs in a relatively long-term horizon (2020-2040), discounting at a yearly rate of 3% in accordance with standard conventions [27]. Therefore, strategies are considered "minimum cost" for the least costly strategy over a 20-year horizon in present-day costs (costs discounted at 3% annually). In all health zones but Mosango, the minimum-cost strategy is the current practice; in Mosango, the current practice is "dominated" because the addition of VC yields cost-savings on a 20-year horizon while averting more DALYs than strategies without VC (Mean AS and Max AS). All strategies designated as "dominated" are more costly and less effective at averting DALYs than another strategy.
A second form of dominance is evident in Budjala, where Max AS is designated "weakly dominated". A strategy is "weakly dominated" if the next more expensive strategy (in absolute terms) has a lower ICER -indicating that although more expensive, the latter is more efficient on a per-DALY-averted basis than the former. In Budjala, the Max AS strategy is weakly dominated by the Max AS & VC strategy, because although the latter is more expensive the additional cost is better justified on a per-DALY basis. For an extended discussion on dominated and weakly dominated strategies, see Section S1.5. In short, both dominated and weakly dominated strategies present inefficiencies.
A strategy is considered cost-effective, though not cost-saving, when it maximizes DALYs averted while having an ICER below an acceptable threshold, considered as the willingness-to-pay (WTP) threshold. The WTP value is denominated in costs-per-DALY averted and we do not prescribe any specific WTP value, as per WHO recommendations.
6 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.   We leave the choice of threshold for the policy-maker, though we note that other analyses have shown that the WTP in DRC for any new strategy is approximately $5-230 per DALY averted [28], and past guidelines have delineated cost-effectiveness threshold as equal to the GDP per capita of the country [27], which was $557 in DRC in 2018. In all places but Kwamouth, the ICER for strategies other than the minimum-cost-strategy are moderately high, indicating low comparative efficiency for settings in DRC. In Kwamouth, the Max AS & VC strategy maximizes DALYs averted, at an ICER of $4,421, which is considered too high by any current or historical metric, but the Mean AS & VC strategy has an ICER of $4, below historical WTP thresholds, so it satisfies both criteria of maximizing DALYs averted and staying under the WTP threshold.

Cost-effectiveness and elimination of transmission by 2030
In Yasa Bonga and in Boma Bungu, EOT is estimated to occur by 2030 with >99% probability, and any more ambitious operations would be to detect existing cases -yielding high ICERs. In Budjala, perhaps because of minimal AS coverage, there is a 36% probability that RS would be deployed in Budjala after cessation (see Table 2), and while VC would nearly halve that prospect, it would come at a steep cost of $2,932 per DALY averted. Lastly, while in Mosango adding VC to the current practice of AS both saves money and increases the probability of EOT by 2030 from 79% to >99%, in Kwamouth the same shift makes EOT by 2030 possible, but at a cost of $4 per DALY averted.

7
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review)  Table 4: Summary of cost-effectiveness, assuming a time horizon of 2020-2040. Cost differences and DALYs averted are relative to the comparator-first strategy listed for each location. Mean DALYs averted and mean cost differences are shown; these estimates are discounted at 3 percent per year in accordance with guidelines. The uncertainty analysis (columns 5-8) shows the probability that a strategy is cost-effective. Strategies highlighted in pink are optimal strategies: the strategies for which the mean net monetary benefit (NMB) is highest. ICER: incremental cost-effectiveness ratio, DALY: disability adjusted life-years. For an extended discussion of these terms, see supplement section S1.5.

Cost-effectiveness in the presence of probabilistic uncertainty
While the ICER provides an intuitive measure of efficiency, it does not integrate uncertainty in decision-making, so we present the ICERs here for comparability with other literature, but we base our conclusions of cost-effectiveness on the results within the net benefits framework, which examines the efficiency of gHAT end-game strategies in the presence of uncertainty (see sections S1.5, S2). We consider a strategy 'optimal' or cost-effective if the net monetary benefit (NMB) is maximum at a given willingness-to-pay value. In Yasa Bonga, Boma Bungu and Budjala, the predicted minimum-cost strategy has a 97-99% probability of EOT by 2030. After accounting for uncertainty, the current practice (Mean AS & VC) in Yasa Bonga is the optimal strategy across all WTP thresholds <$1,000 in more than 78% of simulations. In Budjala, there is even more support (76-91% of simulations) indicating that the least ambitious intervention (Mean AS) is optimal across all WTP values. In Mosango, Mean AS & VC is cost-saving with a 49% probability. In Kwamouth, more than 65% of simulations favoured the addition of VC to Mean AS for WTP values above $250, and in fact, was cost-saving in 42% of simulations. Cost-effectiveness acceptability curves, expressing the same information as in Table 4 but in a more conventional format within the net benefits literature, are shown in Figure S5.

Scenario analyses
In addition to the uncertainty due to probabilistic uncertainty, we perform scenario analyses to understand the robustness of our findings when the uncertainty in assumptions cannot be quantified probabilistically.
In Mosango and Kwamouth, VC would both raise the probability of EOT by 2030 and minimize the probability 8 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted April 8, 2021. ; https://doi.org/10.1101/2020.08.25.20181982 doi: medRxiv preprint of RS. However, these conclusions are contingent on assumptions around VC operations. In Mosango, the potential for cost-savings is contingent on the assumption that an operation needs to be deployed across 100 km of riverbank only. In fact, with a small operation, VC remains cost-saving even if tsetse population reduction is as low as 60% in one year (see Figure S10). However, if an operation is needed in Mosango that is closer in extent to the one in Yasa Bonga (across 210 km of riverbank dotted with 40 targets per kilometer of riverbank) then VC would be optimal if WTP would exceed $1,488-2,051 per DALY averted, depending on the effectiveness of VC to reduce tsetse populations.
In Kwamouth, we found that if target density must be doubled, WTP must exceed $236-343 per DALY averted after considering deviations in VC effectiveness (see Figure S9). Lastly, we performed additional scenario analyses to examine the impact of our default assumptions around time horizons, discounting (SI section S1.7). The results across shorter (2020-2030) and longer (2020-2050) time horizons are very similar (see SI Figures S7-S8 and SI  Tables S25-S26). The results assuming no discounting of costs or effects ( Figure S6) lead to similar recommendations for WTP values <$1000.

Discussion
Applying a decision-analytic framework across five health zones of DRC we have assessed current efforts to control gHAT, as well as the efficiency of alternative strategies to transition beyond suppression to EOT. The transmission model predicted substantial declines in observed gHAT cases and the underlying transmission in all locations using any strategy, but the cumulative burden of disease and the capacity to reach EOT by 2030 varied considerably.
In Mosango and Kwamouth the addition of VC to Mean AS is predicted to expedite EOT and be cost-effective for low WTP (Table 4). However, in Boma Bungu and Budjala current practice appears sufficient to ensure EOT by 2030 despite Boma Bungu's "moderate" incidence and Budjala's low AS coverage (<0.5%) (Tables 1 and 4). Although the expected total costs in each location vary between $480,000 in Boma Bungu to $5.26 million in Kwamouth, the optimal strategies have a expected cost of at most $3.71 million (Table 3). Per person, the optimal strategies would not exceed $1.35 per year per person protected (Table 3) -comparable to many other global health interventions [15]. Notably, while healthcare costs and DALYs can be considered narrow criteria to justify investments in elimination, we chose to conduct an analysis that would allow us to understand the relative efficiency of these strategies by the same rubric employed to assess the cost-effectiveness of other disease programs; we found that local gHAT elimination is possible and efficient at modest economic costs.
Additionally, multiple CEAs across zones could inform allocative efficiency in case of budget constraints, as it shows that, among these five health zones, the most efficient use of scale-up dollars would first be in Mosango and then in Kwamouth.
While a full analysis of past costs is nearly impossible due to the heterogenous nature of funding for gHAT activities, and because national health accounts contain no details on a disease like gHAT, our cost parameters were informed by recent cost studies in Yasa Bonga and Mosango, therefore updating the costs used in the past [16,29]. A decomposition of the costs showed that while there are contributions from varied factors, some broad patterns emerged. First, AS costs and, when applicable, VC costs will play a large role in overall costs everywhere except in Budjala, where no more than PS and minimal AS is recommended. Second, cessation of AS and VC activities means costs are expected to decrease during the 2020's, and in Kwamouth, investments in VC in the early 2020's could be recovered by the late-2030's (Fig-S4). Therefore, while cost-effectiveness analysis cannot take into account all potential elements of implementation, such as the political will to implement certain components of the strategies proposed, analysis is insightful to understand the comparative efficiency within any two strategies or to raise funds, as we have presented the time-scale for a return on investment.
Importantly, we have found that our conclusions are quite robust to the choice of time-horizon as well as discount rate, and our framework has accounted for parameter uncertainty. We have chosen to consider a 20-year horizon with supplemental results for a 10-and 30-year horizon because strategies with long-term consequences like EOT would be under-appreciated in a short horizon, but conclusions do not differ substantially (see tables S25 and S26 and figures S5 and S6).
Simulation results presented here suggest that local EOT of gHAT is epidemiologically and operationally feasible across different risk settings with the current toolbox. As expected, higher-burden health zones expect later EOT, but VC raises the possibility and substantially expedites EOT, as seen across DRC [21]. The model utilised here 9 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted April 8, 2021  did not explicitly include potential hindrances such as asymptomatic human infection or animal reservoirs [30], however the long-tailed distribution used for disease progression captures possible long-term asymptomatic carriage. Skin-only infections or asymptomatic individuals whose infection may self-resolve would require further model modifications 10 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted April 8, 2021. ; https://doi.org/10.1101/2020.08.25.20181982 doi: medRxiv preprint and ongoing modelling work will examine potential risks and impact. Previous modelling analyses indicate that the existence of animal reservoirs may not greatly alter timelines to EOT, especially if VC is implemented [6,23]. This is the first analysis where fexinidazole is the default treatment, but to reflect the caution in the WHO guidelines 65% of simulated patients were assumed to be treated with NECT or on an in-patient basis due to late-stage detection, low body weight, or the young age of the patient (see Table S2 and sections S6.5.14 to S6.5.16). The prospect of a single-dose oral treatment for gHAT (acoziborole), currently in trials, indicated for broad swaths of the population could further bolster the impact of medical interventions if it allows the treatment of gHAT suspects without confirmation; evaluating the impact of a drug that overcomes the limitations of the current treatment arsenal is beyond the scope of this analysis, and a careful treatment of the matter ought to be pursued [13]. Similarly, the impact of more specific diagnostics that may become available could change the conclusions of this cost-effectiveness analyses, as it is comparative efficiency that it informs.
Previous cost-effectiveness analyses have underlined key factors associated with suppression and elimination of gHAT, but this is the first study to confront the question of cost-effectiveness with data from five specific locations [16,31,32]. An analysis by Bessel et al shows that RDTs in DRC are cost-effective when employed within AS and PS [31], and an analysis by Davis et al shows that once-per-year AS campaigns within endemic villages are cost-effective to control disease [32]. Only one other analysis, performed by Sutherland et al, examined the combination of multi-faceted strategies, and they found that elimination strategies are likely cost-effective at a WTP of $400-1500 (in 2013 US$) in high-and moderate-incidence places [16]. One important difference to the analysis by Sutherland and colleagues is that we assumed yearly (rather than biennial) AS campaigns at coverage levels that match historical averages (rather than a fixed value of 80% of the population -a high value relative to historical coverage observed in DRC) [16]. We also assumed that AS would persist for three years after the last detected case in a health zone, irrespective of the initial incidence of the health zone, whereas Sutherland et al assumed no AS in low-incidence health zones and therefore had to recommend VC activities to reach EOT in those contexts. In our lower incidence health zones (Boma Bungu and Budjala), the presence of AS meant that EOT could be reached without VC (Table-4).
Our analyses reinforce previous findings that VC would be both an expedited method of achieving EOT and cost-effective in one moderate-incidence health zone (Mosango) and one high-incidence health zone (Kwamouth). However, determining the amount of VC necessary to reach a desired tsetse population reduction is a complicated task, as it depends on local ecology. For transparency, we have shown cost-effectiveness results in a three-way sensitivity analysis in places where VC might be warranted and was not previously in place (Mosango and Kwamouth; see S10 and S9). In Mosango, a modest VC operation could be cost-saving, but these results are contingent on the geographic distribution of current cases as well as on the sensitivity of the tsetse population to targets (see S10). Because Kwamouth is substantially larger than other health zones, the geographic clustering of cases will be important to determine whether all high-risk areas can be addressed in a cost-effective way (see Figure S9).
Our findings took into account historical improvements in PS, which was made possible by more recent data and novel model calibration; this element of the strategies has difficult-to-quantify impacts that might explain some of the difference between our results and previous findings. The fact that our analysis was more optimistic about current practice in lower-incidence health zones underscores the potential of a well-equipped health system that can serve self-presenting gHAT cases. Future analyses are warranted on the impact of gHAT screening integrated into the primary care system.
Within the last year, DRC has contended with the emergence of COVID-19, which triggered an untimely but temporary disruption to AS activities [33]. But while these interruptions are lamentable, the EOT goal could still be within reach as long as disruptions remain short [34]. While reported gHAT cases now number under 1000 globally, there remain 51 million people in environments that could sustain gHAT transmission [11]. The epidemic potential of gHAT in locations as resource-constrained as DRC underscores the importance of careful deliberation around strategies of gHAT elimination.

Methods
Word count (methods): 1438 words 11 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

Transmission model
We employed a previously published dynamic, transmission "susceptible-exposed-infected-recovered-susceptible (SEIRS)" model. While the model is detailed elsewhere [17,21], briefly, we simulated gHAT infection and vector transmission in a deterministic compartmental model using a set of ordinary differential equations which represent average expected dynamics. Stochastic features of the system -corresponding to probabilistic, or chance events were included for case detections and deaths and were simulated by sampling from the prevalence of infected individuals for reporting to care in fixed facilities or mobile teams. Specifically, new infections and person-years spent in stage 1 and 2 are deterministic outputs, dependent on the parameter inputs, but case reporting and deaths also include random sampling using the beta-binomial distribution. The absolute population that is seen by active screening is fixed, although the population is projected to increase by 3% per year. The key outputs of the dynamic and diagnostic models include mortality in undetected cases, detected cases in stages 1 and 2, and DALYs before and after presenting to care for all interventions. Further details on the development of the model are found in Section S1.2.1.

Strategies
All strategies (see Figure 1) are constituted of a combination of interventions, detailed below: • Active screening (AS): at-risk village populations are screened each year by mobile teams and suspected cases are confirmed before treatment. In two strategies (1 and 3, "Mean AS") the number of people screened annually per health zone is equal to the mean number of individuals screened during 2014-2018 (the last five years for which there are data). In two other strategies (2 and 4, "Max AS"), the number of people screened annually is the maximum coverage achieved in that health zone in any single year between 2000-2018. • Passive surveillance (PS): testing of symptomatic patients who self-present to local health centres with gHAT diagnostics, followed by confirmation and treatment [3]. • Large-scale vector control (VC): twice-yearly deployment of tiny targets to control the population of tsetse, with an assumed 80% reduction after one year, consistent with the lower bound efficacy estimates of field studies [5,6,8,9]. In Yasa Bonga, VC reduced the tsetse population by 90% in one year, so our assumption in that health zone is a bit higher. This intervention is assumed to span the areas with ongoing transmission in the previous five years, in keeping with existing operations. In all strategies, AS (at either 'Mean' or 'Max' coverage) and PS are in place, with supplemental VC deployed in Strategies 3 and 4 (see Figure 1). In Yasa Bonga, VC has taken place since mid-2015, so only the two VC strategies (3 and 4) were considered.
The transition between the suppression and post-elimination phases, not included in the previous projections [21], is simulated as the cessation of both AS and VC after three consecutive years of zero detected cases by any screening modality (AS or PS). Should a new case present to fixed health facilities (through PS), reactive screening (RS) begins the following year: equivalent in form and intensity to the previous AS, RS is simulated until there are two consecutive years of zero case detection by any mechanism (either AS or PS). The availability of screening in health facilities (PS) is assumed to remain constant for the duration of our simulations, even after cessation of AS and VC and EOT.

Feasibility of EOT by 2030
For each realisation we record the year that EOT is achieved using a proxy threshold of less than one new infection per year in the health zone. We presented feasibility with a variety of metrics indicating the timing of EOT as well as the uncertainty around the goal by 2030: • The mean and 95% prediction interval for the year of EOT is calculated by taking the year in which cases fall under 1 in each iteration of the model and calculating the arithmetic mean and 2.5th and 97.5th percentile of the sample. • The probability of EOT by 2030 was the percentage of iterations in which the year of EOT was 2029 or earlier.
• The estimates for year of AS cessation -and if applicable, VC cessation -were calculated by taking the year when AS first stopped occurring, and calculating the mean and 95% prediction intervals as above. • The probability of RS by 2050 was the percentage of iterations in which there was at least one year AS after AS had first ceased (after at least one year of no AS).
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

Cases, deaths, and treatment outcome model
Using a probability tree (Figure 4), we simulated the treatment process separately for stage 1 and stage 2 disease, sorting patients into the WHO-recommended treatment, then into treatment success or failure, correct diagnosis in the event of treatment failure, and progression to rescue treatment or to death (for further details, see S1.2.3) [4]. Disability-adjusted life years (DALYs) were calculated for both detected and undetected cases (which were assumed to all be fatal). DALYs are the sum of the years of life lost to disability (illness) (YLD), weighted by the severity of the disease, and the years of life lost by fatal cases (YLL) [27,35].
Uncertainty arose from two sources. First, the uncertainty in dynamic model parameters fit to historical data, formally drawn by an adaptive MCMC algorithm (for more information, see our previous publication, [17]). Second, uncertainty surrounding treatment outcomes in literature-based estimates, which was formalized by assigning probability distributions according to standard conventions (for details, see Section S6). Uncertainty in final estimates of cases, deaths, and treatment outcomes were therefore presented as means and 95% prediction intervals, calculated by taking the 2.5th and 97.5th percentiles of the simulated samples.

Costs
Based on field experience from Yasa Bonga and Mosango [18,19], we devised a cost function to calculate the costs of detection via AS and PS, prevention via VC, and treatment. Screening costs include diagnostics, mobile teams, and fixed facilities (more details can be found in section S1.4). Costs of screening were a mix of fixed and variable costs: fixed costs calculated according to health facilities or the mobile teams needed, staff training, and variable costs according to the people screened by health facility or mobile team, and people confirmed.
Treatment costs include staging via lumbar puncture (when indicated), and drug administration; these costs were largely based on the literature. VC costs were scaled according to the geographic expanse treated (in terms of km of riverbank) and the density of targets placed along riverbanks (final costs are described in tables S4-S14). Because no published estimates of VC cost exist from DRC, we used published estimates from Uganda.
Costs were from the perspective of the healthcare system collectively (including the government at all levels and donors). Costs were parameterised using values from the literature converted to 2018 US$ values (see Section S6). Start-up costs were annualised. Uncertainty surrounding costs was formalized by assigning probability distributions to each parameter according to standard conventions and drawing Monte Carlo samples (for details, see Section S6). Uncertainty in final estimates of costs were therefore presented as means and 95%, calculated by taking the 2.5th and 97.5th percentiles of the samples.

Cost-effectiveness analysis
We computed incremental cost-effectiveness ratios (ICERs) to perform our initial cost-effectiveness studies using the mean DALYs averted and mean costs for each strategy compared to the current strategy (Mean AS for all but Yasa Bonga, where the current strategy is Mean AS & VC).
To account for parameter uncertainty in the economic evaluation, we adopted the net benefits framework, which expresses the probability that an intervention is optimal at a range of willingness-to-pay thresholds (WTP). We consider a strategy 'optimal' or cost-effective if the net monetary benefit (NMB), or the difference between the monetary value of DALYs averted and the additional costs (versus the comparator) is favourable (maximum) in the following formulation: where WTP is the willingness-to-pay, or the cost-effectiveness threshold. Prediction intervals of ICERs are mathematically problematic [36], so convention dictates that we show the probability that a strategy is optimal at range of WTP values, for which we chose $0-1000 in $250 increments. Although WTP values should be designated by the countries, we note that some values for DRC in the literature are approximately $5-230 per DALY averted [28], or $557 (equivalent to the country's GDP per capita in 2018) [27]. For further elaboration on the framework and our implementation, see section S1.5.
We also performed scenario analyses to examine the impact of our default assumptions around time horizons, discounting, and on the efficacy and cost of VC (see supplemental section S1.7).

13
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) The copyright holder for this preprint this version posted April 8, 2021. ; https://doi.org/10.1101/2020.08. 25 Figure 4: Treatment for diagnosed gHAT patients is modeled as a branching tree process of possible health outcomes including eligibility for novel fexinidazole. Abbreviations: SAE: Serious adverse events, IP: inpatient care, OP: outpatient care.
14 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) The copyright holder for this preprint this version posted April 8, 2021. ; https://doi.org/10.1101/2020.08.25.20181982 doi: medRxiv preprint Analyses were performed using R version 3.6. Computational considerations are detailed in the Open Science Framework repository for this analysis at https://osf.io/xbwte/. For interested readers and policy-makers, we created a project website to showcase results and sensitivity analysis: https://hatmepp.warwick.ac.uk/5HZCEA/v2/.

Role of the funding source
The funders of the study had no role in study design, data collection, data analysis, data interpretation, or writing of the report.

Data Availability
Information about the data used for fitting is described in Crump et al [17]. The data used for clinical outcomes and costs were estimates from the literature. Screening data from the five health zones in this study were used to inform future potential screening coverage and were obtained through the WHO HAT Atlas. Data cannot be shared publicly because they were aggregated from the World Health Organisation's HAT Atlas which is under the stewardship of the WHO. Data are available from the WHO (contact neglected.diseases@who.int or visit https://www.who.int/trypanosomiasis_african/country/foci_AFRO/en/) for researchers who meet the criteria for access to confidential data. Assumptions and estimates were parameterized according to conventions in the economic evaluation literature [37]. For assumptions around intervention, treatment effects and costs, see supplemental sections S3.2.1, S1.3.1, S2.

Code Availability
Full access to the code via Open Science Framework: https://osf.io/xbwte/.
15 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) The copyright holder for this preprint this version posted April 8, 2021. ; https://doi.org/10.1101/2020.08.25.20181982 doi: medRxiv preprint . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) The copyright holder for this preprint this version posted April 8, 2021. ; https://doi.org/10.1101/2020.08.25.20181982 doi: medRxiv preprint