Skip to main content

Canine serological survey and dog culling and its relationship with human visceral leishmaniasis in an endemic urban area



Visceral leishmaniasis is an important but neglected disease that is spreading and is highly lethal when left untreated. This study sought to measure the Leishmania infantum seroprevalence in dogs, the coverage of its control activities (identification of the canine reservoir by serological survey, dog culling and insecticide spraying) and to evaluate its relationship with the occurrence of the disease in humans in the municipalities of Araçatuba and Birigui, state of São Paulo, Brazil.


Information from 2006 to 2015 was georeferenced for each municipality and modeling was performed for the two municipalities together. To do this, latent Gaussian Bayesian models with the incorporation of a spatio-temporal structure and Poisson distribution were used. The Besag-York-Mollie models were applied for random spatial effects, as also were autoregressive models of order 1 for random temporal effects. The modeling was performed using the INLA (Integrated Nested Laplace Approximations) deterministic approach, considering both the numbers of cases as well as the coverage paired year by year and lagged at one and two years.


Control activity coverage was observed to be generally low. The behavior of the temporal tendency in the human disease presented distinct patterns in the two municipalities, however, in both the tendency was to decline. The canine serological survey presented as a protective factor only in the two-year lag model.


The canine serological coverage, even at low intensity, carried out jointly with the culling of the positive dogs, suggested a decreasing effect on the occurrence of the disease in humans, whose effects would be seen two years after it was carried out.

Peer Review reports


Visceral leishmaniasis (VL) is considered a neglected tropical disease. It presents a high fatality rate, which can attain 100% if left untreated [1, 2]. In Brazil, VL control is based on the Visceral Leishmaniasis Control Program (VLCP), the main strategies of which include the diagnosis and timely treatment of human cases, environmental management and chemical control of the vector with residual insecticide spraying, and the canine serological survey and the canine culling of positive dogs in order to reduce the vectors’ sources of infection, since the domestic dog is considered the main reservoir of the disease in urban areas [3].

Despite decades of investment in these control measures in the country, in recent years we have observed the expansion of the transmission areas combined with the emergence of autochthonous cases in previously disease-free areas [4]. In addition, there are also doubts regarding the effectiveness of control measures for disease reduction [5,6,7]. Furthermore, few studies have been carried out with the objective of evaluating the effectiveness of the control measures. Those few studies have methodological problems that have been recognized in their experimental designs [8].

In the state of São Paulo, Brazil, the first autochthonous cases of VL in humans (HVL) occurred in 1999 in Araçatuba and Birigui, contiguous municipalities located in the northwestern region of the state, followed by other municipalities in the state. These two municipalities have been classified as an area of intense transmission of VL [3, 9, 10].

Because cases of HVL arise two years after the occurrence of canine cases [11], it was decided to test the hypothesis that the effects of VL control measures on the incidence of HVL can be observed after a period of 1 to 2 years. Thus, the objective of this study was to measure the coverage of two VL control activities and one VL surveillance measure in an endemic urban area and to evaluate their relationship with the occurrence of VL in humans. The VL control methods evaluated are: i) insecticide spraying for vector control; ii) culling of the seropositive dogs; and the iii) the canine serological survey.

For our knowledge, this is the first study to propose a spatio-temporal method for modeling the variability of VL and the effects of control measure covariates. Other studies related to VL used spatio-temporal modeling in a Bayesian approach, however, evaluating other effects, such as demographic, socioeconomic, climatic and environmental [12, 13].


Study area, type and period

This is an ecological study using data from 2006 to 2015 and the framework we used is the classical one to investigate spatio-temporal variability in an ecological setting [14, 15]. The study area considered was the urban region of the municipalities of Araçatuba (21° 12′ 41″ S; 50° 25′ 34″ W) and Birigui (21° 17′ 19“ S; 50° 20’ 24” W), plus a rural census tract of Birigui, as a connecting area between them (Fig. 1). The study area is located in the northwest of the state of São Paulo. We chose these municipalities as they were the first to present autochthonous HVL in the state of São Paulo and they are adjacent, thus allowing us to consider them as a single study area.

Fig. 1
figure 1

Map of the study area composed of the municipalities of Birigui and Araçatuba. a Urban areas of the municipalities of Araçatuba and Birigui and rural area of the municipality of Birigui, State of São Paulo, Brazil. b Indication of the area of study, represented by the red rectangle, located in the State of São Paulo, Brazil, and South America. Maps were created by the authors with shape files from:

Sources of information and GEOREFERENCING

To apply zoonosis control measures, the municipalities divide the urban area into sectors (denominated as “SUCEN Sectors”) that were chosen as the units of analysis (Fig. 1). Additionally, the Brazilian Institute of Geography and Statistics (IBGE) provides the population census data every 10 years and they divide the area in a different track (nominated as IBGE census track). Thus, digital maps of the area were constructed referencing both SUCEN sectors and IBGE census track, using the latter for the rural area.

The Zoonoses Control Center of the respective municipalities provided information according to the units of analysis and year of the human cases (residential address and year of occurrence), the control activities (number of dogs with seropositive tests for VL culled and number of properties sprayed with insecticide), and the surveillance measure (number of dogs evaluated in the canine serological survey).

