https://orcid.org/0000-0003-2922-6751
Figures
Abstract
Malaria elimination in Latin America is becoming an elusive goal. Malaria cases reached a historical ~1 million in 2017 and 2018, with Venezuela contributing 53% and 51% of those cases, respectively. Historically, malaria incidence in southern Venezuela has accounted for most of the country's total number of cases. The efficient deployment of disease prevention measures and prediction of disease spread to new regions requires an in-depth understanding of spatial heterogeneity on malaria transmission dynamics. Herein, we characterized the spatial epidemiology of malaria in southern Venezuela from 2007 through 2017 and described the extent to which malaria distribution has changed country-wide over the recent years. We found that disease transmission was focal and more prevalent in the southeast region of southern Venezuela where two persistent hotspots of Plasmodium vivax (76%) and P. falciparum (18%) accounted for ~60% of the total number of cases. Such hotspots are linked to deforestation as a consequence of illegal gold mining activities. Incidence has increased nearly tenfold over the last decade, showing an explosive epidemic growth due to a significant lack of disease control programs. Our findings highlight the importance of spatially oriented interventions to contain the ongoing malaria epidemic in Venezuela. This work also provides baseline epidemiological data to assess cross-border malaria dynamics and advocates for innovative control efforts in the Latin American region.
Author summary
Despite major progress in malaria control or elimination in South America in the last years, this disease is still an important public health concern in all the Amazon basin countries, particularly during 2017 and 2018, when a substantial malaria incidence (~ one million) was reported. Most of the rise in cases has been due to increases in malaria transmission in Venezuela. Within this country, populations in the southeastern (Guiana Shield) region have been at the highest risk of Plasmodium vivax and P. falciparum infection. In this study, reported malaria cases between 2007–2017 in this region were analyzed to stratify the spatial risk of disease and identify the malaria surge's main driver over the last years. Malaria risk was highly focal and prevalent in 2 consistent disease hotspots that accounted for most Plasmodium transmission in the whole region. Illegal gold-mining activities seem to drive malaria in these disease pockets and seem critical in malaria's surge throughout the country. Analyses of the Venezuelan surveillance data showed that from 2014 onwards, local malaria transmission has reemerged in new areas of the country, leading to a shift in this disease's epidemiology. Successful control of Venezuela's ongoing malaria epidemic requires hotspot-targeted control at the national level and regional coordination to avoid cross-border malaria spillover.
Citation: Grillet ME, Moreno JE, Hernández-Villena JV, Vincenti-González MF, Noya O, Tami A, et al. (2021) Malaria in Southern Venezuela: The hottest hotspot in Latin America. PLoS Negl Trop Dis 15(1): e0008211. https://doi.org/10.1371/journal.pntd.0008211
Editor: Rhoel Ramos Dinglasan, University of Florida, UNITED STATES
Received: March 11, 2020; Accepted: November 25, 2020; Published: January 25, 2021
This is an open access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication.
Data Availability: Data of this publication are publicly available online from https://figshare.com/account/home, accession number https://doi.org/10.6084/m9.figshare.12749426.v1.
Funding: Funding was provided by the Council for Sciences and Humanities Development (Grant CDCH-PG-0382182011 to MEG and JEM) and by the National Institutes of Health, USA, Grant 2R01AI110112 to JEC and Grant U19 AI089681 to AAE. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist
Introduction
In the Americas, malaria control programs face significant obstacles that threaten the global strategy for malaria elimination set for the region for 2030 [1–3]. Progress made between 2000–2015 towards reducing malaria morbidity and mortality stalled in 2016 and reversed in 2017–2018 [3–5] when cases escalated to ~1 million in each year [2,3]. Since 2014, the Venezuelan economy's collapse, and with it, its entire healthcare system, have been a significant driver for this surge [6,7]. In particular, Venezuela reported a total of 1,255,299 cases between 2015–2018, the highest number in Latin America, and in 2017 exhibited one of the most substantial increases in malaria cases worldwide [1,2,8].
As in other endemic countries in South America [5], Plasmodium vivax accounts for the majority of reported cases (76%) in Venezuela followed by Plasmodium falciparum (17.7%), mixed P.vivax /P. falciparum infections (6%) and Plasmodium malariae <1% [3]. Venezuela successfully achieved malaria elimination in approximately 75% of its territory during the early 1960s [9]. However, low to moderate malaria transmission by P. vivax and P. falciparum persisted in the lowland Amazon rainforests and savannas of the remote Guayana region, mainly in isolated Amerindians communities, South of the Orinoco River [10]. During the 1980s, P. vivax malaria reemerged along the coastal wetlands in the country's northeastern region [11], but its transmission was interrupted twenty years later [12]. However, malaria transmission in Bolívar state, a region in southeastern Venezuela bordering Brazil and Guyana, has persisted. Since 1990, this region has contributed > 60% (1992–1995) to 88% (2000–2014) of Venezuela’s total malaria cases [12,13]. Gold mining has been associated with the high malaria incidence in specific municipalities of this region where gold miners seem to account for an estimated 47% to 80% of malaria cases [4,5,13]. Clearing forest for mining activities is observed to create favorable conditions for the main vector malaria species in this area [5,7,13]. The health situation in Bolívar state has worsened significantly in recent years, and current limitations in malaria surveillance and control are a cause of major concern [7,13]. Additionally, since 2014, local malaria transmission has reemerged in previously endemic and new areas across the country, including Venezuela's northeastern region [7].
At present, the uncontrolled upsurge of malaria, coupled with Venezuelans' mass migration to neighboring countries, poses a serious threat to the broader region [6,7,14], jeopardizing efforts in bordering malaria-endemic countries to achieve their goals for disease control. Human population movement among countries included in the Guiana Shield (Guyana, French Guiana, Suriname, and parts of Colombia, Venezuela, and Brazil) serve to amplify the regional impact of malaria increase in the state of Bolívar. In addition, there is a latent risk that a regional malaria corridor will form [15] from Bolívar state to northern Brazil via the mass migration of displaced individuals who take advantage of the existing roads and trails connecting Venezuela and Brazil (called Route 10). This situation is worsened by the presence of novel mutations linked to artemisinin resistance in Guyana [16], raising concerns about the potential regional spread of mutations linked to drug-resistant and artemisinin delayed parasite clearance in P. falciparum where Venezuela could act as a hub [17].
Considering all these factors, a first step to better understand how to manage this public health crisis is to characterize malaria transmission's spatiotemporal dynamics in southern Venezuela. Heterogeneity and clustering of cases within an endemic focus (an area with suitable environmental and epidemiological characteristics that sustain transmission) is an emerging pattern common to mosquito borne-infectious diseases [18]. As part of such spatial heterogeneity, there could be hotspots (areas where transmission intensity exceeds the regional average) responsible for most cases in a focus. Such areas play a critical role in malaria's resilience to interventions because they frequently maintain persistent transmission, acting as a source of infection for the entire region [18]. Thus, we analyzed reported (passive surveillance) malaria cases data (P. vivax and P. falciparum) from 2007 through 2017 in Bolivar state. The study aimed to identify critical hotspots and determine their transient or persistence dynamics, drivers, and role in the malaria surge observed over the last few years. Particularly, we expected that illegal gold mining activities disproportionally contributed to the transmission of both parasites in the detected hotspots, a socioeconomic malaria risk factor that has been suggested in previous studies [7,13,19]. Finally, we described the spatial epidemiology of malaria at a national level and assess the extent to which landscape epidemiology of this infection has changed throughout the 2014–2017 period.
Method
Study area
Venezuela, located on the northern coast of South America, has a surface area of ~916,445 km2 (Fig 1A). Bolívar state (Fig 1B), the focus of this study, covers an extension of 240,500 km2 and as of the 2018 Venezuelan census [20], has a population of 1,837,485 inhabitants heterogeneously distributed across 11 municipalities and 46 parishes. This region has a low average population density (5.1 inhabitants per km2) and extensive territory. Population densities per municipality are variable and range from 553 (e.g., Caroni municipality) to < 1 inhabitant per km2 (e.g., Gran Sabana municipality). Most of the population resides in the northernmost and easternmost parts of the state (Fig 1B), where economic activities are related to mining (iron, gold, bauxite, and diamonds), production of steel, aluminum, and hydroelectric industries, businesses and services, forestry, cattle, and agricultural development. The state capital is Ciudad Bolivar (Heres municipality) but Ciudad Guayana (Caroni municipality) is the main urban center. Although these northern cities are well-connected with the rest of Venezuela, the most southern and western rural areas are sparsely populated with few or no connecting roads. There is one major paved road, Route 10, that runs along the eastern side of Bolívar and connects Venezuela with Brazil, crossing through the state's main endemic malaria region (Fig 1B). In this region, malaria does not exhibit a seasonal pattern, particularly for P. vivax, which is characterized by a near perennial incidence [12,13]. The primary vectors for all Plasmodium species in the state are Nyssorhynchus darlingi (also referred to as Anopheles darlingi) and Ny. albitarsis s.l. (also referred to as Anopheles albitarsis) [19,21,22,23].
A) Map of Venezuela, South America, and the study area (Bolívar state in green). B) Map of Bolívar state showing the population density by civil parish (people per square kilometer) and other relevant features of the landscape. Maps were created with the Q-GIS software (https://www.qgis.org/es/site/).
Epidemiological and sociodemographic data
The yearly number of P. vivax and P. falciparum cases per parish (2007–2017) in Bolivar state were provided by the Malaria Control Program, a local branch of the Venezuelan Ministry of Health [24]. Time series of annual malaria incidence (2014–2017) per municipality at the country level were obtained from the PAHO Malaria Surveillance Indicators [25]. Malaria incidence rates per 1,000 inhabitants by parish and by Plasmodium species were accounted. These rates were calculated considering the human population growth rate predicted for the studied period according to the demographic data (at-risk population) from the National Statistics Office of Venezuela [20]. In addition, we obtained information (retrospective analysis) from the local (Bolivar state) malaria database [24] of the demographic and epidemiological profiles of all individuals diagnosed with malaria during the studied period (2007–2017). Particularly, we were looking for whether or not miners were well represented in the occupational category of all the Bolívar state patients. We also examined other demographic factors such as age and sex.
Historically, Venezuela's local health services have reported malaria cases detected by passive surveillance on a weekly basis according to geographic origin and symptomatic cases. Case confirmation usually requires clinical and epidemiologic data (e.g., exposure history in endemic areas) plus a positive blood smear. The Venezuelan Ministry of Health (MofH) uninterruptedly has issued weekly and monthly epidemiological reports known as the “Boletín Epidemiológico” reporting on the number of malaria cases, type of infection (parasite species), distribution per state, municipality and parish, distribution per age and sex, occupation of the patient, risk maps, and the number of detected imported cases. All reported malaria cases are notified at the country level via the National Notifiable Diseases Surveillance System (NDSS) by filling out a questionnaire containing relevant epidemiological inquiries. Despite the Venezuela’s health-care crisis, local malaria surveillance and data collection continued to occur and weekly bulletins were made available until 2017 but have been interrupted [7].
Data analyses
We conducted a secondary analysis of passive malaria surveillance data reported by the Venezuelan Minister of Health in Venezuela between 2007 and 2017. Specifically, the annual malaria incidence per parish in Bolivar was mapped to characterize spatio-temporal disease dynamics. Local spatial autocorrelation analysis was performed using Anselin Local Moran's I test to identify malaria's significant spatial patterns [26]. Specifically, we evaluated the likelihood of malaria occurring equally at any location (parish centroid) within Bolivar state vs. the detection of unusual aggregations of malaria incidence. This test evaluates adjacent positions and finds a strong positive spatial autocorrelation when the surrounding incidence of disease over the entire study area has analogous values (called "High-High" areas). A positive value for 'I' indicates that the adjacent parishes are bounded by P. vivax or P. falciparum incidences with similar values. Such a feature is considered a cluster or hotspot. To test whether the observed clustering/dispersing is statistically significant, an estimated Z-score and a 99% level of significance (P < 0.001) were selected.
We also explored occupational risks and demographic factors in hotspot areas by analyzing the proportion of the diagnosed cases during the studied period (2007–2017) that worked primarily in mining, agriculture, forestry, fishing, retail, and other, as well as examining the ratio of male-to-female and age-structure. The chi-squared test allowed us to determine whether there were significant differences between the expected frequencies and the observed frequencies for one or more categories.
To identify illegal mining operations in the study area, we obtained data on the deforested land cover distribution areas in the parishes classified as hotspots across the study. Data on deforestation were obtained from the Global Forest Watch [27] (http://www.globalforestwatch.org) and remote sensing images from Landsat Thematic Mapper imagery [28] (www.http//earthexplorer.usgs.gov/). Illegal gold mining has been the main cause of deforestation in southern Venezuela in recent years [29]. In order to determine if malaria cases exhibited an increasing trend (distinguished from random behavior) along the years or along the temporal cover tree loss, the best fit line (curve fitting approach) to the time series was applied to aid interpretation of data.
Finally, we used kriging, a local geostatistical interpolation method [30], to generate an estimated continuous surface from the scattered set of points (i.e., municipality centroids) with z-values (i.e., annual parasite incidence in each municipality) to better capture the local spatial variation of malaria spread (spatial risk diffusion map) across the country during the 2014–2017 period. The main assumption of this technique is that the distance or direction between sample points reflects a spatial correlation that can be used to explain the variation of the values of the variable, malaria incidence, on the whole surface. To do that, first, it creates a variogram and covariance function to estimate the statistical dependence (autocorrelation) values. The next step is to fit a model-that, is a continuous function or curve- to the points forming the empirical variogram and predict attribute values (malaria incidence) at unsampled locations. Several variogram models were tested using a cross-validation procedure. The best model was adjusted for any directional spatial trend in our data (anisotropy) in the semi-variogram [30].
Maps, hotspots, descriptive analyses, remote sensing image visualization, and kriging analyses were performed in Q-GIS software (version 2.18.9-Las Palmas de G.C., GNU-General Public License, https://www.qgis.org/es/site/) and the R-software (The R-Development Core Team, http://www.r-project.org).
Results
The total number of accumulated cases of malaria in Venezuela during the sampled decade was 1,207,348 (range: 32,037–411,586), with overall malaria incidence rates (cases/1,000 inhabitants-year) increasing from 5.2 (2007) to 28 (2017; Fig 2). During this period malaria incidence increased nearly 10-fold (2007: 41,749 cases; 2017: 411,586) with a sharp and significant increase since 2014 (Exponential model, R2 = 0.92, P < 0.05). Due to malaria infection, the overall mortality increased nearly 20-fold from 16 (2007) to 312 cases in 2017. Bolívar contributed with ~47% of the total cases in Venezuela during 2017; but during the previous years, this region accounted for 60 to 80% of the reported malaria for the whole country (Fig 2). In Bolívar, 20–30% of malaria cases were due to P. falciparum, with the remaining 70–80% attributed to P. vivax, with a mean ratio of P. vivax/P.falciparum of 3.04 (+/- 0.17 SE). Fig 2 also reveals that the malaria burden due to P.vivax in Bolívar State has had a similar and steady 4.5-fold increment (R2 = 0.93, P < 0.05) since 2014 (from 32,791 to 146,885 cases) with an accumulated 496,847 cases for that time period. A more moderate increase (~2.5-fold, R2 = 0.89, P < 0.05) was detected for malaria morbidity due to P. falciparum (18,127 cases in 2014 to 43,152 in 2017), with 152,589 accumulated cases for the same period (Fig 2).
Annual parasite incidence (API: Number of confirmed malaria cases/1,000 inhabitants, blue bars) in Venezuela (2007–2017).
In Bolívar state, annual P. vivax incidence during 2007–2017 (Figs 3A and 3B and S1) was heterogeneously distributed, with most cases spatially concentrated in the mid-eastern and southern parishes of San Isidro, Dalla Costa and Ikabarú (the latter bordering Brazil), followed by northern Pedro Cova, Santa Barbara, and El Callao Parishes. Few or no cases were reported in the remaining parishes. Parishes with the presence of malaria cases had annual P. vivax incidences as high as 4,672 cases/1,000 inhabitants (2017) and mean annual incidences of up to 1,919 cases/1,000 inhabitants-year (e.g., San Isidro). The temporal dynamics of P. vivax cases in San Isidro were representative of the dynamics of malaria overall in the state Bolívar accounting for ~43.4% of the overall malaria morbidity during the study period. As with P.vivax, there was also substantial spatial heterogeneity for P. falciparum (Figs 3C and 3D and S2). This parasite species had the highest annual incidence with 1,549 cases/1,000 inhabitants (2017) and the highest annual mean of 738 cases /1,000 inhabitants (2017) both reported for San Isidro Parish. Further spatial analyses revealed significant local clustering of P. vivax and P. falciparum incidences between the two parishes of San Isidro and Dalla Costa, providing additional evidence that malaria transmission has focused and persisted in these localities supporting their ‘hotspot’ status within the state of Bolívar over the last decade (Figs 4 and S3 and S4).
Annual P. vivax (A, B) and P. falciparum (C, D) incidence (API: Number of confirmed malaria cases/1,000 inhabitants) per parish in Bolívar state (2007 and 2017), south-eastern Venezuela. Maps were created with the Q-GIS software (https://www.qgis.org/es/site/).
Significant clusters (hotspots) of annual Plasmodium vivax (A, B) and P. falciparum (C, D) incidence (red-red) in Bolívar state (south-eastern Venezuela), during 2007 and 2017. Maps were created with the Q-GIS software (https://www.qgis.org/es/site/).
In San Isidro, throughout the study period P. vivax cases predominantly occurred in men (70%: 145,492/206,455, X2 = 34,608.48, 1df, P < 0.0001, Table 1), mainly affecting age-groups between 21–30 (32%), 11–20 (22%) and 31–40 (20%) years (Table 1 and S5 Fig). From an occupational perspective, most of the individuals who were infected with P. vivax during the study period were linked to gold mining activities (miners) compared with other common occupations including housekeeper, machine operator, and student, among others (62%: 127,898/206,455; range: 56%-68%, X2 = 11,792.08, 1df, P < 0.0001, Table 2). A similar significant association was found for P. falciparum: gold miners accounted for ~ 66% of cases (50,013/75,124; range: 63%-70%, X2 = 8,254.48, 1df, P < 0.0001, Table 2).
The concentration of malaria cases in southern San Isidro coincided with deforestation areas as a result of illegal mining activities. During the study period, most villages from Sifontes Municipality that had >1,000 cases (San Isidro and Dalla Costa parishes) were located in this open (deforested) area (Fig 5A), where the percentage of deforestation (tree cover loss) has dramatically increased in recent years (Fig 5B). Since 2007, San Isidro concurrently lost 3,058 Hectares tree cover (~1.02% decrease) while malaria increased by ~746% (Fig 6A). Similarly, Ikabaru Parish (municipality of Gran Sabana), lost 2,934 Hectares of tree cover (~1.06% decrease) and registered a significant and paralleled increase in malaria cases (Fig 6B).
Left: Sifontes Municipality and its three parishes showing the spatial distribution of all localities and some general landscape features. Right: satellite image (source: Landsat 8, spatial resolution: 15m, date: 12/19/2017) of “Las Claritas” (framed in red window in San Isidro Parish); green triangles: localities with > 1,000 malaria cases across the study period in the deforested mining area; red dots are localities with < 1,000 malaria cases. Maps were created with the Q-GIS software (https://www.qgis.org/es/site/).
Percentage of deforestation (tree cover lost) and accumulated malaria (P. vivax + P. falciparum cases) in Sifontes Municipality (A) and Gran Sabana Municipality (B) in Bolívar State.
Fig 7 depicts both the spatial spread of malaria (annual parasite incidence) across Venezuela that has expanded from southern Guayana toward the northern-central-western areas during 2014–2017 (Epidemiological week-EW 21) and the intensification of disease transmission in the South, an endemic area of sustained concern. Particularly, this spatial prediction analysis emphasizes that the primary high-risk malaria areas and potential sources of parasite dispersal within the country are the hotpots present in the state of Bolívar, followed by the southwestern state of Amazonas. Considering the population growth during that period, the national percentage of individuals living in areas at risk of contracting malaria increased from 34.4% (9,907,708 people) to 50% (15,988,534 people) between 2014 and 2017.
Spatial prediction maps (malaria risk) for 2014 (A), 2015 (B), 2016 (C) and 2017 (D) years derived from the ordinary gaussian kriging interpolation model of the annual parasite incidence (API). Note Year 2017, data available only through Epidemiology Week 21. Maps were created with the ArcGIS desktop software (https://www.esri.com).
Discussion
The risk for malaria in southeastern Venezuela varied widely with most cases reported in the mid-eastern and southern parishes of the state of Bolívar, where we identified two persistent hotspots. Local transmission from these infectious disease pockets accounted for highest malaria transmission in the entire region (~ 61%) and country (> 60%) over time (2007–2017). Both hotspots have been a long-standing regional source of P. vivax and P. falciparum transmission, as suggested by earlier studies [7,13,21]. Our results support previous findings from Venezuela, Brazil and Peru showing that P. vivax malaria spatial heterogeneity is characterized by high-risk localities interspersed with others showing low to moderate risk [11,31–33]. Due to their stability and / or persistence over time, these well characterized hotspots could be predictive of prospective malaria incidence in the surrounding areas as has been found in similar studies [34–36].
Our results highlight that disease patterns at larger spatial scales are driven by a sum of factors acting at local scales [37], such as mosquito ecology (especially larval habitats and host-seeking behavior) and at-risk human population dynamics (e.g., density, distribution, and mobility). In particular, our findings support the hypothesis that illegal gold mining is one of the leading local socioeconomic drivers for malaria in southeastern Venezuela and a major contributing factor to the upsurge of malaria in the recent years. First, we found that areas deforested by gold mining had more clusters of Plasmodium cases compared to areas nearby. Secondly, both P. vivax and P. falciparum increased in incidence (4-8-fold) over time in those regions along with a concomitant decrease of vegetation cover (3-6-fold) resulted from such mining activities. Finally, our results reporting age and sex patterns within the region's malaria patients are compatible with the occupation as a factor determining the greater and prolonged exposure to vector bites of young males in the mining fields compared to girls and women.
Illegal gold mining and the associated deforestation have rapidly increased and expanded in southern Venezuela since 2009 [29], particularly in the municipality of Sifontes, state of Bolívar. Changes in land cover patterns resulting from deforestation can promote the emergence of larval habitats for these Nyssorhynchus (also known as Anopheles) vectors, thereby increasing mosquito abundance [38], vector-host contact [39], survival [40] and consequently transmission risk [41], especially in human settlements located near the forest fringe. Earlier studies in our study area have indicated that the most productive breeding site types for Ny. albitarsis s.l. and to a lesser degree for Ny. darlingi, are abandoned open lagoons or mining dug-outs left after clearing vegetation [42,43]. Thus, a hypothesis to be tested in future studies would be how forest fragmentation by mining activities may influence both vector species' ecology and vectorial capacity. This ecological mechanism should be addressed regardless of the increase in malaria in the area due to an increased population attracted by the gold mining activities.
The malaria epidemic in Venezuela has been fueled by financial constraints for the procurement of malaria control commodities (such as insecticides, drugs, diagnostic supplies, and mosquito nets) and surveillance activities, and lack of provision and implementation of services [6,44]. In particular, the economic crisis has ignited the informal sector, with illegal gold mining featuring as one of the country's fastest-growing shadow economic activities [45]. As a result, migration within the country has increased towards southeastern Venezuela (Sifontes Municipality) where gold mining activities are clustered and where malaria hotspots have been identified [44]. Highly mobile human populations migrate from different regions of the country to these mining areas in search of economic opportunities. Some of these internal migrants return back to the region of origin where viable anopheles vector populations exist, reintroducing malaria to areas where autochthonous transmission had been previously eliminated leading to a shift in the epidemiology of this disease [4,5,7]. Concomitantly, the dismantling of the epidemiological surveillance systems at the national level has unintentionally not prevented the reappearance of malaria across the country [7]
The spatial connectedness of malaria foci related to (infected) individuals' travel behavior has been used to identify so-called sources and sinks of Plasmodium within a transmission network approach [46]. Sources (originating and exporting cases) and sinks (receiving imported cases) of malaria parasites mostly occur in epidemiological settings when humans frequently travel between locations characterized by substantial heterogeneity in malaria transmission [15,47]. Consequently, we hypothesize that a source-sink malaria parasite (metapopulation) dynamic accounts for the current spatial distribution and persistence of this infection in Venezuela. As we showed, since 2014, malaria has spread from the South and northeastern coastal regions to the lowland central savannas and west Andes piedmont ecoregions where malaria had previously been eliminated, increasing the population at-risk to around 50% compared to 34.4% in 2010. Our results also suggest the importance of malaria linked to mobile illegal gold miners and their role in its spread among the Guiana Shield countries [48–51]. Understanding the spatiotemporal variability between Plasmodium and human movement (spatial demography) will help design and implement a strategic malaria control program in Venezuela.
The rapidly increasing malaria burden in Venezuela is affecting neighboring countries, particularly Brazil and Colombia [7]. Migration to Brazil occurs via Route 10, the sole highway that links southern Venezuela and northern Brazil. Unfortunately, this road runs through the hotspots of San Isidro and Dalla Costa. Infected people crossing these borders complicate each country's malaria program by the potential reintroduction of the disease in areas where it had previously been reduced or eliminated; enhancing malaria transmission near border areas; promoting case spillover, and possibly spreading drug resistance alleles across Brazil. Another important risk includes the effects of inadequate treatment that can increase malaria-related mortality. According to the Brazilian Ministry of Health, malaria cases imported from Venezuela into neighboring Roraima State, Brazil have increased from 1,538 (2014) to 4,478 (2018), representing up to 85% of the reported malaria cases of that border state [7,25].
Finally, the spatial distribution of mutations linked to drug-resistant P. falciparum in Venezuela remains poorly understood [17]. However, surveillance for resistance markers is now recommended in some Guiana Shield countries such as Guyana and in the mining areas where artemisinin is available for sale and self-treatment is common [52]. Considering that there are multi-drug resistant Plasmodium falciparum lineages against common anti-malarial drugs in Venezuela [17,52,53], and the recent report of novel mutations linked to the delay of artemisinin efficacy in Guyana [16], it is possible that multi-drug resistant parasite lineages, that also include artemisinin-response delayed mutations, could emerge in the southern Venezuela hotspots.
The study's main limitations are related to the nature of passive malaria surveillance data and the lack of local demographic information to standardize the analysis of occupational risk. Sub-reporting of malaria incidence is presumed due to a limited diagnostic capacity of the parasite detection in the region and the whole country as a result of the Venezuelan Health crisis. The lack of completeness of the surveillance systems was also a critical issue when we analyzed the 2017 malaria data. Although these facts underestimate the total number of malaria infections, incidence data analyzed still gave us insights about the upward trend of malaria transmission in Venezuela and the localization of the main source (hotpots) of P. vivax and P. falciparum cases in recent years as suggested by earlier studies [7,13,21]. Similarly, the gold mining as a significant disease social driver and the land-use change (deforestation) as an important local environmental associated proxy of local malaria increase were such strong patterns that they emerged even with this limited data. The current findings support observations of our previous work [5–7,12,13,17,19,21,22,42, 43]; nevertheless, future studies will be needed to further confirm these patterns.
Conclusions
We show evidence that gold-mining activities seem to drive malaria hotspots in Venezuela and those high transmission pockets were critical in the surge of malaria observed from 2014 onwards. These gold-mining areas likely make malaria transmission resilient to interventions. In particular, they not only sustain transmission but also can restore it ("rescue effect") after interventions have reduced malaria locally or even achieved local elimination in other areas, as has been observed in Venezuela [15,18]. Transmission in such settings is considered the Gordian knot for achieving malaria elimination in Latin America [5,54–56].
Hotspot-targeted control has been hypothesized as a practical approach to reducing the malaria burden in areas of heterogeneous malaria transmission [18,56]. Thus, a program focused on rapid diagnosis and timely treatment, vector control, and monitoring for drug/insecticide resistance is urgent and essential in these hotspots. Malaria surveillance and reporting have been particularly affected by Venezuela's healthcare crisis since mid-2017 onwards [7]. Hence, to recover and strengthen this surveillance system is a necessary step to control malaria in Venezuela. Increased enforcement of the malaria control program in the rest of Venezuela is pivotal to lower the risks of reintroduction to vulnerable areas. Given the current context, successful control of the ongoing malaria epidemic in Venezuela requires national and regional coordination, as evidenced by the cross-border malaria spillover. Without coordinated international efforts, the progress achieved toward malaria elimination in Latin America over the past 18 years could be easily reversed.
Supporting information
S1 Fig. Annual P. vivax incidence (Number of confirmed malaria cases/1,000 inhabitants) per parish in Bolívar state (2008–2016), south-eastern Venezuela.
Maps were created with the Q-GIS software (https://www.qgis.org/es/site/)
https://doi.org/10.1371/journal.pntd.0008211.s001
(TIF)
S2 Fig. Annual P. falciparum incidence (Number of confirmed malaria cases/1,000 inhabitants) per parish in Bolívar state (2008–2016), south-eastern Venezuela.
Maps were created with the Q-GIS software (https://www.qgis.org/es/site/).
https://doi.org/10.1371/journal.pntd.0008211.s002
(TIF)
S3 Fig. Significant clusters (hotspots) of annual Plasmodium vivax incidence (red-red) in Bolívar state (south-eastern Venezuela), from 2008–2016.
Maps were created with the Q-GIS software (https://www.qgis.org/es/site/)
https://doi.org/10.1371/journal.pntd.0008211.s003
(TIF)
S4 Fig. Significant clusters (hotspots) of annual Plasmodium falciparum incidence (red-red) in Bolívar state (south-eastern Venezuela), from 2008–2016.
Maps were created with the Q-GIS software (https://www.qgis.org/es/site/).
https://doi.org/10.1371/journal.pntd.0008211.s004
(TIF)
S5 Fig.
Age distribution of P. vivax (upper panel: A, B) and P. falciparum (lower panel: C, D) among malaria patients of San Isidro (left side: A, C) and Dalla Costa (right side: B, D) parishes (northeastern of Bolivar state) across the study period.
https://doi.org/10.1371/journal.pntd.0008211.s005
(TIF)
Acknowledgments
We acknowledge the logistic support (field and laboratory) provided by Angela Martinez, Porfirio Acevedo, Nelson Moncada, Mayida El Souki, Virginia Behm, Manuel Amarista, and Maria A. Oliveira-Miranda. This work is the result of the collaborative work carried out by the research network for vector-borne diseases surveillance & control in Venezuela & Andean region (VeConVen: https://www.vbdvenezuelanetwork.com/).
References
- 1.
WHO 2017. World malaria report 2017. Geneva: World Health Organization; 2017. Licence: CC BY-NC-SA 3.0 IGO. https://www.who.int/malaria/publications/world-malaria-report-2017/report/en/ (accessed 15 January, 2020).
- 2.
WHO 2018. World Malaria Report 2018. Geneva: World Health Organization; 2018. Licence: CC BY-NC-SA 3.0 IGO. https://www.who.int/malaria/publications/world-malaria-report-2018/report/en/(accessed 15 January, 2020).
- 3.
WHO 2019. World malaria report 2019. Geneva: World Health Organization; 2019. Licence: CC BY-NC-SA 3.0 IGO. https://www.who.int/publications-detail/world-malaria-report-2019 (accessed 15 January, 2020).
- 4. Recht J, Siqueira AM, Monteiro WM, Herrera SM, Herrera S, Lacerda MVG. Malaria in Brazil, Colombia, Peru and Venezuela: current challenges in malaria control and elimination. Malar J. 2017; 16:273. pmid:28676055
- 5.
Conn JE, Grillet ME, Correa M, Sallum MAM. Malaria Transmission in South America- Present Status and Prospects for Elimination. In: Manguin S, Dev V, editors. Towards Malaria Elimination—A Leap Forward. London: IntechOpen; 2018. pp. 281–313. https://doi.org/10.5772/intechopen.76964
- 6. Grillet ME, Villegas L, Oletta JF, Tami A, Conn JE. Malaria in Venezuela requires response. Science. 2018; 359:528. pmid:29420282
- 7. Grillet ME, Hernández-Villena JV, Llewellyn MS, Paniz-Mondolfi AE, Tami A, Vincenti-Gonzalez MF, et al. Venezuela's humanitarian crisis, resurgence of vector-borne diseases, and implications for spillover in the region. Lancet Infect Dis. 2019; 19(5):149–161. pmid:30826190
- 8. Battle KE, Lucas TCD, Nguyen M, Howes RE, Nandi AK, Twohig KA, et al. Mapping the global endemicity and clinical burden of Plasmodium vivax, 2000–17: a spatial and temporal modelling study. Lancet. 2019; 394(10195):332–343. pmid:31229233
- 9. Gabaldon A. Malaria eradication in Venezuela: doctrine, practice, achievements after twenty years. Am J Trop Med Hyg. 1983; 32:203–211. pmid:6340536
- 10. Magris M, Rubio-Palis Y, Menares C, Villegas L. Vector bionomics and malaria transmission in the Upper Orinoco River, Southern Venezuela. Mem Inst Oswaldo Cruz. 2007; 102(3):303–311. pmid:17568935
- 11. Barrera R, Grillet ME, Rangel Y, Berti J, Aché A. Temporal and spatial patterns of malaria reinfection in north-eastern Venezuela. Am J Trop Med Hyg. 1999; 61:784–790. pmid:10586913
- 12. Grillet ME, El Souki M, Laguna F, León JR. The periodicity of Plasmodium vivax and Plasmodium falciparum in Venezuela. Acta Trop. 2014; 129:52–60. pmid:24149288
- 13. Moreno JE, Rubio-Palis Y, Martínez Ángela R, Porfirio A. Evolución espacial y temporal de la malaria en el municipio Sifontes del estado Bolívar, Venezuela. 1980–2013. Bol Malariol Salud Ambient. 2014; 54:236–249.
- 14. Jaramillo-Ochoa R, Sippy R, Farrel DF, Cueva-Aponte C, Beltrán-Ayala E, Gonzaga JL, et al. Effects of political instability in Venezuela on Malaria resurgence at Ecuador-Peru Border, 2018. Emerg Infect Dis. 2019; 25(4):834–836. pmid:30698522
- 15. Pacheco MA, Schneider KA, Céspedes N, Herrera S, Arévalo-Herrera M, Escalante AA. Limited differentiation among Plasmodium vivax populations from the northwest and to the south Pacific Coast of Colombia: A malaria corridor? PLoS Negl Trop Dis. 2019; 13(3):e0007310. pmid:30921317
- 16. Chenet SM, Akinyi Okoth S, Huber CS, Chandrabose J, Lucchi NW, Talundzic E, et al. Independent emergence of the Plasmodium falciparum Kelch propeller domain mutant allele C580Y in Guyana. J Infect Dis. 2016; 213:1472–1475. pmid:26690347
- 17. Pacheco MA, Forero-Peña DA, Schneider KA, Chavero M, Gamardo A, Figuera L, et al. Malaria in Venezuela: changes in the complexity of infection reflects the increment in transmission intensity. Malar J. 2020; 19:176. pmid:32380999
- 18. Bousema T, Griffin J, Sauerwein R, Smith D, Churcher T, et al. Hitting hotspots: Spatial targeting of malaria for control and elimination. PLoS Med. 2012; 9:e1001165. pmid:22303287
- 19. Moreno JE, Rubio-Palis Y, Paez E, Perez E, Sanchez V, Vaccari E. Malaria entomological inoculation rates in gold mining areas of Southern Venezuela. Mem Inst Oswaldo Cruz. 2009; 104(5):764–768. pmid:19820839
- 20.
Instituto Nacional de Estadística. INE [Internet]. Ine.gov.ve. [cited 10 January 2020]. Available from: http://www.ine.gov.ve/.
- 21. Moreno JE, Rubio-Palis Y, Paez E, Perez E, Sanchez V. Abundance, biting behaviour and parous rate of anopheline mosquito species in relation to malaria incidence in gold-mining areas of southern Venezuela. Med Vet Entomol. 2007; 21(4):339–349. pmid:18092972
- 22. Abou Orm S, Moreno JE, Carrozza M, Acevedo P, Herrera F. Plasmodium spp. infection rates for some Anopheles spp. from Sifontes Municipality, Bolívar State, Venezuela. Bol Malariol Salud Ambient. 2017; 57:17–25.
- 23. Foster PG, de Oliveira TMP, Bergo ES, Conn JE, Sant’Ana DC, Nagaki SS, et al. Phylogeny of Anophelinae using mitocondrial protein coding genes. R Soc Open Sci. 2017; 4(11):170758. pmid:29291068
- 24.
Dirección de Salud Ambiental del Estado Bolívar [DSAEB] (2007–2017). Servicio de Vigilancia Epidemiológica del Estado Bolívar, Registros Epidemiológicos, Ciudad Bolívar, Venezuela.
- 25.
Pan American Health Organization. PLISA: Plataforma de información en salud para las Américas. [online] (accessed: 15 January 2020) Available from: https://www.paho.org/data/index.php/es/temas/indicadores-malaria.html
- 26. Anselin L. Local indicators of Spatial Association-LISA. Geogr Anal. 2010; 27:93–115.
- 27.
Global Forest Watch. World Resources Institute [online] (accessed: 15 January 2020) Available from: www.globalforestwatch.org.
- 28.
United States Geological Survey (USGS) [online] (accessed: 19 December 2017) Available from: http://earthexplorer.usgs.gov/.
- 29.
Red Amazónica de Información Socioambiental Georreferenciada (RAISG). 2015. Deforestación en la Amazonía (1970–2013). 48 págs. [online] (accessed: 15 January 2020). Available from: https://www.amazoniasocioambiental.org/es/sobre/
- 30.
Dale MRT, Fortin MJ. Spatial Analysis. A Guide for Ecologists. 2nd ed. Cambridge: Cambridge University press; 2014.
- 31. Chowell G, Munayco CV, Escalante AA, McKenzie FE. The spatial and temporal patterns of falciparum and vivax malaria in Perú: 1994–2006. Malar J. 2009;8:142. pmid:19558695
- 32. Soto-Calle V, Rosas-Aguirre A, Llanos-Cuentas A, Abatih E, DeDeken R, Rodriguez H, et al. Spatio-temporal analysis of malaria incidence in the Peruvian Amazon Region between 2002 and 2013. Sci Rep. 2017; 7:40350. pmid:28091560
- 33. Grillet ME, Barrera R, Martinez JE, Berti J, Fortin MJ. Disentangling the effect of local and global spatial variation on a mosquito-borne infection in a neotropical heterogeneous environment. Am J Trop Med Hyg. 2010; 82(2):194–201. pmid:20133991
- 34. Mosha JF, Sturrock HJ, Greenwood B, Sutherland CJ, Gadalla NB, Atwal S, et al. Hot Spot or not: A comparison of spatial statistical methods to predict prospective malaria infections. Malar. J. 2014; 13:53. pmid:24517452
- 35. Noé A, Zaman SI, Rahman M, Saha AK, Aktaruzzaman MM, Maude RJ. Mapping the stability of malaria hotspots in Bangladesh from 2013 to 2016. Malar J. 2018; 17:259. pmid:29996835
- 36. Stresman G, Bousema T, Cook J. Malaria Hotspots: Is there epidemiological evidence for fine-scale spatial targeting of interventions? Trends Parasitol. 2019; 35(10):822–834. pmid:31474558
- 37. Smith DL, Perkins TA, Reiner RC Jr, Barker CM, Niu T, Chaves LF, et al. Recasting the theory of mosquito-borne pathogen transmission dynamics and control. Trans R Soc Trop Med Hyg. 2014; 108:185–197. pmid:24591453
- 38. Vittor AY, Pan W, Gilman RH, Tielsch J, Glass G, Shields T, et al. Linking deforestation to malaria in the Amazon: Characterization of the breeding habitat of the principal malaria vector, Anopheles darlingi. Am J Trop Med Hyg. 2009; 81:5–12. PMID: pmid:19556558.
- 39. Vittor AY, Gilman RH, Tielsch J, Glass G, Shields T, Lozano WS. The effect of deforestation on the human-biting rate of Anopheles darlingi, the primary vector of falciparum malaria in the Peruvian Amazon. Am J Trop Med Hyg. 2006; 74(1):3–11. PMID: pmid:16407338.
- 40. Barros FS, Honorio NA. Deforestation and malaria on the Amazon Frontier: Larval clustering of Anopheles darlingi (Diptera: Culicidae) determines focal distribution of malaria. Am J Trop Med Hyg. 2015; 93:939–953. pmid:26416110
- 41. Kweka EJ, Kimaro EE, Munga S. Effect of deforestation and land use changes on mosquito productivity and development in Western Kenya Highlands: Implication for malaria risk. Front Public Health. 2016; 4:238. pmid:27833907
- 42. Moreno JE, Rubio-Palis Y, Sánchez V, Martínez A. Fluctuación poblacional y hábitat larval de anofelinos en el municipio Sifontes, estado Bolívar, Venezuela. Bol Malariol Salud Ambient. 2015a; 55:52–68.
- 43. Moreno JE, Rubio-Palis Y, Sánchez V, Martínez Á. Caracterización de hábitats larvales de anofelinos en el municipio Sifontes del estado Bolívar, Venezuela. Bol Malariol Salud Ambient. 2015b; 55:117–131.
- 44. Hotez PJ, Basáñez MG, Acosta-Serrano A, Grillet ME. Venezuela and its rising vector-borne neglected diseases. PLoS Negl Trop Dis. 2017; 11(6):e0005423. pmid:28662038
- 45. Pániz-Mondolfi A, Tami A, Grillet ME, Márquez M, Hernández-Villena J, Escalona-Rodríguez MA, et al. Resurgence of Vaccine-Preventable Diseases in Venezuela as a Regional Public Health Threat in the Américas. Emerg Infect Dis. 2019; 25(4):625–632. pmid:30698523
- 46. Wesolowski A, Eagle N, Tatem AJ, Smith DL, Noor AM, Snow RW, et al. Quantifying the impact of human mobility on malaria. Science. 2012; 338(6104):267–270. pmid:23066082
- 47. Ihantamalala F, Herbreteau V, Rakotoarimanana FMJ, Rakotondramanga JM, Cauchemez S, Rahoilijaona B, et al. Estimating sources and sinks of malaria parasites in Madagascar. Nat Commun. 2018;9:3897. pmid:30254280
- 48. Musset L, Pelleau S, Girod R, Ardillon V, Carvalho L, Dusfour I, et al. Malaria on the Guiana Shield: A review of the situation in French Guiana. Men Inst Oswaldo Cruz. 2014; 109(5):525–533. pmid:25184998
- 49. Pommier de Santi V, Dia A, Adde A, Hyvert G, Galant J, Mazevet M, et al. Malaria in French Guiana linked to illegal gold mining. Emerg Infect Dis. 2016; 22:344–346. pmid:26811971
- 50. Douine M, Adenis A, Musset L, Demar M, Djossou F, Brousse P, et al. Malaria and illegal gold mines in a high-income country: first description of the disease burden. AM J Trop Med Hyd. 2015; 93(4):92–92.
- 51. Douine M, Sanna A, Hiwat H, Briolant S, Nacher M, Belleoud D, et al. Investigation of a possible malaria epidemic in an illegal gold mine in French Guiana: an original approach in the remote Amazonian forest. Malar J. 2019; 18, 91. pmid:30902054
- 52. Chenet SM, Schneider KA, Villegas L, Escalante, AA. Local population structure of Plasmodium: impact on malaria control and elimination. Malar J. 2012; 11:412. pmid:23232077
- 53. McCollum AM, Mueller K, Villegas L, Udhayakumar V, Escalante AA. Common origin and fixation of Plasmodium falciparum dhrf and dhps mutations associated with sulfadoxine-pyrimethamine resistance in a low-transmission area in South America. Antimicrob Agents Chemother. 2007; 51(6):2085–2091. pmid:17283199
- 54. Douine M, Musset L, Corlin F, Pelleau S, Pasquier J, Mutricy L, et al. Prevalence of Plasmodium spp. In illegal gold miners in French Guiana in 2015: a hidden but critical malaria reservoir. Malar J. 2016; 15:315. pmid:27277831
- 55. Sanchez JF, Carnero AM, Rivera E, Rosales LA, Baldeviano GC, Asencios JL, et al. Unstable malaria transmission in the southern Peruvian Amazon and its association with gold mining, Madre de Dios, 2001–2012. Am J Trop Med Hyg. 2017; 96(2):304–311. pmid:27879461
- 56. Bousema T, Stresman G, Baidjoe AY, Bradley J, Knight P, Stone W, et al. The Impact of Hotspot-Targeted Interventions on Malaria Transmission in Rachuonyo South District in the Western Kenyan Highlands: A Cluster-Randomized Controlled Trial. PLoS Med. 2016;13:e1001993.doi: https://doi.org/https://doi.org/10.1371/journal.pmed.1001993. pmid:27071072