Skip to main content
Advertisement
  • Loading metrics

Clinical and immunological characteristics of tegumentary leishmaniasis cases in Bolivia

  • Cristina Ballart ,

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing

    cristinaballartferrer@ub.edu (CB); albert.picado@finddx.org (AP)

    Affiliations ISGlobal, Hospital Clínic—Universitat de Barcelona, Barcelona, Spain, Secció de Parasitologia, Departament de Biologia, Sanitat i Medi Ambient, Facultat de Farmàcia i Ciències de l’Alimentació, Universitat de Barcelona, Barcelona, Spain

  • Mary Cruz Torrico,

    Roles Data curation, Investigation, Methodology, Project administration, Supervision, Validation, Writing – review & editing

    Affiliations Fundación CEADES Salud y Medio Ambiente, Cochabamba, Bolivia, Facultad de Medicina, Universidad Mayor de San Simón (U.M.S.S.), Cochabamba, Bolivia

  • Gisela Vidal,

    Roles Data curation, Investigation, Writing – review & editing

    Affiliation Fundación CEADES Salud y Medio Ambiente, Cochabamba, Bolivia

  • Faustino Torrico,

    Roles Investigation, Methodology, Project administration, Resources, Writing – review & editing

    Affiliations Fundación CEADES Salud y Medio Ambiente, Cochabamba, Bolivia, Facultad de Medicina, Universidad Mayor de San Simón (U.M.S.S.), Cochabamba, Bolivia

  • Daniel Lozano,

    Roles Investigation, Methodology, Project administration, Resources, Supervision, Writing – review & editing

    Affiliations Fundación CEADES Salud y Medio Ambiente, Cochabamba, Bolivia, Facultad de Medicina, Universidad Mayor de San Simón (U.M.S.S.), Cochabamba, Bolivia

  • Montserrat Gállego,

    Roles Investigation, Validation, Visualization, Writing – original draft, Writing – review & editing

    Affiliations ISGlobal, Hospital Clínic—Universitat de Barcelona, Barcelona, Spain, Secció de Parasitologia, Departament de Biologia, Sanitat i Medi Ambient, Facultat de Farmàcia i Ciències de l’Alimentació, Universitat de Barcelona, Barcelona, Spain

  • Lilian Pinto,

    Roles Investigation, Writing – review & editing

    Affiliation Fundación CEADES Salud y Medio Ambiente, Cochabamba, Bolivia

  • Ernesto Rojas,

    Roles Investigation, Methodology, Project administration, Validation, Writing – review & editing

    Affiliation Facultad de Medicina, Universidad Mayor de San Simón (U.M.S.S.), Cochabamba, Bolivia

  • Ruth Aguilar,

    Roles Formal analysis, Investigation, Validation, Writing – original draft, Writing – review & editing

    Affiliation ISGlobal, Hospital Clínic—Universitat de Barcelona, Barcelona, Spain

  • Carlota Dobaño,

    Roles Formal analysis, Investigation, Validation, Writing – original draft, Writing – review & editing

    Affiliation ISGlobal, Hospital Clínic—Universitat de Barcelona, Barcelona, Spain

  • Sonia Ares-Gomez,

    Roles Formal analysis, Software, Writing – original draft, Writing – review & editing

    Affiliation ISGlobal, Hospital Clínic—Universitat de Barcelona, Barcelona, Spain

  • Albert Picado

    Roles Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing

    cristinaballartferrer@ub.edu (CB); albert.picado@finddx.org (AP)

    Affiliations ISGlobal, Hospital Clínic—Universitat de Barcelona, Barcelona, Spain, Foundation for Innovative New Diagnostics (FIND), Geneva, Switzerland

Abstract

Background

Tegumentary leishmaniasis (TL) is a parasitic disease that can present a cutaneous or mucocutaneous clinical form (CL and MCL, respectively). The disease is caused by different Leishmania species and transmitted by phlebotomine sand flies. Bolivia has one of the highest incidences of the disease in South America and the diagnosis is done by parasitological techniques. Our aim was to describe the clinical and immunological characteristics of CL and MCL patients attending the leishmaniasis reference center in Cochabamba, Bolivia, in order to gain updated clinical and epidemiological information, to evaluate the diagnostic methods used and to identify biomarkers related to clinical disease and its evolution.

Methodology/Principal findings

The study was conducted from September 2014 to November 2015 and 135 patients with lesions compatible with CL or MCL were included. Epidemiological and clinical data were collected using a semi-structured questionnaire. Two parasitological diagnostic methods were used: Giemsa-stained smears and culture of lesion aspirates. Blood samples obtained from participants were used to measure the concentrations of different cytokines. 59.2% (80/135) were leishmaniasis confirmed cases (CL: 71.3%; MCL: 28.7%). Sixty percent of the confirmed cases were positive by smears and 90.6% were positive by culture. 53.8% were primo-infections. Eotaxin and monokine induced by IFN-γ presented higher serum concentrations in the MCL clinical presentation compared to CL cases and no-cases. None of the cytokines presented different concentrations between primo-infections and secondary infections due to treatment failure.

Conclusions/Significance

In Bolivia, parasitological diagnosis remains the reference standard in diagnosis of leishmaniasis because of its high specificity, whereas the sensitivity varies over a wide range leading to loss of cases. Until more accurate tools are implemented, all patients should be tested by both smears and culture of lesion aspirates to minimize the risk of false negatives. Our results showed higher concentrations of several cytokines in MCL compared to CL, but no differences were observed between CL and no-cases. In addition, none of the cytokines differed between primary and secondary infections. These results highlight the need of further research to identify biomarkers of susceptibility and disease progression, in addition to looking at the local cellular immune responses in the lesions.

