Skip to main content
Advertisement
  • Loading metrics

Leptospirosis in Rio Grande do Sul, Brazil: An Ecosystem Approach in the Animal-Human Interface

  • Maria Cristina Schneider ,

    schneidc@paho.org

    Affiliation Department of Communicable Diseases and Health Analysis, Pan American Health Organization, Washington, D.C., United States of America

  • Patricia Najera,

    Affiliation Department of Communicable Diseases and Health Analysis, Pan American Health Organization, Washington, D.C., United States of America

  • Martha M. Pereira,

    Affiliation Laboratório de Referência Nacional para Leptospirose, Centro Colaborador da Organização Mundial da Saúde, Fundação Oswaldo Cruz, Rio de Janeiro, Brazil

  • Gustavo Machado,

    Affiliation Faculdade de Medicina Veterinária, Departamento de Medicina Veterinária Preventiva, Laboratório de Epidemiologia Veterinária, Universidade Federal do Rio Grande do Sul, Porto Alegre, Rio Grande do Sul, Brazil

  • Celso B. dos Anjos,

    Affiliation Secretaria de Saúde do Rio Grande do Sul, Porto Alegre, Rio Grande do Sul, Brazil

  • Rogério O. Rodrigues,

    Affiliation Instituto de Pesquisas Veterinárias Desidério Finamor (IPVDF), Fundacão Estadual de Pesquisa Agropecuária (FEPAGRO), Eldorado do Sul, Rio Grande do Sul, Brazil

  • Gabriela M. Cavagni,

    Affiliation Secretaria da Agricultura Pecuária do Sul, Porto Alegre, Rio Grande do Sul, Brazil

  • Claudia Muñoz-Zanzi,

    Affiliation Division of Epidemiology and Community Health, School of Public Health, University of Minnesota, Minneapolis, Minnesota, United States of America

  • Luis G. Corbellini,

    Affiliation Faculdade de Medicina Veterinária, Departamento de Medicina Veterinária Preventiva, Laboratório de Epidemiologia Veterinária, Universidade Federal do Rio Grande do Sul, Porto Alegre, Rio Grande do Sul, Brazil

  • Mariana Leone,

    Affiliation Department of Communicable Diseases and Health Analysis, Pan American Health Organization, Washington, D.C., United States of America

  • Daniel F. Buss,

    Affiliation Laboratório de Avaliação e Promoção da Saúde Ambiental, Instituto Oswaldo Cruz, Fundação Oswaldo Cruz, Rio de Janeiro, Brazil

  • Sylvain Aldighieri,

    Affiliation Department of Communicable Diseases and Health Analysis, Pan American Health Organization, Washington, D.C., United States of America

  • Marcos A. Espinal

    Affiliation Department of Communicable Diseases and Health Analysis, Pan American Health Organization, Washington, D.C., United States of America

Abstract

Background

Leptospirosis is an epidemic-prone neglected disease that affects humans and animals, mostly in vulnerable populations. The One Health approach is a recommended strategy to identify drivers of the disease and plan for its prevention and control. In that context, the aim of this study was to analyze the distribution of human cases of leptospirosis in the State of Rio Grande do Sul, Brazil, and to explore possible drivers. Additionally, it sought to provide further evidence to support interventions and to identify hypotheses for new research at the human-animal-ecosystem interface.

Methodology and findings

The risk for human infection was described in relation to environmental, socioeconomic, and livestock variables. This ecological study used aggregated data by municipality (all 496). Data were extracted from secondary, publicly available sources. Thematic maps were constructed and univariate analysis performed for all variables. Negative binomial regression was used for multivariable statistical analysis of leptospirosis cases. An annual average of 428 human cases of leptospirosis was reported in the state from 2008 to 2012. The cumulative incidence in rural populations was eight times higher than in urban populations. Variables significantly associated with leptospirosis cases in the final model were: Parana/Paraiba ecoregion (RR: 2.25; CI95%: 2.03–2.49); Neossolo Litolítico soil (RR: 1.93; CI95%: 1.26–2.96); and, to a lesser extent, the production of tobacco (RR: 1.10; CI95%: 1.09–1.11) and rice (RR: 1.003; CI95%: 1.002–1.04).

Conclusion

Urban cases were concentrated in the capital and rural cases in a specific ecoregion. The major drivers identified in this study were related to environmental and production processes that are permanent features of the state. This study contributes to the basic knowledge on leptospirosis distribution and drivers in the state and encourages a comprehensive approach to address the disease in the animal-human-ecosystem interface.

Author Summary

Leptospirosis is a bacterial disease affecting humans and several animal species, which serve as reservoirs for infection. Contamination happens through exposure to urine of infected animals in water and soil.

A better understanding of the factors that affect the transmission of the disease, incorporating the relationship between humans, animals, and ecosystems within the One Health approach, would allow for tailored prevention and control measures at the local level. The aims of this study were to analyze the distribution of leptospirosis human cases and the possible factors that influence the transmission of the disease in the state of Rio Grande do Sul, as well as provide evidence to support the development of interventions and guide new studies in the region.

Leptospirosis cases and possible environmental, socioeconomic, and livestock factors were analyzed. The study used data by municipality (all 496) obtained via an open access database. Georeferenced maps were constructed and statistical analysis performed for 26 variables. A multivariable regression analysis was carried out.

An annual average of 428 human cases of leptospirosis was reported in the state from 2008 to 2012. The cumulative incidence in rural populations was eight times higher than in urban populations. Variables significantly associated with leptospirosis cases in the final model were: ecoregion, type of soil and to a lesser extent tobacco and rice production.

Results showed that leptospirosis cases were concentrated in the state capital and in specific ecoregions. Environmental and agricultural production processes were identified as possible risk factors and these conditions are probably permanent characteristics of the state. This study contributes to the knowledge on leptospirosis distribution and risk factors in the state, and also highlights the importance of a holistic view and intersectoral approach in preventing the disease.

Introduction

Even though leptospirosis is now recognized as a disease of epidemic potential with a significant health impact in many parts of the world, it remains a neglected disease. Its burden is estimated at 500,000 severe human cases per year worldwide, but a WHO expert group recently put its annual global incidence at 1.03 million people with 58,900 deaths [1,2]. Nevertheless, leptospirosis continues to be a silent disease [3], mainly due to the paucity of data in many countries, including of the Americas [4,5].

The Leptospira bacteria may also affect a wide variety of animal species, both wild and domestic, which serve as reservoirs for human infection [6]. The diversity of animal carriers represents a significant challenge for prevention and control [7]. Exposure to water and soil contaminated by the urine of infected animals is the most common route of transmission to people and domestic animals [6].

