Accelerating Progress Towards the 2030 Neglected Tropical Diseases Targets: How Can Quantitative Modeling Support Programmatic Decisions?

Over the past decade, considerable progress has been made in the control, elimination, and eradication of neglected tropical diseases (NTDs). Despite these advances, most NTD programs have recently experienced important setbacks; for example, NTD interventions were some of the most frequently and severely impacted by service disruptions due to the coronavirus disease 2019 (COVID-19) pandemic. Mathematical modeling can help inform selection of interventions to meet the targets set out in the NTD road map 2021-2030, and such studies should prioritize questions that are relevant for decision-makers, especially those designing, implementing, and evaluating national and subnational programs. In September 2022, the World Health Organization hosted a stakeholder meeting to identify such priority modeling questions across a range of NTDs and to consider how modeling could inform local decision making. Here, we summarize the outputs of the meeting, highlight common themes in the questions being asked, and discuss how quantitative modeling can support programmatic decisions that may accelerate progress towards the 2030 targets.

Understanding the Potential Impact of Different Drug Properties on Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) Transmission and Disease Burden: A Modelling Analysis.

BACKGROUND: The public health impact of the coronavirus disease 2019 (COVID-19) pandemic has motivated a rapid search for potential therapeutics, with some key successes. However, the potential impact of different treatments, and consequently research and procurement priorities, have not been clear. METHODS: Using a mathematical model of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) transmission, COVID-19 disease and clinical care, we explore the public-health impact of different potential therapeutics, under a range of scenarios varying healthcare capacity, epidemic trajectories; and drug efficacy in the absence of supportive care. RESULTS: The impact of drugs like dexamethasone (delivered to the most critically-ill in hospital and whose therapeutic benefit is expected to depend on the availability of supportive care such as oxygen and mechanical ventilation) is likely to be limited in settings where healthcare capacity is lowest or where uncontrolled epidemics result in hospitals being overwhelmed. As such, it may avert 22% of deaths in high-income countries but only 8% in low-income countries (assuming R = 1.35). Therapeutics for different patient populations (those not in hospital, early in the course of infection) and types of benefit (reducing disease severity or infectiousness, preventing hospitalization) could have much greater benefits, particularly in resource-poor settings facing large epidemics. CONCLUSIONS: Advances in the treatment of COVID-19 to date have been focused on hospitalized-patients and predicated on an assumption of adequate access to supportive care. Therapeutics delivered earlier in the course of infection that reduce the need for healthcare or reduce infectiousness could have significant impact, and research into their efficacy and means of delivery should be a priority.

Estimating the health impact of vaccination against ten pathogens in 98 low-income and middle-income countries from 2000 to 2030: a modelling study.

BACKGROUND: The past two decades have seen expansion of childhood vaccination programmes in low-income and middle-income countries (LMICs). We quantify the health impact of these programmes by estimating the deaths and disability-adjusted life-years (DALYs) averted by vaccination against ten pathogens in 98 LMICs between 2000 and 2030. METHODS: 16 independent research groups provided model-based disease burden estimates under a range of vaccination coverage scenarios for ten pathogens: hepatitis B virus, Haemophilus influenzae type B, human papillomavirus, Japanese encephalitis, measles, Neisseria meningitidis serogroup A, Streptococcus pneumoniae, rotavirus, rubella, and yellow fever. Using standardised demographic data and vaccine coverage, the impact of vaccination programmes was determined by comparing model estimates from a no-vaccination counterfactual scenario with those from a reported and projected vaccination scenario. We present deaths and DALYs averted between 2000 and 2030 by calendar year and by annual birth cohort. FINDINGS: We estimate that vaccination of the ten selected pathogens will have averted 69 million (95% credible interval 52-88) deaths between 2000 and 2030, of which 37 million (30-48) were averted between 2000 and 2019. From 2000 to 2019, this represents a 45% (36-58) reduction in deaths compared with the counterfactual scenario of no vaccination. Most of this impact is concentrated in a reduction in mortality among children younger than 5 years (57% reduction [52-66]), most notably from measles. Over the lifetime of birth cohorts born between 2000 and 2030, we predict that 120 million (93-150) deaths will be averted by vaccination, of which 58 million (39-76) are due to measles vaccination and 38 million (25-52) are due to hepatitis B vaccination. We estimate that increases in vaccine coverage and introductions of additional vaccines will result in a 72% (59-81) reduction in lifetime mortality in the 2019 birth cohort. INTERPRETATION: Increases in vaccine coverage and the introduction of new vaccines into LMICs have had a major impact in reducing mortality. These public health gains are predicted to increase in coming decades if progress in increasing coverage is sustained. FUNDING: Gavi, the Vaccine Alliance and the Bill & Melinda Gates Foundation.

