Spatial patterns and climate drivers of malaria in three border areas of Brazil, Venezuela and Guyana, 2016-2018

In 2020, 77% of malaria cases in the Americas were concentrated in Venezuela, Brazil, and Colombia. These countries are characterized by a heterogeneous malaria landscape and malaria hotspots. Furthermore, the political unrest in Venezuela has led to significant cross-border population movement. Hence, the aim of this study was to describe spatial patterns and identify significant climatic drivers of malaria transmission along the Venezuela-Brazil-Guyana border, focusing on Bolivar state, Venezuela and Roraima state, Brazil. Malaria case data, stratified by species from 2016 to 2018, were obtained from the Brazilian Malaria Epidemiology Surveillance Information System, the Guyana Vector Borne Diseases Program, the Venezuelan Ministry of Health, and civil society organizations. Spatial autocorrelation in malaria incidence was explored using Getis-Ord (Gi*) statistics. A Poisson regression model was developed with a conditional autoregressive prior structure and posterior parameters were estimated using the Bayesian Markov chain Monte Carlo simulation with Gibbs sampling. There were 685,498 malaria cases during the study period. Plasmodium vivax was the predominant species (71.7%, 490,861). Malaria hotspots were located in eight municipalities along the Venezuela and Guyana international borders with Brazil. Plasmodium falciparum increased by 2.6% (95% credible interval [CrI] 2.1%, 2.8%) for one meter increase in altitude, decreased by 1.6% (95% CrI 1.5%, 2.3%) and 0.9% (95% CrI 0.7%, 2.4%) per 1 cm increase in 6-month lagged precipitation and each 1 °C increase of minimum temperature without lag. Each 1 °C increase of 1-month lagged maximum temperature increased P. falciparum by 0.6% (95% CrI 0.4%, 1.9%). P. vivax cases increased by 1.5% (95% CrI 1.3%, 1.6%) for one meter increase in altitude and decreased by 1.1% (95% CrI 1.0%, 1.2%) and 7.3% (95% CrI 6.7%, 9.7%) for each 1 cm increase of precipitation lagged at 6-months and 1 °C increase in minimum temperature lagged at 6-months. Each 1°C increase of two-month lagged maximum temperature increased P. vivax by 1.5% (95% CrI 0.6%, 7.1%). There was no significant residual spatial clustering after accounting for climatic covariates. Malaria hotspots were located along the Venezuela and Guyana international border with Roraima state, Brazil. In addition to population movement, climatic variables were important drivers of malaria transmission in these areas.