Leptospirosis is an excellent example for the “One Health” approach, where the relationship between humans, animals and ecosystems is studied to improve knowledge on a disease and to enhance collaborative intersectoral and multidisciplinary control strategies [8]. The One Health working definition states that “it is feasible to integrate human, animal, and environmental health efforts to predict and control certain diseases at the human–animal–ecosystem interface; integrated approaches that consider human, animal, and environmental health components can improve prediction and control of certain diseases” [9]. Yet this approach is rarely used to advance knowledge on leptospirosis transmission, develop evidence-based policies, and create tools to save human lives and reduce the impact on domestic animals.

Leptospirosis is one of the major neglected diseases in Latin America [10]. It has been reported in a variety of settings, from large urban centers, to remote rural areas [11,12,13,14,15]. Socioeconomic drivers include living in dense urban or peri-urban areas with inadequate waste collection and sanitation, which is often associated with vulnerable populations [16]. Leptospirosis has been linked to poverty, lack of water and sanitation, and poor housing conditions [12,17].

Environmental drivers have been identified in previous studies. Heavy rains or floods have been linked to a higher number of cases of leptospirosis [12,16,18,19,20,21]. Alkaline and neutral soil are suspected of promoting a longer survival of the bacteria, especially in young, not-yet-structured soils like those of volcanic origin [12,22]. In addition, soil temperature and proximity to water bodies were also reported as potential enablers for bacterial survival [23]. Leptospirosis is also considered an occupational disease, affecting rice laborers, sewer workers, animal handlers and gold miners [6,24,25].

A better understanding of the drivers for leptospirosis would provide crucial information for decision-makers to be able to target risk areas for priority interventions. Indeed, the current gaps in scientific and technological knowledge hinder the detection of cases and limit surveillance and control programs. Finally, because the instruments needed for control and elimination–such as broad rapid tests and vaccines–are not available at this time, leptospirosis is not considered to be “tool ready” and as a result, it is not targeted through large-scale global initiatives.

Brazil is the fifth most populous country in the world (approximately 200 million people) and has the seventh highest gross domestic product (USD 2,250,673). Even though the economy has been growing steadily, there are still around 20 million people (10% of the population) living in poverty in the country [26]. This population group is the most vulnerable to neglected diseases and other poverty-related infections. Previous studies have demonstrated the impact of neglected tropical diseases in Brazil and the need to develop new tools and technologies to fight them [27,28].

There is no accurate surveillance system in place for leptospirosis globally or in the Americas; however, some countries’ surveillance systems include the disease and estimates of its public health burden are available [5,29]. Notification of human leptospirosis is mandatory in Brazil and an annual average of 3,888 cases with 9.48% fatality are officially reported by the country’s surveillance system [30].

The Brazilian state of Rio Grande do Sul ranks fifth in the incidence rate (4.7 cases per 10,000 population) and presents around 15% of the total number of cases in the country [31]. In a previous study conducted to identify high transmission areas and possible ecological components of leptospirosis transmission in Rio Grande do Sul, the highest incidence rates were found in the coastal sedimentary areas with low altitude and predominantly agricultural land use in the central valley [32]. The state economy is based on agribusiness, including cattle and rice paddies, with an associated increased risk of leptospirosis in some areas that needs to be evaluated.

The objective of this study is to analyze the distribution of human cases of leptospirosis in the state of Rio Grande do Sul and explore possible drivers using the One Health approach. This analysis may orient further studies of the disease in the human-animal-ecosystem interface and the results of this study may be used as evidence to inform decision makers in the state.

Methods

Study design and data collection

An ecological study was carried out using aggregated data by municipality to analyze the situation of leptospirosis in all 496 municipalities (corresponding to the second subnational administrative level) in the state of Rio Grande do Sul, Brazil, between 2008 and 2012.

A geo-coded database was created using different sources. Variables were either downloaded or created from original sources. The data source used for each variable is described in the S1 Supporting Information.

Human leptospirosis case data, de-identified and aggregated at the municipality level, were obtained from the Ministry of Health of Brazil national surveillance system database (acronym in Portuguese SINAN) [33]. All case data were publicly available by open consultation on the government website.

Seven environmental variables (ecoregion, type of soil, temperature, precipitation of the wettest month, altitude, slope of the land (hill incline), and drainage) were gathered from diverse “open-access” data sources and used in the study. When disaggregated, the ecoregion and type of soil variables turned into six and twelve variables, respectively. Altitude and hydrology information used in the background were obtained from the USGS-EROS, HYDRO1k Elevation Derivative Database [34]. Bioclimatic variables were calculated from monthly temperatures and rainfall values [35]. Geo-processing techniques were applied to assign and measure environmental variables for each municipality. Data for the socioeconomic variables (gross domestic product per capita, Gini coefficient, illiteracy rate) were gathered from the Brazilian Institute of Geography and Statistics (acronym in Portuguese IBGE) [36]. As the economy of Rio Grande do Sul is based on agribusiness, variables related to rice and tobacco production were collected from the IBGE database and data on bovine raising were provided by the Department of Agriculture of the state of Rio Grande do Sul and the Federal University of Rio Grande do Sul [37].

Definitions

Leptospirosis cases: According to the Ministry of Health of Brazil, human cases of leptospirosis are those presenting clinical symptoms consistent with the disease and confirmed by laboratory diagnosis either with ELISA–IgM or MAT [38]. These laboratory confirmation techniques are available at the state level at the Central Public Health Laboratory, which is part of the National Public Health Laboratory Network. In Brazil, a case could also be confirmed by clinical-epidemiological criteria (selected symptoms with epidemiological history) [39]. All cases in the SINAN database were considered confirmed. In the original database, the cases were classified according to area of residence (urban, peri-urban or rural) [39].

Cumulative incidence: the number or proportion of a group (cohort) of people who experience the onset of a health-related event during a specific time interval [40]. In this study, it was estimated per 10,000 inhabitants.

Criteria for risk stratification: The criteria used were based on a previous risk stratification study conducted in Nicaragua [12]. In summary, geographical areas were classified into two categories: i) Silent Area (no cases were reported during the study period) and ii) Productive Area (active transmission was reported during the study period) that could be an endemic or critical area (higher quintile).

Study area

Rio Grande do Sul is the southernmost state in Brazil bordering Argentina and Uruguay. The state covers an area of 281,731,445 km², divided into 496 municipalities. In 2010, the population was 10,693,929 inhabitants, 85.1% of which lived in urban areas, among them the Porto Alegre metropolitan area where 15% of the state population resides (1,472,482 inhabitants) [41].