Development, environmental degradation, and disease spread in the Brazilian Amazon.

The Amazon is Brazil's greatest natural resource and invaluable to the rest of the world as a buffer against climate change. The recent election of Brazil's president brought disputes over development plans for the region back into the spotlight. Historically, the development model for the Amazon has focused on exploitation of natural resources, resulting in environmental degradation, particularly deforestation. Although considerable attention has focused on the long-term global cost of "losing the Amazon," too little attention has focused on the emergence and reemergence of vector-borne diseases that directly impact the local population, with spillover effects to other neighboring areas. We discuss the impact of Amazon development models on human health, with a focus on vector-borne disease risk. We outline policy actions that could mitigate these negative impacts while creating opportunities for environmentally sensitive economic activities.

Spatial and temporal invasion dynamics of the 2014-2017 Zika and chikungunya epidemics in Colombia.

Zika virus (ZIKV) and chikungunya virus (CHIKV) were recently introduced into the Americas resulting in significant disease burdens. Understanding their spatial and temporal dynamics at the subnational level is key to informing surveillance and preparedness for future epidemics. We analyzed anonymized line list data on approximately 105,000 Zika virus disease and 412,000 chikungunya fever suspected and laboratory-confirmed cases during the 2014-2017 epidemics. We first determined the week of invasion in each city. Out of 1,122, 288 cities met criteria for epidemic invasion by ZIKV and 338 cities by CHIKV. We analyzed risk factors for invasion using linear and logistic regression models. We also estimated that the geographic origin of both epidemics was located in Barranquilla, north Colombia. We assessed the spatial and temporal invasion dynamics of both viruses to analyze transmission between cities using a suite of (i) gravity models, (ii) Stouffer's rank models, and (iii) radiation models with two types of distance metrics, geographic distance and travel time between cities. Invasion risk was best captured by a gravity model when accounting for geographic distance and intermediate levels of density dependence; Stouffer's rank model with geographic distance performed similarly well. Although a few long-distance invasion events occurred at the beginning of the epidemics, an estimated distance power of 1.7 (95% CrI: 1.5-2.0) from the gravity models suggests that spatial spread was primarily driven by short-distance transmission. Similarities between the epidemics were highlighted by jointly fitted models, which were preferred over individual models when the transmission intensity was allowed to vary across arboviruses. However, ZIKV spread considerably faster than CHIKV.

Spatiotemporal variations in exposure: Chagas disease in Colombia as a case study.

Age-stratified serosurvey data are often used to understand spatiotemporal trends in disease incidence and exposure through estimating the Force-of-Infection (FoI). Typically, median or mean FoI estimates are used as the response variable in predictive models, often overlooking the uncertainty in estimated FoI values when fitting models and evaluating their predictive ability. To assess how this uncertainty impact predictions, we compared three approaches with three levels of uncertainty integration. We propose a performance indicator to assess how predictions reflect initial uncertainty.In Colombia, 76 serosurveys (1980-2014) conducted at municipality level provided age-stratified Chagas disease prevalence data. The yearly FoI was estimated at the serosurvey level using a time-varying catalytic model. Environmental, demographic and entomological predictors were used to fit and predict the FoI at municipality level from 1980 to 2010 across Colombia.A stratified bootstrap method was used to fit the models without temporal autocorrelation at the serosurvey level. The predictive ability of each model was evaluated to select the best-fit models within urban, rural and (Amerindian) indigenous settings. Model averaging, with the 10 best-fit models identified, was used to generate predictions.Our analysis shows a risk of overconfidence in model predictions when median estimates of FoI alone are used to fit and evaluate models, failing to account for uncertainty in FoI estimates. Our proposed methodology fully propagates uncertainty in the estimated FoI onto the generated predictions, providing realistic assessments of both central tendency and current uncertainty surrounding exposure to Chagas disease.

Genomic Epidemiology of Severe Acute Respiratory Syndrome Coronavirus 2, Colombia.

