Exploring Spatial and Temporal Distribution of Cutaneous Leishmaniasis in the Americas, 2001-2011

Abstract

Leishmaniasis is an important health problem in several countries in the Americas and cases notification is limited and underreported. In 2008, the Pan American Health Organization (PAHO/WHO) met with endemic countries to discuss the status and need of improvement of systems region-wide. The objective is to describe the temporal and spatial distribution of cutaneous leishmaniasis (CL) cases reported to PAHO/WHO by the endemic countries between 2001 and 2011 in the Americas.

Methods

Cases reported in the period of 2001-2011 from 14/18 CL endemic countries were included in this study by using two spreadsheet to collect the data. Two indicators were analyzed: CL cases and incidence rate. The local regression method was used to analyze case trends and incidence rates for all the studied period, and for 2011 the spatial distribution of each indicator was analyzed by quartile and stratified into four groups.

Results

From 2001-2011, 636,683 CL cases were reported by 14 countries and with an increase of 30% of the reported cases. The average incidence rate in the Americas was 15.89/100,000 inhabitants. In 2011, 15 countries reported cases in 180 from a total of 292 units of first subnational level. The global incidence rate for all countries was 17.42 cases per 100,000 inhabitants; while in 180 administrative units at the first subnational level, the average incidence rate was 57.52/100,000 inhabitants. Nicaragua and Panama had the highest incidence but more cases occurred in Brazil and Colombia. Spatial distribution was heterogeneous for each indicator, and when analyzed in different administrative level. The results showed different distribution patterns, illustrating the limitation of the use of individual indicators and the need to classify higher-risk areas in order to prioritize the actions. This study shows the epidemiological patterns using secondary data and the importance of using multiple indicators to define and characterize smaller territorial units for surveillance and control of leishmaniasis.

About the Authors:

Ana Nilce Silveira Maia-Elkhoury

* E-mail: [email protected]

Affiliation: Communicable Diseases and Health Analysis, Pan American Health Organization, Duque de Caxias, Rio de Janeiro, Brazil

ORCID http://orcid.org/0000-0001-5369-5860

Zaida E. Yadón

Affiliation: Communicable Diseases and Health Analysis, Pan American Health Organization, Washington, District of Columbia, United States of America

Martha Idali Saboyá Díaz

Affiliation: Communicable Diseases and Health Analysis, Pan American Health Organization, Washington, District of Columbia, United States of America

Francisca de Fátima de Araújo Lucena

Affiliation: Ministry of Social Development and Against Hungry, Brasilia, Federal District, Brazil

Luis Gerardo Castellanos

Affiliation: Communicable Diseases and Health Analysis, Pan American Health Organization, Washington, District of Columbia, United States of America

Manuel J. Sanchez-Vazquez

Affiliation: Pan American Foot-and-Mouth Disease Center-PANAFTOSA/PAHO, Duque de Caxias, Rio de Janeiro, BrazilAbstract

Leishmaniasis is an important health problem in several countries in the Americas and cases notification is limited and underreported. In 2008, the Pan American Health Organization (PAHO/WHO) met with endemic countries to discuss the status and need of improvement of systems region-wide. The objective is to describe the temporal and spatial distribution of cutaneous leishmaniasis (CL) cases reported to PAHO/WHO by the endemic countries between 2001 and 2011 in the Americas.

Methods

Cases reported in the period of 2001-2011 from 14/18 CL endemic countries were included in this study by using two spreadsheet to collect the data. Two indicators were analyzed: CL cases and incidence rate. The local regression method was used to analyze case trends and incidence rates for all the studied period, and for 2011 the spatial distribution of each indicator was analyzed by quartile and stratified into four groups.

Results

From 2001-2011, 636,683 CL cases were reported by 14 countries and with an increase of 30% of the reported cases. The average incidence rate in the Americas was 15.89/100,000 inhabitants. In 2011, 15 countries reported cases in 180 from a total of 292 units of first subnational level. The global incidence rate for all countries was 17.42 cases per 100,000 inhabitants; while in 180 administrative units at the first subnational level, the average incidence rate was 57.52/100,000 inhabitants. Nicaragua and Panama had the highest incidence but more cases occurred in Brazil and Colombia. Spatial distribution was heterogeneous for each indicator, and when analyzed in different administrative level. The results showed different distribution patterns, illustrating the limitation of the use of individual indicators and the need to classify higher-risk areas in order to prioritize the actions. This study shows the epidemiological patterns using secondary data and the importance of using multiple indicators to define and characterize smaller territorial units for surveillance and control of leishmaniasis.

Author Summary

In the Americas, cutaneous leishmaniasis (CL) cases are notified across a wide geographic area, extending from southern United States to northern Argentina. Currently, 70-75% of all estimated cases of CL worldwide occur in ten countries, including four in the Americas (Brazil, Colombia, Nicaragua, and Peru) and six countries from Africa and the Eastern Mediterranean. This study shows the epidemiological situation and geographical distribution of CL cases reported by passive surveillance system of endemic countries of the American Region from 2001 to 2011, collected and consolidated by the Pan American Health Organization (PAHO). Furthermore, it represents a joint effort of the National Programs of leishmaniases. Despite some limitations, the data were analyzed and discussed and this study represents the first step to understand the global epidemiological situation in this Region and it contributes to the improvement of the surveillance of leishmaniases in the Americas.