Author summary

Tegumentary leishmaniasis is a parasitic disease caused by different Leishmania species and transmitted by phlebotomine sand flies. The disease is distributed worldwide and it can present a cutaneous or mucocutaneous clinical form. The clinical outcomes of the disease depend on factors inherent to the parasite, the vector and the host. Bolivia has one of the highest incidences of the disease in South America, where the diagnosis is mainly done by parasitological techniques. The present study was conducted from September 2014 to November 2015 with 135 patients attending the leishmaniasis reference centre in Cochabamba (Bolivia). We aimed to describe the clinical, epidemiological and immunological characteristics of leishmaniasis patients to gain knowledge on tegumentary leishmaniasis and evaluate the suitability of the diagnostic methods used. In South America, Giemsa-stained smears and the culture of lesion aspirates remain the reference standards in diagnosis because of their high specificity, however the sensitivity varies over a wide range resulting in some loss of case detection. Until more accurate tools are implemented, all patients should be screened by both tests to minimize the risk of false negatives. Additionally, for better disease management, identification of biomarkers of susceptibility and disease progression would be desirable.

Introduction

Tegumentary leishmaniasis (TL) has two different clinical presentations: ulcerative skin lesions (cutaneous leishmaniasis [CL]) or a destructive mucosal inflammation (mucocutaneous leishmaniasis [MCL]) that affects oral-nasal-pharyngeal cavities. MCL usually appears after a CL episode but both clinical forms can present together [1]. The clinical outcomes of the disease depend on factors inherent to the parasite, the vector, and the host. CL patients should be diagnosed and treated even if CL lesions may self-resolve in up to 70–80% of cases, as they can cause a high degree of morbidity, social stigmatization, and in some cases evolve to MCL [1].

Bolivia is one of the endemic countries with the highest incidence of CL and MCL in South America [1,2]. The Leishmania (Viannia) braziliensis species, transmitted by phlebotomine sand flies, is the predominant etiological agent causing CL in Bolivia (85% of cases) [1,2,3]. CL is endemic to 70% of Bolivia, with 54.7 cases per 100,000 individuals in 2018 [1,4,5,6].

In Bolivia, diagnosis of leishmaniasis is done by Parasitological Direct Exam (PDE) and in a few cases, by parasite culture method (CM). Patients cover the cost of both diagnostics. TL treatment is free of charge but it is only provided when the case has been confirmed. In Bolivia, first-line TL therapy is meglumine antimoniate (Glucantime), administered by intravenous or intramuscular injection, depending on the clinical form of the disease (CL and MCL, respectively), at a dosage of 20 mg SbV/kg/day over a period of 20 to 30 days. Second-line treatment is amphotericin B deoxycholate (Fungizone) at a dosage of 0.5-1mg/kg/day administered intravenously for 30 to 45 days requiring a prolonged hospitalization [4,7,8]. However, therapy failures and variable efficacy of these drugs have been reported [1,7,9].

This study aimed to describe the clinical and immunological characteristics of CL and MCL patients attending the leishmaniasis reference center in Cochabamba. The study is relevant to have a realistic description of the population affected by TL in this low-income country area, and to identify potential problems associated with disease management, such as gaps in the diagnostic tools currently used. The study also assessed the usefulness of cytokines as biomarkers of clinical presentation and treatment failure, and highlights the need for further research on biomarkers of susceptibility and disease progression in leishmaniasis.

Materials and methods

Ethics statement

The study was approved by the Ethics Committees of the Hospital Clínic de Barcelona, Spain (HCB/2014/0582); CEADES Salud y Medio Ambiente and the Facultad de Medicina UMSS in Bolivia.

All suspected cases of leishmaniasis provided written informed consent (parents or guardians in case of patients were under 18 years old) before enrolling in the study. All leishmaniasis suspects attending LABIMED during the study period, including those not enrolled in the study, were diagnosed for free. CL and MCL cases were referred for treatment.

Study population, data and sample collection

The study was conducted in patients attending the Institute of Biomedical Research, IIBISMED, and Laboratorios de Investigación Biomédica (LABIMED), from the Faculty of Medicine of the Universidad Mayor de San Simón (UMSS) in Cochabamba city (Bolivia). Suspects of leishmaniasis are referred to the IIBISMED from other health centres for diagnosis and treatment.

Patients with lesions compatibles with CL or MCL attending the LABIMED from September 2014 to November 2015 were invited to participate in the study. Subjects from both sexes and all ages were recruited. After obtaining written informed consent, data on the patients’ demographic, epidemiological and clinical characteristics (number, location and lesion type, time of evolution, previous drug treatment, duration, and outcome) were collected using a semi-structured questionnaire (S1 Fig). Smears and aspirates of active lesions were obtained to diagnose CL and MCL cases. More than one sample was collected per patient if multiple lesions were present. Blood samples were also obtained by venipuncture.

Diagnostic methods

The diagnostic methods used were: i) microscopic examination of lesion smears obtained by scrape and stained with Giemsa for the visualization of amastigotes (PDE) and ii) culture of lesion aspirates from areas of apparent parasitic activity in TSTB (Torrico-Solano-Torrico-Bermúdez) biphasic medium (CM) [10]. Cultures were maintained for four weeks and considered either positive or negative if parasite growth was detected or not and considered contaminated if non-specific microorganisms were also present. PDE and CM were performed on samples from three different areas of the same lesion for all the patients who participated in the study.

Definitions

Based on the clinical and epidemiological data and the results of the diagnostic tests, the following groups of patients were defined.

Clinical forms.

