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.

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'.

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.

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.

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.

Spatial distribution and risk factors for human cysticercosis in Colombia.

BACKGROUND: Cysticercosis is a zoonotic neglected tropical disease (NTD) that affects humans and pigs following the ingestion of Taenia solium eggs. Human cysticercosis poses a substantial public health burden in endemic countries. The World Health Organization (WHO) aims to target high-endemicity settings with enhanced interventions in 17 countries by 2030. Between 2008 and 2010, Colombia undertook a national baseline serosurvey of unprecedented scale, which led to an estimated seroprevalence of T. solium cysticercus antibodies among the general population of 8.6%. Here, we use contemporary geostatistical approaches to analyse this unique dataset with the aim of understanding the spatial distribution and risk factors associated with human cysticercosis in Colombia to inform how best to target intervention strategies. METHODS: We used a geostatistical model to estimate individual and household risk factors associated with seropositivity to T. solium cysticercus antibodies from 29,253 people from 133 municipalities in Colombia. We used both independent and spatially structured random effects at neighbourhood/village and municipality levels to account for potential clustering of exposure to T. solium. We present estimates of the distribution and residual correlation of seropositivity at the municipality level. RESULTS: High seroprevalence was identified in municipalities located in the north and south of Colombia, with spatial correlation in seropositivity estimated up to approximately 140 km. Statistically significant risk factors associated with seropositivity to T. solium cysticercus were related to age, sex, educational level, socioeconomic status, use of rainwater, consumption of partially cooked/raw pork meat and possession of dogs. CONCLUSIONS: In Colombia, the distribution of human cysticercosis is influenced by socioeconomic considerations, education and environmental factors related to the spread of T. solium eggs. This information can be used to tailor national intervention strategies, such as targeting spatial hotspots and more highly exposed groups, including displaced people and women. Large-scale seroprevalence surveys accompanied by geospatial mapping are an essential step towards reaching the WHO's 2021‒2030 NTD roadmap targets.

COVID-19 and public transport: an overview and recommendations applicable to Latin America / COVID-19 y transporte público: una revisión y recomendaciones aplicables a América Latina

Abstract The article presents a general overview on COVID-19 transmission in the context of public transport, particularly applicable to decision making in Latin America. Based on recent findings on COVID-19 transmission and the relative importance of each factor (droplets, fomites, and aerosol routes) in such transmission, we seek to update the discussion on the topic that has generally been based on social distance as the only parameter for reducing the risk of transmission and broadens the vision to integrate ventilation, users' behavior (mask and eye protection use, silence while in the transport system) and travel distance. Recommendations to improve mobility conditions reducing the risk of COVID-19 contagion are provided.

Resumen El artículo presenta una revisión de transmisión de COVID-19 en el contexto de transporte público, con aplicación particular para toma de decisiones en América Latina. Con base en los hallazgos recientes sobre transmisión de Covid-19 y la importancia relativa de cada factor (gotículas, fómites y rutas de aerosoles) en dicha transmisión, buscamos actualizar la discusión sobre el tema que generalmente se ha basado en la "distancia social" como parámetro único de reducción de riesgo de transmisión y amplía esta visión para integrar la ventilación, el comportamiento de usuarios (uso de mascarilla, protección ocular, silencio), y la distancia de viaje. Se indican al final recomendaciones para mejorar las condiciones de movilidad en general sin aumentar el riesgo de contagio de Covid-19.

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.

Descriptive analysis of surveillance data for Zika virus disease and Zika virus-associated neurological complications in Colombia, 2015-2017.

Zika virus (ZIKV) is a mosquito-borne pathogen that recently caused a major epidemic in the Americas. Although the majority of ZIKV infections are asymptomatic, the virus has been associated with birth defects in fetuses and newborns of infected mothers as well as neurological complications in adults. We performed a descriptive analysis on approximately 106,000 suspected and laboratory-confirmed cases of Zika virus disease (ZVD) that were reported during the 2015-2017 epidemic in Colombia. We also analyzed a dataset containing patients with neurological complications and recent febrile illness compatible with ZVD. Females had higher cumulative incidence of ZVD than males. Compared to the general population, cases were more likely to be reported in young adults (20 to 39 years of age). We estimated the cumulative incidence of ZVD in pregnant females at 3,120 reported cases per 100,000 population (95% CI: 3,077-3,164), which was considerably higher than the incidence in both males and non-pregnant females. ZVD cases were reported in all 32 departments. Four-hundred and eighteen patients suffered from ZIKV-associated neurological complications, of which 85% were diagnosed with Guillain-Barré syndrome. The median age of ZIKV cases with neurological complications was 12 years older than that of ZVD cases. ZIKV-associated neurological complications increased with age, and the highest incidence was reported among individuals aged 75 and older. Even though neurological complications and deaths due to ZIKV were rare in this epidemic, better risk communication is needed for people living in or traveling to ZIKV-affected areas.

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

Mayaro virus (MAYV) is an arbovirus that is endemic to tropical forests in Central and South America, particularly within the Amazon basin. In recent years, concern has increased regarding MAYV's ability to invade urban areas and cause epidemics across the region. We conducted a systematic literature review to characterise the evolutionary history of MAYV, its transmission potential, and exposure patterns to the virus. We analysed data from the literature on MAYV infection to produce estimates of key epidemiological parameters, including the generation time and the basic reproduction number, R0. We also estimated the force-of-infection (FOI) in epidemic and endemic settings. Seventy-six publications met our inclusion criteria. Evidence of MAYV infection in humans, animals, or vectors was reported in 14 Latin American countries. Nine countries reported evidence of acute infection in humans confirmed by viral isolation or reverse transcription-PCR (RT-PCR). We identified at least five MAYV outbreaks. Seroprevalence from population based cross-sectional studies ranged from 21% to 72%. The estimated mean generation time of MAYV was 15.2 days (95% CrI: 11.7-19.8) with a standard deviation of 6.3 days (95% CrI: 4.2-9.5). The per-capita risk of MAYV infection (FOI) ranged between 0.01 and 0.05 per year. The mean R0 estimates ranged between 2.1 and 2.9 in the Amazon basin areas and between 1.1 and 1.3 in the regions outside of the Amazon basin. Although MAYV has been identified in urban vectors, there is not yet evidence of sustained urban transmission. MAYV's enzootic cycle could become established in forested areas within cities similar to yellow fever virus.

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.

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.

