- Short report
- Open Access
Molecular signatures of neutrophil extracellular traps in human visceral leishmaniasis
Parasites & Vectorsvolume 10, Article number: 285 (2017)
Infections with parasites of the Leishmania donovani complex result in clinical outcomes that range from asymptomatic infection to severe and fatal visceral leishmaniasis (VL). Neutrophils are major players of the immune response against Leishmania, but their contribution to distinct states of infection is unknown. Gene expression data suggest the activation of the NETosis pathway during human visceral leishmaniasis. Thus, we conducted an exploratory study to evaluate NET-related molecules in retrospective sera from VL patients, asymptomatic individuals and uninfected endemic controls.
We demonstrate that VL patients and asymptomatic individuals exhibit differential regulation of molecules associated with neutrophil extracellular traps (NET). These differences were observed at the transcriptional level of genes encoding NET-associated proteins; in quantifications of cell free DNA and metalloproteinase 9; and in enzymatic activity of DNAse and elastase. Moreover, multivariate analysis resulted in class-specific signatures, and ROC curves demonstrate the ability of these molecules in discriminating asymptomatic infection from uninfected controls.
Molecules that are associated with NETs are differentially regulated between distinct states of infection with L. infantum, suggesting that NETs might have distinct roles depending on the clinical status of infection. Although unlikely to be exclusive for VL, these signatures can be useful to better characterize asymptomatic infections in endemic regions of this disease.
Visceral leishmaniasis (VL) is a life-threatening disease caused by infections with the protozoan parasites Leishmania donovani and L. infantum. It is estimated that over 20,000 deaths occur annually due to VL, however, most infections remain asymptomatic . The prognosis of asymptomatic infection and its role in transmission are unknown; furthermore, the mechanisms underlying either active VL or asymptomatic infections have not been fully elucidated. An effective immune response should account for a favorable outcome after infection, thus asymptomatic infections can provide important insights towards understanding host responses during infections with Leishmania . Neutropenia is an independent predictor of risk of dying from VL , whereby systemic inflammation plays a major role. These findings indicate that neutrophil activity may be involved in the distinct states of infection; however, the mechanisms linking neutrophils to active VL or asymptomatic infection are unknown.
Neutrophils rapidly migrate and infiltrate into sites of infection and inflammation. They are important sources of anti-microbial effector molecules involved in host defense, some of which cause tissue damage. Microbicidal activity of neutrophils is exerted by proteolytic enzymes, reactive oxygen species and neutrophil extracellular traps (NETs), which exhibit critical roles during infections with a wide range of pathogens. The release of NETs emerged as an important process of host response to extracellular bacteria and fungi [4, 5]; we demonstrated that Leishmania also trigger this response : L. donovani and L. infantum induce release of NETs by human neutrophils, but they are capable of evading killing through expression of lipophosphoglycan or 3′-nucleotidase/nuclease activity, respectively [7, 8]. In spite of this knowledge, the impact of NETosis during distinct outcomes of human infections with Leishmania remains unexplored.
Here, we evaluated gene expression data from blood of VL patients, asymptomatic individuals and uninfected controls; we observed significant differences between these groups in the modulation of genes encoding NET-associated proteins. We also quantified levels of cell-free DNA and metalloproteinase 9; or measured the enzymatic activity of DNAse, elastase and myeloperoxidase in retrospective sera from distinct states of infection with L. infantum. We identified profiles that suggest a differential regulation of NETosis pathway between states of infection with L. infantum. Importantly, these signatures point towards novel perspectives for characterization of asymptomatic infections in endemic areas for VL.
In this study, we used retrospective serum samples collected from patients with uncomplicated VL at University Hospital, UFS, Aracaju-SE, Brazil (n = 35; male to female ratio = 1.9; mean age ± standard deviation (SD) = 31.5 ± 18.3 years), and stored at -80 °C. All samples were collected before treatment of patients. Diagnosis was confirmed by identification of Leishmania amastigotes in Giemsa-stained smears, positive culture on NNN media of bone marrow aspirate or rK39 positive serology. Retrospective serum samples of healthy individuals living in the same areas were also used. Individuals were tested for delayed type hypersensitivity (reflecting cellular immune response) and IgG reactivity (reflecting humoral immune response) to leishmanial antigens as described previously . Only subjects exhibiting positive reactions to both tests were considered to be asymptomatically infected (n = 28; male to female ratio = 1.1; mean age ± SD = 31.5 ± 15.0 years). Controls included individuals who exhibited negative reactions (n = 25; male to female ratio = 0.9; mean age ± SD = 30.6 ± 18.1 years). The groups did not present significant differences with respect to age (tested with ANOVA) or sex distribution (tested with Chi-square test).
Whole blood transcriptional analyses
Blood collection, RNA extraction, hybridization, scanning and raw data pre-processing pipeline were described previously . Differential gene expressions between distinct groups were identified with linear models and moderated t-statistic. Differentially expressed genes (DEGs) were evaluated with GeneGo MetaCore (Thomson Reuters, NY, USA). P-values were adjusted with Benjamini-Hochberg, and DEGs or pathway over-representation were identified by a FDR or q < 0.05. Gene set enrichment analysis (GSEA) was performed in pre-ranked list mode with 1000 permutations and weighted enrichment statistic.
Measurement of extracellular dsDNA
Serum was diluted 1:4 in TE (Tris-HCl 10 mM, pH 8,0; EDTA 1 mM) and quantified with Quant-iT™ PicoGreen® dsDNA Assay kit (Thermo Fisher Scientific, Waltham, USA) as reported previously .
DNAse activity assay
NET-DNA enriched supernatants were obtained as described previously . Serum was diluted 1:10 and further incubated with 7 μg/ml of NET-DNA enriched supernatants and DNAse cofactors (500 μM CaCl2; 5 mM MgCl2) at 35 °C and 5% CO2. After 5 h, DNAse activity was calculated as the percentage of NET-DNA degradation using Quant-iT™ PicoGreen® dsDNA Assay kit (Thermo Fisher Scientific, Waltham, USA). Purified DNAse-I (2 units; Promega, Fitchburg, USA) was used as a positive control of DNAse activity. Untreated NET-DNA supernatants were used as reference for calculations.
Elastase, myeloperoxidase (MPO) and metalloproteinase 9 (MMP-9) assays
Serum was diluted 1:4 in buffer (50 mM HEPES, 100 mM NaCl and 0.01% Triton X-100) and elastase activity was measured with the fluorogenic substrate N-Methoxysuccinyl-Ala-Ala-Pro-Val-7-amido-4-methylcoumarin (Sigma-Aldrich, St. Louis, USA) according to manufacturer’s instructions. Serum was diluted 1:5 in PBS 0.5% Hexadecyltrimethylammonium bromide (5 mg/ml) and peroxidase activity was measured using 3,3′,5,5′-Tetramethylbenzidine (Sigma) after incubation for 1 h at 35 °C. Absorbance was detected on a SpectraMax® Paradigm® microplate reader at 630 nm (Molecular Devices, California, USA). Quantification of serum MMP-9 was performed with Human MMP-9 Quantikine ELISA kit (R&D Systems, Minneapolis, USA) according to instructions provided by the supplier.
One-way ANOVA with Bonferroni’s or Kruskal-Wallis with Dunn’s multiple comparison tests were applied to evaluate differences among independent groups. Spearman’s rank correlation was used to assess non-parametric associations. P-values less than 0.05 were considered significant. Hierarchical clustering was performed with Euclidian distance for metric calculations and the complete linkage method. Receiver Operating Characteristic (ROC) curves were generated with the online analysis platform MetaboAnalyst 3.0 , with PLS-DA as a classification method and t-statistic as a ranking metric.
Results and discussion
We recently demonstrated that, compared to asymptomatic individuals or controls, microarray data expression from whole blood of VL patients was associated with the perturbation of several immune related pathways annotated in the GeneGO Metacore database . Differentially expressed genes (DEGs, q < 0.05) from VL patients were enriched in the “NETosis in SLE” pathway (Additional file 1: Figure S1). To further investigate this finding, we generated a gene-set termed “NET-associated proteins”, based on genes coding for proteins shown to be enriched in NETs released from human neutrophils with mass spectrometry [5, 12]. GSEA with DEGs from VL patients pre-ranked by expression fold-change resulted in significant associations between the transcriptional profile of VL patients and the “NET-associated proteins” gene-set (Fig. 1a, NES = 2.20, P = 0.004 for VL and uninfected controls; NES = 1.80, P = 0.02 for VL and asymptomatic infection). Indeed, the majority of genes coding for proteins associated with NETs were upregulated in VL patients (Fig. 1a, b). However, compared to uninfected controls, VL patients exhibited more DEGs in the core enrichment than compared to asymptomatic individuals (Fig. 1b). Due the inflammatory nature of NETs, it is unlikely that this activity of the NETosis pathway is exclusive for Leishmania infection. Thus, we evaluated publicly available gene expression profiles (GEO database) of patients infected with Mycobacterium tuberculosis (TB), individuals with latent infection with M. tuberculosis (LTB) and uninfected controls . Of interest, we also identified differential regulation of “NET-associated proteins” gene set among the distinct states of infection with M. tuberculosis (Additional file 2: Figure S2; NES = 2.77, P < 0.0001 for TB and uninfected controls; NES = 2.06, P < 0.0001 for TB and LTB). Taken together, those findings suggest that expression profiles might reflect the activity of the NETosis pathway during infections, and might be differentially regulated between active disease and asymptomatic status as demonstrated for other immune compartments such as antibody profiles in tuberculosis  and VL , and transcriptional activity of type I interferon responses in VL .
To investigate further, we screened sera for known components of NETs in a retrospective cohort of subjects presenting with distinct statuses of infection with L. infantum. We quantified levels of cell free DNA and matrix metalloproteinase 9 (MMP-9), and measured the enzymatic activity of DNAse, elastase and myeloperoxidase (MPO). Compared to uninfected controls, levels of DNA were elevated in sera from VL patients and asymptomatic individuals (Fig. 1c, H = 18.43, df = 2, P < 0.0001). The degradation of NET-DNA, reflecting the enzymatic activity of DNAse, was also elevated in sera from VL patients and asymptomatic controls (Fig. 1c, H = 13.76, df = 2, P = 0.001). In contrast, the levels of MMP-9 were reduced in VL patients compared to both asymptomatic and uninfected controls (Fig. 1c, H = 25.15, df = 2, P < 0.0001), which did not differ between these two groups, as reported . Surprisingly, the enzymatic activity of elastase was elevated in asymptomatic individuals, compared to both VL patients and uninfected controls (Fig. 1c, H = 20.12, df = 2, P < 0.0001), but VL and uninfected controls did not differ. Although a tendency for increased activity of MPO was observed for sera of asymptomatic individuals, there was not significant difference between the groups. Correlation analysis revealed a positive association between levels of DNA and elastase activity in sera from VL patients (Fig. 1d, r = 0.5, P = 0.03). Of note, levels of DNA in sera from asymptomatic individuals correlated negatively with activity of DNAse (Fig. 1d, r = -0.6, P = 0.03). Moreover, activity of DNAse and MPO were negatively correlated for VL patients (r = -0.5, P = 0.04). The lack of correlations between NET-associated molecules should be considered suggestive that the activity of those molecules might have a different origin than neutrophils. However, we speculate that the constitution and/or activity of the NETosis pathway could depend on the state of infection, as different stimulus can impact NET protein composition . Whether the differential regulation between those molecules is indeed resulting from the release of NETs still needs detailed investigation. In this context, recent work demonstrated that VL patients present elevated levels of plasma arginase, MPO and elastase, which return to baseline levels after treatment, however asymptomatic infections were not analyzed . We evaluated chronically infected individuals that should be carrying intracellular forms of Leishmania, raising questions about of the role of extracellular traps in asymptomatically infected and diseased subjects. One explanation could be that when amastigotes emerge from their niches to invade other cells they become susceptible to the activity of NETs. Recently, we reported an intimate association of extracellular amastigotes with NETs in cutaneous leishmaniasis lesions . The spontaneous release of NETs plays an important role in several inflammatory diseases , giving another perspective to NETs in infections with L. infantum: spontaneous release of NETs during chronic VL could induce and/or maintain a systemic inflammatory response, contributing to the pathology of the disease. In a previous study, neutrophils from VL patients exhibited impaired release of NETs during stimulations in vitro . However, the inflammatory environment to which neutrophils are submitted during the course of an infection has profound effects on neutrophils functions , and needs to be accounted for. Our study was limited to retrospective sera samples collected previously , but our data provide a proof of principle that the NETosis pathway can be important and should be investigated with more details in VL and asymptomatic infections. While neutrophil counts might also influence the quantification of the molecules evaluated herein, we highlight the fact that despite reduced proportions of granulocytes, VL patients still exhibited high activity of the NETosis pathway at the transcriptional level .
Next, we conducted principal components analysis (PCA) to identify profiles of expression data of genes coding for NET-associated proteins or of serum indicators of NETosis from the groups of study. We observed a clear clustering pattern in expression data from VL patients and, as expected, those of asymptomatic individuals and uninfected controls were not distinguishable (Fig. 2a). Unfortunately, we were unable to perform all the serum assays for each of the subjects enrolled in the study, however, for a subset of individuals the analysis of principal components with indicators of NETosis revealed a successful segregation of subjects from distinct states of infection (Fig. 2b, left panel). Therefore, to gain power in the analysis by increasing the number of individuals from each group, we excluded measurements from MMP-9 and again found a consistent clustering pattern among subjects of specific groups (Fig. 2b, right panel). We next performed PCA based on peer-comparisons between the groups of study that confirmed distinct profiles of serum indicators of NETosis (Fig. 2c). To characterize the contribution of those molecules to distinct profiles between asymptomatic infection and active VL, we performed an unsupervised hierarchical clustering that retrieved biosignatures with the ability to discriminate between states of infection (Fig. 2d). Of note, asymptomatic individuals present an overall increase in levels of serum DNA and activity of MPO and elastase, while VL patients exhibit a strong signature of the activity of DNAse (Fig. 2d). Clustering of subjects in distinct groups, in addition to their original classification, suggests that the molecular signatures identified herein could also help in the design of prognosis scores, leading to a faster identification of individuals prone to develop VL and facilitate planning of treatments. For that purpose, longitudinal studies following asymptomatic carriers during a substantial period of time after infection might be useful to understand the relationship between molecular signatures of NETs and turnover of clinical status.
More striking was the observation that signatures composed by NET-related molecules distinguish asymptomatic individuals from uninfected controls (Fig. 2e). Most biomarkers of infection with Leishmania and diagnostic assays have been evaluated for their capacity to detect clinical disease, whereas the absence of robust markers for asymptomatic infection remains a problem . Identification of these individuals is time consuming, expensive, and the availability of appropriate leishmanial antigens to performed DTH tests is often limited. The investigation of asymptomatic infections remains largely neglected, but can provide significant information to understand the immunobiology of VL . NETosis is not an exclusive process for infections with Leishmania, and our data do not point for biomarkers of active VL. Moreover, an ideal biomarker exhibits unique associations with a particular phenotype. However, these molecular signatures might contribute to a better characterization of asymptomatic infection. This can be particularly useful in settings where only one method of characterization is available in endemic regions for VL. To support this concept, ROC curves were generated to verify the predictive value of these molecules in distinguishing asymptomatic L. infantum infection from uninfected individuals (Fig. 2f). By comparing models with varying numbers of features, the model with the 4 features exhibit the largest area under the curve (AUC - 0.948) and highest predictive accuracy value of 90.6%. This multivariate approach might reflect better the activity of the NETosis pathway in comparison to univariate analyses. Nevertheless, further investigation with larger numbers of individuals, and longitudinal designs will be necessary to validate these findings.
In summary, our study prompts a more comprehensive evaluation of NETs during distinct states of infection with Leishmania. The signatures identified herein raise the hypothesis that NETs can play distinct roles depending on the progression of the infection. Of importance, the predictive value of these signatures point for novel directions to better characterize asymptomatic infections in endemic areas for VL.
Analysis of variance
Differentially expressed genes
Delayed type hypersensitivity
False discovery rate
Gene set enrichment analysis
Neutrophil extracellular traps
Principal components analysis
Systemic lupus erythematosus
Badaró R, Jones TC, Lorenco R, Cerf BJ, Sampaio D, Carvalho EM, et al. A prospective study of visceral leishmaniasis in an endemic area of Brazil. J Infect Dis. 1986;154:639–49.
Singh OP, Hasker E, Sacks D, Boelaert M, Sundar S. Asymptomatic Leishmania infection: a new challenge for Leishmania control. Clin Infect Dis. 2014;58:1424–9.
Sampaio MJA de Q, Cavalcanti NV, JGB A, MJCF F, Correia JB. Risk factors for death in children with visceral leishmaniasis. PLoS Negl Trop Dis. 2010;4:e877.
Brinkmann V, Reichard U, Goosmann C, Fauler B, Uhlemann Y, Weiss DS, et al. Neutrophil extracellular traps kill bacteria. Science. 2004;303:1532–5.
Urban CF, Ermert D, Schmid M, Abu-Abed U, Goosmann C, Nacken W, et al. Neutrophil extracellular traps contain calprotectin, a cytosolic protein complex involved in host defense against Candida albicans. PLoS Pathog. 2009;5:e1000639.
Guimarães-Costa AB, Nascimento MTC, Froment GS, Soares RPP, Morgado FN, Conceição-Silva F, et al. Leishmania amazonensis promastigotes induce and are killed by neutrophil extracellular traps. Proc Natl Acad Sci U S A. 2009;106:6748–53.
Gabriel C, McMaster WR, Girard D, Descoteaux A. Leishmania donovani promastigotes evade the antimicrobial activity of neutrophil extracellular traps. J Immunol. 2010;185:4319–27.
Guimarães-Costa AB, DeSouza-Vieira TS, Paletta-Silva R, Freitas-Mesquita AL, Meyer-Fernandes JR, Saraiva EM. 3′-nucleotidase/nuclease activity allows Leishmania parasites to escape killing by neutrophil extracellular traps. Infect Immun. 2014;82:1732–40.
Gardinassi LG, Dotz V, Ederveen AH, de Almeida RP, CHN C, Costa DL, et al. Clinical severity of visceral leishmaniasis is associated with changes in immunoglobulin G Fc N-glycosylation. MBio. 2014;5:e01844–14.
Gardinassi LG, Garcia GR, Costa CHN, Silva VC, Santos IKF de M. Blood transcriptional profiling reveals immunological signatures of distinct states of infection of humans with Leishmania infantum. PLoS Negl Trop Dis. 2016;10:e0005123.
Xia J, Sinelnikov IV, Han B, Wishart DS. MetaboAnalyst 3.0 - making metabolomics more meaningful. Nucleic Acids Res. 2015;43:W251–7.
O’Donoghue AJ, Jin Y, Knudsen GM, Perera NC, Jenne DE, Murphy JE, et al. Global substrate profiling of proteases in human neutrophil extracellular traps reveals consensus motif predominantly contributed by elastase. PLoS One. 2013;8:e75141.
Berry MPR, Graham CM, McNab FW, Xu Z, Bloch SAA, Oni T, et al. An interferon-inducible neutrophil-driven blood transcriptional signature in human tuberculosis. Nature. 2010;466:973–7.
Lu LL, Chung AW, Rosebrock TR, Ghebremichael M, Yu WH, Grace PS, et al. A functional role for antibodies in tuberculosis. Cell. 2016;167:433–443.e14.
de Oliveira FA, Vanessa Oliveira Silva C, Damascena NP, Passos RO, Duthie MS, Guderian JA, et al. High levels of soluble CD40 ligand and matrix metalloproteinase-9 in serum are associated with favorable clinical evolution in human visceral leishmaniasis. BMC Infect Dis. 2013;13:331.
Khandpur R, Carmona-Rivera C, Vivekanandan-Giri A, Gizinski A, Yalavarthi S, Knight JS, et al. NETs are a source of citrullinated autoantigens and stimulate inflammatory responses in rheumatoid arthritis. Sci Transl Med. 2013;5:178ra40.
Yizengaw E, Getahun M, Tajebe F, Cruz Cervera E, Adem E, Mesfin G, et al. Visceral leishmaniasis patients display altered composition and maturity of neutrophils as well as impaired neutrophil effector functions. Front Immunol. 2016;7:517.
Morgado FN, Nascimento MTC, Saraiva EM, de Oliveira-Ribeiro C, Madeira M de F, da Costa-Santos M, et al. Are neutrophil extracellular traps playing a role in the parasite control in active American tegumentary leishmaniasis lesions? PLoS One. 2015;10:e0133063.
Leoratti FM de S, Trevelin SC, Cunha FQ, Rocha BC, PAC C, Gravina HD, et al. Neutrophil paralysis in Plasmodium vivax malaria. PLoS Negl Trop Dis. 2012;6:e1710.
The authors are grateful to the participants of the study.
This work was supported by Fundação Carlos Chagas Filho de Amparo à Pesquisa do Estado do Rio de Janeiro (FAPERJ), Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) and Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES). LGG and GRG were supported by scholarships from the Fundação de Amparo à Pesquisa do Estado de São Paulo - FAPESP (2011/23819-0 and 2014/25856-8; 2013/00382-0, respectively).
Availability of data and materials
Probe-level expression data files are accessible at the Gene Expression Omnibus repository under accession number GSE77528 for infection with Leishmania infantum, and GSE19491 for infection with Mycobacterium tuberculosis.
LGG, IKFMS and EMS conceived the study; LGG, TSDV, NOS, GRG, VMB and RNSC performed experiments and data analysis; RPA conducted clinical assessment and follow up; IKFMS and EMS supervised the study; LGG, TSDV, VMB, RPA, IKFMS and EMS wrote the paper. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
This study was approved by the Federal University of Sergipe (CAAE: 4,587,312.2.0000.0058) and written informed consent was obtained from all participants.
“NETosis in SLE” pathway map generated with GeneGO Metacore. Colored asterisks indicate differentially expressed genes of VL patients compared to uninfected controls, where red indicates upregulated genes and blue indicates downregulated genes. See MetaCore website for detailed legend at https://portal.genego.com/legends/MetaCoreQuickReferenceGuide.pdf. (TIFF 2787 kb)
Gene set enrichment analysis (GSEA) showing statistically significant enrichment of differentially expressed genes from patients infected with Mycobacterium tuberculosis in the “NET-associated proteins” gene-set (FDR = false discovery rate). M. tuberculosis patients (TB), uninfected controls (CTRL) and latent M. tuberculosis infection (LTB). (TIFF 785 kb)