The diagnostic tests used to identify infected dogs are the TRdPP®-Bio-Manguinhos test for screening and the ELISA test as a test to confirm the previous results. The TRdPP®-Bio-Manguinhos advertises to have a sensitivity of 92.9–100% and specificity of 87.5–91.7%, compared to parasitological test. The ELISA test presents sensitivity of 94.54% and specificity of 91.76%, compared to indirect immunofluorescence test [16]. The seropositive dogs are culled according to the laws with ethical principles mentioned in the surveillance manual [9]. The insecticide used in spraying operations is Alfacypermetrina with 40 mg of ingredient dose active per m2 and with a loading weight of 50 ml. This weight of the load was calculated for use in a standard spray pump with 10 l of capacity [9].

The human cases were geocoded using their residential address and information on the control measures was joined with the “SUCEN Sectors” map. The populations of the SUCEN sectors were estimated based on the information provided by IBGE (IBGE, 2018) and Foundation for the State System of Data Analysis [17, 18].

The dataset with information of human visceral leishmaniasis incidence, canine seroprevalence and control measures in Araçatuba/SP and Birigui/SP (period 2006–2015), is available according to the official code of municipalities of study in the ‘Additional File 1’ of Supplementary Information.

Dependent variable

Incidence rates of HVL, according to the year of study, were calculated for Araçatuba and Birigui, dividing each year’s number of cases by the respective populations. After the geocoding of HVL cases, the number of cases of HVL per year, per unit of analysis was obtained.


The following covariates were obtained by unit of analysis and represented the control measure’s coverage:

  • Canine serological survey coverage: the ratio between the number of dogs evaluated in the survey, multiplied by 100, and the estimated number of dogs that should be evaluated. The number of dogs was estimated based on information about the human population (ratio of 1 dog: 5 people) [19];

  • Canine seroprevalence for Leishmania infantum: ratio of number of seropositive dogs, multiplied by 100, and the number of dogs evaluated in a serological canine survey;

  • Canine culling coverage: the ratio between the number of dogs that were culled, multiplied by 100, concerning the number of dogs seropositive for VL;

  • Insecticide spraying coverage: the ratio between the number of properties that received the chemical control, multiplied by 100, and the number of properties expected to receive this control. This was estimated based on a block of radius of 200 m around each human case.

  • We also considered the municipality of occurrence of HVL cases (Araçatuba or Birigui) as a covariate.

Data analysis

To evaluate the spatio-temporal variation of HLV and the role of the covariates, the natural framework of analyses is that of ecological regression, as the data are provided as number of cases for each municipality and year [20]. The number of cases of HVL was assumed to follow a Poisson distribution, characterized by an offset given by the population in that area and by a rate specified as a latent Gaussian model with the incorporation of a spatio-temporal structure.

Initially, the number of HVL cases was considered to be a function of an intercept function and random spatial and temporal effects (intercept model). As for random spatial effects, they were considered a component of the Besag-York-Mollié (BYM) type, with random effects representing the spatial dependence of the number of cases of HVL as a function of the geographical location of the units of analysis. In particular, the BYM model is the sum of two random effects, one assuming that neighboring areas are correlated (local dependency) and one which assumes that all the areas have a degree of correlation (global dependency) [21]. To represent the temporal autocorrelation present in the data, an autoregressive model of order 1 of the random walk-1 type (RW1) was adopted [22]. This modeling was performed for the entire study period (2006 to 2015).

For the initial construction of covariate models, an exploratory data analysis was performed. Insecticide spraying coverage, according to the year and the unit of analysis, had a large number of values equal to or approximately zero. The canine culling coverage revealed collinearity with the canine serological survey coverage (correlation coefficient of 0.58) and, for this reason, insecticide spraying was not included in the model. The canine seroprevalence for Leishmania infantum was, moreover, not included because it was a result of the canine serological survey and not a control measure. Therefore, only the municipality and the canine serological survey coverage were considered in the covariate models.

Three models were tested, the first considering the yearly number of human cases and the canine serological survey coverage of the same year (with no lag) and the second and third models with a respective one and two-year lag (lagged models of 1 or 2 years). The first model considered the number of cases occurring in the period from 2006 to 2015, the second from 2007 to 2015, and the third from 2008 to 2015.

To obtain the number of expected cases, the global incidence rate was calculated by dividing the total number of HVL cases throughout the study period by the sum of the populations of the two municipalities of every year. With this global rate, the number of expected cases for each unit of analysis in each year was obtained by relating the expected rate to the population of the unit of analysis for the respective year. For the intercept model and the covariate model with no lag, the global incidence rate was calculated based on cases occurring between 2006 and 2015. For covariate and lagged models, the respective rates were calculated for the periods from 2007 to 2015 and from 2008 to 2015.

By considering the number of cases expected in the modeling, per year and unit of analysis, we can interpret the coefficients used, after exponentiation, as relative risks (RR) concerning the average incidence rate of the entire study period. The coverage of control activities used in the models were standardized by the respective means and standard deviation.

Modeling was performed in a Bayesian context using the INLA (Integrated Nested Laplace Approximation) methodology, which obtains deterministically the posterior distributions of the parameters of interest. This approach, concerning the computational time consumption, is more efficient than the MCMC (Markov Chain Monte Carlo) method, especially when dealing with latent Gaussian Bayesian models [23].