Endemic and Epidemic Human Alphavirus Infections in Eastern Panama: An Analysis of Population-Based Cross-Sectional Surveys.

Madariaga virus (MADV) has recently been associated with severe human disease in Panama, where the closely related Venezuelan equine encephalitis virus (VEEV) also circulates. In June 2017, a fatal MADV infection was confirmed in a community of Darien Province. We conducted a cross-sectional outbreak investigation with human and mosquito collections in July 2017, where sera were tested for alphavirus antibodies and viral RNA. In addition, by applying a catalytic, force-of-infection (FOI) statistical model to two serosurveys from Darien Province in 2012 and 2017, we investigated whether endemic or epidemic alphavirus transmission occurred historically. In 2017, MADV and VEEV IgM seroprevalences were 1.6% and 4.4%, respectively; IgG antibody prevalences were MADV: 13.2%, VEEV: 16.8%, Una virus (UNAV): 16.0%, and Mayaro virus: 1.1%. Active viral circulation was not detected. Evidence of MADV and UNAV infection was found near households, raising questions about its vectors and enzootic transmission cycles. Insomnia was associated with MADV and VEEV infections, depression symptoms were associated with MADV, and dizziness with VEEV and UNAV. Force-of-infection analyses suggest endemic alphavirus transmission historically, with recent increased human exposure to MADV and VEEV in Aruza and Mercadeo, respectively. The lack of additional neurological cases suggests that severe MADV and VEEV infections occur only rarely. Our results indicate that over the past five decades, alphavirus infections have occurred at low levels in eastern Panama, but that MADV and VEEV infections have recently increased-potentially during the past decade. Endemic infections and outbreaks of MADV and VEEV appear to differ spatially in some locations of eastern Panama.

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.

Heterogeneity of Trypanosoma cruzi infection rates in vectors and animal reservoirs in Colombia: a systematic review and meta-analysis.

BACKGROUND: The heterogeneity of Trypanosoma cruzi infection rates among triatomines insects and animal reservoirs has been studied in independent studies, but little information has been systematised to allow pooled and comparative estimates. Unravelling the main patterns of this heterogeneity could contribute to a further understanding of T. cruzi transmission in Colombia. METHODS: A systematic search was conducted in PubMed, Medline, LILACS, Embase, Web of Knowledge, Google Scholar and secondary sources with no filters of language or time and until April 2018. Based on selection criteria, all relevant studies reporting T. cruzi infection rates in reservoirs or triatomines were chosen. For pooled analyses, a random effects model for binomial distribution was used. Heterogeneity among studies is reported as I2. Subgroup analyses included: taxonomic classification, ecotope and diagnostic methods. Publication bias and sensitivity analyses were performed. RESULTS: Overall, 39 studies reporting infection rates in Colombia were found (22 for potential reservoirs and 28 for triatomine insects) for a total sample of 22,838 potential animals and 11,307 triatomines evaluated for T. cruzi infection. We have found evidence of 38/71 different animal species as potential T. cruzi reservoirs and 14/18 species as triatomine vectors for T. cruzi. Among animals, the species with the highest pooled prevalence were opossum (Didelphis marsupialis) with 48.0% (95% CI: 26-71%; I2 = 88%, τ2 = 0.07, P < 0.01) and domestic dog (Canis lupus familiaris) with 22.0% (95% CI: 4-48%; I2 = 96%, τ2 = 0.01, P < 0.01). Among triatomines, the highest prevalence was found for Triatoma maculata in the peridomestic ecotope (68.0%, 95% CI: 62-74%; I2 = 0%, τ2 = 0, P < 0.0001), followed by Rhodnius prolixus (62.0%, 95% CI: 38-84%; I2 = 95%, τ2 = 0.05, P < 0.01) and Rhodnius pallescens (54.0%, 95% CI: 37-71%; I2 = 86%, τ2 = 0.035, P < 0.01) in the sylvatic ecotope. CONCLUSIONS: To our knowledge, this is the first systematic and quantitative analyses of triatomine insects and potential animal reservoirs for T. cruzi infection in Colombia. The results highlight a marked heterogeneity between species and provide initial estimates of infection rates heterogeneity.

Estimating spatiotemporally varying malaria reproduction numbers in a near elimination setting.

In 2016 the World Health Organization identified 21 countries that could eliminate malaria by 2020. Monitoring progress towards this goal requires tracking ongoing transmission. Here we develop methods that estimate individual reproduction numbers and their variation through time and space. Individual reproduction numbers, Rc, describe the state of transmission at a point in time and differ from mean reproduction numbers, which are averages of the number of people infected by a typical case. We assess elimination progress in El Salvador using data for confirmed cases of malaria from 2010 to 2016. Our results demonstrate that whilst the average number of secondary malaria cases was below one (0.61, 95% CI 0.55-0.65), individual reproduction numbers often exceeded one. We estimate a decline in Rc between 2010 and 2016. However we also show that if importation is maintained at the same rate, the country may not achieve malaria elimination by 2020.