In 2010, the state had the 4th highest gross domestic product per capita in the country and the Gini index (0.5472) was lower than the national level [36]. There are three major economic regions: 1) the south with greater land concentration, large cattle raising farms, and mechanized plantation of rice, soybean and wheat. This area also presents higher income inequality; 2) the northeast region, which includes the state capital, with more industries and predominantly small properties; and 3) the northern region, mostly colonized by European immigrants, with higher forest coverage, valleys, and plain areas with small agricultural lands. The most common agrarian structure in the state (90% of properties) is a small family farm covering an area smaller than 100 hectares [42].

The hydrology of the state is basically divided into two main areas: the La Plata and Atlantic East Coast watersheds. The La Plata watershed is located in the northern area bordering with Argentina, and comprises the Uruguay River and mayor tributaries (Caboa, Pelotas, Ibicui and Mirinay). The Atlantic East Coast watershed is mainly shaped by the Guaiba and Litoranea basins. The Guaiba basin includes tributaries of the Dos Patos coastal lagoon, the Jacui and Tacuari rivers, which run in the central area of the state and wash its most populated areas; other tributaries include the Sinos and Cai rivers, which also flow into the Dos Patos lagoon. (S2 Supporting Information). The Litoranea basin encloses the Camaqua and Piratini rivers also flowing into the lagoon. The city of Porto Alegre is located in the Atlantic East Coast watershed.

Data analysis

Two types of analyses were carried out to investigate associations between the risk of human leptospirosis and 26 variables selected as possible environmental, socioeconomic or livestock drivers: 1) a spatial analysis and thematic mapping of possible drivers, showing the municipalities in the higher quintiles of rates over the distribution of the variable also divided by quintiles (range cuts) and 2) a statistical analysis with univariate and multivariable regressions as described below.

The spatial analysis consisted in computing zonal statistics and surface for the environmental variables, in addition to the geocoding of health and socioeconomic data. ArcGIS zonal statistics by municipality (min, mean, max, standard deviation, range) were calculated for the altitude, slope, temperature, and rain variables. Geo-processing geometric intersection of environmental features shaped the municipal surface of ecoregions and soils. Quintile thematic mapping was conducted once the municipal statistics of environmental, health and socioeconomic variables were geo-processed.

For the multivariable regression, the dependent variable was the case count per municipality and the independent variables all 26 possible drivers. The ecoregion and soil variables were dichotomized as follows: they were coded 0 when they were a minority (when land coverage or proportion was less than or equal to 50% of the municipality) and they were coded 1 when they represented a majority (when land coverage or proportion was above 50.01% of the municipality.) The proportion of farms with up to ten animals per farm (small farms) was dichotomized in the same fashion. In addition, two bovine-related variables were created and used as continuous variable: the proportion of farms per km2 and the proportion of bovines per km2. Productive processes such as the cultivation of tobacco and rice were analyzed on a unit of 10,000 tons.

Predictors were analyzed with negative binomial regression (NB) and robust variance was used to estimate the relative risk (RR) and 95% confidence interval (CI) of the estimates [43]. The link function used was the default logit-function and the offset was the natural log-transformed number of population per municipality.

Univariate analyses were first run for each of all 26 variables and 14 were preselected due to P ≤ 0.15. Variance inflation factors (VIF) were estimated to verify the relations among all selected independent variables and check for potential collinearity. When a high VIF was found (VIF>4), the variable with a lower P-value was eliminated and the process was reiterated until only variables with a VIF<4 were left. Confounding effects were investigated by checking changes in the point estimates of the variables that remained in the model. Parameters with changes in their estimates > 25% were considered confounders and were retained during the variable selection process. Finally, two-way interaction terms between environmental variables with biological plausibility were investigated (slope and altitude, altitude and drainage, altitude and precipitation of the wettest month (mm)). We used deviance performance as a goodness of fit test for the overall model. Results of the variables with P value > 0.15 are included in the S3 Supporting Information.

Statistical tests were performed for over-dispersion evaluation and model fit comparison using literature-recommended approaches [44]. For the negative binomial model, the dispersion parameters were tested for difference with chi-squared statistics. To compare goodness of fit between pairs of the proposed regression models, Vuong statistics were calculated [45]. The following models were compared: Poisson regression, zero-inflated negative binomial (ZINB) and a negative binomial regression (NB) [46]. The differences in AIC and Vuong statistics were computed for all pairs of non-nested models (i.e., Poisson vs. ZINB, Poisson vs. NB, NB vs. ZINB). Based on the AIC and Vuong tests, the negative binomial regression model fit the data better in comparison with others. More details about the procedures carried out in order to compare models can be found in the S4 Supporting Information.

Results

Cases, cumulative incidence and risk stratification

During the study period, 2,141 confirmed cases of leptospirosis were reported with an average of 428 cases annually. The yearly incidence remained similar over the period (412–543 cases), except for a lower than average number of cases (277) in 2012 (Table 1).

thumbnail
Table 1. Human cases of leptospirosis by area of residence, Rio Grande do Sul, 2008–2012.

https://doi.org/10.1371/journal.pntd.0004095.t001

A total of 233 municipalities out of 496 reported at least one case (range 1–208) during the study period. Among the 263 municipalities not reporting cases, most had small populations (less than 10,000 people). On average, there were 4.32 cases per municipality.

A high number of urban cases were clustered in the Metropolitan Region of Porto Alegre (50.52% of all urban cases) and groups of predominantly rural cases were located in the Center Oriental Region along the Jacui and Tacuari rivers and in the Southeast Region, near the great lagoons (Fig 1). Forty-seven municipalities concentrate 79.40% of the cases. Most of these municipalities are located in the Atlantic East Coast watershed, specifically in the Jacui-Tacuari and Cai basins (1,116 cases).

thumbnail
Fig 1. Human cases of leptospirosis according to the residence site, by municipality, Rio Grande do Sul, 2008–2012.

https://doi.org/10.1371/journal.pntd.0004095.g001

The cumulative incidence for the entire state was 2 cases per 10,000 inhabitants (range: 0–56.56). The spatial distribution of the municipalities’ cumulative incidence by quintiles, as well as of those not reporting cases, are presented in Fig 2. The cumulative incidence for rural areas (6.16 per 10,000 people) was eight times higher than for urban areas (0.74 per 10,000 people).

thumbnail
Fig 2. Cumulative incidence of leptospirosis (10,000 habitants), by municipality, Rio Grande do Sul, 2008–2012.

https://doi.org/10.1371/journal.pntd.0004095.g002

Risk stratification was carried out for the 496 municipalities in Rio Grande do Sul to inform planning of leptospirosis prevention and control activities. It shows that out of the 496 municipalities in the state, 263 (53.02%) can be considered silent areas. Among the 233 productive areas (46.98%), 58 were found to be endemic and 75 municipalities were considered critical areas (S5 Supporting Information). These 75 municipalities reported 1,766 cases (82.48%) out of the total of 2,141 cases in the entire state.