Individuals who were positive by microscopy and/or culture were considered leishmaniasis-confirmed patients. Patients negative by microscopy and culture were considered no-cases. Patients presenting only cutaneous lesions in any area of the body were classified as CL, and those with mucosal lesions in mouth and/or nostrils, or concomitant cutaneous and mucosal lesions, were classified as MCL.

Primo and secondary infections.

CL patients who had no previous history of leishmaniasis were defined as primo-infection (first episode). MCL patients with no history of CL or visible scars were considered to have primary lesions [11]. The rest were considered secondary infections (second episode).

Clinical cure and treatment failure.

Clinical cure was defined as complete wound closure and re-epithelization without inflammation or filtration [12].

We used historical clinical data to identify patients who had a relapse due to treatment failure. CL and MCL patients who had previously received a complete anti-leishmanial treatment in the past but were not cured, regardless of its duration, were considered treatment failures. Those patients presented an absent or incomplete scarring of lesion(s) and/or a persistence of inflammation around the initial lesion, a clinical regression of a healed lesion or the presence of new mucosal lesion(s) [1,7,12,13]. For those cases with cutaneous lesions to be considered treatment failures, the previous treatment had to be directed to the same lesion(s) identified during our study.

Cytokine, chemokine and growth factor multiplex bead array assay

The Cytokine Human Magnetic 30-Plex Panel from Life Technologies was used to measure the concentrations (pg/mL) of the following cytokines, chemokines and growth factors in sera: epidermal growth factor (EGF), fibroblast growth factor (FGF), granulocyte colony-stimulating factor (G-CSF), granulocyte-macrophage colony-stimulating factor (GM-CSF), hepatocyte growth factor (HGF), vascular endothelial growth factor (VEGF), tumor necrosis factor (TNF), interferon (IFN)-α, IFN-γ, interleukin (IL)-1RA, IL-1β, IL-2, IL-2R, IL-4, IL-5, IL-6, IL-7, IL-8, IL-10, IL-12(p40/p70), IL-13, IL-15, IL-17, IFN-γ induced protein (IP-10), monocyte chemoattractant protein (MCP-1), monokine induced by IFN-γ (MIG), macrophage inflammatory protein (MIP)-1α, MIP-1β, regulated on activation normal T cell expressed and secreted (RANTES) and eotaxin.

Twenty five μL of sera were tested by applying a modification to the manufacturer’s protocol by using half the volume of all reagents including the standards. This modification was previously tested and showed no difference in assay performance compared to the original protocol and has been used in prior studies [14,15]. Each plate included 16 serial dilutions (2-fold) of a standard sample with known concentrations of each analyte, two blanks and three positive controls of high, medium and low concentrations in duplicates, prepared from the standard for quality control purposes. Samples were acquired on a Luminex 100/200 instrument and analyzed with xPONENT software 3.1. Concentrations of analytes were obtained by interpolating the median fluorescent intensity (MFI) (after blank MFI subtraction) to a 5-parameter logistic regression curve calculated by xPONENT. Values below the lower limit of detection (mean of blanks + 2 standard deviations) were assigned half the expected concentration at the lower limit of quantification, and values above the upper limit of quantification were assigned twice the expected concentration at the upper limit of quantification.

Statistical analysis

Clinical and epidemiological variables were all categorical and included: age, sex, occupation, clinical status, probable place of infection, episode of the disease, type and location of lesion, and questions related to treatment (previous treatments, drug used, dosage and resolution), among others. Differences in sex, age, number of episodes, characteristics of lesions, and characteristics of previous treatments between clinical forms, were assessed by Pearson’s Chi-square test at the significance level α = 0.05. The number of patients identified by both diagnostic techniques (PDE and CM) was estimated. The agreement between the diagnostic methods was evaluated using the Cohen’s kappa coefficient (K) calculated after excluding those cultures that were contaminated. Calculations were performed using the function epi.kappa from the epiR library in R version 3.5.2 [16].

Cytokine, chemokine and growth factor concentrations were tested for normality to inform the appropriate statistical analysis, showing a non-normal distribution. Differences in analyte concentrations between CL, MCL and no-cases, and between primo-infections and treatment failures (shown as box plots), were compared through the Kruskal-Wallis Chi-squared test adjusting for multiple comparisons by the Bonferroni correction. The effect of the clinical manifestation (CL or MCL) on analyte concentrations was assessed through univariable and multivariable linear regression models. Those analytes significant in the univariable models were included in a stepwise backward multiple regression selection at α = 0.05. The final selected model was adjusted by age. IBM SPSS Statistics version 24.0.0 (IBM Corp. Released 2016. IBM SPSS Statistics for Windows, Version 24.0.0 Armonk, NY: IBM Corp.) was used to conduct the statistical analyses related to the cytokine, chemokine, and growth factor concentrations.

Results

From the 141 suspected leishmaniasis patients who reported to IIBISMED during the study period, 135 (96%) provided informed consent and were included in the study. From those, 78% (105/135) and 22% (30/135) were men and women, respectively, from one to 87 years old (≤10: 6%; [1120]: 13%; [2140]: 47%; [41–60]: 21%; >60: 13%). Eighty (59%) patients were leishmaniasis-confirmed cases. Excluding the patients with CM contamination (n = 5) we estimated that only 60% of the cases (45/75; 33 CL and 12 MCL cases) would have been diagnosed by PDE alone, and 90.6% (68/75; 51 CL and 17 MCL cases) if only CM was conducted. Thirty-eight out of 75 cases (50.6%) were positive for both PDE and CM (28 CL and 10 MCL cases), showing a moderate agreement between both diagnostic techniques (K = 0.43; 95% CI: 0.29–0.58).