We used the TerraView software [24] to geocode the LV case address. The software QGIS [25] was used for spatial analysis tools and R software [26] for statistical analysis. The Figs. 1 and 3 were created by the authors, through QGIS Software Development Team (QGIS 2.18.20), available at <> and using the layers obtained from the Brazilian Institute of Geography and Statistics (IBGE), available at <>. The Fig. 2 was also created by the authors using Microsoft Excel, available at <> and Fig. 4 was created by the authors using R software, available at <>.

Fig. 2
figure 2

Human Visceral Leishmaniasis incidence rates and control measures, Araçatuba and Birigui, from 2006 to 2015. a Incidence rates of human visceral leishmaniasis according to the municipalities of Bririgui and Araçatuba, state of São Paulo, Brazil. b Canine serological survey coverage and canine culling coverage in both municipalities. c Insecticide spraying coverage and canine seroprevalence according to the municipality. Graphics were created by the authors


The HVL incidence in Araçatuba initially decreased and then stabilized in 2010 (Fig. 2a), while in Birigui, incidence continually decreased until 2009, followed by an increase in 2011, and then a subsequent decrease again. The decrease and stabilization of the incidence rates observed in Araçatuba were accompanied by canine serological survey coverage ranging from 0 to 10% (Fig. 2b). In Birigui, the increase in incidence in 2010 and 2011 was preceded by a period in which control activities were suspended (2006 to 2009).

Considering only the years in which canine serological surveys were performed, the HVL incidence curve in Araçatuba (2007 and 2015) behaves similarly to that of canine seroprevalence (Fig. 2c). In Birigui, the decrease in the incidence of HVL (2010 to 2015) presents a one year lag in canine seroprevalence. With regards to the coverage of the control measures, it was observed that the coverage of canine culling was highest, reaching values up to 70% in years in which canine serological surveys were performed (Fig. 2b). On the other hand, insecticide spraying coverage for vector control was very low in Araçatuba, the highest percentage being 15% in 2013 and 0% in Birigui (Fig. 2c). In general, control measures were performed with low coverage, though with greater intensity and consistency in Araçatuba.

The relative risks (RR) for HVL in the SUCEN sectors for the intercept model, from 2006 to 2015, revealed higher values in Birigui than in Araçatuba (Fig. 3). While 60% (12/20) of Birigui’s SUCEN sectors (excluding the rural census track) had RR between the third quintile and the maximum value of RR, in Araçatuba the percentage of sectors lying between these two values was 27.8% (10/36). The annual values of the RR for HVL for the intercept model (Fig. 4) revealed a decreasing temporal tendency, with a peak in 2011.

Fig. 3
figure 3

Relative risk posterior means for the human visceral leishmaniasis occurrence for the intercept model. Maps were created by the authors with shape files from:

Fig. 4
figure 4

Annual relative risk posterior means for the human visceral leishmaniasis occurrence for the intercept model. Graphic was created by the authors

Figure 3 shows the risk relative for the intercept model (intercept and spatial and temporal random effects) for the human visceral leishmaniasis occurrence in municipalities of Araçatuba and Birigui, state of São Paulo, Brazil, from 2006 to 2015. The graphic presented in Fig. 4 shows the annual relative risk for the intercept model (intercept and spatial and temporal random effects) for the occurrence of human visceral leishmaniasis municipalities of Araçatuba and Birigui, state of São Paulo, Brazil, from 2006 to 2015.

The municipality covariate was considered, in the three covariate models analyzed, significantly and positively associated with HVL incidence. The risk of occurrence of HVL in Birigui was between 1 and 2 times greater in Araçatuba (RR between 1.96 and 2.76). The effect of the canine serological survey coverage and subsequent canine culling of the positive dogs was an important result when they were lagged two years of the incidence of HVL. In this situation, the coverage of the canine serological survey was negatively associated with the occurrence of human cases and seen to have a protective effect on the incidence (RR = 0.76). The increase of one standard deviation in this coverage corresponded to a decrease of 34% in the risk of occurrence of HVL, Table 1.

Table 1 Relative risk posterior means for the Human Visceral Leishmaniasis occurrence for three covariate models

The table shows the relative risks and their respective 95% confidence intervals (CI 95%) for three covariate models (intercept, covariates and spatial and temporal random effects) for human visceral leishmaniasis incidence. The first one considered no delay between the disease incidence and the control measures. The second and third models considered a delay between the disease incidence and the control measures, respectively, of one and two years. These results are related to the municipalities of Araçatuba and Birigui, state of São Paulo, Brazil, from 2006 to 2015. CSSC, Canine Serological Survey Coverage.


The finding of this study goes in the same direction as the theoretical study of Costa et al. [27]. These authors suggested, using mathematical models, possible impact on canine VL seroprevalence after the long-term continuous program, targeting both asymptomatic and symptomatic dogs, for lower and moderate transmission of HVL rate scenarios. However, it is worth mentioning that this study did not evaluated the impacts on LVH. In our study, areas are classified as having intense transmission, the intervention coverage was low and, even that, the finding suggested a decreasing effect on the occurrence of the disease in humans, after two years of canine serological coverage, jointly with the culling of asymptomatic and symptomatic positives dogs.

