- 1WHO Collaborating Centre for Leishmaniasis, Spanish National Center for Microbiology, Instituto de Salud Carlos III, Madrid, Spain
- 2Centro de Investigación Biomédica en Red de Enfermedades Infecciosas (CIBERINFEC), Instituto de Salud Carlos III, Madrid, Spain
- 3Chronic Disease Program (UFIEC), Instituto de Salud Carlos III, Madrid, Spain
- 4Centro de Investigación Biomédica en Red de Fragilidad y Envejecimiento Saludable (CIBERFES), Instituto de Salud Carlos III, Madrid, Spain
Introduction: Visceral leishmaniasis (VL) occurs more frequently in immunosuppressed individuals, especially those undergoing immunosuppressive drug therapy for an autoimmune disease. In those receiving TNF antagonist therapy (anti-TNF), the course of VL is more severe and the response to traditional leishmanicidal treatments, such as antimonials (Sb), is often reduced. This effect of anti-TNF treatment is observed in our immunosuppressed-mouse model of VL. In this model, we compared anti-TNF immunosuppression with no immunosuppression before and after VL treatment with Sb.
Methods: Serum-derived extracellular vesicles (EVs) were analyzed through label-free quantitative proteomics to identify proteins involved in both VL severity and the impact of anti-TNF immunosuppression on treatment outcome.
Results: In total, 223 dysregulated proteins were found in the pre-treatment groups, the majority of which, such as vitronectin, haemopexin or caveolin-1, were downregulated in the anti-TNF samples. In contrast, 173 proteins were identified in the Sb-treatment groups, most of which were found enriched in the anti-TNF plus treatment samples (anti-TNF+Sb) including fibronectin, transferrin, vitronectin and dipeptidyl peptidase-4. These differentially-expressed proteins were associated with pathways related to the immune system, liver regeneration, and ion transport.
Conclusion: Our findings have useful implications for the clinical management of VL patients under anti-TNF immunosuppression.
1 Introduction
Leishmaniasis is a neglected vector-borne tropical disease, whose most severe form –visceral leishmaniasis (VL)– is fatal if left untreated (1). Recently, an increasing number of cases of VL have been reported among individuals receiving immunosuppressive therapy to treat autoimmune diseases such as psoriasis or rheumatoid arthritis (2). Anti-TNF therapies have improved the quality of life of these patients by reducing the inflammatory effects of the TNF cytokine. However, by reducing the protective capacity of the immune system, the patient becomes more susceptible to opportunistic pathogens such as those causing VL. This occurs because TNF plays an important role in the activation and differentiation of immune cells such as macrophages (3, 4). In addition, immunosuppression, particularly with anti-TNF therapies, is responsible for approximately 50% of VL cases among immunocompromised individuals, which compromise the efficacy of VL treatments and, consequently, increase the risk of VL relapse (5, 6). Thus the challenges faced by these immunosuppressed VL patients are exacerbated, especially given the currently limited availability of antileishmania chemotherapeutic agents.
Considering that immunosuppression is the main individual risk factor for a person to develop VL, there is a real need to understand the mechanisms driving the increased disease severity and reduced efficacy of leishmanicidal treatments observed in these immunocompromised patients. Few studies have addressed this issue as most of the pertinent literature has focused on describing the severe symptoms and complications experienced by these patients. We early reported in a mouse model of VL, that the administration of anti-TNF antibodies modulates the natural course of Leishmania infantum infection, drastically increasing parasite load in the liver and suppressing the Th1-protective effect of CD4+ and CD8+ T cells when compared to what happens in immunocompetent mice (7). Factors found to contribute to parasite persistence and disease severity were a lack of the production of pro-inflammatory cytokines such as IFN-γ, along with an increase in IL-10-producing regulatory T cells (Treg) and in the cell exhaustion marker PD-1 (8–10). We also described that the lowest response to leishmanicidal treatment with pentavalent antimonials (Sb) in anti-TNF mice, with respect to the immunocompetent group, was mainly the consequence of reduced activation of the immune system’s defence capacity (11, 12). Under this specific immunosuppression, Sb therapy failed to promote the recruitment of dendritic cell populations, and, instead, increased the frequencies of IL-10 and PD-1-producing B cells. This diminished antigen presentation and the reduced activation of T lymphocytes could explain the greater risk of relapse observed in clinical cases of VL receiving anti-TNF immunosuppression therapy.
As these cell mechanisms affected by anti-TNF immunosuppression become clearer, more research is needed to improve the clinical management of these patients in endemic areas of leishmaniasis. Clinically useful biomarkers are urgently needed to monitor the efficacy of leishmanicidal treatments in anti-TNF immunosuppressed patients (13). In recent years, extracellular vesicles (EVs) have revolutionised the field of biomarker research due to their roles in host-parasite communication, the immune response and drug resistance, among others (14–17). EVs are small particles enclosed in a lipid bilayer membrane (18) that contain a wide range of biomolecules, such as proteins, that may serve as biomarkers (19, 20), particularly those derived from readily-available samples like blood plasma or serum.
To gain insight into the underlying mechanisms of VL and its treatment, the primary aim of this study was to compare the protein contents of EVs using a proteomic approach, in serum samples from Leishmania-infected non-immunosuppressed and anti-TNF- immunosuppressed mice before and after Sb treatment. Subsequently, the most prominent proteins were evaluated by ELISA in serum to assess whether the candidate proteins identified in EVs could also be detectable and relevant in an easily accessible in clinical matrix for clinical purposes. Our ultimate goal was to identify proteins that could be associated with the prognosis and treatment response to VL under anti-TNF immunosuppression.
2 Materials and methods
2.1 Sample collection
Serum samples were obtained from a VL model in anti-TNF immunosuppressed BALB/c mice intravenously infected with 1 x107 L. infantum promastigotes (12). Mice were assigned to the groups (n=6 each): control, PBS 1X (Gibco, USA); anti-TNF, 20 mg/kg (Leinco Technologies, USA); control+Sb, PBS + Glucantime® (Sanofi, France) and anti-TNF+Sb, anti-TNF + Glucantime® (Figure 1). The control and anti-TNF groups were administered PBS or anti-TNF twice per week and blood samples collected after six weeks of infection. The control+Sb and anti-TNF+Sb groups were administered PBS or anti-TNF twice per week, and Sb was given in the last 21 days of the experiment prior to blood collection after nine weeks of infection. Serum was isolated from blood samples (200 µL) collected from the submaxillary vein by centrifugation at 17,000 g for 15 min and kept at -20°C until use.