The prevalence (71.3% CL and 28.7% MCL) and demographic characteristics of leishmaniasis confirmed cases are shown in Table 1. Among the MCL cases, 19 presented a scar from an old cutaneous lesion and 4 presented active concomitant cutaneous and mucosal lesions. The majority of leishmaniasis cases were found in men (83.8%), but taking into account the relationship between the number of confirmed and suspected cases for each sex, no statistical difference was observed (p = 0.4) (men: 67 confirmed cases/105 suspected cases (p = 0.6): women: 13 confirmed cases/30 suspected cases (p = 0.4)). The age range most affected was between 21–40 years old (58.8%; p = 0.004) with a population median age being 32.5 years old. All patients ≤ 20 years old presented a CL form. Forty-three cases (53.8%) were primo-infections, and two MCL cases without history of cutaneous lesions were classified as primary lesions (2.4%) (Table 1).

thumbnail
Table 1. Patient characteristics by clinical manifestation of leishmaniasis disease attending the IIBISMED Service in Cochabamba Department between 2014 and 2015.

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

The location, number, and type of lesions, as well as time of evolution are shown in Table 2. Among CL cases, infected ulcerative lesions were the most common (57.9%). Single lesions were the most frequent (61.4%) (Median: 1; Maximum number of lesions: 8), predominantly in lower members (43.9%) and with less than 6 months of evolution (82.5%). MCL patients presented mainly mucosal inflammation.

thumbnail
Table 2. Characteristics of clinical forms of leishmaniasis disease in patients attending the IIBISMED Service in Cochabamba Department between 2014 and 2015.

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

Information about the previous treatment in patients presenting a second episode (same lesion or other lesions in different locations) is summarized in Table 3. Fourteen patients were considered as treatment failure (considering those who had received a complete treatment): including ten patients with CL and four patients with MCL. Glucantime (meglumine antimoniate) was the most frequently used drug (85%), followed by amphotericine B (5%). No information on the treatment and its duration was obtained for the rest of patients (10%). The large majority of patients receiving a second treatment completed it in the past (14/20; 70%) and, from those, 42.9% (n = 6) showed complete recovery for the past leishmaniasis episode. One primo-infection patient died during the treatment after an apparent recovery, however the cause of death was not recorded, although he was reported to have severe malnutrition.

thumbnail
Table 3. Information about the previous treatment in the 35 patients presenting a second episode of the disease.

https://doi.org/10.1371/journal.pntd.0009223.t003

Cytokine profiles differ between CL and MCL clinical forms

Higher concentrations of several cytokines were found in the MCL clinical cases compared to CL cases (Fig 1): the growth factors VEGF (P = 0.008) and HGF (P = 0.018); the chemokines eotaxin (P = 0.005), MIP-1α (P = 0.025), MIG (P<0.001) and IP-10 (P = 0.006); the pro-inflammatory cytokine IL-6 (P = 0.003) and the TH2 cytokine IL-13 (P = 0.042). However, after correcting for multiple comparisons only MIG remained significantly higher (P = 0.006). Accordingly, higher concentrations of most of these cytokines were associated with the MCL phenotype in the univariable regression models, however, after the stepwise multiple regression selection and adjusting by age, only eotaxin and MIG remained significant (Table 4).

thumbnail
Fig 1. Differences in cytokine concentrations between CL, MCL and no-cases.

Box plots representing the median and interquartile range of each marker concentration (pg/mL). Levels between both groups were compared by the Kruskal-Wallis Chi-squared test adjusting for multiple comparisons by using the Bonferroni correction (P-values in brackets). Significant P-values (<0.05) are in bold.

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

thumbnail
Table 4. Univariate and multivariable linear regression models to assess the effect of the MCL clinical manifestation on cytokine concentrations.

Linear regression models to assess the effect of MCL on cytokines concentrations.

https://doi.org/10.1371/journal.pntd.0009223.t004

Cytokine profiles do not differ between primo and secondary infections

None of the cytokines presented different concentrations between primo-infections and secondary infections due to treatment failure. Only the growth factor EGF showed a higher trend in the secondary infections, but the difference was borderline (P = 0.055), and not statistically significant after multiple comparisons (P = 1) (S2 Fig).

Discussion

In South America, parasitological diagnosis remains the reference standard in diagnosis of CL because of its high specificity, whereas the sensitivity varies over a wide range, and is much lower in MCL cases due to the scarcity of amastigotes [17]. In this study, nearly 60% of the 135 CL and MCL suspects tested by microscopy and culture were confirmed to be leishmaniasis cases. This result is consistent with the prevalence previously found in the tropic of Cochabamba, where 64% of cutaneous lesions proved to be leishmaniasis, with the rest being from other etiological origins [18]. However, it must be pointed out that some cases are missed due to the low sensitivity of the two diagnostic methods [17]. Until more accurate tools are implemented (e.g. PCR) all patients should be screened by both tests to minimize the risk of false negatives. As found in the current study, by only using PDE which is a more subjective technique, 40% of the cases would have been missed due to its low sensitivity that ranges from 15 to 70% [19]. On the contrary, by using CM alone, only 9.4% of cases would have been overlooked. Five patients could be diagnosed solely by PDE due to culture contamination (6.3%), which were mainly MCL cases. In our study, the rate of culture contamination was low (9/135; 6.7%), and stemmed only from fungal microorganisms. The low frequency of contamination was most likely due to the high concentration of antibiotic in the culture medium, the proper cleaning of the lesions before sampling and the expertise of the laboratory staff [10]. The over-infection of lesions was frequently due to lack of hygiene and/or the use of domestic treatments, causing a decrease in the sensitivity of the diagnostic techniques if not solved [20].