It is worth noting that, in the descriptive part of this study, the result was pointed out that Araçatuba has a peak in the incidence of LVH in 2007, followed by a valley in 2010. After this period, the incidence presented a stabilization. In Birigui, there were two peaks, one in 2006 (about 21 cases / 100 thousand inhabitants-year) and the other in 2011(about 16 cases / 100 thousand inhabitants-year), followed also by stabilization. According to data from the Notifiable Diseases Information System [28] this stabilization continues until 2019. Thus, the incidence of human VL remains below 8 cases / 100 thousand inhabitants per year for at least ten years in Araçatuba and, at least, six years in Birigui.

In Araçatuba and Birigui, due to its territorial characteristics and the intensity of HVL transmission, the visceral leishmaniasis control program (VLCP) recommends that census surveys be carried out to identify naturally infected dogs [3]. It was, therefore, expected that the coverage of the canine serological survey would be 100%. However, it remained below 10% throughout the entire study period. The canine serological surveys only covered some of the sectors of these two municipalities, which was the reason these activities were not carried out as had been expected. This situation occurs in other municipalities elsewhere in the country and can be explained by public resistance to the VLCP, by structural difficulties from municipal administration [29] or by underfunding and lack of political support for the VLCP [8].

Results indicate that even in more populous municipalities, control measures targeting the canine reservoir in high priority areas may be effective at reducing human incidence. The effectiveness of the canine serological survey in reducing the risk for HVL, observed in this study (less than 40%), suggests that increasing the coverage of the serological survey, to the percentage recommended by the VLCP and may result in an even greater risk reduction. However, it is noteworthy that, as pointed out by other authors, given the operational challenges for the complete implementation of the control activities recommended in the VLCP, there are doubts about the feasibility of the suggested coverage levels [8, 29,30,31].

The coverage of positives dog culling, which is dependent on the canine serological survey, was high where the latter activity was performed. It is already known that when the infected dogs are eliminated the source of infection for the vector in that locality is reduced [32, 33]. Moreover, some studies suggest that canine culling has a lower cost compared to vector control [34, 35].

Some authors have found results similar to those of this study regarding the effectiveness of canine control. Costa et al. [36] demonstrated that canine culling reduced the incidence of HVL by 80% when compared to areas that only received Insecticide spraying in Teresina, PI, Brazil. Nunes et al. [37], moreover, revealed, in a two-year lag analysis, that the reduction in the incidence of HVL was statistically correlated with the increase in canine culling. Werneck et al. [38] showed that, as compared to insecticide spraying, only canine culling was effective in reducing the incidence of HVL, presenting an inverse relationship between these measures. Although those authors considered the positivity of the result to be biased due to the selective loss during follow-up, our results reinforce their findings.

Nevertheless, there is no consensus in the literature about the effectiveness of canine reservoir elimination. Dietze et al. [39] showed that canine culling did not result in a significant decrease in human incidence, pointing to the possibility of humans being reservoirs. Ashford et al. [40], although they suggested that eliminating the majority of the positive dogs could reduce the human incidence, concluded that the attempt to remove them was insufficient to guarantee the elimination of the disease. Thus, they suggested that other reservoirs could be acting in the transmission, while also questioning the serological tests for the identification of the seropositive dogs. In addition, some authors point to problems related to canine cullings such as the lack of precision of the serological tests, the long delay between the identification and the canine culling, the refusal of the owners to cooperate and the replacement of the dogs that were culled [31, 38, 41, 42].

In addition to the technical issues, the most controversial aspect of canine reservoir elimination has been the ethical question considering the emotional involvement of both the owners and the subjects involved in the culling procedure [33,34,35,36,37,38,39,40,41,42,43,44,45,46]. However, it is noteworthy that HVL is a disease with high lethality if untreated, and it is in the process of geographical expansion [3, 47]. Thus, to stop canine reservoir control measures, it is first necessary that there is an efffective alternativepublic health programs.

Among these alternatives are the canine treatment, canine vaccination and the use of the insecticide-impregnated collar, available so far through the initiative of the owners. To date, there is canine treatment that was legally approved by the Ministries of Agriculture and Health in 2016, highlighting the fact that parasitological cure has not yet been proven. This treatment, in addition to being costly, requires the use of the insecticide-impregnated collar and repellent against the vector, as well as a follow-up of the canine patient by a veterinarian every 4 months, which makes it difficult to apply to the population in general [48, 49]. However, there are published recommendations of preserving the right to choose between responsible treatment and euthanasia for owners of infected dogs [50].

Dog vaccination, whose efficacy ranges from 58.1 to 80.8%, has been considered a good alternative [51], however, there is a lack of randomized studies that might demonstrate its efficacy for general use as a measure of VL prevention [42, 52]. In addition to the cost, one of the difficulties for large-scale dog vaccination is the need for three initial doses, with annual boosters, a strategy that is impractical for public health programs.