Figure 1. Schematic representation of the experimental design. BALB/c mice were randomly divided into two groups and received intraperitoneal (i.p.) administrations of either PBS (control group) or anti-TNF at 20 mg/kg twice weekly. These regimens were maintained throughout the duration of the experiment. One week after the initiation of immunosuppressive treatment (week 0), mice were intravenously infected with L. infantum promastigotes. At six weeks post-infection (W6), blood samples were collected for serum isolation prior to VL treatment. At this point, two additional groups of mice –previously immunosuppressed and infected- began a 21-day course of Glucantime treatment (20 mg/kg 7day, i.p.), forming the control+Sb and anti-TNF+Sb groups. A second blood collection was performed at the end of the treatment period (W9).
BALB/c mice used in this study were specific pathogen-free (SPF) and housed under controlled environmental conditions in individually ventilated cages. Food and water were provided ad libitum, both having been previously autoclaved to ensure sterility. All procedures were approved by the Committee on Ethics and Animal Welfare of the Instituto de Salud Carlos III (CBA 04_2018, PROEX 072/18) and animals handled according to Spanish legislation for the protection of animals used for scientific purposes (Royal Decree 53/203, law 32/2007).
2.2 Extracellular vesicle isolation
For EV isolation, a homogeneous 1 mL pool of serum was prepared for each study group. These samples were centrifuged at 300 g for 10 min and the supernatant diluted with an equal volume of filtered PBS to reduce viscosity. Next, the samples were centrifuged at 2,000 g for 30 min and at 12,000 g for another 30 min to pellet larger vesicles. EVs were than isolated from this final supernatant by size exclusion chromatography (SEC) followed by ultracentrifugation (UC) as previously described (21). Briefly, serum samples were first overlaid on the 70nm/qEV size exclusion column (Izon Science, New Zealand) and the flow-through was collected in 500 µL fractions. Next, fractions 7 to 9 were pooled and subjected to two ultracentrifugation steps at 100,000 g and 4°C for 2 h and 30 min, in a Beckman Coulter Optima XPN-100 ultracentrifuge with a swinging-bucket rotor (SW60Ti; Beckman Coulter Inc, USA). Finally, the pellet was resuspended in a final volume of 400 µL filtered PBS.
2.3 EV sample characterization
The protein contents of the serum-derived EVs were measured with the BCA Protein Assay Kit (Thermo Fisher Scientific, USA) following the manufacturer’s instructions. Briefly, 10-µL aliquots were used in the reaction mixture followed by incubation for 30 min at 37°C and absorbance measurement at 562 nm. Protein concentration was calculated based on a standard curve prepared with bovine serum albumin (BSA) and fitted using a four-parameter logistic model.
Particle size and concentration measurements of the isolated EVs were made in each sample using a NanoSight NS300 instrument (Malvern, Worcestershire, UK) and analysed using NTA 3.2 software. Samples were diluted 1:50 in filtered PBS and run with default settings. Mean particle size and concentration were calculated based on three independent records (60-s each) obtained for each sample.
Finally, EV samples were examined by transmission electron microscopy (TEM) through negative staining on glow-discharged carbon-coated copper grids. Preparations were fixed in 2% paraformaldehyde for 5 min and then washed two times with MiliQ water and negatively stained with 2% aqueous uranyl acetate for 1 min. Particles were visualized using a FEI Tecnai 12 electron microscope equipped with a LaB6 filament operating at 120 kV. Images were captured using a FEI Ceta digital camera at a nominal magnification of 28,000x (Microscopy Service of the National Centre for Microbiology, Instituto de Salud Carlos III).
2.4 Label free quantification proteomics
To assess protein dysregulation associated with anti-TNF immunosuppression and Sb treatment, we used a label free quantification (LFQ) proteomics approach. To this end, 2.8 µg of each EV sample were resuspended in lysis buffer (RIPA, Sigma-Aldrich, USA), made up to a final volume of 150 µL, and incubated 5 min on ice and 5 min at 95°C (x5) to ensure complete lysis of the EVs. The protein extracts were then reduced in 1:10 diluted 100 mM tris (2-carboxyethyl) phosphine (TCEP, Sigma-Aldrich) for 45 min at 37°C and 1,000 rpm, and alkylated with 1:10 diluted 400 mM chloroacetamide (Sigma-Aldrich) for 30 min at room temperature (RT) and 1,000 rpm in the dark. For protein anchoring, the reduced extracts were incubated with 100 µL of SeraMag magnetic beads mix (2.5 µL of hydrophilic beads-2.5 µL of hydrophobic beads per sample, Cytiva, UK) and 200 μL of 100% acetonitrile (ACN) for 35 min at RT and 1,000 rpm. Next, the magnetic beads were washed twice with ethanol 70% and once with ACN 100%, and, finally, beads were incubated overnight at 37°C with 0.2 μg of porcine trypsin (Thermo Fisher Scientific) in 100 μL of 50 mM ammonium bicarbonate, pH 8.0. To recover peptides, the samples were sonicated twice and supernatants containing the digested proteins collected, dried under vacuum, and stored at -80°C until analysis in an Orbitrap Exploris 480 mass spectrometer equipped with a FAIMS pro Duo interface (22).
2.5 LC-MS/MS analysis
Peptide samples were resuspended in 11 μL of 0.1% formic acid (FA) H2O and 2 μL of each sample (800 ng) were injected four times using the Vanquish Neo UHPLC System (Thermo Fisher Scientific). For liquid chromatography (LC), samples were loaded onto a precolumn PepMap 100 C18 3 µm, 75 µm × 2 cm Nanoviper Trap 1200BA (Thermo Fisher Scientific) and eluted in an Easy-Spray PepMap RSLC C18 2 µm, 75 µm × 50 cm (Thermo Fisher Scientific) heated to 50°C. The mobile phase flow rate was 300 nL/min using 0.1% FA H2O (buffer A) and 0.1% FA in 80% ACN (buffer B). The 2-h elution gradient was: 2% buffer B for 5 min, 2-20% buffer B for 100 min, 20-42% buffer B for 10 min, 42–95% buffer B for 1 min, and 95% buffer B for 10 min.
For tandem mass spectrometry (MS/MS) analysis, 2300 V of liquid junction voltage and 280°C capillary temperature were used for ionization. The full scan method employed a m/z 350–1400 mass selection, an Orbitrap resolution of 60,000 (at m/z 200), an automatic gain control (AGC) value of 300%, and a maximum injection time (IT) of 25 ms. For MS2, the 12 most intense precursor ions were selected for fragmentation with a normalized collision energy of 32. MS2 scans were acquired with a 100 m/z first mass, an AGC target of 200%, a resolution of 15,000 (at m/z 200), an intensity threshold of 5×104, an isolation window of 1.3 m/z, and a maximum IT of 22 ms. Charge state screening was enabled to reject unassigned, singly charged, and greater than or equal to seven protonated ions. A dynamic exclusion time of 30 s was used to discriminate against previously selected ions. For FAIMS, a gas flow of 4 L/min and CVs = -45 V and -60 V were used.
2.6 MS data analysis and Statistical analysis
MS data were analyzed with MaxQuant (version 2.1.3) using standardized workflows. Mass spectra *.raw files were searched against the Uniprot UP000000589_10090.fasta Mus musculus (mouse), accessed in October 2023 (17,114 protein entries) through standard procedures. Trypsin/P was specified as cleavage enzyme, and precursor and reporter mass tolerances were set to 4.5 ppm and 0.003 Da, respectively, allowing 2 missed cleavages. Carbamidomethylation of cysteines was set as a fixed modification, and methionine oxidation, N-terminal acetylation, and Ser, Thr, and Tyr phosphorylation were set as variable modifications. Label-free quantification and Fast LFQ were selected, with a LFQ minimum number of neighbours of 3, and a LFQ average number of neighbours of 6. Unique and razor peptides were considered for quantification. Minimal peptide length and maximal peptide mass were fixed to 7 amino acids and 4600 Da, respectively. Identified peptides were filtered by their precursor intensity fraction with a false discovery rate (FDR) threshold of 0.01. Proteins identified with at least one unique peptide and an ion score above 99% were considered for evaluation, whereas proteins identified as potential contaminants were excluded from the analysis. Protein sequence coverage was estimated as the percentage of matching amino acids from the identified peptides having a confidence level greater than or equal to 95% in the specific proteins divided by the total number of amino acids in the sequence.
Raw proteomics data obtained with the Orbitrap Exploris 480 mass spectrometer equipped with FAINS pro DUO interface were deposited to the ProteomeXchange Consortium via the PRIDE partner repository with the dataset identifier PXD060935.
Sample loading normalization was performed with R Studio (version 4.1.1) according to established protocols (https://github.com/pwilmart, accessed on 2 November 2022), using the “tidyverse”, “psych”, “gridExtra”, “scales”, and “ggplot2” packages (version 4.1.1). Finally, statistical analysis was performed using an empirical Bayes-moderated t-statistics method in R Studio (version 4.1.1) using the packages “limma”, “dplyr”, “tidyverse”, “ggplot2”, and “rstatix”, according to previously described procedures (22–24). Only proteins identified in at least 60% of samples analyzed in each comparison were considered for the analysis, and missing values imputed by random draws from a gaussian using the “imputeLCMD” R package.
2.7 Biological analysis
Functional annotations of the identified proteins in the EV preparations were obtained using the Database for Annotation, Visualization and Integrated Discovery (DAVID 2021) (25, 26). A gene ontology (GO) enrichment analysis was performed for cellular component and biological processes. Pathways with a p-value ≤ 0.05 and FDR ≤ 0.01 were taken into consideration. Results were created using Graphpad Prism software version 9.0 (GraphPad Software, USA), Functional Enrichment Analysis tool (FunRich) (27), and Flaski tool box (version 3.16.22) (28).
2.8 Enzyme-linked immunosorbent assay on serum
100 µL of both individual serum samples and standard curve solutions were added to pre-coated plates specific for each protein. The proteins assayed were Cav1 (serum dilution 1:2; EM0904, Fine Test), Dpp4 (1:100; orb391051, Biorbyt), Fn (1:500; EM0079, Fine Test), Hmgb-1 (1:10; orb1807923, Biorbyt), Hpx (1:100,000; EM2003, Fine Test), Tf (1:1000; EM1426, Fine Test) and Vtn (1:100; orb565383, Biorbyt). Dilutions of the serum samples were optimized according to the manufacturer’s instructions.
In brief, after an incubation of 90 min at 37°C, the plates were washed and biotin-conjugated detection antibody added for 60 min at 37°C. Next, streptavidin was added and the plates incubated for 30 to 45 min, depending on the specific protein. Following an additional wash, the plates were incubated with 3,3’,5,5’-tetramethylbenzidine (TMB) substrate solution and the reaction stopped with the stop solution provided in each ELISA kit. Absorbance was measured at 450 nm in a Multiskan FC microplate reader (Thermo Fisher Scientific, USA). Data analysis was performed using a four-parameter logistic curve implemented in GraphPad Prism version 9.0. Normally distributed data were analysed via ANOVA, followed by a Tukey’s post hoc test for multiple comparison. Significance was set up at p ≤ 0.05.
3 Results
3.1 Characterization and protein profiles of serum-derived EVs
After purification of EVs by SEC and UC, samples were first characterized in terms of particle size and concentration by NTA, with a distribution profile comparable across the groups despite differences in concentrations (Figure 2A). In all groups, average EV size was in the range 100 to 200 nm, representing 73.80% (control), 67.91% (anti-TNF), 60.60% (control+Sb), and 72.13% (anti-TNF+Sb) of the total number of particles isolated (Figure 2B). Additionally, mean particle size was similar among groups, with no significant differences between them: 175.5 ± 1.9 nm (control), 172.9 ± 3.8 nm (control+Sb), and 173.3 ± 5.1 nm (anti-TNF+Sb). EVs recovered from the anti-TNF group showed the lowest mean size (137.3 ± 3 nm). In addition, NTA analysis confirmed an efficient isolation process with a particle concentration above 109 in all extracts: 7.06 x 109 ± 2.36 x 108 (control), 8.29 x 109 ± 4.55 x 108 (anti-TNF), 1.34 x 109 ± 1.34 x 108 (control+Sb) and 5.92 x 109 ± 2.96 x 108 (anti-TNF+Sb) particles/mL. Next, the presence of EVs in all the preparations was confirmed by TEM (Figure 2C). In all images, it was possible to observe lipid bilayer-enveloped structures with a characteristic cup-shaped appearance, with maximum diameters in the range 100 to 200 nm.

Figure 2. Characterization of EVs from control and anti-TNF samples before and after Sb treatment. (A) NTA was used to determine particle sizes and concentrations to identify the total number of isolated EVs. Although differences in concentrations were found, the distribution profile was similar across the groups. (B) Particle size distribution profile shown as the percentage of total particles detected within each sample. (C) EV morphology and size were determined by TEM. Scale bar set at 100 nm.
Next, a label-free proteomics analysis served to identify and quantify a total of 223 and 173 proteins in samples before (control and anti-TNF) and after (control+Sb and anti-TNF+Sb) antimonial treatment, respectively (Supplementary Tables 1, 2). Before data analysis, we performed a principal component analysis (PCA) to ensure the distribution of samples and the reproducibility of replicates (Figure 3A). We observed different and distinguishable clusters between the four groups, with pre-treatment samples clustering separately from post-treatment samples, thus allowing us to compare the protein contents of EVs between samples. Further analysis to map these proteins against the Mus musculus genome as background was performed according to Gene Ontology (GO) terms, using the Database for Annotation, Visualization and Integrated Discovery (DAVID) and Uniprot database. An initial Cellular Component (CC) analysis was performed to classify proteins within the GO terms related to EVs and plasma such as Extracellular exosome (GO: 0070062), Extracellular region (GO: 0005576), Extracellular space (GO: 0005615), and Plasma membrane (GO: 0005886) (Figure 3B). The percentages of all these GO terms were similar among all samples, the term Extracellular space being the one with the greater proportion of proteins mapped (54.94% to 60.51%). 40.17 to 51.15% of these proteins were also associated with the term Plasma membrane and 18.78 to 27.41% to Extracellular region. For the GO term Extracellular exosome (14.43 to 17.71%), 17 proteins were found in all the samples (p= 1.73E-14, FDR= 8.44E-13 and p= 3.70E-16, FDR= 2.42E-14 when comparing control and anti-TNF and control+Sb and anti-TNF+Sb, respectively) (Figure 3C).

Figure 3. Proteomic analysis of anti-TNF and control EVs before and after antimonial treatment. (A) Principal component analysis (PCA) using the Flasky tool of four replicates of proteomic data obtained in each group. (B) Using the DAVID database, proteins identified in all samples were classified according to the Cellular Component term related to EVs and plasma. The FunRich tool was used to represent the results. (C) Summary of the proteins detected in all the preparations mapped to the GO term Extracellular exosome.
According to MISEV guidelines (18), an evaluation of the protein content was performed to confirm that all samples were enriched in EVs (Supplementary Table 3). In line with recommendations, at least one protein from each of the three major categories was present in the proteomics analysis of all samples. Identification of category 1 and 2 proteins indicates the presence of EVs by characterizing their lipid-bilayer membrane. Within category 1, we found tetraspanins CD82 and CD9 (Uniprot: P40237 and P40240), CD47 multi-pass membrane protein (Uniprot: Q61735), and many integrins, alpha and beta among others. A total of 16 proteins from category 1 were found in control and anti-TNF samples, while 15 were found in control+Sb and anti-TNF+Sb samples. In category 2, the proteins identified (21 in each of the four groups) were caveolin-1 (Uniprot: P49817), heat shock protein Hsp70 (Uniprot: P63017), together with annexins, tubulins, and actins. Finally, in category 3, which describes non-EV structures isolated with EVs, we identified albumin (Uniprot: P07724), apolipoproteins and immunoglobulins, totalling 16 in control and anti-TNF samples, and 12 in samples after Sb treatment.
3.2 Immunosuppression with anti-TNF leads to the reduced expression in EVs of proteins involved in biological processes associated with the response to L. infantum infection
To explore the proteins involved in the increased severity of L. infantum infection under pharmacological immunosuppression with anti-TNF, as observed in the highest parasite loads in the liver of this group (7, 12), protein EV contents were compared between immunosuppressed infected and non-immunosuppressed infected mice (anti-TNF and control, respectively) (Figure 4; Supplementary Table 1). Of the 223 common proteins identified in the proteomics analysis, 76 showed a fold change (FC) ratio ≥ 1.5 or ≤ 0.67 and significant expression level differences between groups (p value ≤ 0.05 and FDR ≤ 0.05) (Figure 4). Compared to the control group (no immunosuppression treatment), anti-TNF immunosuppression led to downregulation (FC ≤ 0.67) of 60 proteins, whereas 16 were upregulated (FC ≥ 1.5) (Figure 4A). Next, to identify the biological processes in which these proteins were involved, a further analysis using the DAVID database was conducted to determine the GO term Biological Process (BP) (Figure 4B). Anti-TNF immunosuppression led to increases in many BP related to the immune system such as negative regulation of cell migration (GO:0030336), positive regulation of IL-8 production (GO:0032757), negative regulation of macrophage differentiation (GO:0045650), negative regulation of T cell differentiation (GO:0045581), and positive regulation of IL-10 production (GO:0032733). In addition, other BPs enriched under this condition were negative regulation of DNA biosynthetic process (GO:2000279), negative regulation of intracellular protein transport (GO:0090317), and negative regulation of phagocytosis (GO:0050765). In contrast, the 60 downregulated proteins in EVs in the anti-TNF group were actively involved in, at least, 10 BPs related to the immune system. Among these BPs, we identified reductions in T cell activation (GO:0042110), regulation of T cell mediated immunity (GO:0002709), liver regeneration (GO:0097421), cytokine-mediated signalling pathway (GO:0019221), type II (GO:0060333) and (GO:0060337) IFN-mediated signalling pathways, and leukocyte migration involved in the inflammatory response (GO:0002523).

Figure 4. Differential expression of EV-derived proteins linked to immunosuppression by the TNF antagonist (anti-TNF vs control group). (A) Volcano plot of proteins found modified in the anti-TNF compared to control group. Coloured dots indicate differentially upregulated (red) or downregulated (blue) proteins in anti-TNF (FDR≤ 0.05; indicated by a horizontal dotted line). Vertical dotted lines indicate log2 fold changes (± 1.5) in expression. (B) Main biological processes altered by anti-TNF immunosuppression according to the DAVID and Uniprot databases. The FunRich tool was used for this analysis.
Among the proteins found to be upregulated in the anti-TNF group were the apolipoproteins Apoc1 (P34928), Apoc4 (Q61268) and Apoe (P08226), which are non-EV co-isolated structures, along with the Adipoq (Q60994) protein (Table 1). We also identified other upregulated proteins such as Epb41 (P48193), Vcp (Q01853), Prdx2 (Q61171) or Hmgb1 (P63158). Among the 60-downregulated proteins detected in the anti-TNF samples, we identified Dpp4 (P28843), Cav1 (P49817), Vtn (P29788), Hpx (Q1X72), Ezr (P26040), Serpina1a (P07758), or Anxa6 (P14824), among others (Table 1).
3.3 Anti-TNF therapy dysregulates VL-related proteins after treatment with antimonials
Our next goal was to determine whether any proteins could be responsible for a reduced response to VL treatment under this type of immunosuppression, as anti-TNF altered the immune response towards the Th2-type profile (12). To do this, we compared protein EV-contents in anti-TNF immunosuppressed-infected mice and control animals after VL treatment with pentavalent antimony (anti-TNF+Sb and control+Sb groups, respectively) (Figure 5). This proteomics analysis identified 173 proteins (Supplementary Table 2), of which 43 showed a FC ratio ≥ 1.5 or ≤ 0.67 and significantly different expression levels between the groups (p value ≤ 0.05 and FDR ≤ 0.05). Anti-TNF therapy led to the upregulation of 30 proteins in infected mice after antimonial treatment (Figure 5A). These proteins were associated with increases in the BPs regulation of iron ion transport (GO:0034756), T cell activation (GO:0042110), liver regeneration (GO:0097421), metal ion transport (GO:0030001), DNA repair (GO:0006281), positive regulation of IL-6 production (GO:0032755), and negative regulation of IL-2 production (GO:0032703) (Figure 5B). In contrast, only 13 proteins were downregulated in the anti-TNF+Sb group, involved in T cell migration (GO:2000406), leukocyte migration (GO:0002687), dendritic cell differentiation (GO:0097028), response to toxic substance (GO:0009636), IL-5 mediated signalling pathway (GO:0038043), and apoptotic cell clearance (GO:0043277).

Figure 5. Differential expression of EV-derived proteins linked to anti-TNF immunosuppression after antimonial treatment (anti-TNF+Sb vs control+Sb group). (A) Volcano plot of proteins found modified in the anti-TNF+Sb compared to control+Sb group. Coloured dots indicate the differentially upregulated (red) or downregulated (blue) proteins in anti-TNF+Sb (FDR ≤ 0.05; indicated by a horizontal dotted line). Vertical dotted lines indicate log2 fold changes (± 1.5) in expression. (B) Main biological processes altered by anti-TNF immunosuppression after Sb treatment according to the DAVID and Uniprot databases. The FunRich tool was used for this analysis.
The following proteins (Table 2) were identified among the 30 upregulated proteins in the anti-TNF+Sb samples: Tf (Q921I1), Dpp4 (P28843), Ezr (P26040), Fn1 (P11276), Cfh (P06909), Serpina3k (P07759), and Cd5l (Q9QWK4). In contrast, Ehd4 (Q9EQP2), Itgb3 (O54890), Itga2b (Q9QUMO), and Slc4a1 (P04919), were identified among the most downregulated proteins in the anti-TNF+Sb group.

Table 2. List of proteins found dysregulated in the anti-TNF+Sb group compared to the Control+Sb group.
3.4 Analysis of dysregulated proteins in serum samples
To determine whether dysregulated proteins found in EVs could be detected in clinical samples and therefore be of use to identify situations of anti-TNF immunosuppression during L.infantum infection, we evaluated seven proteins (Fn, Vtn, Tf, Hpx, Cav, Dpp4 and Hmgb1) in mouse serum samples through ELISA.
Serum expression levels of Fn, Vtn and Tf (Figures 6A–C) showed no significant variation across the different groups. Expression levels of Hpx (103 µg/mL) and Cav (1.41 ng/mL) were significantly higher in the serum samples in the anti-TNF group compared to the control group (p= 0.0088 and p= 0.049 respectively) (Figures 6D, E). In addition, there is a tendency to an increase of the protein levels of Dpp4 after Sb treatment, whereas a partial reduction was observed in Hmgb1 in the post-treatment groups, with no statistical changes observed (Figures 6F, G). The differences in serum expression levels observed between groups did not align with the findings from the proteomic analysis of EVs.

Figure 6. Serum concentrations of selected proteins. Serum samples were obtained after blood collection at two time points post-infection, week 6 (control and anti-TNF groups) and week 9 (control+Sb and anti-TNF+Sb groups). Seven proteins were evaluated by sandwich ELISA: (A) Fibronectin, (B) Vitronectin, (C) Transferrin, (D) Hemopexin, (E) Caveolin, (F) Dipeptidyl pptidase-4 and (G) High mobility group box 1. The concentration of each selected protein is shown individually. Also provided are the means and standard errors for each study group. Comparisons were made between groups sharing the same immunosuppression status (control vs. control + Sb and anti-TNF vs. anti-TNF + Sb), as well as between groups before (control vs. anti-TNF) and after (control + Sb vs. anti-TNF + Sb) VL treatment. Significant differences were assessed using an ANOVA followed by a Tukey’s post hoc test and are indicated as *p ≤ 0.05; **p ≤ 0.01.
4 Discussion
Quantitative proteomics analysis of EVs provides useful information on proteins associated with disease and the biological processes in which they are involved. However, while most EV studies have focused on cancer patients (29, 30), data are now starting to emerge on EV proteomic profiles in infectious diseases, such as those caused by parasites (31, 32). In the context of leishmaniasis, the few studies performed to date have mainly investigated how EVs released by Leishmania parasites modulate cell communication and the host immune response in vitro (33, 34). Although proteomics approaches have been employed in EVs recovered from canine and human VL patients (35) (Torres et al. Front Immunol, in review), these studies have not addressed the effects and complexity of this disease in conditions of immunosuppression. Such studies are essential to better understand the changes that occur during VL and its treatment in this subset of patients showing a higher risk of relapse due to anti-TNF therapy. In the present study, we compared the protein contents of plasma-derived EVs from L. infantum-infected mice undergoing, or not, anti-TNF therapy both before and after Sb treatment. Our ultimate goal was to identify proteins linked to disease progression and to antileishmania agent efficacy in conditions of anti-TNF immunosuppression.
Our study reveals that anti-TNF immunosuppression significantly impacts the natural progression of VL. Blocking TNF prevents granuloma formation and leads to parasite dissemination throughout the liver, causing histological lesions and hepatic necrosis (36, 37). Compared to the situation in control animals, one of the biological processes found to be negatively regulated, as reflected by the EV contents of our anti-TNF group, was liver regeneration. Accordingly, vitronectin (Vtn), a glycoprotein derived from hepatocytes and then secreted into the circulation (38), was significantly depleted in these animals. The authors of studies conducted in Vtn-/- animal models have reported that deletion of the Vtn gene gives rise to delayed wound healing and tissue repair processes (39, 40). Another protein found downregulated in this anti-TNF group was haemopexin (Hpx), which is usually expressed in the liver following an inflammatory event. This protein has a high affinity for haem, which not only protects the cell from oxidative stress but also inhibits the growth of infectious agents like Plasmodium falciparum, thereby preventing their pathogenesis (41–43). The fact that Vtn and Hpx were depleted in the EV-contents of our anti-TNF group may explain the higher liver parasite loads found in these immunocompromised animals (9.82 x 104 vs 3.71 x 104 total parasites in liver in the anti-TNF vs control groups; p=0.0099 (12)). However, this was not observed in the serum expression levels of both proteins, as Hpx was more expressed in the anti-TNF group, whereas serum levels of Vtn remained comparable among all groups. Given that serum is not an EV-enriched matrix and may contain co-isolated proteins or degradation products, discrepancies between EV proteomic profiles and ELISA-based quantification are expected, particularly when targeting low-abundance or EV-specific proteins. Moreover, we also found that peroxiredoxin 2 (Prdx2) protein was upregulated in the EVs of the anti-TNF group compared to control animals. Peroxiredoxins (Prdx) are known to be determinants of the survival and virulence of parasites such as Leishmania, by protecting them from the oxidative damage caused by the host immune system (44). Other studies have shown that P. falciparum can import host Prdx to increase its infectivity by boosting its antioxidant defence capacity (45, 46). Addressing the question of whether Leishmania parasites could also sequester Prdx from the host would provide valuable insight, as this mechanism might also be involved in the increased presence of the parasite detected here in the livers of anti-TNF immunosuppressed mice.
In addition, several immune-mediated pathways were found modified in the anti-TNF immunosuppressed mice such as those of negative regulation of T cell differentiation and positive regulation of both IL-10 and IL-8 production. IL-8 and IL-10 cytokines play a critical role in VL progression as they are strongly correlated with parasite persistence (47). Enrichment in the BP IL-8 production in our anti-TNF serum-derived EV samples could be mainly related to the increased cytokine production that takes place during VL progression due to neutrophil infiltration (48). One of the proteins found upregulated in the anti-TNF group was high mobility group box 1 (Hmgb1). This protein is generally secreted in a damage-associated molecular pattern by activated macrophages and neutrophils, thus triggering a pro-inflammatory response to control infection (49, 50). However, in the context of leishmaniasis, the oxidative environment generated during intracellular infection via reactive oxygen species, changes the conformational state of the protein such that it is unable to interact with its receptors (51). This status shifts to an anti-inflammatory function, increasing the production of T regulatory cells (Tregs) and IL-10 (52). The resultant M2 phenotype leads to reduced nitric oxide production (53) and diminishes the immune system response by inhibiting the frequency of effector T cells, which contributes to parasite persistence and activation of the VL-associated pathology (54, 55) under anti-TNF immunosuppression. In addition, the downregulation of caveolin-1 (Cav1) protein further exacerbates this immune deficiency, by reducing lymphocyte activation and weakening the immune response, thereby diminishing the host’s ability to deal with Leishmania infection (56). This situation leads to an increased parasite burden and pathogenesis, as previously observed in our anti-TNF immunosuppressed mice at six weeks post-infection (7).
Following treatment with Sb, fibronectin (Fn) was enriched in the anti-TNF group compared to the control mice. Leishmania promastigotes interact with this protein to impair the activation of parasite-infected macrophages (57). Thus, Fn overexpression could lead to parasite persistence by binding to parasite receptors, promoting macrophage invasion and parasite spreading (57, 58). Another mechanism of Leishmania persistence is the acquisition of iron (59) by fusion of the host transferrin (Tf) receptor to the parasitophorus vacuole (60). Upregulation of this main iron-carrier protein in the post-Sb treatment anti-TNF group could play a role in parasite survival and proliferation. We also found an increase in Vtn, a protein associated with liver regeneration, in line with the partial reduction of parasite load observed after Sb treatment in immunosuppressed mice (12). Alternatively, Vtn can also bind pathogens promoting their internalization and evading a complement response via the integrin complex (61, 62). In effect, the affinity of the Vtn protein for L. donovani promastigotes has been reported to lead to increased circulating levels of parasites in post kala-azar dermal leishmaniasis (PKDL) patients (63, 64). Further, high levels of Vtn have been associated with the regulatory cytokine TGF-β, which plays an important role during both PKDL and VL progression, as it suppresses iNOS and IFN-γ Th1-type responses (64–66). Thus, the poorer immune response in these anti-TNF mice, together with Fn, Tf and Vtn enrichment could explain why our anti-TNF immunosuppressed mice were unable to effectively eliminate parasites despite receiving leishmanicidal treatment. Moreover, dipeptidyl peptidase-4 (Dpp4) upregulation also exacerbates the reduced capacity of anti-TNF immunosuppressed mice to eliminate the infection, as its increased presence has been observed in unresolved cases of leishmaniasis compared to levels in cured patients and controls (67). We also observed downregulation of the band 3 anion protein transport (Slc4a1) and of myosin-9 (Myh9) in our anti-TNF immunosuppressed mice after Sb treatment. These cure-related proteins have been recently described in plasma-derived EVs obtained in immunocompetent patients when compared to levels observed in patients with active VL treated with Ambisome® (Torres et al. Front Immunol, in review). The impaired enrichment of these proteins after Sb treatment supports the idea of a limited response to VL treatment under conditions of anti-TNF immunosuppression and may be related to the higher risk of VL relapse seen in clinical cases.
In summary, our study provides the first integrative evidence that anti-TNF immunosuppressive therapy not only reshapes host immune responses during Leishmania infection but also profoundly alters the extracellular vesicle proteome, revealing a dysregulation of key biological pathways—particularly those involved in liver regeneration and iron metabolism—that may contribute to increased disease severity and reduced responsiveness to antimonial treatment. These findings underscore the need to tailor clinical management strategies for VL in individuals receiving anti-TNF therapy and highlight the potential of EV proteomics to uncover mechanistic insights and identify novel prognostic biomarkers.
Data availability statement
The mass spectrometry proteomics data have been deposited at the ProteomeXchange Consortium via the PRIDE partner repository with the dataset identifier PXD060935.
Ethics statement
The animal study was approved by Committee on Ethics and Animal Welfare of the Instituto de Salud Carlos III. The study was conducted in accordance with the local legislation and institutional requirements.
Author contributions
LB: Conceptualization, Methodology, Writing – review & editing, Data curation, Investigation, Writing – original draft, Formal analysis. AM-C: Writing – original draft, Formal analysis, Data curation, Investigation, Methodology. JS: Formal analysis, Writing – review & editing, Investigation, Writing – original draft, Conceptualization, Data curation. ML-R: Methodology, Formal analysis, Investigation, Writing – original draft. AT: Formal analysis, Writing – original draft, Methodology, Investigation. CS: Methodology, Formal analysis, Investigation, Writing – original draft. RB: Methodology, Writing – review & editing. JM: Funding acquisiton, Writing – review & editing. EC: Conceptualization, Funding acquisition, Writing – review & editing, Writing – original draft.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. This study was funded by the Instituto de Salud Carlos III through ISCIII-AES projects (PI21CIII/00005, PI22/00009, and PI24CIII/00023). JCS was supported by a contract from CIBERINFEC (CB21/13/00018).
Acknowledgments
The authors thank Felix Docando, María Carmen Terrón, and Daniel Luque from the Microscopy Unit (National Centre for Microbiology, Instituto de Salud Carlos III, Majadahonda, Spain) for their help with the TEM imaging.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2025.1634080/full#supplementary-material
References
1. Gradoni L, Lopez-Velez R, and Mokni M. Manual on case management and surveillance of the leishmaniases in the WHO European Region. Geneva, Switzerland: World Health Organisation (2017).
2. Kurizky PS, Marianelli FF, Cesetti MV, Damiani G, Sampaio RNR, Goncalves LMT, et al. A comprehensive systematic review of leishmaniasis in patients undergoing drug-induced immunosuppression for the treatment of dermatological, rheumatological and gastroenterological diseases. Rev Inst Med Trop Sao Paulo. (2020) 62:e28. doi: 10.1590/s1678-9946202062028
3. Bosch-Nicolau P, Ubals M, Salvador F, Sanchez-Montalva A, Aparicio G, Erra A, et al. Leishmaniasis and tumor necrosis factor alpha antagonists in the Mediterranean basin. A switch in clinical expression. PloS Negl Trop Dis. (2019) 13:e0007708. doi: 10.1371/journal.pntd.0007708
4. Neumayr AL, Morizot G, Visser LG, Lockwood DN, Beck BR, Schneider S, et al. Clinical aspects and management of cutaneous leishmaniasis in rheumatoid patients treated with TNF-alpha antagonists. Travel Med Infect Dis. (2013) 11:412–20. doi: 10.1016/j.tmaid.2013.05.003
5. van Griensven J, Carrillo E, López-Vélez R, Lynen L, and Moreno J. Leishmaniasis in immunosuppressed individuals. Clin Microbiol Infection. (2014) 20:286–99. doi: 10.1111/1469-0691.12556
6. Botana L, Ibarra-Meneses AV, Sanchez C, Matia B, San Martin JV, Moreno J, et al. Leishmaniasis: A new method for confirming cure and detecting asymptomatic infection in patients receiving immunosuppressive treatment for autoimmune disease. PloS Negl Trop Dis. (2021) 15:e0009662. doi: 10.1371/journal.pntd.0009662
7. Bernardo L, Solana JC, Romero-Kauss A, Sanchez C, Carrillo E, and Moreno J. Effect of immunosuppressants on the parasite load developed in, and immune response to, visceral leishmaniasis: A comparative study in a mouse model. PloS Negl Trop Dis. (2021) 15:e0009126. doi: 10.1371/journal.pntd.0009126
8. Rodrigues LS, Barreto AS, Bomfim LGS, Gomes MC, Ferreira NLC, da Cruz GS, et al. Multifunctional, TNF-alpha and IFN-gamma-secreting CD4 and CD8 T cells and CD8(High) T cells are associated with the cure of human visceral leishmaniasis. Front Immunol. (2021) 12:773983. doi: 10.3389/fimmu.2021.773983
9. Jafarzadeh A, Kumar S, Bodhale N, Jafarzadeh S, Nemati M, Sharifi I, et al. The expression of PD-1 and its ligands increases in Leishmania infection and its blockade reduces the parasite burden. Cytokine. (2022) 153:155839. doi: 10.1016/j.cyto.2022.155839
10. Yasmin H, Adhikary A, Al-Ahdal MN, Roy S, and Kishore U. Host–pathogen interaction in leishmaniasis: immune response and vaccination strategies. Immuno. (2022) 2:218–54. doi: 10.3390/immuno2010015
11. Murray HW. Clinical and experimental advances in treatment of visceral leishmaniasis. Antimicrob Agents Chemother. (2001) 45:2185–97. doi: 10.1128/AAC.45.8.2185-2197.2001
12. Bernardo L, Solana JC, Sanchez C, Torres A, Reyes-Cruz EY, Carrillo E, et al. Immunosuppressants alter the immune response associated with Glucantime((R)) treatment for Leishmania infantum infection in a mouse model. Front Immunol. (2023) 14:1285943. doi: 10.3389/fimmu.2023.1285943
13. Carrillo E and Moreno J. Editorial: biomarkers in leishmaniasis. Front Cell Infect Microbiol. (2019) 9:388. doi: 10.3389/fcimb.2019.00388
14. Douanne N, Dong G, Douanne M, Olivier M, and Fernandez-Prada C. Unravelling the proteomic signature of extracellular vesicles released by drug-resistant Leishmania infantum parasites. PloS Negl Trop Dis. (2020) 14:e0008439. doi: 10.1371/journal.pntd.0008439
15. Douanne N, Dong G, Amin A, Bernardo L, Blanchette M, Langlais D, et al. Leishmania parasites exchange drug-resistance genes through extracellular vesicles. Cell Rep. (2022) 40:111121. doi: 10.1016/j.celrep.2022.111121
16. Wu L, Zhou L, An J, Shao X, Zhang H, Wang C, et al. Comprehensive profiling of extracellular vesicles in uveitis and scleritis enables biomarker discovery and mechanism exploration. J Transl Med. (2023) 21:388. doi: 10.1186/s12967-023-04228-x
17. Marcilla A, Martin-Jaular L, Trelis M, de Menezes-Neto A, Osuna A, Bernal D, et al. Extracellular vesicles in parasitic diseases. J Extracell Vesicles. (2014) 3:25040. doi: 10.3402/jev.v3.25040
18. Welsh JA, Goberdhan DCI, O'Driscoll L, Buzas EI, Blenkiron C, Bussolati B, et al. Minimal information for studies of extracellular vesicles (MISEV2023): From basic to advanced approaches. J Extracell Vesicles. (2024) 13:e12404. doi: 10.1002/jev2.12404
19. Fernandez-Becerra C, Xander P, Alfandari D, Dong G, Aparici-Herraiz I, Rosenhek-Goldian I, et al. Guidelines for the purification and characterization of extracellular vesicles of parasites. J Extracellular Biol. (2023) 2:e117. doi: 10.1002/jex2.v2.10
20. Doyle LM and Wang MZ. Overview of extracellular vesicles, their origin, composition, purpose, and methods for exosome isolation and analysis. Cells. (2019) 8(7):727. doi: 10.3390/cells8070727
21. Torres A, Bernardo L, Sanchez C, Morato E, Solana JC, and Carrillo E. Comparing the proteomic profiles of extracellular vesicles isolated using different methods from long-term stored plasma samples. Biol Proced Online. (2024) 26:18. doi: 10.1186/s12575-024-00243-4
22. Montero-Calle A, Garranzo-Asensio M, Rejas-Gonzalez R, Feliu J, Mendiola M, Pelaez-Garcia A, et al. Benefits of FAIMS to improve the proteome coverage of deteriorated and/or cross-linked TMT 10-plex FFPE tissue and plasma-derived exosomes samples. Proteomes. (2023) 11(4):35. doi: 10.3390/proteomes11040035
23. Montero-Calle A, Jimenez de Ocana S, Benavente-Naranjo R, Rejas-Gonzalez R, Bartolome RA, Martinez-Useros J, et al. Functional proteomics characterization of the role of SPRYD7 in colorectal cancer progression and metastasis. Cells. (2023) 12(21):2548. doi: 10.3390/cells12212548
24. Montero-Calle A, Coronel R, Garranzo-Asensio M, Solis-Fernandez G, Rabano A, de Los Rios V, et al. Proteomics analysis of prefrontal cortex of Alzheimer's disease patients revealed dysregulated proteins in the disease and novel proteins associated with amyloid-beta pathology. Cell Mol Life Sci. (2023) 80:141. doi: 10.1007/s00018-023-04791-y
25. Sherman BT, Hao M, Qiu J, Jiao X, Baseler MW, Lane HC, et al. DAVID: a web server for functional enrichment analysis and functional annotation of gene lists (2021 update). Nucleic Acids Res. (2022) 50:W216–W21. doi: 10.1093/nar/gkac194
26. Huang da W, Sherman BT, and Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. (2009) 4:44–57. doi: 10.1038/nprot.2008.211
27. Fonseka P, Pathan M, Chitti SV, Kang T, and Mathivanan S. FunRich enables enrichment analysis of OMICs datasets. J Mol Biol. (2021) 433:166747. doi: 10.1016/j.jmb.2020.166747
28. Iqbal A, Duitama C, Metge F, Rosskopp D, and Boucas J. Flaski. Geneva, Switzerland: Zenodo (2021). doi: 10.5281/zenodo.4849516
29. Xu R, Rai A, Chen M, Suwakulsiri W, Greening DW, and Simpson RJ. Extracellular vesicles in cancer - implications for future improvements in cancer care. Nat Rev Clin Oncol. (2018) 15:617–38. doi: 10.1038/s41571-018-0036-9
30. Tomiyama E, Matsuzaki K, Fujita K, Shiromizu T, Narumi R, Jingushi K, et al. Proteomic analysis of urinary and tissue-exudative extracellular vesicles to discover novel bladder cancer biomarkers. Cancer Sci. (2021) 112:2033–45. doi: 10.1111/cas.v112.5
31. Marti M and Johnson PJ. Emerging roles for extracellular vesicles in parasitic infections. Curr Opin Microbiol. (2016) 32:66–70. doi: 10.1016/j.mib.2016.04.008
32. Pinheiro AAS, Torrecilhas AC, Souza BSF, Cruz FF, Guedes HLM, Ramos TD, et al. Potential of extracellular vesicles in the pathogenesis, diagnosis and therapy for parasitic diseases. J Extracell Vesicles. (2024) 13:e12496. doi: 10.1002/jev2.12496
33. Dong G, Filho AL, and Olivier M. Modulation of host-pathogen communication by extracellular vesicles (EVs) of the protozoan parasite leishmania. Front Cell Infect Microbiol. (2019) 9:100. doi: 10.3389/fcimb.2019.00100
34. Dong G, Wagner V, Minguez-Menendez A, Fernandez-Prada C, and Olivier M. Extracellular vesicles and leishmaniasis: Current knowledge and promising avenues for future development. Mol Immunol. (2021) 135:73–83. doi: 10.1016/j.molimm.2021.04.003
35. Esteves S, Lima C, Costa I, Osorio H, Fernandez-Becerra C, Santarem N, et al. Characterization and proteomic analysis of plasma EVs recovered from healthy and diseased dogs with canine leishmaniosis. Int J Mol Sci. (2023) 24(6):5490. doi: 10.3390/ijms24065490
36. Kaye PM, Svensson M, Ato M, Maroof A, Polley R, Stager S, et al. The immunopathology of experimental visceral leishmaniasis. Immunol Rev. (2004) 201:239–53. doi: 10.1111/j.0105-2896.2004.00188.x
37. Murray HW, Jungbluth A, Ritter E, Montelibano C, and Marino MW. Visceral leishmaniasis in mice devoid of tumor necrosis factor and response to treatment. Infect Immun. (2000) 68:6289–93. doi: 10.1128/IAI.68.11.6289-6293.2000
38. Ruzha Y, Ni J, Quan Z, Li H, and Qing H. Role of vitronectin and its receptors in neuronal function and neurodegenerative diseases. Int J Mol Sci. (2022) 23(20):12387. doi: 10.3390/ijms232012387
39. Jang YC, Tsou R, Gibran NS, and Isik FF. Vitronectin deficiency is associated with increased wound fibrinolysis and decreased microvascular angiogenesis in mice. Surgery. (2000) 127:696–704. doi: 10.1067/msy.2000.105858
40. Sano K, Asanuma-Date K, Arisaka F, Hattori S, and Ogawa H. Changes in glycosylation of vitronectin modulate multimerization and collagen binding during liver regeneration. Glycobiology. (2007) 17:784–94. doi: 10.1093/glycob/cwm031
41. Delanghe JR and Langlois MR. Hemopexin: a review of biological aspects and the role in laboratory medicine. Clin Chim Acta. (2001) 312:13–23. doi: 10.1016/S0009-8981(01)00586-1
42. Tolosano E and Altruda F. Hemopexin: structure, function, and regulation. DNA Cell Biol. (2002) 21:297–306. doi: 10.1089/104454902753759717
43. Ramos S, Jeney V, Figueiredo A, Paixao T, Sambo MR, Quinhentos V, et al. Targeting circulating labile heme as a defense strategy against malaria. Life Sci Alliance. (2024) 7(4). doi: 10.26508/lsa.202302276
44. Castro H, Rocha MI, Duarte M, Vilurbina J, Gomes-Alves AG, Leao T, et al. The cytosolic hyperoxidation-sensitive and -robust Leishmania peroxiredoxins cPRX1 and cPRX2 are both dispensable for parasite infectivity. Redox Biol. (2024) 71:103122. doi: 10.1016/j.redox.2024.103122
45. Sadowska-Bartosz I and Bartosz G. Peroxiredoxin 2: an important element of the antioxidant defense of the erythrocyte. Antioxidants (Basel). (2023) 12(5):1012. doi: 10.3390/antiox12051012
46. Koncarevic S, Rohrbach P, Deponte M, Krohne G, Prieto JH, Yates J 3rd, et al. The malarial parasite Plasmodium falciparum imports the human protein peroxiredoxin 2 for peroxide detoxification. Proc Natl Acad Sci U S A. (2009) 106:13323–8. doi: 10.1073/pnas.0905387106
47. Kumar R, Bumb RA, and Salotra P. Evaluation of localized and systemic immune responses in cutaneous leishmaniasis caused by Leishmania tropica: interleukin-8, monocyte chemotactic protein-1 and nitric oxide are major regulatory factors. Immunology. (2010) 130:193–201. doi: 10.1111/j.1365-2567.2009.03223.x
48. Samant M, Sahu U, Pandey SC, and Khare P. Role of cytokines in experimental and human visceral leishmaniasis. Front Cell Infect Microbiol. (2021) 11:624009. doi: 10.3389/fcimb.2021.624009
49. Pugliese M, Sfacteria A, Oliva G, Falcone A, Gizzarelli M, and Passantino A. Clinical significance of ROMs, OXY, SHp and HMGB-1 in canine leishmaniosis. Anim (Basel). (2021) 11(3):754. doi: 10.3390/ani11030754
50. Dejean E, Foisseau M, Lagarrigue F, Lamant L, Prade N, Marfak A, et al. ALK+ALCLs induce cutaneous, HMGB-1-dependent IL-8/CXCL8 production by keratinocytes through NF-kappaB activation. Blood. (2012) 119:4698–707. doi: 10.1182/blood-2011-10-386011
51. Hernandez-Pando R, Barrios-Payan J, Mata-Espinosa D, Marquina-Castillo B, Hernandez-Ramirez D, Bottasso OA, et al. The role of high mobility group box 1 protein (HMGB1) in the immunopathology of experimental pulmonary tuberculosis. PloS One. (2015) 10:e0133200. doi: 10.1371/journal.pone.0133200
52. Wild CA, Bergmann C, Fritz G, Schuler P, Hoffmann TK, Lotfi R, et al. HMGB1 conveys immunosuppressive characteristics on regulatory and conventional T cells. Int Immunol. (2012) 24:485–94. doi: 10.1093/intimm/dxs051
53. Muxel SM, Aoki JI, Fernandes JCR, Laranjeira-Silva MF, Zampieri RA, Acuna SM, et al. Arginine and polyamines fate in leishmania infection. Front Microbiol. (2017) 8:2682. doi: 10.3389/fmicb.2017.02682
54. Rai AK, Thakur CP, Singh A, Seth T, Srivastava SK, Singh P, et al. Regulatory T cells suppress T cell activation at the pathologic site of human visceral leishmaniasis. PloS One. (2012) 7:e31551. doi: 10.1371/journal.pone.0031551
55. Jawed JJ, Dutta S, and Majumdar S. Functional aspects of T cell diversity in visceral leishmaniasis. BioMed Pharmacother. (2019) 117:109098. doi: 10.1016/j.biopha.2019.109098
56. Schaffer AM and Minguet S. Caveolin-1, tetraspanin CD81 and flotillins in lymphocyte cell membrane organization, signaling and immunopathology. Biochem Soc Trans. (2020) 48:2387–97. doi: 10.1042/BST20190387
57. Kulkarni MM, Jones EA, McMaster WR, and McGwire BS. Fibronectin binding and proteolytic degradation by Leishmania and effects on macrophage activation. Infect Immun. (2008) 76:1738–47. doi: 10.1128/IAI.01274-07
58. Speziale P, Arciola CR, and Pietrocola G. Fibronectin and its role in human infective diseases. Cells. (2019) 8(12):1516. doi: 10.3390/cells8121516
59. Nairz M, Haschka D, Demetz E, and Weiss G. Iron at the interface of immunity and infection. Front Pharmacol. (2014) 5:152. doi: 10.3389/fphar.2014.00152
60. Reyes-Lopez M, Pina-Vazquez C, and Serrano-Luna J. Transferrin: endocytosis and cell signaling in parasitic protozoa. BioMed Res Int. (2015) 2015:641392. doi: 10.1155/2015/641392
61. Wang T, Shen H, Deng H, Pan H, He Q, Ni H, et al. Quantitative proteomic analysis of human plasma using tandem mass tags to identify novel biomarkers for herpes zoster. J Proteomics. (2020) 225:103879. doi: 10.1016/j.jprot.2020.103879
62. Singh B, Su YC, and Riesbeck K. Vitronectin in bacterial pathogenesis: a host protein used in complement escape and cellular invasion. Mol Microbiol. (2010) 78:545–60. doi: 10.1111/j.1365-2958.2010.07373.x
63. Fatoux-Ardore M, Peysselon F, Weiss A, Bastien P, Pratlong F, and Ricard-Blum S. Large-scale investigation of Leishmania interaction networks with host extracellular matrix by surface plasmon resonance imaging. Infect Immun. (2014) 82:594–606. doi: 10.1128/IAI.01146-13
64. Jaiswal P, Ghosh M, Patra G, Saha B, and Mukhopadhyay S. Clinical proteomics profiling for biomarker identification among patients suffering with Indian post kala azar dermal leishmaniasis. Front Cell Infect Microbiol. (2020) 10:251. doi: 10.3389/fcimb.2020.00251
65. Torres A, Younis BM, Alamin M, Tesema S, Bernardo L, Solana JC, et al. Differences in the Cellular Immune Response during and after Treatment of Sudanese Patients with Post-kala-azar Dermal Leishmaniasis, and Possible Implications for Outcome. J Epidemiol Glob Health. (2024) 14:1167–79. doi: 10.1007/s44197-024-00270-0
66. Saha S, Mondal S, Ravindran R, Bhowmick S, Modak D, Mallick S, et al. IL-10- and TGF-beta-mediated susceptibility in kala-azar and post-kala-azar dermal leishmaniasis: the significance of amphotericin B in the control of Leishmania donovani infection in India. J Immunol. (2007) 179:5592–603. doi: 10.4049/jimmunol.179.8.5592
Keywords: extracellular vesicles, visceral leishmaniasis, immunosuppression, TNF antagonist, antimonials, quantitative proteomics, LFQ proteomics analyses, biomarkers
Citation: Bernardo L, Montero-Calle A, Solana JC, Lozano-Rendal M, Torres A, Sánchez C, Barderas R, Moreno J and Carrillo E (2025) Protein dysregulation during Leishmania infantum infection in anti-TNF immunosuppressed mice revealed through quantitative proteomics analysis of extracellular vesicles. Front. Immunol. 16:1634080. doi: 10.3389/fimmu.2025.1634080
Received: 23 May 2025; Accepted: 14 July 2025;
Published: 30 July 2025.
Edited by:
Fabrizio Bruschi, University of Pisa, ItalyReviewed by:
Peter Epeh Kima, University of Florida, United StatesCarolina Verónica Poncini, National Scientific and Technical Research Council (CONICET), Argentina
Rafael Pedro Madeira, Hospital Nove De Julho, Brazil
Copyright © 2025 Bernardo, Montero-Calle, Solana, Lozano-Rendal, Torres, Sánchez, Barderas, Moreno and Carrillo. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Lorena Bernardo, bG9yZW5hLmJlcm5hcmRvQGlzY2lpaS5lcw==; Jose Carlos Solana, amMuc29sYW5hQGV4dGVybm9zLmlzY2lpaS5lcw==
†These authors share first authorship
‡These authors share last authorship