In this study, 18% of cases were considered to be treatment failures according to the criteria described above [1,7,12,13]. Llanos-Cuentas and colleagues found that the risk factors associated with treatment failure were young age, presence of an additional lesion and infection with L. braziliensis [12]. Patients not included in the main groups (primo-infection or treatment failure) had different characteristics and were difficult to identify. The majority did not finish treatment, and others used different plants or solvents to dry the wound or simply burn it. Occasionally, the lesion self-resolved leaving a scar, but this is not a guarantee of parasitological cure, and MCL could still appear after several years [1,20]. Mucosal lesions required hospitalization to be treated and may be fatal in some cases [1].

Since the study only included patients arriving passively to be diagnosed, a considerable number of undiagnosed might have been missed [1,4,21]. Access barriers to case management were associated with the absence of diagnostic methods and trained staff in rural areas, their cost, limited accessibility to health centres (scarce and expensive transport, poor road conditions, long travel distances), and adverse climate conditions [4,10,20,21,22]. The local languages (Quechua and Aymara) and the stigmatization of the disease also hampers communication with the patients [4,22].

Locations where patients may have acquired the infection were deduced from interviews, but in some cases this information was not available. As found in other studies, most infections appeared to be acquired in the tropical region of Chapare province (56%), in the department of Cochabamba [22,23]. The risk of transmission in this area is the highest in the region, as shown by entomological studies [24,25,26]. The traditional sylvatic pattern of transmission is the most prevalent in the area [3,22,23,25], but migrations to endemic areas and leishmaniasis urbanization are also potential means of infection as seen in other countries [22,27,28]. A large percentage of our population with therapeutic failure had a history of visiting endemic areas for only short periods, therefore lacking previous exposure to parasite antigens and sand fly saliva that confer protective immunity [5].

In the present study, the age range most affected was between 21–40 years old and infections probably occurred during work activities, with more frequent lesions found in the exposed areas of the body, mainly lower members [1,22,28,29]. These patients described field activities as farming, fruit and vegetable agriculture, or coca cultivation. Occasionally, women and children accompanied men to their job activities [23,29]. However, it is worth noting that the number of women self-presenting was three times lower (30/135) than men. This may be related to a genuine lower prevalence among women and/or to differences in access to health services associated to sex as described before [28,29,30]. In addition, leishmaniasis shows sex-related differences regarding physiology, immune responses, drug metabolism and progression of the disease [31].

Cellular immunity is essential for controlling intracellular pathogens like Leishmania, and of primary importance in the outcome of this infection. In the mouse model of L. major, resistance and susceptibility are determined by functionally distinct T cell subsets [32,33], with a predominant T helper (Th)1 response leading to the control of the infection, and a Th2 response increasing susceptibility and disease progression [34,35]. However, an exaggerated Th1 response has been associated with tissue damage [36,37]. Our results showing higher concentrations of cytokines in MCL agree with previous studies reporting that lymphocyte activation, proliferation and cytokine secretion upon stimulation with Leishmania antigens are higher in MCL compared with CL [37,38], probably related to the severity of tissue destruction in MCL. Among the analytes increased in MCL, we found mainly growth factors (VEGF and HGF) and chemokines (eotaxin, IP-10, MIG and MIP-1α), as well as the pro-inflammatory cytokine IL-6 and the Th2 cytokine IL-13. Higher concentrations of growth factors in MCL may be related to the more severe and extensive tissue destruction and the need of a more active tissue repair. VEGF is associated with changes in the vasculature during wound healing, and has been suggested to be important in the development and resolution of leishmanial lesions [39]. Similarly, HGF is a potent angiogenic factor having a major role in wound healing, acting upon epithelial and endothelial cells, and on haemopoietic progenitor cells and T cells [40].

The higher concentration of several chemokines in MCL could also be related to the more extensive wound healing, as they stimulate the migration of multiple cell types to the wound site, particularly inflammatory cells [41], and contribute to the regulation of re-epithelialization, tissue remodelling, and angiogenesis [42]. In addition, chemokines may have antimicrobial activity involved in parasite clearance. Eotaxin may chemoattract eosinophils to the infection site for Leishmania clearance [43]. Similarly, IP-10, MIG and MIP-1α, three Th1-mobilizing chemokines produced by leishmaniasis lesion cells, recruit monocytes, macrophages, and activated T cells to sites of infection, contributing to wound healing and parasite elimination [44,45]. Neutrophil derived MIP-1α has also been shown to be essential for the recruitment of dendritic cells to the infection site that will further direct the development of an adaptive immune response [46].

The higher concentration of IL-6 in MCL agrees with a previous association with severity in visceral leishmaniasis [47,48]. This detrimental effect of IL-6 has been related to its suppressive effect on TNF production and Th1 responses, and its induction of Th2 responses, which is suggested to down-regulate macrophage microbicidal activity, facilitating Leishmania proliferation [47]. Finally, the higher concentration of IL-13 in MCL may have a dual role (mainly Th2 but also Th1 responses) that is influenced by parasite species and the host’s genetic background [48,49].

Among all these biomarkers, only eotaxin and MIG showed significantly higher concentrations in MCL versus CL after adjustment, indicating their relevant role in the immune response against the MCL clinical form. Their contribution in recruiting cells to the infected tissues may be favoring macrophage activation and parasite clearance; however, excessive production of these molecules could potentially lead to uncontrolled inflammation and tissue destruction [45,50]. Thus, it is difficult to disentangle whether the higher concentrations of these chemokines in MCL versus CL are a positive response to a more extensive infection, or are in fact contributing to the more generalized inflammation and tissue damage observed in MCL cases. The absence of association between the levels of cytokines and primo-infections or secondary infections suggests similar cellular immune responses among them independently of the type of infection.