Exploratory analysis of possible drivers

Environmental variables.

The state of Rio Grande do Sul has six ecoregions. The largest extension is the Uruguayan savanna, covering 230,291.4 km² or 64.41% of the state territory and located in the southwest region bordering with Argentina and Uruguay, followed by the Parana-Paraiba interior forests with 63,504.9 km² (17.76%). More towards the center, the Brazilian Araucaria moist forest covers 57,758.2 km² (16.5%) and borders the State of Santa Catarina (north), and small areas of the Brazilian Atlantic coast (4,181.4 km² (1.17%)) and the Serra do Mar coastal forests (1,678.60 km² (0.47%)) (Fig 3). The Mesopotamian savannas bordering Argentina along the Uruguay River represent a very small area (113.1 km² (0.03%)). A visual exploration of the distribution of leptospirosis rates by quintiles under the ecoregions suggests that municipalities in the higher quintiles are concentrated within the Parana-Paraiba interior forests. A summary of the distribution of these variables is shown in Table 2.

thumbnail
Fig 3. Cumulative incidence for leptospirosis and ecoregions, Rio Grande do Sul, 2008–2012.

https://doi.org/10.1371/journal.pntd.0004095.g003

thumbnail
Table 2. Possible drivers and leptospirosis cases per municipalities (P<0.15 in the univariable and multivariable analysis and P<0.05), Rio Grande do Sul, 2008–2012.

https://doi.org/10.1371/journal.pntd.0004095.t002

A closer look at the bioclimatic variables yielded medians and ranges of the monthly temperature in Celsius degrees and precipitation of the wettest month in millimeters of 18.85°C (14.53–20.92) and 173.1 mm (114.6–203.2), respectively. The average altitude in the state is 352 meters above sea level (range: 2.61–1,169.28) with 25% of the municipalities at an altitude above 157 meters. The interaction between altitude and slope showed to be significantly related to higher number of cases of leptospirosis in the univariate analysis (Table 2).

According to the Brazilian agency for soils classification [47], there are 12 types of soil in the state. The local Neossolo Litolitico (RL) soil prevails, covering 106.6 km² (23.40%), followed by Latosol Vermelho (91.29 km²) and Cambisol Haplico (52.88 km²). Neossolos Litoliticos are distributed near river valley walls surrounding the deeper and well developed Latosols soils. When analyzing the predominant types of soils individually by quintiles, the first visual impression is that a higher human leptospirosis cumulative incidence overlaps with Neossolo Litolitico (Fig 4).

thumbnail
Fig 4. Cumulative incidence for leptospirosis type of soil and neossole type, by municipality Rio Grande do Sul, 2008–2012.

https://doi.org/10.1371/journal.pntd.0004095.g004

Socioeconomic variable.

The average gross domestic product (GDP) per capita in the Rio Grande do Sul State was 19,646.61 Brazilian Real (range: 7,020.39–224,004.4, median 17,121.95). The Gini index measuring inequality, in this case related to the domestic income per capita by municipality, was 0.5472 (range: 0.2841–0.73248). The illiteracy rate in the state was 6.56% (range: 0.9–19.1%, median 6.1%). The univariate analysis suggested some relation between the Gini index and the number of cases (P≤0.15), however this relationship was not observed in the final model.

Productive process.

During 2010, the state of Rio Grande do Sul produced 6,875,077 tons of rice, with a mean of 1,3861.04 and median of 6,5 tons (range: 0–507,788). There were rice paddies in 337 of the 496 municipalities (67.94%). A large portion of the state, the entire southern half, includes municipalities in the highest quintile of rice production (Fig 5). When comparing thematic maps on rice production and leptospirosis cumulative incidence, the overlapping areas suggest a potential relationship between the two variables.

thumbnail
Fig 5. Cumulative incidence for leptospirosis, rice paddy plantation, and tobacco plantation, by municipality, Rio Grande do Sul, 2008–2012.

https://doi.org/10.1371/journal.pntd.0004095.g005

The other agricultural production analyzed was tobacco, which is grown in 319 municipalities (64.31%) with a total of 343,682 tons per year (mean: 692.91, median: 16 and range: 0–21,080). The municipalities in the highest quintile for tobacco production are located in the central and northern areas of the state, especially at lower altitudes (Fig 5). The thematic map overlapping tobacco production and leptospirosis also suggests some correlation.

All municipalities in the state have livestock, in this case bovine, with a total of 13,541,230 animals (mean: 27,300.86, median: 9,482.5, and range: 91–633,381). A concentration of municipalities in the highest quintile of bovine per property was identified in the southwest along the border with Uruguay and Argentina (Fig 6). In the visual analysis of the number of bovine by property, a possible overlap appeared between a high number of properties with less than 10 bovines and leptospirosis cumulative incidence in the highest quintile (Fig 6). The univariate analysis suggested some relation between the number of properties with less than 10 bovines and leptospirosis cases (P≤0.15), however this relationship was not observed in the final model.

thumbnail
Fig 6. Cumulative incidence for leptospirosis and percentage of properties with up to 10 bovines, by municipality, Rio Grande do Sul, 2008–2012.

https://doi.org/10.1371/journal.pntd.0004095.g006

Multivariable analysis

The negative binomial regression modeling identified fourteen variables that suggested some relation with the number of cases per municipalities in the univariate analysis (P≤0.15) (Table 2). The 12 variables that were analyzed as possible drivers (P>0.15) are presented in the S3 Supporting Information.

The final model identified four variables as significantly associated with the leptospirosis case count (P ≤ 0.05) (Table 2): the Parana/Paraiba interior forest ecoregion (RR = 2.25; IC95% = 2.03–2.49; P < 0.001); Neossolo Litolítico soil (RR = 1.93; IC95% = 1.26–2.96; P = 0.006); to a weaker extent, production of tobacco per 10,000 tons (RR = 1.10; IC95% = 1.09–1.11; P <0.001); and finally, a borderline association with production of rice per 10,000 tons (RR = 1.003; IC95% = 1.002–1.04; P <0.001). The goodness-of-fit model was tested by a deviance chi-squared test and was found not to be significant (P > 0.05).

Discussion

Leptospirosis is an important public health problem for the state of Rio Grande do Sul, with an average of 428 cases reported annually. The actual incidence is probably much higher, particularly in vulnerable populations, because leptospirosis is commonly misdiagnosed and unreported.