Coronavirus disease (COVID-19) in Colombia was first diagnosed in a traveler arriving from Italy on February 26, 2020. However, limited data are available on the origins and number of introductions of COVID-19 into the country. We sequenced the causative agent of COVID-19, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), from 43 clinical samples we collected, along with another 79 genome sequences available from Colombia. We investigated the emergence and importation routes for SARS-CoV-2 into Colombia by using epidemiologic, historical air travel, and phylogenetic observations. Our study provides evidence of multiple introductions, mostly from Europe, and documents >12 lineages. Phylogenetic findings validate the lineage diversity, support multiple importation events, and demonstrate the evolutionary relationship of epidemiologically linked transmission chains. Our results reconstruct the early evolutionary history of SARS-CoV-2 in Colombia and highlight the advantages of genome sequencing to complement COVID-19 outbreak investigations.

Early Insights from Statistical and Mathematical Modeling of Key Epidemiologic Parameters of COVID-19.

We report key epidemiologic parameter estimates for coronavirus disease identified in peer-reviewed publications, preprint articles, and online reports. Range estimates for incubation period were 1.8-6.9 days, serial interval 4.0-7.5 days, and doubling time 2.3-7.4 days. The effective reproductive number varied widely, with reductions attributable to interventions. Case burden and infection fatality ratios increased with patient age. Implementation of combined interventions could reduce cases and delay epidemic peak up to 1 month. These parameters for transmission, disease severity, and intervention effectiveness are critical for guiding policy decisions. Estimates will likely change as new information becomes available.

Complementary Paths to Chagas Disease Elimination: The Impact of Combining Vector Control With Etiological Treatment.

Background: The World Health Organization's 2020 goals for Chagas disease are (1) interrupting vector-borne intradomiciliary transmission and (2) having all infected people under care in endemic countries. Insecticide spraying has proved efficacious for reaching the first goal, but active transmission remains in several regions. For the second, treatment has mostly been restricted to recently infected patients, who comprise only a small proportion of all infected individuals. Methods: We extended our previous dynamic transmission model to simulate a domestic Chagas disease transmission cycle and examined the effects of both vector control and etiological treatment on achieving the operational criterion proposed by the Pan American Health Organization for intradomiciliary, vectorial transmission interruption (ie, <2% seroprevalence in children <5 years of age). Results: Depending on endemicity, an antivectorial intervention that decreases vector density by 90% annually would achieve the transmission interruption criterion in 2-3 years (low endemicity) to >30 years (high endemicity). When this strategy is combined with annual etiological treatment in 10% of the infected human population, the seroprevalence criterion would be achieved, respectively, in 1 and 11 years. Conclusions: Combining highly effective vector control with etiological (trypanocidal) treatment in humans would substantially reduce time to transmission interruption as well as infection incidence and prevalence. However, the success of vector control may depend on prevailing vector species. It will be crucial to improve the coverage of screening programs, the performance of diagnostic tests, the proportion of people treated, and the efficacy of trypanocidal drugs. While screening and access can be incremented as part of strengthening the health systems response, improving diagnostics performance and drug efficacy will require further research.

Risk factors associated with Chagas disease in pregnant women in Santander, a highly endemic Colombian area.

OBJECTIVE: To determine the prevalence and risk factors associated with Chagas disease in pregnant women in an endemic area of Santander, Colombia. METHODS: Cross-sectional study included 23 municipalities of Santander, Colombia. Serological IFAT and ELISA tests were undertaken to detect IgG anti- Trypanosoma cruzi. A questionnaire was conducted for assessing the risk factors of each participant. Newborns were evaluated at birth and followed up to 1 year of age to determine congenital infection. RESULTS: An overall prevalence of 3.2% (95% CI 2.4-4.2) among 1518 pregnant women was detected. Prevalences by provinces were as follows: Guanentina: 6.0% (95% CI 4.1-8.5), García Rovira: 2.9% (95% CI: 1.5-4.8) and Comunera: 0.4% (0.4-2.3). The main risk factors identified were age >32 years old (OR: 2.1; 95% CI: 1.1-3.9); currently having a thatched roof (OR: 11.8; CI95% 2.2-63.2) and a thatched roof during childhood (OR: 3.0; 95% CI: 1.4-6.6); having below primary school education level (OR: 4.6; 95% CI: 2.2-9.5); and a history of a close contact with the vector (triatomine bugs) at least once during their lifetime (OR: 6.9; 95% CI: 3.7-12.9). No congenital cases were detected by parasitological or serological techniques. CONCLUSIONS: Prevalence of Chagas disease in pregnant women is a potential source of infection in this Colombian endemic area. The main risk factors associated with seropositivity were related to conditions favouring the contact with the vector. The results show that it is necessary to continue an active surveillance in order to offer diagnosis and treatment to mothers and their newborns in addition to screening to pregnant women from endemic areas.