One limitation of the study is that the analyses did not consider the Leishmania species causing the CL and MCL lesions as those were not identified. However, based on previous studies in the same area, we would expect that over 85% of the cases were caused by L. braziliensis [2,30]. The intrinsic variability of the circulating Leishmania species influences treatment efficacy, thus it is important to be considered for the control of the disease [2]. In addition, the related lack of data on resistance to antimonials is also important for disease control, but its incorporation into the diagnostic routine is not realistic in the Bolivian context. The study does not include healthy individuals as control group for the comparison of cytokine concentrations with leishmanial cases. We only had sera from patients with lesions compatible with leishmaniasis but not confirmed as cases, however some cases may have been missed due to the low sensitivity of the current diagnostic methods [17]. Finally, quality and quantity of inflammatory cells in the lesions may differ between clinical phenotypes and patients, affecting the local cytokine profiles, and this may not be reflected in the serum cytokine profiles. Future studies should address the local cellular immune response in the lesions to better understand its association to wound healing and evolution of treatment.

In conclusion, the results of the present study provide additional information on the epidemiological, clinical and immunological situation of TL in Bolivia. Although parasitological diagnosis remains the reference standard for leishmaniasis the diagnostic techniques evaluated in the present study evidence that more accurate tools must be implemented in the future to avoid false negatives. Also, further studies must be done in a larger number of patients relating clinical manifestations with species characterization considering that Leishmania species identification might be helpful for better disease management. Regarding the immunological analysis these results also highlight that further research would certainly present a further step towards improving knowledge of potential biomarkers of susceptibility and disease progression.

Supporting information

S1 Fig. Semi-structured questionnaire.

Demographic, epidemiological and clinical data collected from the patients recruited for the study.

https://doi.org/10.1371/journal.pntd.0009223.s001

(TIF)

S2 Fig. Differences in cytokine concentrations between primo and secondary infections.

Box plots representing the median and interquartile range of each marker concentration (pg/mL). Levels between both groups were compared by the Kruskal-Wallis Chi-squared test adjusting for multiple comparisons by using the Bonferroni correction (P-values in brackets). Significant P-values (<0.05) are in bold.

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

(TIF)

Acknowledgments

The contents of this publication are the sole responsibility of the authors. The authors would like to thank all the staff of LABIMED for their collaboration, and IIBISMED to the support and helpful scientific discussions in the study. We thank Robert A. Mitchell proofreading. ISGlobal receives support from the Tropical Disease Cooperative Research Network (RICET) (RD12/0018/0010) and the Spanish Ministry of Science, Innovation and Universities through the “Centro de Excelencia Severo Ochoa 2019–2023” Program (CEX2018-000806-S). ISGlobal is member of the CERCA Program, Generalitat de Catalunya.