The study results show that rural populations of the state have an approximately eight times higher risk of contracting leptospirosis than their urban counterparts, even though the number of reported cases was high in both areas. Rural cases are concentrated mostly in two zones, characterized by higher tobacco production (central state) and higher rice production (south). This result suggests that preventive strategies in these regions require strong collaboration between the Health and the Agriculture sectors. Preventing cases among rural workers will be the main goal for these areas.

In order to contribute to the One Health discussion and support possible intersectoral collaboration, suggestions for components of a leptospirosis plan and activities in the state of Rio Grande do Sul are presented in the S6 Supporting Information; using both the working definition of One Health and the results of this study as a basic scenario [9].

The highest concentration of urban cases was found in the metropolitan areas of cities with low altitude, including the state capital, which also suffer from frequent flooding. At the center of the state, high rate areas were found within the Atlantic East Coast watershed, which receive water from the Jacui, Taquari, Cai, Sinos and others rivers as mayor tributaries, all well known for their periodic flooding [48]. This suggests that leptospirosis prevention should also be considered in natural disaster plans and that the main goals for these areas should be to reduce the number of severe cases and save lives during an outbreak.

Other possible drivers for leptospirosis identified in the final regression model were the Parana-Paraiba ecoregion, Neossolo Litolítico soil, and the production of tobacco. The Parana/Paraiba interior forest ecoregion showed an RR of 2.25 when compared to other ecoregions. This ecoregion is part of the Atlantic semi-deciduous forest biome that covers central and northern areas of the state and extends over numerous areas of Brazil and neighboring countries. It ranges from river plains to middle-level highlands. The climate is subtropical [49] so humidity may act as an important factor in the leptospirosis bacteria survival. The Parana-Paraiba ecoregion includes a diversity of soil types that vary from very fertile to impoverished. Neossolo Litolítico, a very young and not yet structured type of soil, is one of those that correlate geographically with high rates of human leptospirosis. Its RR is 1.93 when compared with municipalities with other predominant soil type. It has been suggested that its pH [50] and/or its low drainage capacity [51] may favor the bacteria survival.

With each additional 10.000 tons of tobacco produced by the municipality, there was an increased RR in the leptospirosis cases of 1.10. There is spatial overlapping of municipalities with high production of tobacco and higher rates of the disease in the Jacui-Tacuari valley. In the central area, several conditions concur with higher leptospirosis cases, in particular tobacco (and rice) production in a subtropical forest environment, over slope and inclined terrain with shallow weathered soils.

We did not find information in the literature suggesting an association between tobacco plantation and leptospirosis cases. However, tobacco growth requires a pH of 5–6.5 [52], which matches the leptospirosis bacterial soil requirement of 6.2 for survival over 7 weeks in the environment [22,53]. PH survival requirements vary among serovars, and as an example L. interrogans serovar hardjo prefers a pH of 6.5–6.8 [54]. Further studies analyzing which serovars are circulating in humans and in animals in these high risk areas will provide additional knowledge on the cycle of leptospirosis transmission.

The agro-industrial production and export of tobacco are concentrated in three states in the south of Brazil. Rio Grande do Sul and the neighboring state of Santa Catarina are responsible for 80% of the production. Brazilian tobacco is grown by approximately 186,000 farming families in properties with an average area of 16 hectares [55]. Every step in the plantation of tobacco in Southern Brazil is carried out manually, mostly by family members from several generations of small farmers descended from Europeans immigrants [55]. Most of the farms have a small number of bovines, mostly dairy cows, although swine are also very frequent [37]. In this region, the most common feed for livestock is corn, mostly stocked on the farm where it attracts rodents. Insects and rodents are responsible for an estimated loss of 15% of all corn storage on small farm barns [47]. This region also includes large areas of preserved forests, which probably increases contact with wild rodents. Further studies about the presence of rodents and programs for rodent control would provide important information on their true impact.

In the final model, the association with rice production was weak (around 1 but statically significant), but there was some spatial overlapping of municipalities with high leptospirosis cumulative incidence and rice plantations in the Lagoa dos Patos shoreline. The state of Rio Grande do Sul produces 63% of all rice in Brazil [56]. The association between rice paddy workers and leptospirosis has already been described in the literature [6,57]. Working conditions in rice fields have been studied in Peru and results highlight occupational hazards linked to leptospirosis infection, such as long period of exposure to water, lack of use of any personal protection equipment, and presence of skin wounds [57]. It is very important to raise awareness about the risk of this disease in critical areas and promote personal protection equipment (PPE).

In the state of Rio Grande do Sul, where temperatures exceed 24°Celsius (75°Fahrenheit) for half of the year, educational campaigns for the adequate use of PPE (such as boots) among workers are unlikely to be effective, especially during long shifts. The most important prevention tool for leptospirosis in high risk areas and especially under these conditions is vaccination, however, a human vaccine is not available in the majority of affected countries.

In the central areas of the state, smaller livestock producers (less than 10 bovines by property) are predominant [37], but this variable only presented relevance in the univariate analysis. Similarly, the proportion of bovine farms per km2 has shown some relation with the number of leptospirosis cases. The proximity with animals such as dogs increases the risk for leptospirosis [58]. A previous study found that leptospirosis titers were present in 39% of bovines in the state [59]. Further studies are required to better understand the epidemiological situation in the animal sector.

This study provides evidence of the importance of the environmental component and the One Health approach for leptospirosis. As highlighted here, where the ecoregions were suitable for specific plantations (such as tobacco in Parana-Paraiba) human occupation of these territories has led to transformations of the natural environment and increased the use of natural resources for economic purposes [60]. Further research is required, using primary data, to further explore this finding, as well as the potential relationship between the type of ecoregion and economic processes, and their relationship with the animal-human interface. For other diseases, such as rabies transmitted by vampire bats, outbreaks have been reported rapidly in the Amazon region due to some productive processes such as gold prospection [61].

Rio Grande do Sul has one of the highest GDP per capita in the country [62]. For this reason, the study includes the Gini index that reflects possible inequalities related to the distribution of income among the inhabitants of the municipalities. However, the Gini index presented relevance only in the univariate analysis and no association was found in the final model.

This results of this study provide potentially useful information for government decision-makers in the elaboration of a possible intersectoral plan. Evidence-based policy-making requires advanced social planning capacity, including the capacity to monitor and evaluate programs in the different ministries and agencies that provide services to the most vulnerable populations. Risk stratification generates evidence for decision-making and the critical areas identified by this study could be used as pilots for future local intersectoral studies. In the same way, areas identified as silent require attention from the health authorities to confirm the absence of disease in humans or the need to improve surveillance. The prediction, detection, prevention, and response to outbreaks of leptospirosis will be improved and better defined through knowledge generated by an integrated approach within the animal-human-ecosystem interface.