Citation: Maia-Elkhoury ANS, E. Yadón Z, Idali Saboyá Díaz M, de Fátima de Araújo Lucena F, Gerardo Castellanos L, J. Sanchez-Vazquez M (2016) Exploring Spatial and Temporal Distribution of Cutaneous Leishmaniasis in the Americas, 2001-2011. PLoS Negl Trop Dis 10(11): e0005086. doi:10.1371/journal.pntd.0005086

Editor: Anne-Laure Bañuls, Institut de Recherche pour le Développement, FRANCE

Received: May 9, 2016; Accepted: September 30, 2016; Published: November 8, 2016

Copyright: © 2016 Maia-Elkhoury et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Data Availability: All relevant data are within the paper.

Funding: The authors received no specific funding for this work.

Competing interests: The authors have declared that no competing interests exist.

Introduction

Leishmaniasis is an important health problem for several countries in the Americas. It is caused by different species of protozoans of the genus Leishmania and transmitted to humans and animals by insect vectors of the Psycodidae family infected during blood feeding on vertebrate reservoirs and hosts [1-3]. Various factors including climate, economic, environmental, social, and others, have an influence in the leishmaniasis transmission cycle and thus determining the epidemiological patterns of the disease [3-5].

The disease occurs in different clinical forms (cutaneous, mucosal, and visceral), according to the parasite. Aspects such as the immunogenetic profile of the affected population, malnutrition, and access to health systems, among others, have an influence in the prognosis of the disease, which can progress from mild to severe, causing deformities, mutilation, and even death [3, 5].

A recent study shows that leishmaniasis had an increase of 8.8% (-6.4 to 25.3%) in the disability-adjusted life year (DALYs) for all ages when comparing 2013 to 2005, and a significantly higher increment for the cutaneous and mucocutaneous forms (35.9%-23.7 to 49.0%) [6].

Globally, Cutaneous leishmaniasis (CL) is broadly distributed occurring in 83 countries, with an estimated of 0.7-1.2 million new cases per year [7]. Currently, 70-75% of all estimated cases of CL worldwide occur in ten countries: four in the Americas (Brazil, Colombia, Nicaragua, and Peru) and six in Africa and Eastern Mediterranean (Afghanistan, Algeria, Iran, Syria, Ethiopia, and Sudan) [3, 7]. In the Americas, CL cases have been notified across a wide geographic area extending from southern United States to northern Argentina [3].

In the past decades it has become apparent that CL is more prevalent in many areas of the Latin America than previously thought. This uncertainty about its occurrence is partially because the notification in many endemic countries was often underreported and the collected information remained limited [7, 8]. A comprehensive understanding of the transmission of leishmaniasis, incorporating the knowledge of specific spatial and temporal distribution of the disease, is essential to make decisions in order to direct and prioritize actions for surveillance and control strategies for cutaneous leishmaniasis [9]. At this juncture, the Pan American Health Organization/World Health Organization (PAHO/WHO) Regional Leishmaniasis Program, brought together (in Colombia 2008), the coordinators of the National Leishmaniasis Control Programs of the American Region to discuss the status of the leishmaniasis reporting system and agreed upon the development of a regional information and surveillance systems region-wide coordinated by PAHO/WHO. The meeting helped jumpstart the debate and definition of criteria for an improved epidemiological surveillance system to fully understand the epidemiological status of leishmaniasis in the Americas.

The objective of this study is to describe the temporal and spatial distribution of CL cases between 2001 and 2011 utilizing the data reported to PAHO/WHO by the endemic countries of the Americas, as a first step to determine the time and spatial distribution of the disease in the subregions and region.

Materials and Methods

Data collection

Over 2011 and 2012, the Regional Leishmaniasis Program of PAHO/WHO coordinated the data collection among endemic countries, requesting the annual data of CL from 2001 to 2011 through two spreadsheet templates. One spreadsheet to collect national data standardized by country, year, and clinical form of the disease for the whole study period; and a second spreadsheet to collect only data from 2011 disaggregated at the first subnational administrative level (department, state, or province, depending on the country). Additionally, consolidated, aggregated, and analyzed all the information provided by the endemic countries using the regional database [10].

There were many differences between countries in their leishmaniasis surveillance systems, including the process of reporting. Thus, while case notification for leishmaniasis is mandatory region-wide, reporting can be done on individual or aggregate basis, depending on the specific criteria established by the national epidemiological surveillance system. The definition of cases, notification process, and variables for each data collected were also different among the countries due to the different clinical forms, internal system and information flow established in each country. This scenario of uneven data structure defined the limitations of the potential analyses to be carried out in this study.