Epidemiological findings, estimates of the instantaneous reproduction number, and control strategies of the first Mpox outbreak in Latin America.

BACKGROUND: The 2022-2023 period marked the largest global Mpox outbreak, with Latin America's situation notably underexplored. This study aims to estimate Mpox's instantaneous reproduction number (R(t)), analyze epidemiological trends, and map vaccination efforts in six Latin American countries. METHODS: Utilizing Pan American Health Organization Mpox surveillance data, we examined demographic characteristics, cumulative incidence rates, and epidemic curves, calculated R(t) with weekly sliding windows for each country, alongside a review of vaccination initiatives. RESULTS: From 2022 to 2023, 25,503 Mpox cases and 71 deaths were reported across Argentina, Brazil, Chile, Colombia, Mexico and Peru, with a significant majority (91.8%-98.5%) affecting men, with a mean age of 32-35 years. Maximum R(t) values varied across countries: Argentina (2.63; 0.85 to 5.39), Brazil (3.13; 2.61 to 3.69), Chile (2.91; 1.55 to 4.70), Colombia (3.15; 2.07 to 4.44), Mexico (2.28; 1.18 to 3.75), and Peru (2.84; 2.33 to 3.40). The epidemic's peak occurred between August and September 2022 with R(t) values subsequently dropping below 1. From November 2022, and as of February 2024, only Chile, Peru, and Brazil had initiated Mpox vaccination campaigns, with Colombia launching a Clinical Trial. CONCLUSION: The peak of the Mpox epidemic in the studied countries occurred before the commencement of vaccination programs. This trend may be then partly attributed to a combination of behavioral modifications in key affected communities and contact tracing local programs. Therefore, the proportion of the at-risk population that remains susceptible is still uncertain, highlighting the need for continued surveillance and evaluation of vaccination strategies.

Correction: The epidemiology of Mayaro virus in the Americas: A systematic review and key parameter estimates for outbreak modelling.

[This corrects the article DOI: 10.1371/journal.pntd.0009418.].

From serological surveys to disease burden: a modelling pipeline for Chagas disease.

In 2012, the World Health Organization (WHO) set the elimination of Chagas disease intradomiciliary vectorial transmission as a goal by 2020. After a decade, some progress has been made, but the new 2021-2030 WHO roadmap has set even more ambitious targets. Innovative and robust modelling methods are required to monitor progress towards these goals. We present a modelling pipeline using local seroprevalence data to obtain national disease burden estimates by disease stage. Firstly, local seroprevalence information is used to estimate spatio-temporal trends in the Force-of-Infection (FoI). FoI estimates are then used to predict such trends across larger and fine-scale geographical areas. Finally, predicted FoI values are used to estimate disease burden based on a disease progression model. Using Colombia as a case study, we estimated that the number of infected people would reach 506 000 (95% credible interval (CrI) = 395 000-648 000) in 2020 with a 1.0% (95%CrI = 0.8-1.3%) prevalence in the general population and 2400 (95%CrI = 1900-3400) deaths (approx. 0.5% of those infected). The interplay between a decrease in infection exposure (FoI and relative proportion of acute cases) was overcompensated by a large increase in population size and gradual population ageing, leading to an increase in the absolute number of Chagas disease cases over time. This article is part of the theme issue 'Challenges and opportunities in the fight against neglected tropical diseases: a decade from the London Declaration on NTDs'.

Designing school reopening in the COVID-19 pre-vaccination period in Bogotá, Colombia: A modeling study.