This study clearly highlights the importance of a multidisciplinary and intersectoral approach, which was at the root of the creation of initiatives such as the Global Leptospirosis Environmental Action Network (GLEAN). GLEAN brings together experts from different disciplines and sectors with the objective of orienting knowledge regarding the prediction, prevention, detection, and response to leptospirosis; and transforming relevant research findings into operational tools for affected communities and countries, as well as for organizations [63].

Although the use of aggregated data by municipality is a possible limitation of this study, the geographical approach helped to gather, unify and shape all factors (environmental, socioeconomic, and health) into the same area and detect spatial patterns. Ecological studies are an inexpensive way to use secondary available data; however they are commonly associated with the ecological fallacy [21]. Another limitation possibly resides in the large number of municipalities with small populations (under 10,000 people). However, the possible effect of this data characteristic was considered in the statistical analysis through the use of negative binomial regression, by testing the effects using strata with and without small population. Also, no significant difference was found by analyzing the same approach using only municipalities with at least one positive case.

This study contributes to the basic knowledge on leptospirosis distribution and drivers in the state of Rio Grande do Sul and points to the direction of tool-readiness for the disease. It also encourages a holistic approach for the disease, taking into account human, animal, and ecosystem interactions and supporting research in the development and validation of a rapid diagnostic test for humans and animals, and eventually production of a vaccine that is still desperately needed.

Supporting Information

S2 Supporting Information. Map of hydrology and leptospirosis cumulative incidence.

https://doi.org/10.1371/journal.pntd.0004095.s002

(DOCX)

S3 Supporting Information. Table of possible drivers and leptospirosis rates (P>0.15).

https://doi.org/10.1371/journal.pntd.0004095.s003

(DOCX)

S4 Supporting Information. Details about modeling comparison and verifications.

https://doi.org/10.1371/journal.pntd.0004095.s004

(DOCX)

S6 Supporting Information. Components and activities for a leptospirosis plan in the state of Rio Grande do Sul.

https://doi.org/10.1371/journal.pntd.0004095.s006

(DOCX)

Acknowledgments

The authors thank Diego Viali dos Santos for his vital participation in the initial phase of this study and for supporting the original proposal, Deise Galan for assisting in the elaboration of the proposal with the same title drafted in 2013 and for bibliographical review; the Rio Grande do Sul Agriculture authorities and Health authorities, specially Sirlei Famer from CEVS/SES-RS, and colleagues from both institutions. We would also like to thank Roberto Bazzani from the IDRC for his call for ‘Policy-oriented Research on Agriculture, Environment, and Health in Latin America and the Caribbean’; Eduardo Caldas from SVS/MS for supporting leptospirosis surveillance and control as part of the Rio Grande do Sul State Department and currently as part of the Ministry of Health, and for calling the attention about the Taquari valley river as a higher risk area; Paulo Buss for supporting the study and the One Health approach in the Rede Saúde-RS; Eric Bertherat from WHO for drawing the attention to this important disease; for Ottorino Cosivi and colleagues from PANAFTOSA for their constant collaboration; and Anne Roca for editing and reviewing this article.

Authors’ Disclaimer

The authors are staff members from different organizations, however they alone are responsible for the views expressed in this article and they do not necessarily represent the decisions or the stated policy of their organization.

Author Contributions

Conceived and designed the experiments: MCS PN MMP. Performed the experiments: PN MCS GM ML. Analyzed the data: MCS PN GM LGC CMZ. Contributed reagents/materials/analysis tools: PN GM. Wrote the paper: MCS PN MMP GM CBdA ROR GMC CMZ LGC ML DFB SA MAE.