The use of insecticide-impregnated collars has been recommended as the most effective measure among those cited above [43, 53,54,55]. Silva et al. [56] demonstrated a reduction of 15% (p = 0.004) in Montes Claros, MG, Brazil, and 60% (p < 0.001) in Fortaleza, CE, Brazil, in the vector population of captured. Kazimoto et al. [57] observed that there was a 53–59% reduction in the VL incidence in dogs resulting from the use of these collars. However, this measure requires continuous use in a large percentage of the canine population. Also, the collars must be changed every six months and replaced when lost. As a result, the authors warn of their high cost when used on a large scale [56, 58, 59].

One of the limitations of our study was the fact that insecticide spraying had zero or very low coverage in the municipalities studied, making it impossible to evaluate its effectiveness in reducing HVL incidence. This low coverage is indicative of the operational difficulties faced in the implementation, including the residents’ refusal to apply insecticide in their homes, deficiencies in the quality control of the handling of the insecticide, a brief residual effect of the insecticide spraying and a lack of human and financial resources [7, 36, 60]. Another limitation is the fact that both dogs and humans move and there is no way to be sure that transmission of Leishmania occurred in the SUCEN sector where cases were assigned.

Another issue that must be considered is the lack of accuracy of available serological tests, especially for asymptomatic dogs. Grimaldi et al. [61] highlight the TRdPP®-Bio-Manguinhos presents high sensitivity (98%) for diagnosed dogs with symptoms, but low sensitivity to identify dogs without signs and symptoms (47%). Moreover, these authors discuss that the ELISA test, despite having reliable sensitivity (ranging from 93 to 100%) for symptomatic dogs, also presents less sensitivity for dogs without symptoms. It is noteworthy that, even using tests with low sensitivity for asymptomatic dogs, our results suggested an inverse relationship between HVL incidence and canine survey.

Moreover, we were unable to consider controlling the model for HIV/AIDS incidence during the study period. HIV/ AIDS coinfection is known associated with higher risk of VL human symptomatic infection and has important epidemiological and clinical implications. In endemic areas, it is common for cases of VL to be asymptomatic. However, HIV coinfection increases the risk of developing symptomatic VL by between 100 and 2320 times [62]. Since it was not possible to adjust our models for the incidence of HIV / AIDS, it is not possible to state what effect this variable would have on the relationship between the incidence of HVL and the canine serological survey coverage, an issue that needs to be investigated in new studies.

Another limitation is the use of secondary data. Nevertheless, these data were indispensable and appropriate for the development of the study. It is worth noting that it is the health teams which make the decisions about the planning of the surveillance actions and control of VL in their respective municipalities based on this information. This information is essential for all involved in VL control efforts.

Strengths of the study include the consideration of the spatial and temporal dependence of the phenomenon studied and, consequently, more accurate estimates were obtained than in studies that did not take these dimensions into account. The proximity of Araçatuba to Birigui is also worth mentioning as it made it possible to evaluate the relationship between HVL and control measures more comprehensively and with greater study power, given the greater number of units of analysis.


Our study revealed that the canine serological survey in the municipalities of Araçatuba and Birigui, undertaken concomitantly with the canine culling of positive dogs, although with incomplete coverage, suggests a reduction in the risk for the occurrence of HVL of 34%, with the effects observed after two years of the implementation of these control measures. Therefore, the study suggest that, even when below the recommended levels, the control measures directed at the canine reservoir were effective for the control of HVL.

Thus, our study showed the significant role played by the domestic dog in the occurrence of human disease. It suggests the use of measures for surveying and eliminating these reservoirs should continue until new control measures are proven effective as a large scale alternative reduce the risk of visceral leishmaniasis in humans. Others studies with different methodological approaches must be carried out to affirm the efficacy of this association.

Availability of data and materials

All relevant data are within the Additional File 1.



Confidence intervals


Canine Serological Survey Coverage


Human Visceral Leishmaniasis


Brazilian Institute of Geography and Statistics




Integrated Nested Laplace Approximation


Markov Chain Monte Carlo


Relative risks


Random walk-1 type


Visceral leishmaniasis