The COVID-19 pandemic has affected millions of people around the world. In Colombia, 1.65 million cases and 43,495 deaths were reported in 2020. Schools were closed in many places around the world to slow down the spread of SARS-CoV-2. In Bogotá, Colombia, most of the public schools were closed from March 2020 until the end of the year. School closures can exacerbate poverty, particularly in low- and middle-income countries. To reconcile these two priorities in health and fighting poverty, we estimated the impact of school reopening for in-person instruction in 2021. We used an agent-based model of SARS-CoV-2 transmission calibrated to the daily number of deaths. The model includes schools that represent private and public schools in terms of age, enrollment, location, and size. We simulated school reopening at different capacities, assuming a high level of face-mask use, and evaluated the impact on the number of deaths in the city. We also evaluated the impact of reopening schools based on grade and multidimensional poverty index. We found that school at 35% capacity, assuming face-mask adherence at 75% in>8 years of age, had a small impact on the number of deaths reported in the city during a third wave. The increase in deaths was smallest when only pre-kinder was opened, and largest when secondary school was opened. At larger capacities, the impact on the number of deaths of opening pre-kinder was below 10%. In contrast, reopening other grades above 50% capacity substantially increased the number of deaths. Reopening schools based on their multidimensional poverty index resulted in a similar impact, irrespective of the level of poverty of the schools that were reopened. The impact of schools reopening was lower for pre-kinder grades and the magnitude of additional deaths associated with school reopening can be minimized by adjusting capacity in older grades.

Estimating Zika virus attack rates and risk of Zika virus-associated neurological complications in Colombian capital cities with a Bayesian model.

Zika virus (ZIKV) is a mosquito-borne pathogen that caused a major epidemic in the Americas in 2015-2017. Although the majority of ZIKV infections are asymptomatic, the virus has been associated with congenital birth defects and neurological complications (NC) in adults. We combined multiple data sources to improve estimates of ZIKV infection attack rates (IARs), reporting rates of Zika virus disease (ZVD) and the risk of ZIKV-associated NC for 28 capital cities in Colombia. ZVD surveillance data were combined with post-epidemic seroprevalence data and a dataset on ZIKV-associated NC in a Bayesian hierarchical model. We found substantial heterogeneity in ZIKV IARs across cities. The overall estimated ZIKV IAR across the 28 cities was 0.38 (95% CrI: 0.17-0.92). The estimated ZVD reporting rate was 0.013 (95% CrI: 0.004-0.024), and 0.51 (95% CrI: 0.17-0.92) cases of ZIKV-associated NC were estimated to be reported per 10 000 ZIKV infections. When we assumed the same ZIKV IAR across sex or age group, we found important spatial heterogeneities in ZVD reporting rates and the risk of being reported as a ZVD case with NC. Our results highlight how additional data sources can be used to overcome biases in surveillance data and estimate key epidemiological parameters.

Linear and Machine Learning modelling for spatiotemporal disease predictions: Force-of-Infection of Chagas disease.

BACKGROUND: Chagas disease is a long-lasting disease with a prolonged asymptomatic period. Cumulative indices of infection such as prevalence do not shed light on the current epidemiological situation, as they integrate infection over long periods. Instead, metrics such as the Force-of-Infection (FoI) provide information about the rate at which susceptible people become infected and permit sharper inference about temporal changes in infection rates. FoI is estimated by fitting (catalytic) models to available age-stratified serological (ground-truth) data. Predictive FoI modelling frameworks are then used to understand spatial and temporal trends indicative of heterogeneity in transmission and changes effected by control interventions. Ideally, these frameworks should be able to propagate uncertainty and handle spatiotemporal issues. METHODOLOGY/PRINCIPAL FINDINGS: We compare three methods in their ability to propagate uncertainty and provide reliable estimates of FoI for Chagas disease in Colombia as a case study: two Machine Learning (ML) methods (Boosted Regression Trees (BRT) and Random Forest (RF)), and a Linear Model (LM) framework that we had developed previously. Our analyses show consistent results between the three modelling methods under scrutiny. The predictors (explanatory variables) selected, as well as the location of the most uncertain FoI values, were coherent across frameworks. RF was faster than BRT and LM, and provided estimates with fewer extreme values when extrapolating to areas where no ground-truth data were available. However, BRT and RF were less efficient at propagating uncertainty. CONCLUSIONS/SIGNIFICANCE: The choice of FoI predictive models will depend on the objectives of the analysis. ML methods will help characterise the mean behaviour of the estimates, while LM will provide insight into the uncertainty surrounding such estimates. Our approach can be extended to the modelling of FoI patterns in other Chagas disease-endemic countries and to other infectious diseases for which serosurveys are regularly conducted for surveillance.

Global variation in force-of-infection trends for human Taenia solium taeniasis/cysticercosis.