References

  1. 1. World Health Organization. Report of the First Meeting of the Leptospirosis Burden Epidemiology Reference Group; 2009 Dec 2–4; Geneva Switzerland. Geneva: WHO; 2010.
  2. 2. Hagan JE, Costa F, Calcagno J, Kane M, Torgerson P, Martinez-Silveira MS et al. Global Morbidity and Mortality of Leptospirosis: A Systematic Review. Proceedings of the 8th International Leptospirosis Society Scientific Meeting; 2013 Oct 8–11; Fukuoka, Japan. Fukuoka: ILSS; 2013.
  3. 3. Schneider MC, Jancloes M, Buss DF, Aldighieri S, Bertherat E, Najera P et al. Leptospirosis: A Silent epidemic disease. Int J Environ Res Public Health. 2013; (10):7229–7234. pmid:24351743
  4. 4. Hotez P, Woc-Colburn L, Bottazzi ME. Neglected tropical diseases in Central America and Panama: Review of their prevalence, populations at risk and impact on regional development. Int J Parasitol. 2014 Aug; 44(9):597–603. pmid:24846528
  5. 5. Costa F, Martinez-Silveira MS, Hagan JE, Hartskeerl RA, Dos Reis MG, Ko AI. Surveillance for leptospirosis in the Americas, 1996–2005: a review of data from ministries of health. Rev Panam Salud Publica. 2012 Sep; 32(3):169–77. pmid:23183556
  6. 6. World Health Organization, International Leptospirosis Society. Human Leptospirosis: Guidance for Diagnosis, Surveillance and Control. Geneva, Switzerland: WHO; 2003.
  7. 7. Petrakovsky J, Bianchi A, Fisun H, Nájera-Aguilar P, Pereira MA. Animal Leptospirosis in Latin America and the Caribbean countries: Reported outbreaks and literature review (2002–2014). Int J Environ Res Public Health. 2014;11(10): 10770–89. pmid:25325360
  8. 8. Schneider MC, Aldighieri S, Najera P, Galan DI, Fisun H, Cosivi O et al. O conceito de “Uma Saúde” e sua aplicação. Rev CFMV. 2014; 62:18–21. Portuguese.
  9. 9. Rabinowitz PM, Kock R, Kachani M, Kunkel R, Thomas J, Gilbert J et al. Toward proof of concept of One Health approach to disease prediction and control. Emerg Infect Dis. 2013 Dec; 19(12).
  10. 10. Hotez PJ, Bottazzi ME, Franco-Paredes C, Ault SK, Roses Periago M. The neglected tropical diseases of Latin America and the Caribbean: a review of disease burden and distribution and a roadmap for control and elimination. PLoS Negl Trop Dis. 2008 Sept; 2(9):e300. pmid:18820747
  11. 11. Ko AI, Galvaõ RM, Ribeiro Dourado CM, Johnson WD, Riley LW. Urban epidemic of severe leptospirosis in Brazil. Salvador Leptospirosis Study Group. Lancet. 1999 Sep 4; 354(9181):820–5. pmid:10485724
  12. 12. Schneider MC, Nájera P, Aldighieri S, Bacallao J, Soto A, Marquiño W et al. Leptospirosis outbreaks in Nicaragua: Identifying critical areas and exploring drivers for evidence-based planning. Int J Environ Res Public Health. 2012; 9(11):3883–3910. pmid:23202822
  13. 13. Maciel EAP, de Carvalho ALF, Nascimento SF, de Matos RB, Gouveia EL, Reis MG et al. Household transmission of leptospira infection in urban slum communities. PLoS Negl Trop Dis. 2008 Jan; 2(1):e154. pmid:18357340
  14. 14. Lau CL, Clements AC, Skelly C, Dobson AJ, Smythe LD, Weinstein P. Leptospirosis in American Samoa—Estimating and mapping risk using environmental data. PLoS Negl Trop Dis. 2012; 6(5):e1669. pmid:22666516
  15. 15. Jancloes M, Bertherat E, Schneider MC, Belmain S, Munoz-Zanzi C, Hartskeerl R et al. Towards a “One Health” Strategy against Leptospirosis. Planet@Risk [Internet]. 2014; 2(3):204–206. Available from: http://planet-risk.org/index.php/pr
  16. 16. Barcellos C, Sabroza PC. The place behind the case: Leptospirosis risks and associated environmental conditions in a flood-related outbreak in Rio de Janeiro. Cad Saude Publica. 2001; 17 Suppl: 59–67. pmid:11426266
  17. 17. Bacallão J, Schneider MC, Najera P, Aldighieri S, Soto A, Marquiño W et al. Socioeconomic factors and vulnerability to outbreaks of leptospirosis in Nicaragua. Int J Environ Res Public Health. 2014 Aug 15; 11(8):8301–18. pmid:25153463
  18. 18. Sehgal SC, Sugunan AP, Vijayachari P. Outbreak of leptospirosis after the cyclone in Orissa. Natl Med J India. 2002 Jan-Feb; 15(1):22–3. pmid:11855587
  19. 19. Liverpool J, Francis S, Liverpool CE, Dean GT, Mendez DD. Leptospirosis: case reports of an outbreak in Guyana. Ann Trop Med Parasitol. 2008 Apr; 102(3):239–45. pmid:18348778
  20. 20. Amilasan AS, Ujiie M, Suzuki M, Salva E, Belo MC, Koizumi N et al. Outbreak of leptospirosis after flood, the Philippines, 2009. Emerg Infect Dis. 2012 Jan; 18(1):91–4. pmid:22257492
  21. 21. Weinberger D, Baroux N, Grangeon JP, Ko AI, Gorant C. El Niño Southern oscillation and leptospirosis outbreaks in New Caledonia. PLoS Negl Trop Dis. 2014 Apr 17; 8(4):e2798. pmid:24743322
  22. 22. Acha PN, Szyfres B. Zoonoses and Communicable Diseases Common to Man and Animals. In Vol I: Bacterioses and Mycoses. 3rd ed. Washington, DC: Pan American Health Organization; 2003. p. 157–168.
  23. 23. Henry RA, Johnson RC. Distribution of genus Leptospira in soil and water. Appl Environ Microbiol. 1978 Mar; 35(3):492–499. pmid:637546
  24. 24. Benschop J, Heuer C, Jaros P, Collins-Emerson J, Midwinter A, Wilson P. Sero-prevalence of leptospirosis in workers at a New Zealand slaughterhouse. N Z Med J. 2009 Dec 11; 122(1307):39–47. pmid:20148043
  25. 25. Bertherat E, Renaut A, Nabias R, Dubreuil G, Georges-Courbot MC. Leptospirosis and Ebola virus infection in five gold-panning villages in northeastern Gabon. Am J Trop Med Hyg. 1999 Apr; 60(4):610–5. pmid:10348236
  26. 26. Poverty & Equity, Country Dashboard, Brazil [Internet]. Washington, DC: The World Bank. c 2015—[cited 2015 Mar 3]. Available from: http://povertydata.worldbank.org/poverty/country/BRA.
  27. 27. Hotez PJ, Fujiwara RT. Brazil's neglected tropical diseases: an overview and a report card. Microbes Infect. 2014 Aug; 16(8):601–6. pmid:25088506
  28. 28. Schneider MC, Aguilera XP, Barbosa da Silva J Junior, Ault SK, Nareja P, Martinez J et al. Elimination of neglected diseases in Latin America and the Caribbean: A mapping of selected diseases. PLoS Negl Trop Dis. 2011 Feb 15; 5(2):e964. pmid:21358810
  29. 29. World Health Organization. Report of the Second Meeting of the Leptospirosis Burden Epidemiology Reference Group; 2010 Sep 22–23; Geneva, Switzerland. Geneva: WHO; 2011.
  30. 30. Caldas, E. Leptospirose no Brasil, Ministério da Saúde. 2013[Internet]. Apresentação durante a reunião anual da rede global de ação ambiental contra a leptospirose. Brasília. Brasil. Washington DC: WHO- [about 1 screen]. Available from: http://www.paho.org/hq/index.php?option=com_topics&view=rdmore&cid=4957&Itemid=40934&lang=en.
  31. 31. Caldas, E. Presentations in meetings supported by PAHO. 2012 [Internet]. Presentations of the International Meeting of Leptospirosis of Nicaragua. Managua, Nicaragua. Washington DC: WHO- [about 1 screen]. Available from: http://www.paho.org/hq/index.php?option=com_topics&view=readall&cid=5576&Itemid=40934&lang=en.
  32. 32. Barcellos C, Lammerhirt CB, de Almeida MA, dos Santos E. Distribuição espacial da leptospirose no Rio Grande do Sul, Brasil: recuperando a ecologia dos estudos ecológicos. Cad Saude Publica. 2003 Sep-Oct; 19(5):1283–92. Portuguese.
  33. 33. Sistema de informação de agravos de notificação [Internet]. Brasilia, Brasil: Ministério da Saúde. c 2014—[cited 2015 Feb 24]. Available from: http://dtr2004.saude.gov.br/sinanweb.
  34. 34. Science for a changing world (USGS) [Internet] [cited 2015 Mar 14]. Available from: https://lta.cr.usgs.gov/HYDRO1K.
  35. 35. WorldClim—Global Cimate data. [Internet]. California, United States of America. c 2015—[cited 2015 Feb 10]. Available from: http://www.worldclim.org/bioclim.
  36. 36. Informações de saúde. Datasus [Internet]. Brasilia, Brasil: Ministério da Saúde. c 2010—[cited 2015 Feb 24]. Available from: http://tabnet.datasus.gov.br/cgi/ibge/censo/cnv/ginirs.def.
  37. 37. Silva GS, Costa E, Bernardo FA, Groff FH, Todeschini B, Viali SD et al. [Cattle Rearing in Rio Grande do Sul, Brazil]. Acta Sci Vet. 2014; 42:1215. Portuguese.
  38. 38. Ministério da Saúde do Brasil, Secretaria de Vigilância em Saúde. Vigilância em Saúde: Zoonosis, 2009 [Internet]. Brasilia, Brazil: Ministério da Saúde; 2009 [cited 2015 Jan 15]. Available from: http://189.28.128.100/dab/docs/publicacoes/cadernos_ab/abcad22.pdf. Portuguese.
  39. 39. Ministério da Saúde do Brasil, Secretaria de Vigilância em Saúde. Leptospirose diagnóstico e manejo clínico. Brasilia, Brasil: Ministério da Saúde; 2014. Portuguese.
  40. 40. Porta M. A dictionary of epidemiology. 5th ed. New York: Oxford University Press.
  41. 41. Fundação Instituto Brasileiro de Estatística-IBGE. 2014. [cited 2015 Jan 15]. Available from: from: http://downloads.ibge.gov.br/downloads_estatisticas.htm.
  42. 42. Brumer A. Genero e agricultura: a situacao da mulher na agricultura do Rio Grande do Sul. Universidade Federal do Rio Grande do Sul. Estudos Feministas, Florianolopis, 12(1):360, 2004. http://www.scielo.br/pdf/ref/v12n1/21699.
  43. 43. Deddens JA, Petersen MR.Approaches for estimating prevalence ratios. Occ Envir Med. 2008 65: 501–506.
  44. 44. Barry SC, Welsh AH. Generalized additive modelling and zero inflated count data. Eco Mod. 2002 Nov 3; 157(2–3)179–88.
  45. 45. Vuong QH. Likelihood ratio tests for model selection and non-nested hypotheses. Econo. 1989 Mar; 57(2):307–33
  46. 46. Hu MC, Pavlicova M, Nunes EV. Zero-inflated and Hurdle Models of Count Data with Extra Zeros: Examples from an HIV-Risk Reduction Intervention Trial. Am J Drug Alcohol Abuse. 2011 37(5): 367–375. pmid:21854279
  47. 47. Empresa Brasileira de Pesquisa Agropecuária—EMBRAPA. Milho e Sorgo, cultivo do milho, sistemas de produção, 2 ISSN 1679-012X, 5 edição completar.
  48. 48. Both GC, Haetinger C, Ferreira ER, Diedrich VL, Fay de Azambuja JL. Uso da modelagem matemática para a previsão de enchentes no Vale do Taquari–RS. Anais do Simpósio Brasileiro de Engenharia Ambiental [CD-ROM]; 2008.
  49. 49. FAO GeoNetwork-geo-spatial data [Internet]. [cited 2014 Dec 10]. Available from: http://www.fao.org/geonetwork/srv/en/main.home.
  50. 50. Gordon CE, Turner LH. The effect of pH on the survival of leptospirosis in water. Bull. World Health Organ. 1961; 24(1):35–43. pmid:20604084
  51. 51. Prado H. Pedologia fácil. Enquete #16 –Fatores de formação dos solos. [Internet]. [cited 2015 Jan 24]. Available from: http://www.pedologiafacil.com.br/enquetes/enq16.php. Portuguese.
  52. 52. Crop water information: Tobacco [Internet]. Rome, Italy: Food and Agriculture Organization, Land and water division. c 2013 –[cited 2015 Jan 15]. Available from: http://www.fao.org/nr/water/cropinfo_tobacco.html
  53. 53. Smith DJW, Self HRM. Observations on the survival of Leptospira australis A in soil and water. J. Hyg. 1955; 53(4):436–44. pmid:13278528
  54. 54. Khairani-Bejo S, Bahaman AR, Zamri-Saad M, Mutalib AR. The Survival of Leptospirosis interrogans Serovar Hardjo in the Malaysian Environment. J An Vet Adv. 2004; 3(3):123–9.
  55. 55. Silveira RLL, Dornelles M. A produção de tabaco no sul do Brasil: um olhar geográfico sobre a sua dinâmica produtiva e organização territorial. In: Silveira RLL (Org.). Tabaco, sociedade e território: relações e contradições no sul do Brasil. Santa Cruz do Sul: EDUNISC, 2013. p. 17–38.
  56. 56. Companhia nacional de abastecimento-CONAB. [Internet]. Porto Alegre, Brasil. [cited 2015 Jan 11]. Available from: www.conab.gov.br
  57. 57. Caminiti Tejada R, Romani Romani F, Wong Chero P, Alarcon Villaverde J. Prácticas laborales de riesgo en cultivadores de arroz del valle del Alto Mayo, Región San Martín, Perú. Rev Peru Epidemiol [Internet]. 2011 Apr [cited 2015 Mar 12];15(1):[6 p.]. Available from: http://rpe.epiredperu.net/rpe_ediciones/2011_V15_N01/8AO_Vol15_No1_2011_Practicas_laborales_arroceros.pdf. Spanish.
  58. 58. Lelu M, Muñoz-Zanzi C, Higgins B, Galloway R. Seroepidemiology of leptospirosis in dogs from rural and slum communities of Los Rios Region, Chile. BMC Vet Res. 2015, 11:31. pmid:25880871
  59. 59. Hermann GP, Oliveira Rodrigues R, Machado G, Moreira EC, Lage A, Cerqueira Leite R. Soroprevalencia de leptospirose em bovinos nas mesorregiões sudeste e sudoeste do estado Rio Grande do Sul, Brasil. Ciência Animal Brasileira. 2012;13(1):131–138. Portuguese.
  60. 60. Astudillo VM, Dora FJ & Silva AJM. (1986). Ecosistemas y estrategias regionales de control de la fiebre aftosa. Aplicación al caso de Rio Grande do Sul, Brasil. Bol Cent Panamerican Fiebre aftosa. 1986; 52:47–61. Spanish.
  61. 61. Schneider MC, Santos-Burgoa C. Algunas consideraciones sobre la rabia humana transmitida por murcielago. Salud Publica Mex. 1995; 37 (4): 354–62.
  62. 62. Data, GDP ranking [Internet]. Washington, DC: The World Bank. c 2015 –[cited 2015 Mar 3]. Available from: http://data.worldbank.org/data-catalog/GDP-ranking-table
  63. 63. Durski KJ, Jancloes M, Chowdhary T, Bertherat E. A Global, Multi-Disciplinary, Multi-Sectorial Initiative to Combat Leptospirosis: Global Leptospirosis Environmental Action Network (GLEAN). Int J Environ Res Public Health. 2014; 11(6), 6000–6008. pmid:24905245