The data collection included cases reported from 2001-2011 from 14 of the 18 endemic American countries (Argentina, Brazil, Bolivia, Colombia, Costa Rica, El Salvador, Ecuador, Guatemala, Guyana, Honduras, Nicaragua, Panama, Paraguay, and Peru). French Guiana, Mexico, Suriname, and Venezuela, did not complete the spreadsheets templates provided by PAHO/WHO to report data in the time requested. For the 2011 spatial analysis, 15 countries participated (the 14 countries listed above plus Mexico) and provided information on the following variables: total number of CL cases, total population of the administrative unit (national or first subnational administrative department, state, or province) where the cases occurred.

Data analysis

CL case and incidence rate.

In order to analyze and present basic health indicators [11], PAHO/WHO divides the American Region into six subregions: 1) Andean area (Bolivia, Colombia, Ecuador, and Peru); 2) Central America (Costa Rica, El Salvador, Guatemala, Honduras, Nicaragua, and Panama); 3) Brazil; 4) Southern Cone (Argentina and Paraguay); 5) the non-Latin Caribbean (Guyana); and 6) Mexico. The same subregions were used for data analysis in this study.

Fourteen countries were analyzed for the period of 2001-2001 for two indicators: CL cases and CL incidence. The incidence rate was calculated by dividing the number of confirmed cases by the country population at each year-end and multiplied by 100,000. The population data were obtained from the official census of each country. In 2011, the incidence rate was calculated at the national and first subnational administrative level using the method described above.

Microsoft Office Excel 2007 (Microsoft Corporation, Redmond, Washington, United States) was used to compile and analyze the data.

Temporal distribution of CL cases and incidence rates.

The temporal distribution was evaluated at three levels of geographical units: i.e. region, subregion and country. The local regression method, known by the acronym LOESS (locally weighted scatterplot smoothing) [12, 13], was used to identify the trend line for the number of cases and incidence. LOESS is an exploratory statistical method that allows a least-squares regression curve to be represented through a series of data points which defines two parameters: 1) the value of the smoothing parameter that defines the subset of data used (either the nearest or the furthest), expressed as "[alpha] " and 2) the degree of the local polynomial. To adjust the curve, several combinations of parameters were attempted, minimizing an excessive weight in the smoothing (to avoid over-adjustment). Based on the results of these attempts, an [alpha] value of 0.66 was utilized with a linear weighting (a first-degree local polynomial). Graphs depicting the data for the subregions, were represented with data points and the fitted trend line. The analyses were conducted with the R statistical software environment [14], specifically the "stats" library.

Spatial distribution of CL cases and incidence rates for 2011.

For the spatial analysis, the number of cases and incidence rates for CL reported in 2011 were linked to the national and first subnational administrative data represented by polygons. Subsequently, the data were categorized in quartiles, which allowed distributing the values of each indicator into four groups containing equal number of units. These categories were utilized to develop choropleth maps (scaling the color according to the value) with the distribution of the two CL indicators, both at the national and first subnational administrative levels across the Region.