References

  1. 1. García AL, Parrado R, Rojas E, Delgado R, Dujardin JC, Reithinger R. Leishmaniases in Bolivia: comprehensive review and current status. Am J Trop Med Hyg. 2009;80(5):704–711. pmid:19407110
  2. 2. Bilbao-Ramos P, Dea-Ayuela MA, Cardenas-Alegría O, Salamanca E, Santalla-Vargas JA, Benito C, et al. Leishmaniasis in the major endemic region of Plurinational State of Bolivia: Species identification, phylogeography and drug susceptibility implications. Acta Trop. 2017;176:150–161. pmid:28751163
  3. 3. Rojas E, Parrado R, Delgado R, Reithinger R, Garcia, AL. Leishmaniasis in Chaparé, Bolivia. Emerg Infect Dis. 2009;15(4):678–680. pmid:19331774
  4. 4. Alvar J, Vélez ID, Bern C, Herrero M, Desjeux P, Cano J, et al. World Health Organization, WHO. Leishmaniasis Control Team. Leishmaniasis worldwide and global estimates of its incidence. PLoS One. 2012;7(5):e35671. pmid:22693548
  5. 5. Davies CR, Reithinger R, Campbell-Lendrum D, Feliciangeli D, Borges R, Rodriguez N. The epidemiology and control of leishmaniasis in Andean countries. Cad Saude Publica. 2000;16(4):925–950. pmid:11175518
  6. 6. Organización Panamericana de la Salud. Leishmaniasis: informe epidemiológico de las Américas. Washington, D.C.: OPS; 2019. http://iris.paho.org/xmlui/handle/123456789/505505.
  7. 7. Bermúdez H, Rojas E, García L, Desjeux P, Dujardin JC, Boelaert M, et al. Generic sodium stibogluconate is as safe and effective as branded meglumine antimoniate for the treatment of tegumentary leishmaniasis in Isiboro Secure Park, Bolivia. Ann Trop Med Parasitol. 2006;100:591–600. pmid:16989685
  8. 8. David CV, Craft N. Cutaneous and mucocutaneous leishmaniasis. Dermatol Ther. 2009;22(6):491–502. pmid:19889134
  9. 9. Rojas R, Valderrama L, Valderrama M, Varona MX, Ouellette M, Saravia NG. Resistance to antimony and treatment failure in human Leishmania (Viannia) infection. J Infect Dis. 2006;193:1375–1383. pmid:16619185
  10. 10. Bermúdez H, Solano M, Torrico MC, Carballo M, Lafuente O, Lara M, et al. Diagnóstico de leishmaniasis utilizando medio de cultivo TSTB en pacientes de trópico de Cochabamba. Gac Med. 2005;31–35.
  11. 11. Rojas E. Leishmaniasis mucosa laríngea sin puerta de entrada evidente: presentación de un caso clínico. Gac Med Bol. 2011;34(2):91–93.
  12. 12. Llanos-Cuentas A, Tulliano G, Araujo-Castillo R, Miranda-Verastegui C, Santamaria-Castrellon G, Ramirez L, et al. Clinical and parasite species risk factors for pentavalent antimonial treatment failure in cutaneous leishmaniasis in Peru. Clin Infect Dis. 2008;46(2):223–231. pmid:18171254
  13. 13. Arevalo J, Ramirez L, Adaui V, Zimic M, Tulliano G, Miranda-Verástegui C, et al. Influence of Leishmania (Viannia) species on the response to antimonial treatment in patients with American tegumentary leishmaniasis. J Infect Dis. 2007;195(12):1846–1851. pmid:17492601
  14. 14. Aguilar R, Campo JJ, Chicuecue S, Cisteró P, Català A, Luis L, Ubillos I, Galatas B, Aide P, Guinovart C, Moncunill G, Dobaño C. Changing plasma cytokine, chemokine and growth factor profiles upon differing malaria transmission intensities. Malar J. 2019;18(1):406. pmid:31806027
  15. 15. Natama HM, Moncunill G, Rovira-Vallbona E, Sanz H, Sorgho H, Aguilar R, Coulibaly-Traoré M, Somé MA, Scott S, Valéa I, Mens PF, Schallig HDFH, Kestens L, Tinto H, Dobaño C, Rosanas-Urgell A. Modulation of innate immune responses at birth by prenatal malaria exposure and association with malaria risk during the first year of life. BMC Med. 2018;16(1):198. pmid:30384846
  16. 16. R Core Team (2018). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. URL. Available from: https://www.R-project.org/
  17. 17. World Health Organization (WHO). Report of a meeting of the WHO Expert Committee on the Control of Leishmaniasis.WHO Technical Report Series, 2010; 949.
  18. 18. Bermúdez H, Zapata ME. Etiología y características epidemiológicas y clínicas de las úlceras cutáneas del trópico de Cochabamba. UMMS, DUI, OMS. Informe final, Curso de Diplomado en Investigación en Salud; 1997.
  19. 19. Torrico MC, Zubieta MS. Manual de Normas y Procedimientos Técnicos de Laboratorio para Leishmaniasis; Murillo, La Paz, Bolivia. Ministerio de Salud y Deportes; 2010. ISBN 978-99954-50-13-7. 100pp.
  20. 20. Gómez H, Janco A, Aguilar D. Leishmaniasis tegumentaria en Bolivia. Artículos de revisión. CUMETROP, Facultad de Medicina de la UMMS. 2003:68–70.
  21. 21. Eid D, Guzman-Rivero M, Rojas E, Goicolea I, Hurtig AK, Illanes D, San Sebastian M. Assessment of a Leishmaniasis Reporting System in Tropical Bolivia Using the Capture-Recapture Method. Am J Trop Med Hyg. 2017. pmid:29141751
  22. 22. Eid D, Guzman-Rivero M, Rojas E, Goicolea I, Hurtig AK, Illanes D, San Sebastian M. Risk factors for cutaneous leishmaniasis in the rainforest of Bolivia: a cross-sectional study. Trop Med Health. 2018;46:9. pmid:29692654
  23. 23. Maita García X, Miranda Gutierrez C, Marañón Mendoza LC, Carvajal Yáñez N, Santander López A. Características epidemiológicas de la Leishmaniasis en el departamento de Cochabamba durante el periodo 2002–2010. Rev Cient Cienc Med. 2011;14(2):8–11.
  24. 24. Ballart C, Vidal G, Picado A, Cortez MR, Torrico F, Torrico MC, et al. Intradomiciliary and peridomiciliary captures of sand flies (Diptera: Psychodidae) in the leishmaniasis endemic area of Chapare province, tropic of Cochabamba, Bolivia. Acta Trop. 2016;154:121–124. pmid:26608724
  25. 25. Bustamante M, Diaz M, Espinoza J, Parrado R, Reithinger R, García AL. Sand fly fauna in Chapare, Bolivia: an endemic focus of Leishmania (Viannia) braziliensis. J Med Entomol. 2012;49(5):1159–1162. pmid:23025199
  26. 26. García AL, Tellez T, Parrado R, Rojas E, Bermudez H, Dujardin JC. Epidemiological monitoring of American tegumentary leishmaniasis: molecular characterization of a peridomestic transmission cycle in the Amazonian lowlands of Bolivia. Trans R Soc Trop Med Hyg. 2007;101(12):1208–1213. pmid:17945322
  27. 27. Campbell-Lendrum D, Dujardin JP, Martinez E, Feliciangeli MD, Perez JE, Silans LN, et al. Domestic and peridomestic transmission of American cutaneous leishmaniasis: changing epidemiological patterns present new control opportunities. Mem Inst Oswaldo Cruz. 2001;96(2):159–162. pmid:11285490
  28. 28. Tedesqui VL, Calleja GN, Parra R, Pabón JP, Bóia MN, Carvalho-Costa FA. Active surveillance of American tegumentary leishmaniasis in endemic areas in rural Bolivia. Rev Soc Bras Med Trop. 2012;45(1):30–34. pmid:22370825
  29. 29. Sosa-Estani S, Segura EL, Gomez A, Salomón OD, Peralta M, Coutada V, et al. Cutaneous leishmaniasis in Northern Argentina: identification of risk factors in a case-cohort study of three municipalities in Salta. Rev Soc Bras Med Trop. 2001;34(6):511–517. pmid:11813056
  30. 30. García Bustos MF, González-Prieto G, Ramos F, Mora MC, Hashiguchi Y, Parodi C, et al. Clinical and epidemiological features of leishmaniasis in northwestern-Argentina through a retrospective analysis of recent cases. Acta Trop. 2016;154:125–132. pmid:26611809
  31. 31. Lockard RD, Wilson ME, Rodríguez NE. Sex-Related Differences in Immune Response and Symptomatic Manifestations to Infection with Leishmania Species. J Immunol Res. 2019;10;2019:4103819. pmid:30756088
  32. 32. Gumy A, Louis JA, Launois P. The murine model of infection with Leishmania major and its importance for the deciphering of mechanisms underlying differences in Th cell differentiation in mice from different genetic backgrounds. Int J Parasitol. 2004;34:433–444. pmid:15013733
  33. 33. Reiner SL, Locksley RM. The regulation of immunity to Leishmania major. Annu Rev Immunol. 1995;13:151–177. pmid:7612219
  34. 34. Coutinho SG, Oliveira MP, Da Cruz AM, De Luca PM, Mendonça SC, Bertho AL, et al. T-cell responsiveness of American cutaneous leishmaniasis patients to purified Leishmania pifanoi amastigote antigens and Leishmania braziliensis promastigote antigens: immunologic patterns associated with cure. Exp Parasitol. 1996;84:144–155. pmid:8932764
  35. 35. Von Stebut E. Cutaneous Leishmania infection: progress in pathogenesis research and experimental therapy. Exp Dermatol. 2007;16:340–346. pmid:17359341
  36. 36. Bacellar O, Lessa H, Schriefer A, Machado P, Ribeiro de Jesus A, Dutra WO, et al. Up-regulation of Th1-type responses in mucosal leishmaniasis patients. Infect Immun. 2002;70:6734–6740. pmid:12438348
  37. 37. Carvalho EM, Johnson WD, Barreto E, Marsden PD, Costa JL, Reed S, et al. Cell mediated immunity in American cutaneous and mucosal leishmaniasis. J Immunol. 1985;135:4144–4148. pmid:4067312
  38. 38. Ribeiro-de-Jesus A, Almeida RP, Lessa H, Bacellar O, Carvalho EM. Cytokine profile and pathology in human leishmaniasis. Braz J Med Biol Res. 1998;31:143–148. pmid:9686192
  39. 39. Weinkopff T, Konradt C, Christian DA, Discher DE, Hunter CA, Scott P. Leishmania major infection-induced VEGF-A/VEGFR-2 signaling promotes lymphangiogenesis that controls disease. J Immunol. 2016;197(5):1823–1831. pmid:27474074
  40. 40. Bussolino F, Di Renzo MF, Ziche M, Bocchietto E, Olivero M, Naldini L, et al. Hepatocyte growth factor is a potent angiogenic factor which stimulates endothelial cell motility and growth. J Cell Biol. 1992;119(3):629–641. pmid:1383237
  41. 41. Barrientos S, Stojadinovic O, Golinko MS, Brem H, Tomic-Canic M. Growth factors and cytokines in wound healing. Wound Repair Regen. 2008;16(5):585–601. pmid:19128254
  42. 42. Raja , Sivamani K, Garcia MS, Isseroff RR. Wound reepithelialization: modulating keratinocyte migration in wound healing. Front Biosci. 2007;12:2849–2868. pmid:17485264
  43. 43. Rodríguez NE, Wilson ME. Eosinophils and mast cells in leishmaniasis. Immunol Res. 2014;59(1–3):129–141. pmid:24838146
  44. 44. Dufour JH, Dziejman M, Liu MT, Leung JH, Lane TE, Luster AD. IFN-inducible protein 10 (IP-10; CXCL10)-deficient mice reveal a role for IP-10 in effector T cell generation and trafficking. J Immunol. 2002;168:3195–3204. pmid:11907072
  45. 45. Ritter U, Korner H. Divergent expression of inflammatory dermal chemokines in cutaneous leishmaniasis. Parasite Immunol. 2002;24:295–301. pmid:12102714
  46. 46. Charmoy M, Brunner-Agten S, Aebischer D, Auderset F, Launois P, Milon G, et al. Neutrophil-derived CCL3 is essential for the rapid recruitment of dendritic cells to the site of Leishmania major inoculation in resistant mice. PLoS Pathog. 2010;6(2):e1000755. pmid:20140197
  47. 47. Dos Santos PL, de Oliveira FA, Santos ML, Cunha LC, Lino MT, de Oliveira MF, et al. The Severity of Visceral Leishmaniasis Correlates with Elevated Levels of Serum IL-6, IL-27 and sCD14. PLoS Negl Trop Dis. 2016;10(1):e0004375. pmid:26814478
  48. 48. Dayakar A, Chandrasekaran S, Kuchipudi SV, Kalangi SK. Cytokines: Key Determinants of Resistance or Disease Progression in Visceral Leishmaniasis: Opportunities for Novel Diagnostics and Immunotherapy. Front Immunol. 2019; 10:670. eCollection 2019. pmid:31024534
  49. 49. Hurdayal R, Brombacher F. The role of IL-4 and IL-13 in cutaneous Leishmaniasis. Immunol Lett. 2014;161(2):179–83. pmid:24412597
  50. 50. Amato VS, de Andrade HF, Duarte MI. Mucosal leishmaniasis: in situ characterization of the host inflammatory response, before and after treatment. Acta Trop. 2003:85:39–49. pmid:12505182