Infection by Taenia solium poses a major burden across endemic countries. The World Health Organization (WHO) 2021-2030 Neglected Tropical Diseases roadmap has proposed that 30% of endemic countries achieve intensified T. solium control in hyperendemic areas by 2030. Understanding geographical variation in age-prevalence profiles and force-of-infection (FoI) estimates will inform intervention designs across settings. Human taeniasis (HTT) and human cysticercosis (HCC) age-prevalence data from 16 studies in Latin America, Africa, and Asia were extracted through a systematic review. Catalytic models, incorporating diagnostic performance uncertainty, were fitted to the data using Bayesian methods, to estimate rates of antibody (Ab)-seroconversion, infection acquisition and Ab-seroreversion or infection loss. HCC FoI and Ab-seroreversion rates were also estimated across 23 departments in Colombia from 28,100 individuals. Across settings, there was extensive variation in all-ages seroprevalence. Evidence for Ab-seroreversion or infection loss was found in most settings for both HTT and HCC and for HCC Ab-seroreversion in Colombia. The average duration until humans became Ab-seropositive/infected decreased as all-age (sero)prevalence increased. There was no clear relationship between the average duration humans remain Ab-seropositive and all-age seroprevalence. Marked geographical heterogeneity in T. solium transmission rates indicate the need for setting-specific intervention strategies to achieve the WHO goals.

Estimating the COVID-19 infection fatality ratio accounting for seroreversion using statistical modelling.

Background: The infection fatality ratio (IFR) is a key statistic for estimating the burden of coronavirus disease 2019 (COVID-19) and has been continuously debated throughout the COVID-19 pandemic. The age-specific IFR can be quantified using antibody surveys to estimate total infections, but requires consideration of delay-distributions from time from infection to seroconversion, time to death, and time to seroreversion (i.e. antibody waning) alongside serologic test sensitivity and specificity. Previous IFR estimates have not fully propagated uncertainty or accounted for these potential biases, particularly seroreversion. Methods: We built a Bayesian statistical model that incorporates these factors and applied this model to simulated data and 10 serologic studies from different countries. Results: We demonstrate that seroreversion becomes a crucial factor as time accrues but is less important during first-wave, short-term dynamics. We additionally show that disaggregating surveys by regions with higher versus lower disease burden can inform serologic test specificity estimates. The overall IFR in each setting was estimated at 0.49-2.53%. Conclusion: We developed a robust statistical framework to account for full uncertainties in the parameters determining IFR. We provide code for others to apply these methods to further datasets and future epidemics.

Impact of contact tracing on COVID-19 mortality: An impact evaluation using surveillance data from Colombia.

BACKGROUND: Contact tracing is a crucial part of the public health surveillance toolkit. However, it is labor-intensive and costly to carry it out. Some countries have faced challenges implementing contact tracing, and no impact evaluations using empirical data have assessed its impact on COVID-19 mortality. This study assesses the impact of contact tracing in a middle-income country, providing data to support the expansion and optimization of contact tracing strategies to improve infection control. METHODS: We obtained publicly available data on all confirmed COVID-19 cases in Colombia between March 2 and June 16, 2020. (N = 54,931 cases over 135 days of observation). As suggested by WHO guidelines, we proxied contact tracing performance as the proportion of cases identified through contact tracing out of all cases identified. We calculated the daily proportion of cases identified through contact tracing across 37 geographical units (32 departments and five districts). Further, we used a sequential log-log fixed-effects model to estimate the 21-days, 28-days, 42-days, and 56-days lagged impact of the proportion of cases identified through contact tracing on daily COVID-19 mortality. Both the proportion of cases identified through contact tracing and the daily number of COVID-19 deaths are smoothed using 7-day moving averages. Models control for the prevalence of active cases, second-degree polynomials, and mobility indices. Robustness checks to include supply-side variables were performed. RESULTS: We found that a 10 percent increase in the proportion of cases identified through contact tracing is related to COVID-19 mortality reductions between 0.8% and 3.4%. Our models explain between 47%-70% of the variance in mortality. Results are robust to changes of specification and inclusion of supply-side variables. CONCLUSION: Contact tracing is instrumental in containing infectious diseases. Its prioritization as a surveillance strategy will substantially impact reducing deaths while minimizing the impact on the fragile economic systems of lower and middle-income countries. This study provides lessons for other LMIC.