The maps were generated with the ArcGIS program (ESRI, Redlands, California, United States). The administrative boundaries at the country level were obtained from the Global Administrative Areas Project (http://gadm.org/) and the boundaries of the first subnational administrative level were obtained from the Second Administrative Level Boundaries data set (SALB).

Ethical Considerations

In this study the legal principles and ethical aspects were considered and respected. These data refer to the routine of a control program of leishmaniasis and surveillance services of the Health Ministries of endemic countries, which maintained confidential the identification of patients.

Results

CL cases: 2001-2011 period

In the period of 2001-2011, 14 countries reported a total of 636,683 CL cases, with an annual average of 57,923 cases (ranging from 47,286 to 67,949). A total of 270,572 of the reported cases (42.5%) were from Brazil; 256,261 (40.2%) were from countries of the Andean subregion; 100,475 (15.8%) were from Central American countries; and 9,375 (1.5%) were reported by Argentina, Guyana, and Paraguay, Table 1.

Download:

*

PPT

PowerPoint slide

*

PNG

larger image

*

TIFF

original image

Table 1. Number of cases and incidence (Incid) rate of cutaneous leishmaniasis per 100.000 inhabitants, reported per year for the 14 selected endemic countries and by subregion, Americas region, 2001-2011.

Source: PAHO/WHO. Data reported by the leishmaniasis control programs of the countries to PAHO/WHO through May 2012.

http://dx.doi.org/10.1371/journal.pntd.0005086.t001

In the first half of the studied period (2001-2005) there was a continuous increase of the total number of CL cases reported per year, which increased 44% in 2005 when compared to 2001. In the second half of the period (2006-2011) there was a decline of 11% of the CL in this region, but resulting in an overall increase of 30% for all study period, Table 1, Fig 1.

Download:

*

PPT

PowerPoint slide

*

PNG

larger image

*

TIFF

original image

Fig 1. Number of cutaneous leishmaniasis cases reported by 14 endemic countries (overall), and stratified by sub region, Americas region, 2001-2011.

The Plotted points correspond to actual counts of cases. The Lines consist in short term trends (fitted with LOESS). Source: PAHO/WHO. Data reported by the leishmaniasis control programs of the countries to PAHO/WHO through May 2012.

http://dx.doi.org/10.1371/journal.pntd.0005086.g001

The number of cases reported per year by subregion for the period of 2001-2011, showed an overall increase of 68% and 39% in the Andean countries and in Central America, respectively. For the same period, in both Brazil and the Southern Cone there was a decline of the total number of cases. However, in the Southern Cone the number of CL cases reported tripled between 2001 and 2002, Fig 1.

Two countries in the Andean subregion (Colombia and Peru), reported an increase of cases during the overall studied period. In Colombia, the reported cases increased 4.36-fold in 2005 versus 2001. In Peru, there was a continuous increase of cases over the period, with the number of cases increasing 2.13-fold (Table 1). Of the total cases reported between 2001 and 2011 from the Andean subregion, Colombia accounted for 50% of the cases, followed by Peru with 33%, except for 2011, when Peru reported 47.7% of the total CL of the subregions.

In the Central America subregion, the countries that contributed to the overall increase in the occurrence of CL for this period were: Costa Rica (4.7 times higher); Honduras (1.8 times higher); Panama (1.3 times higher) and Nicaragua (2.0 times higher) (Table 1). In the cumulative total of cases for the studied period, Nicaragua reported 36.74%, followed by Panama, which reported 26.82%. This situation was similar every year except in 2006, when Panama reported more cases than Nicaragua, Table 1.

CL incidence rate: period of 2001-2011

The annual CL incidence for the 14 selected endemic countries collectively and by subregion is shown in Table 1. The average incidence of CL was 15.89 cases per 100,000 inhabitants. The incidence increased 35% in the first period (2001 to 2005), but in the second half (2006 to 2011) there was a 21% reduction (Fig 2). In the Andean and Central America subregions, the average incidence rates were higher than the average for all 14 countries, with an incidence of 24.50 and 23.70 cases per 100,000 inhabitants respectively.

Download:

*

PPT

PowerPoint slide

*

PNG

larger image

*

TIFF

original image

Fig 2. Annual incidence of cutaneous leishmaniasis per 100,000 inhabitants reported by 14 endemic countries (overall), and stratified by sub region, Americas region, 2001-2011.

The Plotted points correspond to actual CL incidences. The Lines consist in short term trends (fitted with LOESS). Source: PAHO/WHO. Data reported by the leishmaniasis control programs of the countries to PAHO/WHO through May 2012.

http://dx.doi.org/10.1371/journal.pntd.0005086.g002

Colombia and Peru showed the highest incidence rates for the Andean subregion (Table 1). On the other hand, the incidence rate in the Southern Cone declined (in both Argentina and Paraguay), Table 1.

Central America countries presented the highest CL incidence rates during the period, and among these Panama had the highest number, followed by Nicaragua and Costa Rica, which had higher incidence in 2009 (Table 1).

Spatial distribution of CL cases and incidence rates for 2011

Within the 15 analyzed countries, 180 (61.6%) out of 292 administrative units at the first subnational level (departments, states, or provinces) reported cases of CL, with an average of 309.64 cases (ranging from 1 to 3,667). Eight countries (Bolivia, Brazil, Colombia, Costa Rica, Nicaragua, Panama, Paraguay, and Peru) reported cases in more than 75% of their first-subnational-level units. The overall incidence rate for all 15 countries in 2011 was 17.42 cases per 100,000 inhabitants, while for the 180 affected administrative units at the first subnational level, the average incidence rate was 57.52 per 100,000 inhabitants (ranging from 0.03 to 1,229.01).

Stratification at the national level by quartile.

The 2011 results at national level were distributed by quartile, and the 15 countries were stratified into four groups for each indicator, Fig 3A and 3B. The distribution of the countries by quartile was heterogeneous for each indicator. Thus, for case distribution, Brazil, Colombia, and Peru were in the top quartile (with totals ranging from 3,235 to 21,306 cases), whereas for incidence rate distribution, Nicaragua, Panama, and Peru were in the top quartile (with totals ranging from 29.96 to 86.50 cases per 100,000 inhabitants).

Download:

*

PPT

PowerPoint slide

*

PNG

larger image

*

TIFF

original image

Fig 3.

Spatial distribution of cutaneous leishmaniasis cases (A) and incidence rates per 100,000 inhabitants (B), stratified at national level by quartile, Americas region, 2011. (3A) Group 1- countries with ≥ 3 235 cases; Group 2-->1 376 to ≤ 3 234 cases; Group 3-->184 to ≤ 1 375 cases; and Group 4--≤183 cases. (3B) Group 1- Incidence rate ≥ 29.96/100 000 inhabitants; Group 2--incidence rate >11.07 to ≤ 29.95 cases/100 000 inhabitants; Group 3--incidence rate >1.90 ≤ 11.06 cases/100 000 inhabitants; and Group 4--≤ 1.989 cases/100 000 inhabitants. Source: PAHO/WHO. Data reported by the leishmaniasis control programs of the countries to PAHO/WHO through May 2012.

http://dx.doi.org/10.1371/journal.pntd.0005086.g003

Stratification at the first subnational administrative level by quartile.

The data for each indicator on case distribution and incidence rates at the first subnational administrative level for 2011 were stratified by quartile for each indicator, and by administrative units at the first subnational level, Fig 4A and 4B.

Download:

*

PPT

PowerPoint slide

*

PNG

larger image

*

TIFF

original image

Fig 4.

Spatial distribution of cutaneous leishmaniasis stratified by quartile for cases (A) and incidence rate of cases/100,000 inhabitants (B) at the first subnational administrative level, Americas region, 2011. (4A) Group 1--Units at the first subnational level with ≥299; Group 2--> 100 to ≤ 298 cases; Group 3--> 16 to ≤ 99 cases; and Group 4--≤ cases. (4B) Group 1- Units at the subnational administrative level with incidence rates ≥44.81 per 100 000 inhabitants; Group 2--incidence rates > 11.30 to ≤44.80 cases per 100 000 inhabitants; Group 3--incidence rates > 3.43 ≤ 11.29 cases per 100 000 inhabitants; and Group 4--incidence rates ≤ 3.43 cases per 100 000 inhabitants. Source: PAHO/WHO. Data reported by the leishmaniasis control programs of the countries to PAHO/WHO through May 2012.

http://dx.doi.org/10.1371/journal.pntd.0005086.g004

In 2011, spatial distribution by quartile was also heterogeneous for each of the two indicators in the investigation at the first subnational administrative. For example, Guainía in Colombia, Amazonas in Brazil, and Beni in Bolivia were in the top quartile for incidence rates (44.81 to 1,229 cases per 100,000 inhabitants). Paraná in Brazil and Lima in Peru were in the top quartile for case distribution (299 to 3,667 cases), but in the bottom quartile for incidence (0.03 to 3.42 and 3.43 to 11.29 cases per 100,000 inhabitants respectively).

Discussion

For the first time the temporal and geographical distribution of CL cases is described in the American Region. The study find that from 2001 to 2011 a total of 636,683 CL cases were reported by the passive surveillance system from 14 endemic countries representing a 30% increase in the period with an average incidence rate of 15.89/100,000 inhabitants.

In the Americas, most of the CL endemic countries use passive surveillance as a routine reporting system; while active surveillance is mainly conducted during outbreaks or to assist specific inhabitants living in risk areas and of difficult access [15-18]. This study shows the importance of maintenance of adequate strategies for passive surveillance of CL and the need of a system that combines the information from different countries of the Americas, in order to investigate epidemiological features such as time and geographical distribution.

Our findings indicate a 44% increase in the number of CL cases, during the first five years of investigation, also reflected in the incidence rate. This growth could be attributed to an increase of reported cases from the Andean and Central America countries. The reasons for this increase are unknown, but could be explained by a variety of factors, including a true change in incidence or the improvement of clinical and laboratory diagnosis, and enhancement of the CL surveillance system in the endemic countries [15-19] (e.g. Argentina, Brazil, Colombia, Honduras, Paraguay, and Peru [20-23]). Then again, this increase might be even greater if we consider underreporting, which could arise due to lack of access to health services by patients, lack of awareness of the leishmaniasis reporting system by the staff, and the absence of clinical diagnosis, since subclinical or benign forms that cure spontaneously can occur [2, 3, 7, 8, 17, 19].

Furthermore, during the period of investigation, the Region faced several risk factors that contributed to environment and demographic changes (e.g., deforestation, modifications in the land use, increased migration, urbanization, climate change, etc.) that might have led to an increase of contact of human with vectors, therefore affecting the CL case trend [4, 20-29]. Studies conducted in Argentina, Bolivia, Brazil, Colombia, Mexico, and Peru have shown that frequency of leishmaniasis cases are often associated to climate changes and their impacts on the environment and economic activities, including agriculture [4, 19-30]. In Colombia, the Chaparral County reported an outbreak that lasted 5 years (2003-2008) and was considered the largest CL outbreak ever recorded in the country [21]. Also, Argentina recorded outbreaks in the Bella Vista, Formosa, and Salta Counties [31, 32]. In addition, the El Niño-Southern Oscillation (ENSO) has caused a growth of deforestation and of the Phlebotomine population, which has led to the rise of occurrence of CL [33-36]. Studies conducted in Colombia showed an association between increase of CL incidence rates during the El Niño and decline during the La Niña [4, 23]. On the contrary, in Bolivia the CL incidence rates increased during the La Niña and declined during the El Niño [37]. Therefore, climate forecast might be used to predict the leishmaniasis transmission risk in endemic countries affected by both phenomena [38].

The study also shows a different geographical distribution when applied to different countries and subnational levels, illustrating how limiting the use of each indicator can be when higher-risk areas need to be identified to prioritize actions. For example, Brazil has the highest number of cases reported, and at the same time one of the lowest incidence rates due to the size of the population, which emphasizes the need of new analysis strategies, including composite indicators such as the one used in this study. Moreover, it is important to use, when available, complementary information on the vectors, parasites, environment and social characteristic in order to understand the dynamic and traits of this disease [39].

The exploratory analytical tools used in this study to investigate the spatial and temporal distribution of CL in the Americas helped enhance knowledge of the epidemiology of CL and added value to the data reported by the countries. These exploratory methods and their graphic representations are an essential part for the understanding of the complexity of the data correlated over time [40]. This is particularly true in this study where our "sample" is the whole population under surveillance, and for more practical purposes it could be assumed that cases represented approximately the occurrence of CL in the population; in which statistical inference is not as useful as the descriptive method [41].

This study demonstrates that it is possible to provide, retrospectively, an overall picture of the epidemiological patterns of the endemic countries using secondary data, underlining the importance of passive surveillance as a key tool for management of the national control program and the use of the data analysis results for evidence-based programmatic actions and policy formulation. Our results should be interpreted in light of potential limitations posed by the use of data collected under passive surveillance, as showed and discussed in others studies [7, 8, 15-18]. The differences between reporting systems of different countries might affect the comparability of data. The countries included in this study have different types of denominator to estimate the CL incidence rate. Some countries apply the total population, others only the population within rural and urban areas, or only rural counties with CL transmission. The fact that this study uses a standard population as the denominator allowed the assessment and comparison of the incidence rate of CL between the 14 countries included. The results, however, may differ from data available in the countries or from other sources.

The results presented here have shown important advances in the implementation of leishmaniasis surveillance system in the Region, globally demonstrating the epidemiological status of CL in the Americas, but also supporting managers and international organizations in defining, directing and strengthening the actions of surveillance and control of leishmaniasis, in order to reduce cases, deformities and deaths caused by the disease.

Improvement of the surveillance system will allow a better understanding of the ongoing and changing dynamics of the CL transmission. Further operational research is required to design and implement preventive and control strategies [39]. For instance, operational research to determine the advantages in using multiple combined indicators, such as the two indicators used in this study, and to asses complementary indicators (e.g. biological, environmental, social, and other indicators for each subnational administrative entity) in order to generate a more detailed picture of the occurrence of leishmaniasis in the Americas [39, 42]. In addition, a more detailed and integrated analysis of the indicators along with geographical data, disaggregated at the second (or lower) subnational administrative level (counties), would also help to further clarify the epidemiological status of CL at the regional and country level.

Overall, surveillance activities should also be strengthened to provide disaggregated leishmaniasis data to the countries, ideally at the local level. It is also recommended the characterization of smaller territorial units according to their own unique human and physiographic traits (climatological, geomorphological, hydrographical, etc.), which can be highly complex and dynamic, and have the capacity to establish and/or influence an epidemiological situation. Generally, the characteristics of the area (concerning CL occurrence) are the result of interaction among different social groups that share the same space but may have different lifestyles, economies and environment [9, 20, 42].

In conclusion, this investigation indicates changes in the epidemiological patterns of CL in the American Region for unknown reasons, which could be due to changes in the risk factors as well as improvement of the diagnosis and notification. Further operational research and improvement of the surveillance system is needed, combined with activities to strengthen the capacity of the program. The joint use of these actions would lead to a culture of evidence-based in the leishmaniasis program management.

Acknowledgments

The authors thank all of the professionals working in the national leishmaniasis programs and in epidemiological surveillance in the endemic countries of the Americas; Lia Puppim Buzanovsky, for her support with mapping.

Author Contributions

1. Conceptualization: ANSME ZEY.

2. Formal analysis: MJSV FdFdAL.

3. Methodology: ANSME ZEY MJSV FdFdAL MISD.

4. Software: MJSV FdFdAL.

5. Visualization: ANSME ZEY MJSV LGC.

6. Writing - review & editing: ANSME ZEY MJSV MISD.

References

Lainson R, Shaw JJ. Epidemiology and ecology of leishmaniasis in Latin America. Nature. 1978; 273: 595-600. pmid:351409

Desjeux P. Leishmaniasis: current situation and new perspectives. Comp Immunol Microb. 2004; 27: 305-318

World Health Organization. Control of the leishmaniases. Report of a meeting of the WHO Expert Committee on the Control of Leishmaniases, Geneva, 22-26 March 2010. Geneva: World Health Organization; 2010. (WHO technical report series no. 949). Available from: http://apps.who.int/iris/bitstream/10665/44412/1/WHO_TRS_949_eng.pdf. Accessed 22 May 2014

Cárdenas R, Sandoval CM, Rodriguez-Morales AJ, Vivas P. Zoonoses and climate variability: the example of leishmaniasis in southern departments of Colombia. Ann NY Acad Sci. 2008; 1149: 326-30. doi: 10.1196/annals.1428.094. pmid:19120241

Reithinger R, Dujardin JC, Louzir H, Pirmez C, Alexander B, Brooker S. Cutaneous leishmaniasis. Lancet Infect Dis. 2007; 7: 581-596. doi: 10.1016/S1473-3099(07)70209-8. pmid:17714672

GBD 2013 DALYs and HALE Collaborators, et al. Global, regional, and national disability-adjusted life years (DALYs) for 306 diseases and injuries and healthy life expectancy (HALE) for 188 countries, 1990-2013: quantifying the epidemiological transition. Lancet. 2015; 386: 2145-2191. doi: 10.1016/S0140-6736(15)61340-X. pmid:26321261

Alvar J, Vélez ID, Bern C, Herrero M, Desjeux P, Cano J, et al. Leishmaniasis worldwide and global estimates of its incidence. PLoS One. 2012; 7:e35671. doi: 10.1371/journal.pone.0035671. pmid:22693548

Yadón ZE, Quigley MA, Davies CR, Rodrigues LC, Segura EL. Assessment of leishmaniasis notification system in Santiago del Estero, Argentina, 1990-1993. Am J Trop Med Hyg. 2001; 65: 27-30. pmid:11504403

Kawa H, Sabroza PC, Oliveira RM, Barcellos C. A produção do lugar de transmissão da leishmaniose tegumentar: o caso da localidade pau da fome na cidade do Rio de Janeiro, Brasil. Cad Saúde Publica. 2010; 26: 1495-1507. Portuguese. pmid:21229209

Pan American Health Organization. Leishmaniases: Epidemiological Report of the Americas: Washington: Pan American Health Organization; 2013. Available from: http://http://www.paho.org/hq/index.php?option=com_topics&view=readall&cid=6722&Itemid=40754&lang=es. Accessed 04 April 2014. Spanish

Pan American Health Organization. Health situation in the Americas: basic indicators 2012. Washington: Pan American Health Organization; 2012. Available from: http://ais.paho.org/chi/brochures/2012/BI_2012_ENG.pdf. Accessed 10 May 2014

Cleveland WS. Robust locally weighted regression and smoothing scatterplots. J Am Stat Assoc. 1979; 74: 829-836

Cleveland WS, Devlin SJ. Locallyweighted regression: an approach to regression analysis by local fitting. J Am Stat Assoc. 1998; 83: 596-610

R Core Team. R: A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2014. Available from: http://www.Rproject.org/. Accessed 10 December 2014

Copeland H.W, Arana BA, Navin TR. Comparison of Active And Passive Case Detection Of Cutaneous Leishmaniasis In Guatemala. Am J Trop Med Hyg. 1990; 43(3): 257-259 (90-149). pmid:2221220

Rojas CA, Weigle KA, Barrera L, Collazos C. Surveillance and screening of American cutaneous leishmaniasis by Colombian primary health care workers using a clinical prediction rule. Trans R Soc Trop Med Hyg. 2002; 96(4): 405-410. pmid:12497977

Mansuri FA, Borhany T, Kalar M. Factors Responsible for Under Reporting of Notifiable Infectious Diseases By General Practitioners: A Veiled Reality. Biomedica. 2014; 30 (2): 126-129

Seneviratue SL, Gunatilake SB, Silva HJ. Reporting notifiable diseases: methods for improvement, attitudes and community outcome. Trans R Soc Trop Med Hyg. 1997; 91: 135-137. pmid:9196750

Miranda A1, Carrasco R, Paz H, Pascale JM, Samudio F , Saldaña A, Santamaría G, Mendoz Y, Calzada JE. Molecular epidemiology of American tegumentary leishmaniasis in Panama. Am J Trop Med Hyg. 2009; 81(4): 565-571. doi: 10.4269/ajtmh.2009.08-0265. pmid:19815867

Castellanos PL. Epidemiologia, saúde pública, situação de saúde e condições de vida. Considerações conceituais. In: Barata RB, editor. Condições de vida e situação de saúde. Rio de Janeiro: ABRASCO; 1997. p. 31-75. (Série saúde e movimento no. 4). Portuguese

Valderrama-Ardila C, Alexander N, Ferro C, Cadena H, Marín D, Holford TR, et al. Environmental risk factors for the incidence of American cutaneous leishmaniasis in a sub-Andean zone of Colombia (Chaparral, Tolima). Am J Trop Med Hyg. 2010; 82: 243-250. doi: 10.4269/ajtmh.2010.09-0218. pmid:20134000

Salomón OD, Sosa-Estani S, Ramos K, Orellano PW, Sanguesa G, Fernández G, et al. Tegumentary leishmaniasis outbreak in Bella Vista City, Corrientes, Argentina during 2003. Mem Inst Oswaldo Cruz. 2006; 101: 767-774 pmid:17160285

Cárdenas R, Sandoval CM, Rodríguez-Morales AJ, FrancoParedes C. Impact of climate variability in the occurrence of leishmaniasis in Bolivia. Am J Trop Med Hyg. 2006; 75: 273-237. pmid:16896132

Chaves LF, Pascual M. Climate cycles and forecasts of cutaneous leishmaniasis, a nonstationary vector-borne disease. PLoS Med. 2006; 3:e295. doi: 10.1371/journal.pmed.0030295. pmid:16903778

Arias J, Beltrán F, Desjeux P, Walton B. Epidemiología y control de la leishmaniasis en las Américas, por país o territorio. Washington: Pan American Health Organization; 1996. 52 p. (PAHO technical paper no. 44). [cited 2015 May 05]. Available from: http://www.bvsde.paho.org/texcom/cd045364/epi-y-control.pdf. Accessed 05 May 2015. Spanish

Visintin AM, Beranek MD, Amieva MJ, Rosa JR, Almirón WR, Salomón OD. Spread of Phlebotominae in temperate climates: province of Córdoba, Argentina. Mem Inst Oswaldo Cruz. 2016; 111 (1): 75-78. doi: 10.1590/0074-02760150381. pmid:26814647

Davies CR, Reithinger R, CampbellLendrum D, Feliciangeli D, Borges R, Rodriguez N. The epidemiology and control of leishmaniasis in Andean countries. Cad Sáude Pública. 2000; 16: 925-950. pmid:11175518

Desjeux P. Worldwide increasing risk factors for leishmaniasis. Med Microbiol Immuno. 2001; 190: 77-79

Barcellos C, Monteiro AMV, Corvalán C, Gurgel HC, Carvalho MS, Artaxo P, et al. Climate, enviromental changes and infectious diseases: scenarios and uncertains for the Brazilian health surveillance. 7ª Expoepi, 2007. Available from: http://www.researchgate.net/profile/Christovam_Barcellos/publication/242309674_Climate_Environmental_Changes_and_Infectious_Diseases_Scenarios_and_Uncertainties_for_the_Brazilian_Health_Surveillance/links/0c96052b8be1bb14c3000000.pdf. Accessed 03 May 2015

Salomón OD, Acardi SA, Liotta DJ, Fernández MS, Lestani E, López D, et al. Epidemiological aspects of cutaneous leishmaniasis in the Iguazú falls area of Argentina. Acta Trop. 2009; 109: 5-11. doi: 10.1016/j.actatropica.2008.08.002. pmid:18983809

Salomon OD, Sosa-Estanil S, Ramos K, Orellano PW, Sanguesa G, Fernández G, et al.Tegumentary leishmaniasis outbreak in Bella Vista City, Corrientes, Argentina during 2003. Mem. Inst. Oswaldo Cruz. 2006; 101 (7): 767-774. pmid:17160285

Salomón OD, Sosa Estani S, Canini L, Cordoba Lanus E. Tegumentary leishmaniasis in area with epidemic levels of transmission, Salta, Argentina, 1998. Medicina (B Aires). 2001; 61(3): 284-90. Spanish

Zeilhofer P, Kummer OP, Santos ES, Ribeiro AL, Missawa NA. Spatial modelling of Lutzomyia (Nyssomyia) whitmani s.l. (Antunes & Coutinho, 1939) (Diptera: Psychodidae: Phlebotominae) habitat suitability in the state of Mato Grosso, Brazil. Mem Inst Oswaldo Cruz. 2008; 103: 653-60. pmid:19057814

Rodríguez-Morales AJ, Echezuria L, Risquez A. Impact of climate change on health and disease in Latin America, Climate Change and Variability. In: Simard S, editor. ISBN: 978-953-307-144-2, InTech. 2010. Available from: http://www.intechopen.com/books/climate-change-and-variability/impact-of-climate-change-on-health-and-disease-in-latin-america. Accessed 05 May 2015

Salomón OD, Orellano PW, Lamfri M, Scavuzzo M, Dri L, Farace MI, et al. Phlebotominae spatial distribution associated with a focus of tegumentary leishmaniasis in Las Lomitas, Formosa, Argentina, 2002. Mem Inst Oswaldo Cruz. 2006; 101: 295-9. pmid:16862326

Salomón OD, Quintana MG, Mastrángelo AV, Fernández MS. Leishmaniasis and climate change--case study: Argentina. J Trop Med. 2012; 2012: 601242. Epub 2012 May 20

Gomez C, Rodríguez-Morales AJ, Franco-Paredes C. Impact of Climate Variability in the Occurrence of Leishmaniasis in Bolivia. Am J Trop Med Hyg. 2006; 75 (5 Suppl): 42

Franke CR, Ziller M, Staubach C, Latif M. Impact of the El Niño / Southern Oscillation on visceral leishmaniasis, Brazil. Emerg Infect Dis. 2002; 8: 914-917. doi: 10.3201/eid0809.010523. pmid:12194766

Negrão GN, Ferreira MMC. Considerações sobre a Dispersão da Leishmaniose Tegumentar Americana nas Américas revista Percurso--NEMO Maringá. 2009; 1 (1): 85-103. Portuguese

Diggle PJ. Time series: a biostatistical introduction.1st ed. Oxford: Clarendon Press; 1990

Sanchez-Vazquez MJ, Nielen M, Gunn GJ, Lewis FI. Using seasonal-trend decomposition based on loess (STL) to explore temporal patterns of pneumonic lesions in finishing pigs slaughtered in England, 2005-2011. Prev Vet Med. 2012; 104: 65-73

Santos M. A natureza do espaço: técnica e tempo, razão e emoção. 4.ed. 2 reimp.São Paulo:Editora Universidade de São Paulo; 2006. Portuguese