Visceral leishmaniasis control program


  1. CDC Centers for Disease Control and Prevention : Parasites – Leishmaniasis. (c2013). Accessed 01 jun 2019.

  2. WHO World Health Organization: Visceral Leishmaniasis. (c2018). Acessed 2018 Jul 2018.

  3. Ministry of Health (BR), Secretariat of Health Surveillance. Health Surveillance Guide. Brasília (DF). (2016). Accessed 27 Jul 2018. Portuguese.

  4. SBMT Brazilian Society of Tropical Medicine: the second autochthonous case of calazar is confirmed in Florianópolis (c2017) Accessed 10 Jun 2019. Portuguese.

  5. Werneck GL. Geographic expansion of visceral leishmaniasis in Brazil. Cad Saude Publica. 2010;26:644–5 Portuguese.

    Article  Google Scholar 

  6. Werneck GL. Visceral leishmaniasis in Brazil: rationale and concerns related to reservoir control. Rev Saude Publica. 2014;48:851–6.

    Article  Google Scholar 

  7. Salomón OD, Feliciangeli MD, Quintana MG, Afonso MM, Rangel EF. Lutzomyia longipalpis urbanisation and control. Mem Inst Oswaldo Cruz. 2015;110:831–46.

    Article  Google Scholar 

  8. Romero GA, Boelaert M. Control of visceral leishmaniasis in latin america-a systematic review. PLoS Negl Trop Dis. 2010;4:e584.

    Article  Google Scholar 

  9. Camargo-Neves VLF, Glasser CM, Cruz LL, Almeida RG. Manual of Surveillance and Control of American Visceral Leishmaniasis of the State of São Paulo. (2006). Accessed 27 jan 2019. Portuguese.

  10. SINAN Notification of injury information system: Birigui and Araçatuba (SP). Confirmed cases. (c2013-2015). Accessed 10 fev 2019. Portuguese.

  11. Rosales JC, Yang HM. Mathematical model to describe the transmission of leishmaniasis. TEMA tend. Mat. Apl. Comput 2006; Portuguese.

  12. Karagiannis-Voules DA, Scholte RG, Guimaraes LH, Utzinger J, Vounatsou P. Bayesian geostatistical modeling of leishmaniasis incidence in Brazil. PLoS Negl Trop Dis. 2013;7:e2213.

    Article  Google Scholar 

  13. Godana AA, Mwalili SM, Orwa GO. Dynamic spatiotemporal modeling of the infected rate of visceral leishmaniasis in human in an endemic area of Amhara regional state, Ethiopia. PLoS One. 2019;14:e0212934.

    Article  CAS  Google Scholar 

  14. Spatial and spatio-temporal models with R-INLA, Wiley; Elliott P, Wakefield J, Best N, Briggs D. Spatial Epidemiology. Oxford University Press, Oxford, 2000.

  15. Blangiardo M, Cameletti M. Spatial and Spatio-temporal Bayesian models with R-INLA. New York: John Wiley and Sons Inc.; 2015.

    Book  Google Scholar 

  16. Ribeiro VM, Miranda JB, Marcelino AP, De Andrade HM. Performance of different serological tests in the diagnosis of natural infection by Leishmania infantum in dogs. Vet Parasitol. 2018;274:108920.

    Article  Google Scholar 

  17. RIPSA Interagency information network for health: Birigui and Araçatuba (SP). Resident population. (c2006-2009) accessed 10 Jun 2019. Portuguese.

  18. SEADE Foundation state system of data analysis: Birigui and Araçatuba (SP) Seade system of population projections. (c2011-2015). Accessed 25 Jun 2019. Portuguese.

  19. Andrade AM, Queiroz LH, SHV P, Nunes CM. Descriptive study of the canine population structure of the urban area of Araçatuba, São Paulo, Brazil, from 1994 to 2004. Cad Saude Publica. 2008;24:927–32 Portuguese.

    Article  Google Scholar 

  20. Elliott P, Wakefield J, Best N, Briggs D, eds. Spatial Epidemiology: Methods and Applications. Oxford University Press; 2001. Accessed February 19, 2020.

  21. Besag J, York J, Mollie A. Bayesian image restoration, with two applications in spatial statistics. Ann Inst Stat Math. 1991;43:1–59.

    Article  Google Scholar 

  22. Blangiardo M, Cameletti M, Baio G, Rue H. Spatial and spatio-temporal models with R-INLA. Spat Spatiotemporal Epidemiol. 2013;7:39–55.

    Article  Google Scholar 

  23. Rue H, Martino S, Chopin N. Approximate Bayesian inference for latent Gaussian models by using integrated nested Laplace approximations. J R Stat Soc Ser B Stat Methodol. 2009;71:319–92.

    Article  Google Scholar 

  24. INPE National Institute for Space Research. TerraView. Version 4.2.2 [software]. Accessed 15 Feb 2013.

  25. QGIS Development Team. QGIS Geographic Information System. Open Source Geospatial Foundation Project. Version 2.18 [software]. Accessed 23 jun 2017.

  26. R Development Core Team. R Foundation for Statistical Computing. Version 3.2.3 [software]. Accessed 10 dec 2018.

  27. Costa DN, Codeço CT, Silva MA, Werneck GL. Culling dogs in scenarios of imperfect control: realistic impact on the prevalence of canine visceral leishmaniasis. PLoS Negl Trop Dis. 2013;7:e2355.

    Article  Google Scholar 

  28. SINAN. Notifiable Diseases Information System. Confirmed cases of Visceral Leishmaniasis according to LPI and year of notification, State of São Paulo, 2014 to 2019. Accessed 10 Feb 2020.

  29. Zuben APB, Donalisio MR. Difficulties in implementing the guidelines of the Brazilian visceral Leishmaniasis control program in large cities. Cad Saude Publica. 2016.

  30. Profeta da Luz ZM, Pimenta DN, Cabral AL, Fiúza VO, Rabello A. Leishmaniasis urbanization and low diagnosis capacity in the Metropolitan Region of Belo Horizonte. Rev Sob Bras Med Trop. 2001. Portuguese.

  31. Nunes CM, Lima VM, Paula HB, Perri SH, Andrade AM, Dias FE, et al. Dog culling and replacement in an area endemic for visceral leishmaniasis in Brazil. Vet Parasitol. 2008;153:19–23.

    Article  Google Scholar 

  32. Palatnik-de-Sousa CB, dos Santos WR, França-Silva JC, da Costa RT, Reis AB, Palatnik M, et al. Impact of canine control on the epidemiology of canine and human visceral leishmaniasis in Brazil. Am J Trop Med Hyg. 2001;65:510–7.

    Article  CAS  Google Scholar 

  33. Palatnik-de-Sousa CB, Day MJ. One health: the global challenge of epidemic and endemic leishmaniasis. Parasit Vectors. 2011;4:197.

    Article  Google Scholar 

  34. Neves VLFC. Epidemiologic aspects and evaluation of the control methods American visceral leishmaniasis in São Paulo state, Brazil [thesis] (2004). Accessed 18 Jul 2019. Portuguese.

  35. Shimozako HJ. Optimization model for the control of leishmaniasis: epidemiological and economical analysis. [thesis]. 2015; doi: Portuguese.

  36. Costa CH, Tapety CM, Werneck GL. Control of visceral leishmaniasis in urban areas: randomized factorial intervention trial. Rev Soc Bras Med Trop. 2007;40:415–9.

    Article  Google Scholar 

  37. Nunes CM, Pires MM, da Silva KM, Assis FD, Gonçalves Filho J, Perri SH. Relationship between dog culling and incidence of human visceral leishmaniasis in an endemic area. Vet Parasitol. 2010;170:131–3.

    Article  Google Scholar 

  38. Werneck GL, Costa CH, de Carvalho FA, Pires e Cruz MS, Maguire JH, Castro MC. Effectiveness of insecticide spraying and culling of dogs on the incidence of Leishmania infantum infection in humans: a cluster randomized trial in Teresina, Brazil. PLoS Negl Trop Dis. 2014;8:e3172.

    Article  Google Scholar 

  39. Dietze R, Barros GB, Teixeira L, Harris J, Michelson K, Falqueto A, et al. Effect of eliminating seropositive canines on the transmission of visceral leishmaniasis in Brazil. Clin Infect Dis. 1997;25:1240–2.

    Article  CAS  Google Scholar 

  40. Ashford DA, David JR, Freire M, David R, Sherlock I, Eulalio MC, et al. Studies on control of visceral leishmaniasis: impact of dog control on canine and human visceral leishmaniasis in Jacobina, Bahia, Brazil. Am J Trop Med Hyg. 1998;59:53–7.

    Article  CAS  Google Scholar 

  41. Costa CH. How effective is dog culling in controlling zoonotic visceral leishmaniasis? A critical evaluation of the science, politics and ethics behind this public health policy. Rev Soc Bras Med Trop. 2011;44:232–42.

    Article  Google Scholar 

  42. Otranto D, Dantas-Torres F. The prevention of canine leishmaniasis and its impact on public health. Trends Parasitol. 2013;29:339–45.

    Article  Google Scholar 

  43. Sevá AP, Ovallos FG, Amaku M, Carrillo E, Moreno J, Galati EA, et al. Canine-based strategies for prevention and control of visceral Leishmaniasis in Brazil. PLoS One. 2016;11:e0160058.

    Article  Google Scholar 

  44. Shimozako HJ, Wu J, Massad E. Mathematical modelling for zoonotic visceral Leishmaniasis dynamics: a new analysis considering updated parameters and notified human Brazilian data. Infect Dis Model. 2017;2:143–60.

    PubMed  PubMed Central  Google Scholar 

  45. Machado CJS, Silva EG, Vilani RM. The use of a controversial public health policy instrument: the euthanasia of dogs contaminated by leishmaniasis in Brazil. Saúde Soc. 2016;25:247–58 Portuguese.

    Google Scholar 

  46. Savio TPS, Marques LFV, Lamounier KCC, Castro JMD, Borja-Cabrera GP. Human visceral leishmaniasis: ethical and legal reflections on the canine reservoir control in Brazil. Rev bioét derecho. 2017. Portuguese.

  47. Sevá AD, Mao L, Galvis-Ovallos F, Tucker Lima JM, Valle D. Risk analysis and prediction of visceral leishmaniasis dispersion in São Paulo state, Brazil. PLoS Negl Trop Dis. 2017;11:e0005353.

    Article  Google Scholar 

  48. Ministério da Agricultura, Pecuária e Abastecimento (BR). Technical Note n° 001/2016 mapa/ms. Provides authorization for product registration Milteforan owned by the company vibrac animal health, indicated for the treatment of visceral leishmaniasis of dogs. (2016). Accessed 01 set 2018. Portuguese.

  49. Vibrac animal health :Milteforan™. (c2018). Accessed 20 jun 2019.

  50. Ribeiro VM, Da Silva SM, Menz I, Tabanez P, Nogueira FDS, Werkhaüser M, et al. Control of visceral leishmaniasis in Brazil: Recommendations from Brasileish. Parasites Vectors. 2013;6:8.

  51. Regina-Silva S, Feres AM, França-Silva JC, Dias ES, Michalsky É, de Andrade HM, et al. Field randomized trial to evaluate the efficacy of the Leish-Tec® vaccine against canine visceral leishmaniasis in an endemic area of Brazil. Vaccine. 2016;34:2233–9.

    Article  CAS  Google Scholar 

  52. Wylie CE, Carbonell-Antoñanzas M, Aiassa E, Dhollander S, Zagmutt FJ, Brodbelt DC, et al. A systematic review of the efficacy of prophylactic control measures for naturally-occurring canine leishmaniosis, part I: vaccinations. Prev Vet Med. 2014;117:7–18.

    Article  CAS  Google Scholar 

  53. Brianti E, Napoli E, Gaglio G, Falsone L, Giannetto S, Solari Basano F, et al. Field evaluation of two different treatment approaches and their ability to control fleas and prevent canine Leishmaniosis in a highly endemic area. PLoS Negl Trop Dis. 2016;10:e0004987.

    Article  Google Scholar 

  54. Lopes EG, Sevá AP, Ferreira F, Nunes CM, Keid LB, Hiramoto RM, et al. Vaccine effectiveness and use of collar impregnated with insecticide for reducing incidence of Leishmania infection in dogs in an endemic region for visceral leishmaniasis, in Brazil. Epidemiol Infect. 2018;146:401–6.

    Article  CAS  Google Scholar 

  55. David JR, Stamm LM, Bezerra HS, Souza RN, Killick-Kendrick R, Lima JW. Deltamethrin-impregnated dog collars have a potent anti-feeding and insecticidal effect on Lutzomyia longipalpis and Lutzomyia migonei. Mem Inst Oswaldo Cruz. 2001;96:839–47.

    Article  CAS  Google Scholar 

  56. Silva RAE, Andrade AJ, Quint BB, Raffoul GES, Werneck GL, Rangel EF, et al. Effectiveness of dog collars impregnated with 4% deltamethrin in controlling visceral leishmaniasis in Lutzomyia longipalpis (Diptera: Psychodidade: Phlebotominae) populations. Mem Inst Oswaldo Cruz. 2018;113:e170377.

    PubMed  PubMed Central  Google Scholar 

  57. Kazimoto TA, Amora SSA, Figueiredo FB, Magalhães JME, Freitas YBN, Sousa MLR, et al. Impact of 4% Deltamethrin-impregnated dog collars on the prevalence and incidence of canine visceral Leishmaniasis. Vector Borne Zoonotic Dis. 2018;18:356–63.

    Article  Google Scholar 

  58. Maroli M, Mizzon V, Siragusa C, D'Oorazi A, Gradoni L. Evidence for an impact on the incidence of canine leishmaniasis by the mass use of deltamethrin-impregnated dog collars in southern Italy. Med Vet Entomol. 2001;15:358–63.

    Article  CAS  Google Scholar 

  59. Gomez SA, Picado A. Systemic insecticides used in dogs: potential candidates for phlebotomine vector control? Tropical Med Int Health. 2017;22:755–64.

    Article  Google Scholar 

  60. Romero GAS. The control of visceral leishmaniasis in Brazil: transforming is necessary. Cad Saude Publica. 2016;32:eCO010616–eCO Portuguese.

    Article  Google Scholar 

  61. Grimaldi G Jr, Teva A, Ferreira AL, dos Santos CB, Pinto I, de-Azevedo CT, Falqueto A. Evaluation of a novel chromatographic immunoassay based on Dual-Path Platform technology (DPP(R) CVL rapid test) for the serodiagnosis of canine visceral leishmaniasis. Trans R Soc Trop Med Hyg. 2012;106:54–9.

    Article  CAS  Google Scholar 

  62. WHO. Leishmaniasis and HIV coinfection. Accessed 10 Feb 2020.

