Neglected tropical diseases risk correlates with poverty and early ecosystem destruction
Infectious Diseases of Poverty volume 12, Article number: 32 (2023)
Neglected tropical diseases affect the most vulnerable populations and cause chronic and debilitating disorders. Socioeconomic vulnerability is a well-known and important determinant of neglected tropical diseases. For example, poverty and sanitation could influence parasite transmission. Nevertheless, the quantitative impact of socioeconomic conditions on disease transmission risk remains poorly explored.
This study investigated the role of socioeconomic variables in the predictive capacity of risk models of neglected tropical zoonoses using a decade of epidemiological data (2007–2018) from Brazil. Vector-borne diseases investigated in this study included dengue, malaria, Chagas disease, leishmaniasis, and Brazilian spotted fever, while directly-transmitted zoonotic diseases included schistosomiasis, leptospirosis, and hantaviruses. Environmental and socioeconomic predictors were combined with infectious disease data to build environmental and socioenvironmental sets of ecological niche models and their performances were compared.
Socioeconomic variables were found to be as important as environmental variables in influencing the estimated likelihood of disease transmission across large spatial scales. The combination of socioeconomic and environmental variables improved overall model accuracy (or predictive power) by 10% on average (P < 0.01), reaching a maximum of 18% in the case of dengue fever. Gross domestic product was the most important socioeconomic variable (37% relative variable importance, all individual models exhibited P < 0.00), showing a decreasing relationship with disease indicating poverty as a major factor for disease transmission. Loss of natural vegetation cover between 2008 and 2018 was the most important environmental variable (42% relative variable importance, P < 0.05) among environmental models, exhibiting a decreasing relationship with disease probability, showing that these diseases are especially prevalent in areas where natural ecosystem destruction is on its initial stages and lower when ecosystem destruction is on more advanced stages.
Destruction of natural ecosystems coupled with low income explain macro-scale neglected tropical and zoonotic disease probability in Brazil. Addition of socioeconomic variables improves transmission risk forecasts on tandem with environmental variables. Our results highlight that to efficiently address neglected tropical diseases, public health strategies must target both reduction of poverty and cessation of destruction of natural forests and savannas.
Neglected tropical diseases impacted at least 1.74 billion people globally in 2019 and are associated with significant morbidity and public health burden [1, 2]. One remarkable example is dengue fever, which causes a significant economic impact on governments and households. For instance, a patient may spend 14–18 days in hospital at a cost of USD 514–1500, a severe or unbearable cost for families in poverty without adequate social security [1, 3]. To reduce the risk of occurrence and emergence of infectious diseases, several approaches have been developed . Recent strategies include interdisciplinary and multisectoral collaboration between public health and governmental institutions. For instance, the One Health approach proposes the monitoring of disease transmission between people, animal vectors, and reservoirs, as well as the ecosystem, and encompasses the agriculture, health, and environment sectors of society . The Global Health Security Agenda  seeks to strengthen countries’ capacity for the prevention and surveillance of infectious diseases. These strategies are consistent with the third Sustainable Development Goal of the United Nations which is to ensure healthy lives and well-being (https://www.undp.org/sustainable-development-goals). Global health security is a goal that is conditional upon adequate infectious disease monitoring and spatial risk assessment [6, 7]. Hence, prevention and monitoring are necessary strategies to ensure global health security, and are a common focus for sustainable development achievement efforts [5, 6, 8].
An effective approach of disease surveillance is to investigate the underlying determinants of disease emergence and recurrence in Latin America . The natural history of a disease system is a key aspect to consider when aiming to reconstruct and predict disease outcomes. Heterogeneities and particularities in vector-borne disease cycles are influenced by the environment . For instance, the availability of mosquito breeding sites is a significant factor in the spread of dengue fever and malaria. This includes the former requiring artificial water containers or ponds, as a premise for the spread of Aedes mosquitoes, and the latter requiring forest related water bodies in shadowed areas, Anopheles mosquitoes most common breeding sites [10, 11]. Dengue virus encompasses serotypes 1 to 4, while malaria is caused by protozoans of the Plasmodium genus, mainly by P. falciparum and P. vivax . Furthermore, Chagas disease can be spread by contact with hematophagous triatomine (kissing bug) feces infected by Trypanosoma protozoans [12, 13]. Usually, contamination occurs trough kissing bug bites, or through contaminated food [12, 13]. The presence of the triatomine bugs is associated with poor household construction . Cutaneous, and visceral leishmaniasis are also transmitted by the bites of arthropod vectors (female phlebotomine sandflies) and are caused by protozoans of genus Leishmania [1, 13]. Leishmaniasis transmission is associated with forest fragmentation . Moreover, Brazilian spotted fever is caused by bacteria of the genus Rickettsia, and is transmitted by tick bites. Transmission usually encompasses an animal reservoir such horses, cattle or capibaras  (See Fig. 1). Despite significant progress in understanding the natural history of vector-borne and neglected tropical diseases, the relationship between disease risk and socioeconomic and environmental factors is still poorly understood.
In the case of directly or environmentally transmitted diseases such as schistosomiasis and leptospirosis which are mainly transmitted by contaminated bodies of water, heterogeneities in the environment become direct drivers of transmission [16,17,18]. Intestinal schistosomiasis is caused by the helminth Schistosoma mansoni and its transmission cycle encompasses Biomphalaria snails as reservoir hosts . Leptospirosis is caused by bacteria of Leptospira genus and its transmission cycle develops in the presence of urine of infected domestic rodents . Likewise, hantavirus is transmitted by inhalation of particles contaminated with the feces of infected wild rodents . Usually, it more heavily affects human in rural environments who work with crops, specifically sugar cane . Hantavirus manifests as hantavirus pulmonary syndrome and can be fatal for infected individuals.
As the environment plays a key role in disease transmission and maintenance, models linking environmental covariates with disease occurrence are often used in modern disease transmission risk mapping [21,22,23,24,25,26]. Neglected tropical diseases, for example, have strong ties to environmental parameters such as temperature and land cover, which correspond to the biological responses of vector survival, abundance, and pathogen transmission . Precipitation and temperature have been frequently used to forecast the spatiotemporal occurrence of vector-borne diseases and zoonotic diseases like hantavirus [21, 27, 28]. The predictive ability of the climate on disease spatial distribution is linked to the importance of optimal temperatures and water availability for parasite reproduction and survival . Environmental changes, measured as landscape change or ecosystem loss, have been used as predictor variables of multiple aspects of disease in spatiotemporal analysis . When ecosystem loss occur (hereafter referred to as ‘ecosystem destruction’) it may promote disease outbreaks by facilitating the distributional overlap of pathogens, vectors, and reservoirs with domestic animals and humans . In addition, ecosystem destruction can change the abundance of reservoirs and vectors , thus augmenting the stability of different host–pathogen systems. In Brazil, current evidence shows that deforestation promotes outbreaks of malaria  and others diseases with high morbidity and mortality .
Socioeconomic factors also influence how human diseases emerge and persist, since they affect how individuals live, their quality of life, and their interaction with the environment [32, 33]. Neglected tropical diseases, such as Chagas disease, schistosomiasis, dengue, and leishmaniasis are linked to poverty. Neglected tropical diseases disproportionately affect lower income households, resulting in low productivity, a high epidemiological burden (i.e., mortality and morbidity), disability, and social stigma [32, 34]. Poverty, measured as low proportional income, is a key predictor of disease presence and prevalence at individual or spatial level, therefore, it may help explain the spatial distribution of infectious disease . Urbanization facilitates disease transmission, as demonstrated by dengue fever where mosquitoes thrive in human modified environments that offer abundant resources, such as breeding sites . Similarly, a lack of access to clean water and sanitation has been shown to increase the risk of disease transmission for diseases such as schistosomiasis . Interaction between infectious disease and poverty creates a difficult-to-break cycle, or poverty traps, which can last for generations if proper health interventions are not applied, creating an increasing call to properly surveil infectious diseases [6, 32]. Nevertheless, this challenge is exacerbated in low-income regions where these data are most needed, as data are typically incomplete, sparse, or rarely available . As a result, until recently, the use of the socioeconomic components in disease modeling was mostly overlooked, and was not often included in the ecological niche modeling framework of disease risk mapping .
The ecological niche modeling framework can be used to delimit the realized ecological niche of parasite species. Realized ecological niches provide proxies of the influence of physiological tolerance, biological interactions, and dispersal on the ecology of organisms (e.g., pathogens, vectors) and are measured as a hypervolume along environmental variables occupied by the species . Here, to elucidate the effects of social variables in forecasting the geographic distribution of disease transmission, we combined socioeconomic and environmental variables to address these gaps and improve the predictability of neglected tropical diseases throughout all Brazilian regions. In addition, we sought to investigate the main socioeconomic and environmental drivers affecting the disease transmission risk. We built ecological niche models using different sets of variables: socioeconomic (e.g., mean income, mean income inequality Gini index, housing quality) and environmental (e.g., temperature, precipitation, ecosystem destruction) variables in combined (composite) and environmental only (simplified) models, using disease occurrence as the response variable. Furthermore, we assessed model predictability (capacity to predict disease occurrence) using different model configurations to tease apart which variables best predict disease occurrence. We modeled five vector-borne (Chagas disease, dengue fever, malaria, leishmaniasis, and Brazilian spotted fever) and three directly (environmentally) transmitted zoonotic (hantavirus, schistosomiasis, and leptospirosis) neglected tropical diseases to determine if model predictability and relevant predictor variable were consistent across infectious diseases of different natural history. We hypothesized that incorporating socioeconomic variables into such predictive models would add more precision in limiting disease risk in areas with high socioeconomic inequality and widespread poverty. Moreover, integrating relevant socioeconomic and environmental variables could be a significant step forward aligning the Global Health Security agenda within the United Nation’s Sustainable Development goals (https://www.undp.org/sustainable-development-goals).
Disease occurrence data
The disease occurrence data (i.e., locations where disease was reported) came from the Brazilian Ministry of Health database (DATASUS ). This database contains all the reported cases of diseases that require mandatory notification, such as those of epidemic potential. We collected data of infectious disease from 2007 to 2019, which is the span of time with lab-confirmed occurrences available. Only new occurrences, (i.e., incidence records) were regarded as a result. The chosen diseases were mediated by at least one species of vectors, hosts, or reservoirs, hence being a subject to environmental influence. To reduce uncertainty only laboratory confirmed cases were considered in the study. To avoid interference from vaccination that could potentially bias the models, we excluded diseases that were preventable with high vaccination coverage, such as yellow fever. Eight diseases met these requirements, namely the five vector-borne diseases dengue fever, malaria, Chagas disease, cutaneous and visceral leishmaniasis, and Brazilian spotted fever, as well as three zoonotic diseases including schistosomiasis, leptospirosis, and hantavirus pulmonary syndrome. It is important to keep in mind that the DATASUS disease-occurrence data has several limitations. For instance, the unequal distribution of healthcare coverage in Brazil, especially the diagnostic capacity, could potentially bias the distributions of disease occurrence. Nevertheless, the Brazilian healthcare system is one of the largest integrated healthcare systems in the world and has a reasonable reach capacity, being able to attend isolated populations throughout its territory . While this does not solve biased reporting, it does mitigate potential spatial biases in disease occurrence. Additional file 1: Table S1 expands and explores these and other potential limitations of our method and provides more details on the disease occurrence data (Additional file 1).
The disease occurrences were assigned to the municipality where the infection occurred, as declared in the DATASUS database. Only locality of disease presence was regarded, not the number or prevalence of cases. Exact locations of disease were not available, so they were assigned to their corresponding municipality administrative center (available from IBGE, The Brazilian Institute of Geography and Statistics ). Municipalities were unequal in size, so we standardized occurrences inside an 18 km grid, which resembles the cell size of predictor variables (18 km). To mitigates spatial autocorrelation, we used only one presence point per cell of the spatial extent of the study area. Potential consequences of this approach are discussed in Additional file 1.
We employed temperature-derived and precipitation-derived variables from WorldClim  with a spatial resolution of 10 arc-minutes (approx.18 km) averaged for the period 1970–2000. Worldclim bioclimatic variables are biologically associated to distributional limits of plant and animal species and are often used in ecological niche modeling to estimate environmental limits . For additional landscape-change variables, we used natural habitat cover from the year 2016 (the sum of forest and savannah formation from each municipality in 2016). We used the values from 2016 since it was the median year between occurrence cases. Furthermore, we used ecosystem destruction between 2009 and 2018 (the sum of natural forest and grassland cover that were transformed into urban, pasture or plantation areas) from MapBiomas project to the municipality level (mapbiomas.org) .
As socioeconomic variables, we explored proxies for sanitation, demography, and income. As proxies of sanitation, we used percentages of households with piped water, percentages of households with toilets derived from the 2010 Brazilian socioeconomic census . We also used estimates of population density , the mean human development index for each municipality , and the Gini index from IBGE , which is a coefficient of income disparity. All of these socioeconomic variables are associated with disease vulnerability . We also used the per capita gross domestic product (GDP) as a proxy of socioeconomic status and economic development [43, 44]. All socioeconomic and environmental variables were resampled to approx. 18 km (10’) to match the spatial resolution of occurrence and climatic variables, including variables addressed to municipalities’ bounds. GDP and the Worldclim variables values were extracted at the municipality administrative center coordinate. More details on environmental and socioeconomic variables such as original resolution and source are described at Additional file 1: Table S2.
Variable selection and correlation
To test whether socioeconomic variables improve model performance for each disease, we made two sets of models: one set of models with the environment and socioeconomic variables (combined model) and one set of models with environment-only variables (simple model). In building the models, we started by checking the correlation between all the predictor variables. We allowed the fewest number of variables in each type of model as possible in accordance to the Occam’s razor principle, which prioritizes the simplest explanation. We removed variables with > 70% correlation, aiming to retain the ones with biological linkages with overall disease transmission (Additional file 2). The final set of variables for the combined models included mean annual temperature and annual precipitation, which are linked to vector or reservoir tolerance and suitability [22, 24, 45]; environmental destruction between 2009 and 2018 which was highly correlated with 2016 native forest/vegetation cover (74%), and often associated with pathogen spillover [31, 46, 47]; GDP (highly correlated with human population, 78%), Gini coefficient (proxy for inequality), and proportion of households without a toilet (highly correlated with human development index, − 83%) which are often regarded as proxies of population vulnerability [13, 48]. We used a similar set of variables when building the environmental-only set of variables for the composite model, keeping only the environment-related (temperature, precipitation, loss of natural habitat). Additionally, we used the ‘vifstep’ function from the sdm R package  to choose a less correlated environmental variable to add to the model, so that both models have a similar number of variables and can be compared (see Additional file 2).
Ecological niche modeling
Models for each disease were built using four correlative methods (i.e., generalized linear models, Maxent, random forest, and support vector machines) and one climatic envelop method (Bioclim). These algorithms are more conservative algorithms since they prioritize interpolation over extrapolation [50, 51]. Reduced extrapolation is desirable in disease risk mapping to mitigate overprediction in environmental conditions beyond the observed values . We used cross-validation for model calibration using the h-block strategy that generates nine independent models for each algorithm, for a total of 45 models per disease . Then, using the ‘bin model()’ function from the ntbox R package, we transformed each model using a 10% presence threshold to binarize a continuous probability map into presence and absence maps . After binarization, we generated ensemble models by summarizing the algorithm models together for each disease. Each model ensemble represented consensus among algorithms on the presence of conditions suitable for long-term disease transmission as a proxy for disease transmission risk .
Model evaluation and test
After building two ensemble models for each disease, we evaluated their quality using a partial receiver operating characteristic metric (pROC) that overcomes the limitation of classical ROC and area under the curve (AUC) approaches . ROC and AUC are classic model evaluating metrics that partially ignore the goodness of fit of the models and are usually biased to favor some algorithms over others. The AUC/ROC metric assume algorithms that span broader predicted areas regardless of commission error as accurate. In other words, the extent of the background area affects its outcomes [55, 56]. The pROC values were developed using the ntbox R package  and exhibit a range from zero to two, where values above one represent predictions better than random expectations . Finally, to test whether socioeconomic-environmental models performed overall better than environmental-only models in all diseases, we used a simple t-test for paired samples using pROC values.
Relative variable importance and response curves
To explore the relationship between environmental and socioeconomic variables and disease-case occurrence, we used response curves generated by the evaluation strip method . The evaluation strip method addresses the visualization of predicted responses of a species (in this case, a disease, or pathogen) to environmental variables. Response curves and relative variable importance analysis were performed using the sdm package . The relative variable importance identifies the most important variables in the model, and response curves inform variable effects in relation to disease occurrences. Both response curves and variable importance were generated based on the ensemble model (averaged by all algorithms) for each variable set (environmental only and socioenvironmental variable set).
We modeled the spatial distribution of disease-transmission risk for nine neglected tropical diseases in Brazil from 2007 to 2019. In total we assessed 723,109 disease occurrence records. The number of confirmed cases ranged from 219 for Chagas disease to 429,052 cases for dengue fever. Diseases such as leptospirosis, dengue fever, malaria, and leishmaniasis were widespread across the Brazilian territory. Overall, the composite models including both socioeconomic and environmental predictors performed 10% better (based on partial-ROC estimates, α = 0.05; Table 1) than the simple model using only environmental predictors.
After incorporating socioeconomic variables into the modeling effort, disease risk areas became more defined (Figs. 2 and 3), particularly for dengue fever, malaria, and cutaneous leishmaniasis, which shifted from mild suitability values covering the entire country to more delimited risk in southeast, central, and northeast Brazil, as well as along the Amazon River basin(Figs. 2 and 3). Model performance for dengue fever, malaria, cutaneous leishmaniasis, and leptospirosis, improved the most after adding socioeconomic variables (Table 1, Fig. 4A). The addition of socioeconomic variables had no significant effect on the models of Brazilian spotted fever, visceral leishmaniasis, and Chagas disease (Table 1, Fig. 4A), and actually slightly decreased the performance of the spotted fever model (change in pROC = −1.1%).
Poverty and disease risk
In the composite models, The most significant predictive variables for all diseases were GDP (relative variable importance = 37 ± 13% standard error), followed by annual amount of precipitation (29 ± 7%) and annual mean temperature (24 ± 8%). Low GDP (at or below the poverty level) was associated with higher probability of disease presence, with the strongest effect being found for hantavirus and Chagas disease (Fig. 5). A low GDP was linked to a higher likelihood of sickness, with the exception of dengue fever, which had a positive correlation with GDP (Figs. 4B, 5, Additional file 1: Fig. S3). GDP was the most important socioeconomic variable in the combined models for seven of the nine studied illnesses: schistosomiasis, leptospirosis, malaria, dengue fever, visceral and cutaneous leishmaniasis, and Chagas disease. For hantavirus and Brazilian spotted fever, theproportion of households without toilets was the most significant disease risk predictor (42% and 27% of relative variable importance respectively; Fig. 4B).
Ecosystem destruction and disease risk
In the simple models, the most important environmental predictor for most diseases (except for Brazilian spotted fever) was ecosystem destruction between the years of 2008 and 2019 (mean relative variable importance of 42%), followed by annual mean temperature (32%, Fig. 4). Ecosystem destruction was the most important variable of environmental models for six of the nine studied diseases: leptospirosis, cutaneous leishmaniasis, dengue fever, schistosomiasis, hantavirus, and Chagas disease. Overall, ecosystem destruction had 31% more importance in determining disease presence than the next variable which was mean annual temperature. The relationship between ecosystem destruction and disease likelihood exhibited a negative trend, which meant ecosystem destruction increased disease likelihood when its levels were low (Fig. 5).
Understanding the environmental and socioeconomic predictors which determine disease probability and vulnerability in human populations is a core question in disease surveillance and prevention. Untangling the effects of these variables can help predict transmission events and spatiotemporal emergence of outbreaks, which allow us to forecast the origins and paths of epidemic spread. Here, we show how ecosystem destruction and poverty correlates with the presence of several life cycle zoonoses and neglected tropical diseases on large geographic areas (Fig. 3). Disease probability was relatively high when GDP and natural ecosystem destruction were low, corresponding to its early stages. Furthermore, adding socioeconomic variables in disease risk modeling improved overall model performance and accuracy by up to 18%. Likewise, socioeconomic variables significantly contribute to habitat suitability modeling and disease risk mapping for several neglected tropical and zoonotic diseases.
Disease modeling approaches supported by ecological niche theory have additional advantages and applications concerning traditional spatial epidemiology models, including a biological understanding of the drivers of disease occurrence . Climate-only ecological niche models, however, may end up identifying areas with the environmental conditions suitable for the pathogen, vector, or reservoir, without accounting for the presence of vulnerable human populations or densely populated cities. We found that dengue fever, malaria, Chagas disease, and hantavirus models had a more delimited distribution when susceptible populations were considered as a force facilitating the long-term maintenance of transmission (Fig. 2). Our findings also demonstrate that the absence of socioeconomic variables in disease transmission models may result in a misleading disease transmission risk, which in turn can have a negative influence on public health policy aiming to allocate resources for disease prevention and treatment . More precision is a desirable modeling feature when aiming to conduct data-driven surveillance or deployment of vaccines or treatment to at-risk populations outside the traditional disease-control plans.
Almost all diseases had GDP as a crucial socioeconomic variable, with a tendency to have a higher probability of disease at low-income levels, (i.e., high poverty). This pattern was demonstrated more strongly for Chagas disease and hantavirus (Fig. 5). Overall, these results suggests that poorer populations may have higher probability of disease transmission risk. In fact, poverty is usually associated with infectious disease risk and susceptibility [32, 33]. Low GDP has been linked to malaria, dengue, and a variety of other infectious diseases, where it may act as both a cause and a consequence of disease risk [31, 42]. This relationship between poverty and neglected tropical diseases may hinder populations and communities from escaping disease risk, a process known as poverty trap . Poverty traps may lock communities into a cycle of diseases and poverty that may last generations if interventions with specific policies are not implemented . As a result, public health policies must be combined with social policies that reduce poverty.
In contrast, GDP was positively associated with dengue fever risk (Figs. 2 and 5), where municipalities with higher income also had higher risk of infection. This could be associated with the capacity of dengue transmission to be better sustained is urban areas with high population densities . An analogous result was obtained with Brazilian spotted fever, which had a negative relationship with the proportion of households without a toilet (Additional file 1, Fig. S6). One explanation for this pattern could be that cities with high proportion of houses without toilet are in the Amazon region where Brazilian spotted fever is less prevalent . Brazilian spotted fever occurs in areas with more intense agriculture and cattle economy, towards the central part of the country. These unexpected results encompassing Brazilian spotted and dengue fever demonstrates how different socioeconomic status proxies can denote complex transmission-risk gradients.
We also observed that disease transmission risk for dengue, cutaneous leishmaniasis, Chagas disease, schistosomiasis, leptospirosis, and hantavirus increases under moderate to low levels of ecosystem degradation, in agreement with similar findings on malaria (Fig. 5) . The early stages of ecosystem degradation could alter parasite transmission cycles surrounding pristine forests and savannas facilitating disease emergence . In the case of malaria, the mosquito vector Anopheles cruzii is found in great abundance in tree canopies, but during deforestation the mosquito shifts to ground level increasing malaria prevalence in both simians and humans . This driver of disease emergence is also related to spillover events, which occur when a parasite from a natural cycle of transmission infects a different species, in this case humans , which has been observed during deforestation events . The link between land-use change and outbreaks has been demonstrated recently by epidemics such as Ebola  and COVID-19 . The gradient of recent deforestation or native vegetation loss in Brazil suggests that active deforestation was widespread in municipalities across the country  (see correlation of variables at Additional file 2). Thus, municipalities that went through ecosystem degradation in the past, but do not show ongoing landscape change, may demonstrate different transmission patterns and lower disease risk. Active deforestation and landscape fragmentation may lead to peaks in disease transmission. Our study shows insights on how disease transmission risk varies along ecosystem destruction gradients and how the age of landscape destruction can be used to model vulnerability to disease across large-scale study areas like Brazil.
Our findings also support the argument that disease-transmission risk is driven by the confluence of poverty and ecosystem destruction. Thus, human population health security cannot be achieved without addressing associated sectors such as ecosystem health. Our results indicate that it may be necessary to halt deforestation, even in its early stages to prevent spillover events and new outbreaks. This concept is specifically highlighted in the One Health approach  and it is coherent with the United Nations Sustainable Development Goals and the Global Health Security Agenda . Additionally, our results highlight the importance of multiple sustainable development goals and provide an example of how they are interrelated: to promote, prevent, and surveil good health and well-being it is necessary to extinguish poverty in its different forms and to ensure the protection of terrestrial ecosystems. These are central themes attuned by United Nations by the first, third, and fifteenth sustainable development goals. Unfortunately, actions towards poverty eradication, health and well-being promotion, and the protection of terrestrial land had been reduced by the SARS-CoV-2 ongoing pandemic, which made these themes even more urgent. Thus, to prevent neglected tropical diseases that have strong relationship to poverty and environment destruction, is necessary encompassing common solutions among the divergent development goals.
It is worth noting that our analyses could have potential limitations. As aforementioned, potential diagnostic bias, particularly in terms of regional variations in the capacity to detect diseases, can produce low specificity and sensibility in the health surveillance system, especially in the Amazon region . However, the good reach capacity of Brazilian healthcare  and our presence-only approach, which is more conservative, can mitigate this problem. Furthermore, we used 2010 socioeconomic parameters, and it is possible that supplying outdated variables to the model will affect its accuracy and applicability or bias the models. However, we observed that socioeconomic variables in Brazil, at a large scale, are highly correlated and retain their spatial patterns regardless of timeframe. While the current Brazilian government has been slow to release socioeconomic census data , we believe that it was the best decision to use the variables from the 2010 socioeconomic census to not further hamper the importance of this disease surveillance effort. Thus, we believe our analysis and interpretation of the results are sound. These and other limitations are discussed in greater detail in Additional file 1.
Public health strategies in Brazil should be directed to economically disadvantageous populations adjacent to areas undergoing fast ecosystem destruction. We found that the high socioeconomic inequalities and disparities widespread in Brazil are expected to exacerbate the burden of neglected tropical diseases. This is particularly remarkable given the recent country’s increasing scenario of severe and rapid deforestation, which is influenced by ongoing deliberate decisions from the Brazilian government . For example, during the administration of president Jair Bolsonaro (2019–2022), Brazil’s environmental agency spent less than half of its budget to protect the environment in 2021 despite the record-breaking deforestation . Also in 2021, Brazil was in a severe economic crisis, making populations in poverty even more vulnerable to disease outbreaks. Strikingly, scientists in Brazil now have limited access to updated data on both the socioeconomic-demographic census and infectious disease occurrences at the national level implying that critical decisions are being conducted with incomplete information. To achieve satisfactory public health results and sustainable development, Brazil must recover disease surveillance and sociodemographic census on a country scale and implement policies aiming to reduce socioeconomic imbalances and destruction of natural habitats.
Neglected tropical and zoonotic diseases are main explained, at macro-scales by poverty and early ecosystem destruction in Brazil. To disrupt the cycle of disease transmission, it is vital to ensure that strategies for public health are aligned with policies of poverty diminishment and ecosystem conservation policies. Seeking equitable public health interventions and socioeconomic rescue well aide in disrupting poverty cycles, balanced with coordinated country-specific strategies on deforestation suppression, should be a continental priority to control the neglected tropical diseases and prevent poverty traps caused by a combination of environmental and social factors that facilitate disease emergence and transmission.
Availability of data and materials
The data and code to reproduce the analyses are available in the Additional files1, 2 and in in the GitHub repository: https://github.com/arthurama/poverty-and-habitat-loss-are-predictors-of-NTD-in-Brazil.git.
Informatics Department of the Unified Brazilian Health System
Brazilian institute of geography and statistics
Gross domestic product
Partial receiver operating characteristic
Area under the curve
World Health Organization. A global brief on vector-borne diseases. World Heal Organ. 2014.
The Global Health Observatory. Reported number of people requiring interventions against NTDs. World Heal. Organ. 2021. https://www.who.int/data/gho/data/indicators/indicator-details/GHO/reported-number-of-people-requiring-interventions-against-ntds. Accessed 1 Mar 2021.
Teich V, Arinelli R, Fahham L. Aedes aegypti e sociedade: o impacto econômico das arboviroses no Brasil. J Bras Econ da Saúde. 2017;9:267–76.
Pike J, Bogich T, Elwood S, Finnoff DC, Daszak P. Economic optimization of a global strategy to address the pandemic threat. Proc Natl Acad Sci USA. 2014;111:18519–23.
Berthe FCJ, Bouley T, Karesh WB, Le Gall FG, Machalaba CC, Planté CA, et al. One Health—Operational framework for strengthening human, animal and environmental public health systems at their interface. Washington, D.C.; 2018.
Global Health Security Agenda Steering Group. Global Health Security Agenda (GHSA) 2024 framework. 2018. p. 35. https://ghsagenda.org/wp-content/uploads/2020/06/ghsa2024-framework.pdf. Accessed 25 Oct 2021.
Worsley-Tonks KEL, Bender JB, Deem SL, Ferguson AW, Fèvre EM, Martins DJ, et al. Strengthening global health security by improving disease surveillance in remote rural areas of low-income and middle-income countries. Lancet Glob Heal. 2022;10:e579–84.
Di Marco M, Baker ML, Daszak P, De Barro P, Eskew EA, Godde CM, et al. Opinion: sustainable development must account for pandemic risk. Proc Natl Acad Sci USA. 2020;117:3888–92. https://doi.org/10.1073/pnas.2001655117.
Alexander KA, Lewis BL, Marathe M, Eubank S, Blackburn JK. Modeling of wildlife-associated zoonoses: applications and caveats. Vector-Borne Zoonotic Dis. 2012;12:1005–18.
Carlos BC, Rona LDP, Christophides GK, Souza-Neto JA. A comprehensive analysis of malaria transmission in Brazil. Pathog Glob Health. 2019;113:1–13.
Chersich MF, Wright CY, Venter F, Rees H, Scorgie F, Erasmus B. Impacts of climate change on health and Wellbeing in South Africa. Int J Environ Res Public Health. 2018;15:1884.
Coura JR. The main sceneries of Chagas disease transmission. The vectors, blood and oral transmissions—a comprehensive review. Mem Inst Oswaldo Cruz. 2015;110:277–82.
Codeço CT, Dal’Asta AP, Rorato AC, Lana RM, Neves TC, Andreazzi CS, et al. Epidemiology, biodiversity, and technological trajectories in the Brazilian Amazon: from malaria to COVID-19. Front Public Heal. 2021;9.
de Curi NHA, de Paschoal AMO, Massara RL, Marcelino AP, Ribeiro AA, Passamani M, et al. Factors associated with the seroprevalence of leishmaniasis in dogs living around Atlantic forest fragments. PLoS ONE. 2014;9:e104003. https://doi.org/10.1371/journal.pone.0104003.
de Oliveira SV, Guimarães JN, Reckziegel GC, Neves BM da C, Araújo-Vilges KM de, Fonseca LX, et al. An update on the epidemiological situation of spotted fever in Brazil. J Venom Anim Toxins Incl Trop Dis. 2016;22.
Noya O, Katz N, Pointier JP, Theron A, de Noya BA. Schistosomiasis in America. Neglected Tropical Dis—Lat Am Caribb. Vienna: Springer; 2015. p. 11–43.
Grimes JET, Croll D, Harrison WE, Utzinger J, Freeman MC, Templeton MR. The relationship between water, sanitation and schistosomiasis: a systematic review and meta-analysis. PLoS Negl Trop Dis. 2014;8:e3296. https://doi.org/10.1371/journal.pntd.0003296.
Mwachui MA, Crump L, Hartskeerl R, Zinsstag J, Hattendorf J. Environmental and behavioural determinants of leptospirosis transmission: a systematic review. PLoS Negl Trop Dis. 2015;9:e0003843. https://doi.org/10.1371/journal.pntd.0003843.
Haake DA. Molecular epidemiology of leptospirosis in the Amazon. PLoS Med. 2006;3:1214–5.
Chandy S, Mathai D. Globally emerging hantaviruses: an overview. Indian J Med Microbiol. 2017;35:165–75.
Prist PR, Uriarte M, Tambosi LR, Prado A, Pardini R, D’Andrea PS, et al. Landscape, environmental and social predictors of Hantavirus Risk in São Paulo, Brazil. PLoS ONE. 2016;11:e0163459. https://doi.org/10.1371/journal.pone.0163459.
Escobar LE, Craft ME. Advances and limitations of disease biogeography using ecological niche modeling. Front Microbiol. 2016;7:1–21.
Chipperfield JD, Benito BM, O’Hara R, Telford RJ, Carlson CJ. On the inadequacy of species distribution models for modelling the spread of SARS-CoV-2: response to Araújo and Naimi. EcoEvoRxiv; 2020. https://osf.io/mr6pn.
Costa J, Peterson AT. Ecological niche modeling as a tool for understanding distributions and interactions of vectors, hosts, and etiologic agents of Chagas disease. Adv Exp Med Biol. 2012;710:59–70.
Conquest L. Analysing ecological data. J Am Stat Assoc. 2008. https://doi.org/10.1198/016214508000000715.
Fuller TL, Gilbert M, Martin V, Cappelle J, Hosseini P, Njabo KY, et al. Predicting hotspots for influenza virus reassortment. Emerg Infect Dis. 2013;19:581–8.
Pascual M, Ahumada JA, Chaves LF, Rodó X, Bouma M. Malaria resurgence in the East African highlands: temperature trends revisited. Proc Natl Acad Sci USA. 2006;103:5829–34.
Tian H, Yu P, Cazelles B, Xu L, Tan H, Yang J, et al. Interannual cycles of Hantaan virus outbreaks at the human-animal interface in Central China are controlled by temperature and rainfall. Proc Natl Acad Sci U S A. 2017;114:8041–6. https://doi.org/10.1073/pnas.1701777114.
Myers SS, Gaffikin L, Golden CD, Ostfeld RS, Redford KH, Ricketts TH, et al. Human health impacts of ecosystem alteration. Proc Natl Acad Sci USA. 2013;110:18753–60.
Faust C, McCallum HI, Bloomfield L, Nicole G, Gillespie T, Dobson A, et al. Pathogen spillover during land conversion. Ecol Lett. 2018;21:471–83. https://doi.org/10.1111/ele.12904.
Macdonald AJ, Mordecai EA. Amazon deforestation drives malaria transmission, and malaria burden reduces forest clearing. Proc Natl Acad Sci USA. 2019;117:20335.
Bhutta ZA, Sommerfeld J, Lassi ZS, Salam RA, Das JK. Global burden, distribution, and interventions for infectious diseases of poverty. Infect Dis Poverty. 2014;3:1–7.
Tatem AJ, Adamo S, Bharti N, Burgert CR, Castro M, Dorelien A, et al. Mapping populations at risk: improving spatial demographic data for infectious disease modeling and metric derivation. Popul Health Metr. 2012;10:1–14.
Aagaard-Hansen J, Chaignat CL. Neglected tropical diseases: equity and social determinants. In: Blas E, Sivasankara Kurup A, editors. Equity, social determinants and public health programmes Geneva: World Health Organization; 2010. p. 135–57. https://www.who.int/neglected_diseases/Social_determinants_NTD.pdf.
Soberon J, Peterson AT. Interpretation of models of fundamental ecological niches and species’ distributional areas. Biodivers Informatics. 2005;2.
Brazillian Ministry of Health. DATASUS. 2019. http://www2.datasus.gov.br/.
Castro MC, Massuda A, Almeida G, Menezes-Filho NA, Andrade MV, de Souza Noronha KVM, et al. Brazil’s unified health system: the first 30 years and prospects for the future. Lancet. 2019;394:345–56.
IBGE. Censo Demográfico. 2010. p. 12. : http://www.sidra.ibge.gov.br/bda/.
Fick SE, Hijmans RJ. WorldClim 2: new 1-km spatial resolution climate surfaces for global land areas. Int J Climatol. 2017;37:4302–15.
Souza CM, Shimbo JZ, Rosa MR, Parente LL, Alencar AA, Rudorff BFT, et al. Reconstructing three decades of land use and land cover changes in brazilian biomes with landsat archive and earth engine. Remote Sens. 2020;12:2735.
Tatem AJ. WorldPop, open data for spatial demography. Sci Data. 2017;4:2–5.
Sun TT, Tao R, Su CW, Umar M. How do economic fluctuations affect the mortality of infectious diseases? Front Public Heal. 2021;9:1–13.
Kummu M, Taka M, Guillaume JHA. Gridded global datasets for Gross Domestic Product and Human Development Index over 1990–2015. Sci Data. 2018;5:1–15.
Astrom C, Rocklöv J, Hales S, Béguin A, Louis V, Sauerborn R. Potential distribution of dengue fever under scenarios of climate change and economic development. EcoHealth. 2012;9:448–54.
Li R, Xu L, Bjørnstad ON, Liu K, Song T, Chen A, et al. Climate-driven variation in mosquito density predicts the spatiotemporal dynamics of dengue. Proc Natl Acad Sci. 2019;116:3624–9. https://doi.org/10.1073/pnas.1806094116.
Dos SCVB, Sevá ADP, Werneck GL, Struchiner CJ. Does deforestation drive visceral leishmaniasis transmission? A causal analysis. Proc R Soc B. 2021;288:1–9.
Olivero J, Fa JE, Real R, Márquez AL, Farfán MA, Vargas JM, et al. Recent loss of closed forests is associated with Ebola virus disease outbreaks. Sci Rep. 2017;7:1–9.
Fernández MDP, Gaspe MS, Gürtler RE. Inequalities in the social determinants of health and Chagas disease transmission risk in indigenous and creole households in the Argentine Chaco. Parasites Vectors. 2019;12:1–18. https://doi.org/10.1186/s13071-019-3444-5.
Naimi B, Araújo MB. Sdm: a reproducible and extensible R platform for species distribution modelling. Ecography (Cop). 2016;39:368–75.
Qiao H, Soberón J, Peterson AT. No silver bullets in correlative ecological niche modelling: insights from testing among many potential algorithms for niche estimation. Methods Ecol Evol. 2015;6:1126–36.
Escobar LE, Qiao H, Cabello J, Peterson AT. Ecological niche modeling re-examined: a case study with the Darwin’s fox. Ecol Evol. 2018;8:4757–70.
Muscarella R, Galante PJ, Soley-Guardia M, Boria RA, Kass JM, Uriarte M, et al. ENMeval: an R package for conducting spatially independent evaluations and estimating optimal model complexity for Maxent ecological niche models. Methods Ecol Evol. 2014;5:1198–205.
Osorio-Olvera L, Lira-Noriega A, Soberón J, Peterson AT, Falconi M, Contreras-Díaz RG, et al. ntbox: an r package with graphical user interface for modelling and evaluating multidimensional ecological niches. Methods Ecol Evol. 2020;11:1199–206.
Peterson AT. Mapping disease transmission risk: enriching models using biogeography and ecology. Baltimore: Johns Hopkins University Press; 2014. p. 2014.
Peterson AT, Papeş M, Soberón J. Rethinking receiver operating characteristic analysis applications in ecological niche modeling. Ecol Modell. 2008;213:63–72.
Lobo JM, Jiménez-valverde A, Real R. AUC: a misleading measure of the performance of predictive distribution models. Glob Ecol Biogeogr. 2008;17:145–51.
Elith J, Ferrier S, Huettmann F, Leathwick J. The evaluation strip: a new and robust method for plotting predicted responses from species distribution models. Ecol Modell. 2005;186:280–9.
Bonds MH, Keenan DC, Rohani P, Sachs JD. Poverty trap formed by the ecology of infectious diseases. Proc R Soc B. 2010;277:1185–92.
Multini LC, de Souza AL da S, Marrelli MT, Wilke ABB. The influence of anthropogenic habitat fragmentation on the genetic structure and diversity of the malaria vector Anopheles cruzii (Diptera: Culicidae). Sci Rep. 2020;10.
Johnson EE, Escobar LE, Zambrana-Torrelio C. An ecological framework for modeling the geography of disease transmission. Trends Ecol Evol. 2019;34:655–68. https://doi.org/10.1016/j.tree.2019.03.004.
Rulli MC, Santini M, Hayman DTS, D’Odorico P. The nexus between forest fragmentation in Africa and Ebola virus disease outbreaks. Sci Rep. 2017;7.
Holmes EC. COVID-19—lessons for zoonotic disease. Science. 2022;375:1114–5.
Saúde OP-A da. Relatório 30 anos de SUS, que SUS para 2030? Relatório 30 anos SUS, que SUS para 2030? 2018.
Raphael Tsavkko Garcia. Census Delays Threaten Latin America’s Vulnerable. Bloomberg. 2021; https://www.bloomberg.com/news/articles/2021-07-21/census-delays-threaten-latin-america-s-vulnerable.
Silva Junior CHL, Pessôa ACM, Carvalho NS, Reis JBC, Anderson LO, Aragão LEOC. The Brazilian Amazon deforestation rate in 2020 is the greatest of the decade. Nat Ecol Evol. 2021;5:144–5. https://doi.org/10.1038/s41559-020-01368-x.
Spring J. Brazil spent less than half its 2021 environmental enforcement budget. Reuters. 2022: https://www.reuters.com/world/americas/brazil-spent-less-than-half-its-2021-environmental-enforcement-budget-2022-02-01/. Accessed in 4 Feb 2022.
Fornace KM, Diaz A V., Lines J, Drakeley CJ. Achieving global malaria eradication in changing landscapes. Malar J. 2021;20.
Pandit PS, Doyle MM, Smart KM, Young CCW, Drape GW, Johnson CK. Predicting wildlife reservoirs and global vulnerability to zoonotic Flaviviruses. Nat Commun. 2018;9.
We would like to thank A. Townsend Peterson for providing a free online and open course that provided the fundamental tools for the first author to run the analysis, especially in an underfunded Global South country. We would like Luke Perry, Paula R. Prist and Dr. Rodrigo Gurgel-Gonçalves for advice in an early version of this manuscript. Mauricio Vancine helped providing R code and advice in spatial analysis.
ARM is supported by CAPES Grant—Finance Code 88887.633837/2021-00. LEE was supported by the National Science Foundation under Grant HEGS 2116748 and Institute for Critical Technology and Applied Science, Virginia Tech: ICTAS-JFP-2022–2023. CTC was supported by CNPQ-MCTIC- SINBIOSE 442357/2019-2. JCS considers this work a contribution to his VILLUM Investigator project “Biodiversity Dynamics in a Changing World” funded by VILLUM FONDEN (Grant 16549) and JCS’ Center for Ecological Dynamics in a Novel Biosphere (ECONOVO), funded by Danish National Research Foundation.
Ethics approval and consent to participate
Consent for publication
Authors declare that they have no competing interests.
This material comprises several components, including an assessment of potential limitations, response curves for all diseases analyzed, both in composite and simple models, an exhaustive list of the variables employed in the analysis, and a comprehensive model reproducibility checklist containing essential details about disease occurrence and data processing.
This file includes a comprehensive tutorial on the construction of niche models and maps, utilizing the R programming language with instructive guidance.
About this article
Cite this article
Magalhães, A.R., Codeço, C.T., Svenning, JC. et al. Neglected tropical diseases risk correlates with poverty and early ecosystem destruction. Infect Dis Poverty 12, 32 (2023). https://doi.org/10.1186/s40249-023-01084-1