Download references


Not applicable.


This study was supported by: The São Paulo Research Foundation ( PMMB, 2016/25878–7; FCN, 2014/17764–6; DNCCC, 2014/ 06518–4; and FCN, 2017/09671–6. And by The National Council for Scientific and Technological Development ( CMN and FCN, Research Fellows. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Author information

Authors and Affiliations



Study conception and design: PMMB, DNCCC, CMN, JET, RMH, LACR, RSC, MB, FCN. Acquisition of data: DNCCC. Data Curation: PMMB, FCN. Supervision: FCN. Writing Original Draft: PMMB, DNCCC, CMN, JET, RMH, FCN. Formal Analysis: PMMB, MB, FCN. All authors made the critical reading, contributed to revision of the manuscript and approved the final version of the manuscript.

Corresponding author

Correspondence to Patricia Marques Moralejo Bermudi.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the Ethics in Research Committee - COEP - of the School of Public Health of the USP (CAAE: 38170514.4.0000.5421, December 11, 2014). No informed consent was used because we only worked with secondary data provided by the surveillance system. All data analyzed were anonymized.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Table S1:

Dataset with information of human visceral leishmaniasis incidence and control measures. The dataset contains information on the incidence rates of human visceral leishmaniasis, control measure coverages, and canine seroprevalence for Leishmania infantum, according official code of municipalities of study, Araçatuba/SP and Birigui/SP and period 2006–2015. The data sources were the Zoonoses Control Center of Araçatuba/SP and Birigui/SP.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Bermudi, P.M.M., Costa, D.N.C.C., Nunes, C.M. et al. Canine serological survey and dog culling and its relationship with human visceral leishmaniasis in an endemic urban area. BMC Infect Dis 20, 401 (2020).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Visceral leishmaniasis
  • Ecological study
  • Control measures
  • Brazil