Introduction
============

Leishmaniasis refers to a variety of diseases caused by more than 20 species of intracellular protozoan parasites belonging to the genus *Leishmania*. The clinical spectrum of the disease ranges from simple self-limiting cutaneous ulcers to severe disfiguring mucocutaneous, and even to a fatal visceral disease known as Kala Azar. About 98 tropical and subtropical countries are known to be endemic for this disease, with 350 million people at risk and overall prevalence of 12 million worldwide ([@B1]). Iran is endemic to both cutaneous and visceral leishmaniasis. Cutaneous Leishmaniasis (CL) is commonly caused by *L. major* and *L. tropica* ([@B2]). Visceral leishmaniasis (VL), the most life threatening form, is caused by *L. infantum* and in very rare occasions by *L. tropica* ([@B2], [@B3]).

In VL fever and hepato-splenomegaly are the main clinical signs in which *Leishmania* parasite is dispersed to the internal viscera like spleen, liver and bone marrow ([@B4]). Based on the leishmaniasis clinical symptoms, it is evident that the host immunity factors, *Leishmania* species, and in some cases the *Leishmania* strain, determines the measure of pathogenecity ([@B5]). *Leishmania* spp. has about 8000 genes among only 78 genes are restricted to individual species ([@B6]). In spite of a few species parasite genes implicated in pathogenesis and clinical presentation, the parasite gene expression rates differ greatly among species ([@B6]).

In leishmaniasis, parasites are challenged by the host immune conditions throughout their life cycle such as temperature increase of visceral tissues (liver, spleen or bone marrow). Such challenges causes *Leishmania* experience biochemical changes in which post transcriptional modification are activated and may eventuate into the emergence of the leishmaniasis pathogencity ([@B7]--[@B14]). Proteomics is an invaluable tool for systematic analysis of the proteome. Analysis of proteome is most commonly performed by a combination of 2-DE and mass spectrometry (MS). 2-DE method could separate proteins in first and second dimensions according to their isoelectric and molecular weight points. With the help of 2-DE and the MS, a variable mixture of proteins is separated, visualized and then identified ([@B15]--[@B16]).

In this preliminary study, we compared the proteome mapping, in three Iranian isolates of *Leishmania* species including *L. tropica, L. major* and *L. infantum,* with immobilized pH gradient stripes with linear pH 4--7. Moreover, Liquid Chromatography (LC) - mass spectrometry was used for identification of a number of differentially expressed proteins among the three species.

Materials and Methods
=====================

Leishmania isolates and cell culture
------------------------------------

The proteome of three *Leishmania* species including *L. tropica* (GenBank accession nos. **EF653267**, *L. major* (**JN860745**) and *L. infantum* (**JX289853**) compared and were analyzed. promastigote forms recovered from the "Iranian *Leishmania* parasite bank" located in Leishmaniasis lab, School of Public Heath, Tehran University of Medical Sciences (TUMS). The identity of these strains was already obtained by other molecular DNA based methods ([@B2], [@B17]).

***Cell culture --*** Promastigotes recovered from liquid nitrogen (−196 °C), were mass cultured in RPMI1640 medium (Gibco, Life technologies GmbH, Frankfurt, Germany) supplemented with 15% heat inactivated fetal bovine serum (Gibco, Germany) and 100U/ml penicillin and 100ug/ml streptomycin (Gibco, Germany) and incubated at 24°C. Promastigotes harvested in the stationary phase.

Parasites were harvested washed in sterile Phosphate Buffered Saline (PBS, pH: 7.2--7.4) and were used for protein extraction.

Protein Extraction
------------------

Proteomics analysis was performed on *L. tropica, L. major* and *L. infantum*, three species at the time of study. Promastigotes were harvested by centrifugation at 3000rpm, 4° C and 20 minutes and washed three times in sterile PBS (pH: 7.2--7.4) in the same condition for 10 minutes. The cells were resuspended in 5 mM Tris--HCl, pH 7.8, containing 1 mM phenylmethylsulfonyl fluoride (PMSF (Merck, Germany)). Proteins were precipitated by 10% (w/v) trichloroacetic acid (Merck, Germany) in acetone (Merck) with 0.07 %(w/v) dithiothreitol (DTT) (Merck) for 1 hour at −20°C. The samples were then centrifuged at 17500 g (Hettich, Germany) for 15 minutes at 4° C and the pellets were washed with ice-cold acetone containing 0.07% DTT, incubated at −20°C for 1 h and centrifuged at 4 °C. The samples were then solubilized in lysis buffer (9.5 M urea (Merck), 2% (w/v) CHAPS (Merck), 0.8% (w/v) Ampholyte (Bio-Rad, USA) pH 3--10, 1% (w/v) DTT) ([@B18], [@B19]). The concentration of protein was measured by the Bradford assay with bovine serum albumin BSA (Merck) as the standard ([@B20]).

Two-Dimensional electrophoresis (2-DE)
--------------------------------------

For analytical and preparative gels, 120 μg and 1.2 mg of extracted promastigotes proteins were loaded respectively. IEF was carried out on the 18 cm immobilized pH gradient (IPG) strips (pH 4--7) (Bio-Rad, USA). IPG strips were rehydrated overnight by loading the samples diluted with rehydration buffer containing 8 M urea, 4% CHAPS, 2% ampholyte, 50 mM DTT, and traces of bromophenol blue (Merck). Isoelectric focusing was conducted at 20 °C with Mutiphor II and a DryStrip kit (GE Healthcare, Germany). The running condition was as follows: 300 V for 90 minute, followed by 500 V for 90 min, 1000 V for 3 h and finally 3500 V for 16 h. The focused strips were equilibrated twice in equilibration solution. The first equilibration was performed in a solution containing 6 M urea, 20% (w/v) glycerol, 2% (w/v) SDS (Merck), 1% (w/v) DTT, and 50 mM Tris-HCl (Merck) buffer, pH 8.8. The second equilibration was performed in a solution with 2.5% (w/v) iodoacetamide (Merck). Separation in the second dimension was performed by SDS-PAGE in a vertical slab of acrylamide (Merck) (12% total monomer, with 2.6% cross-linker) using a PROTEAN II Multi Cell (BioRad). The protein spots in analytical and preparative gels were visualized by silver nitrate (Merck, Germany) and Coomassie brilliant blue CBB/ G-250 (Sigma, Germany) respectively (19-21-22).

Gel image Analysis
------------------

GS-800 densitometer (Bio-Rad) was used for scanning of silver stain gels. Gels were analyzed using the Melanie 6 software (GeneBio, Geneva, Switzerland). The molecular masses of protein on gels were determined by co electrophoresis of standard protein markers (GE Heathcare) and *pI* of the proteins were determined by migration of the protein spots on 18 cm IPG (pH 4--7, linear) strips. 2-DE per sample (each species) was run for three biologically independent replicates, percent volume of each spot was estimated and analyzed by one-way analysis of variance (ANOVA) SAS software, and means were compared by the LSD test at P ≤ 0.01. Spots were only considered to be significantly different in abundance at least between two *Leishmania* species when/at *P* ≤ 0.01.

Peptide extraction and mass analysis
------------------------------------

The protein spots of interest were excised from coomassie brilliant blue (CBB) stained gels and analyzed using an Amazon ion trab MS/MS (Bruker Daltonics) Mass spectrometer. Briefly, peptides were solubilized in 0.5 % formic acid and fractionated on a nano flow uHPLC system (Thermo RSLCnano) before online analysis by electrospray ionisation (ESI) mass spectrometry on an Amazon ion trap MS/MS (Bruker Daltonics). Peptide separation was performed on a Pepmap C18 reversed phase column (LC Packings), using a 5 -- 85% v/v acetonitrile gradient (in 0.5% v/v formic acid) run over 45 min. at a flow rate of 0.2 μl / min. Mass spectrometric (MS) analysis was performed using a continuous duty cycle of survey MS scan followed by up to ten MS/MS analyses of the most abundant peptides, choosing the most intense multiply charged ions with dynamic exclusion for 120s. MS data was processed using Data Analysis software (Bruker) and the automated Matrix Science Mascot Daemon server (v2.1.06) ([@B23]). Protein identifications were assigned using the Mascot search engine to interrogate in house databases of protein sequences for *L. major*.

Results
=======

Comparison of proteome patterns
-------------------------------

Protein extracts from the three *Leishmania* species including *L. tropica, L. major* and *L. infantum* were separated by 2-DE gel electrophoresis. Multiple gels from the three independent replications were run to ensure reproducibility of the protein homogenates on the 2-DE gels. The analytical gels were visualized by silver staining. [Fig. 1](#F1){ref-type="fig"} show a representative example of the proteins separated/detected on a 2-DE gel in *L. tropica*, *L. major* and *L. infantum*, where a total weight of 120 μg of proteins had been applied. The gel images were analyzed by Melanie software, and the percent volume of the spots was estimated and compared across the gels. We succeeded in detecting 600 ± 100 spots on the 2-DE gels, from which 638, 590, and 546 spots were statistically analyzed across the replicates in *L. tropica*, *L. major*, and *L. infantum* respectively. A number of 478 spots could be paired in all of the three species (Supplementary Table 1).

![2-DE gels analysis of proteins extracted from different *Leishmania* species: *L. tropica*, *L. major*, and *L. infantum* . In first dimension (IEF), protein was loaded on a 18-cm IPG strip with a linear gradient of pH 4--7. In the second dimension, 12% SDS-PAGE gels were used, with a well for molecular weight standards. Proteins were visualized by silver staining. Arrows represent spots identified by LC Mass spectrometry. All spot numbers correspond to the spot numbers of [Table 1](#T1){ref-type="table"}](IJPA-10-530-g001){#F1}

###### Proteins identified in all three Iranian isolates of *Leishmania* species using LC/MS analysis

  **Spot No. [^a^](#TFN1){ref-type="table-fn"}**   **Protein name & biological process**                                                                                         **Access. No**   **(PI/MW) Exp. [^c^](#TFN3){ref-type="table-fn"}**   **(PI/MW) Theo. [^d^](#TFN4){ref-type="table-fn"}**   **% Volume**          **Score. [^e^](#TFN5){ref-type="table-fn"}**   
  ------------------------------------------------ ----------------------------------------------------------------------------------------------------------------------------- ---------------- ---------------------------------------------------- ----------------------------------------------------- -------------- ------ ---------------------------------------------- -----
                                                   **Cell motility and cytoskeleton**                                                                                                                                                                                                                                                                                             
  382                                              *beta tubulin*                                                                                                                A4HTR1           4.79/45                                              4.74/49                                               0.05           0.02   0.02                                           407
  457                                              *hypothetical protein, conserved (Contains Myosin_tail_1 domain) 1 HECT (E6AP-type E3 ubiquitin-protein ligase) domain.)*     Q4QCC9           6.08/14                                              4.95/347                                              0.02           ND     ND                                             43
  510                                              *ADF/Cofilin*                                                                                                                 A4I4A3           6.35/17                                              6.61/15                                               0.17           0.14   0.10                                           540
                                                   **Cell signaling and vesicular trafficking**                                                                                                                                                                                                                                                                                   
  309                                              *hypothetical protein, conserved (Contains Ras_like_GTPase superfamily)*                                                      A4I8E6           6.46/42                                              5.92/24                                               0.02           0.07   0.05                                           820
  501                                              *calmodulin, putative*                                                                                                        A4HU13           4.46/15                                              4.10/16                                               0.05           0.03   0.15                                           773
                                                   **Intracellular survival / proteolysis**                                                                                                                                                                                                                                                                                       
  31                                               *calpain-like cysteine peptidase, putative,cysteine peptidase, Clan CA, family C2, putative*                                  E9AD27           6.14/15                                              5.14/700                                              0.06           0.12   0.07                                           83
  35                                               *calpain-like cysteine peptidase, putative,cysteine peptidase, Clan CA, family C2, putative*                                  E9AD27           5.96/14                                              5.14/700                                              0.04           ND     ND                                             57
  507                                              *calpain-like cysteine peptidase, putative,cysteine peptidase, Clan CA, family C2, putative*                                  E9AD27           6.54/14                                              5.14/700                                              ND             0.10   ND                                             42
  774                                              *calpain-like cysteine peptidase, putative,cysteine peptidase, Clan CA, family C2, putative*                                  E9AD27           5.68/66                                              5.14/700                                              ND             0.09   ND                                             45
  121                                              *hs1vu complex proteolytic subunit-like,hs1vu complex proteolytic subunit-like, threonine peptidase, Clan T(1), family T1B*   Q8I496           5.38/42                                              6.05/24                                               0.07           ND     ND                                             501
                                                   **Lipid metabolism**                                                                                                                                                                                                                                                                                                           
  70                                               *hypothetical protein, conserved (Contains Acyl transferase domain ;Acyl transferase domain. S-malonyltransferase* ***)***    E9AFH1           6.68/19                                              6.3/59                                                0.12           0.34   0.25                                           41
                                                   **microtubule motor movement**                                                                                                                                                                                                                                                                                                 
  53                                               *hypothetical protein, conserved (Contains Ankarin repeat domain)*                                                            A4IC52           5.42/16                                              5.63/38                                               0.52           0.48   0.28                                           587
                                                   **Nucleotid metabolism**                                                                                                                                                                                                                                                                                                       
  369                                              *hypothetical protein, conserved (Contains Nucleoside 2-deoxyribosyltransferase domain)*                                      E9AH20           5.1/17                                               4.95/16                                               0.03           0.04   0.30                                           351
                                                   **oxidative stress defense**                                                                                                                                                                                                                                                                                                   
  54                                               *tryparedoxin*                                                                                                                Q6RYT3           5.22/18                                              5.20/16                                               0.35           0.06   0.01                                           322
                                                   **Protein synthesis**                                                                                                                                                                                                                                                                                                          
  142                                              *translation elongation factor 1-beta, putative*                                                                              A4IDB2           4.88/51                                              4.61/23                                               0.16           0.01   0.06                                           742
  503                                              *60S acidic ribosomal protein P2*                                                                                             E9AGM1           4.61/16                                              4.18/10                                               ND             0.15   0.18                                           408
                                                   **Amino-acid biosynthesis**                                                                                                                                                                                                                                                                                                    
  527                                              *pyrroline-5-carboxylate reductase*                                                                                           E9AGK4           6.32/43                                              6.21/28                                               0.12           0.03   0.14                                           529
                                                   **Protein ubiquitination / proteolysis**                                                                                                                                                                                                                                                                                       
  214                                              *hypothetical protein, conserved (Contains 1 HECT (E6AP-type E3 ubiquitin-protein ligase) domain.)*                           Q4Q4M9           5.88/19                                              6.12/454                                              0.03           0.51   ND                                             45
                                                   **Transport**                                                                                                                                                                                                                                                                                                                  
  438                                              *calcium channel protein, putative,ion transporter, putative*                                                                 Q4Q3G0           5.36/16                                              8.6/291                                               0.06           0.34   0.05                                           50
                                                   **Stress related proteins/Protein folding**                                                                                                                                                                                                                                                                                    
  6                                                *fk506-binding protein 1-like protein*                                                                                        Q4QBK4           5.02/14                                                                                                    0.23           0.35   ND                                             204
  222                                              *hypothetical protein, conserved (Contains alphacrystallin-Hsps_p23-like Superfamily)*                                        A4IBY7           4.48/39                                              4.15/21                                               0.52           0.07   0.17                                           441
                                                   **Unknown**                                                                                                                                                                                                                                                                                                                    
  402                                              *hypothetical protein, conserved*                                                                                             Q4Q8Y7           4.52/54                                              6.27/95                                               0.02           0.06   0.02                                           45
  419                                              *hypothetical protein, conserved*                                                                                             Q4Q2Z8           5.42/72                                              9.93/129                                              0.02           0.10   0.01                                           32
  475                                              *hypothetical protein*                                                                                                        Q4QE59           4.87/26                                              10.88/58                                              0.12           0.01   ND                                             35

**LC/MS analysis./**

The numbering corresponds to the 2-DE gel in [Figure 1](#F1){ref-type="fig"}./

Accession number in Swiss-Prot./

Experimental p*I* and molecular weight./

Theoretical p*I* and molecular weight./

Mascot score./ **ND**: Not Detected (Spot)

However 34, 33, and 4 protein spots of *L. tropica*, *L. major*, and *L. infantum* were identified as *Leishmania* species-specific spots ([Fig. 2](#F2){ref-type="fig"}). Altogether, 265 protein spots exhibited reproducible quantitative (*P* ≤ 0.01) changes across the three samples in *Leishmania* species (Supplementary Table 2). Among them, 35, 22 and 11 protein spots were different between the *L. tropica* and *L. major* (presented by **a**), *L. tropica and L. infantum* (presented by **b**), *L. major* and *L. infantum* (presented by **c**) respectively. Seven protein spots were different between all of the three species (presented by **abc**).

![Venn diagram showing shared and unique protein expression between different *Leishmania* species: *L. tropica, L. major*, and *L. infantum*. Regions of overlap between circles indicate gene expression which paired in all the three species/two species. Regions that do not overlap between circles indicate unique gene expression of particular/each species as species-specific spots](IJPA-10-530-g002){#F2}

Protein identification
----------------------

Within differentially expressed proteins, which were detected on the analytical gels, we could reliably detect and excite a total number of 28 protein spots on CBB-stained preparative gels. Due to the lack of the protein amount, the remaining proteins could not be detected. The excised protein spots were then analyzed by LC/MS leading to the identification of 24 proteins ([Table 1](#T1){ref-type="table"}). The identified protein spots by mass spectrometry have been labeled in [Fig. 1](#F1){ref-type="fig"}. These proteins were classified in multiple categories according to their species, functions, and biological processes: Cell motility and cytoskeleton, Cell signaling and vesicular trafficking, Intracellular survival / Nucleotid metabolism, Protein synthesis, oxidative stress defense, microtubule motor movement proteolysis, Lipid metabolism, Amino-acid biosynthesis, Protein ubiquitination / proteolysis, Transport, Stress related proteins/Protein folding and hypothetical proteins (Unknown) ([Table 1](#T1){ref-type="table"}). It is worth noting that some of these proteins are hypothetical and their functions in *Leishmania* still remain to be elucidated.

The clustering of protein expression pattern of differentially expressed proteins in *L. tropica, L. major* and *L. infantum* is presented in [Fig. 3](#F3){ref-type="fig"}. All quantitative information is showed using a color scale in which the color ranges from green for the highest down-regulation to red for the highest up-regulation. Black color indicates no changes in expression pattern of the 3 *Leishmania* species. A comparison among the 3 species revealed that the changes in expression pattern were more pronounced in *L. infantum* compared with *L. tropica* and *L. major.* In addition, the number of up regulated proteins was higher than that of down-regulated proteins. The functional annotation of the 3 species identified proteins in *L. tropica, L. major* and *L. infantum* classified by biological function and processes described in [Table 1](#T1){ref-type="table"}. The biological function pie charts of the cutaneous species (*L. tropica* and *L. major*) are highly similar whereas, the pie chart for the biological function of the visceral species (*L. infantum*) is different from those of cutaneous species ([Fig. 4](#F4){ref-type="fig"}). Among the cutaneous species, most of the proteins functionalities are observed in the following categories: intracellular survival / nucleotid metabolism and cell motility/ cytoskeleton (20--22%). Moreover, the least functionalities are observed among cell signaling/ vesicular trafficking and lipid metabolism (5%). Moreover, in the *L. infantum* species, most functionality are observed among protein synthesis and cell motility/ cytoskeleton (15%) and least of them are observed among amino-acid biosynthesis and oxidative stress defense (7%).

![Clustering pattern of differential proteins expressed in *L. tropica, L. major* and *L. infantum.* All quantitative information is showed using a color scale in which the color ranges from green for the highest down-regulation to red for the highest up-regulation. Black color indicates no changes in expression pattern of 3 *Leishmania* species](IJPA-10-530-g003){#F3}

![Functional annotation of the identified proteins of 3 species *L. tropica, L. major* and *L. infantum* classified by biological function and processes in described in [Table 1](#T1){ref-type="table"}](IJPA-10-530-g004){#F4}

Discussion
==========

The aim of this preliminary study was to apply 2-DE, which is a valuable method in the proteomics arena to analyze the protein profile patterns of the three Iranian cutaneous and visceral *leishmania* species including *L. tropica, L. major* and *L. infantum* to search for species-specific *Leishmania* proteins. In recent years by completion of genome sequencing *Leishmania* parasites coupled with protein separation techniques analysis has given us an insight in better understanding the mechanisms of pathogenesis of *Leishmania* species. Proteomics approaches at the protein expression level provide additional information for further analysis with a biological function. Moreover, comparative proteomics has been successful in determining the virulence biomarkers ([@B24]).

Overall, the proteome 2-DE maps of *L. tropica, L. major* and *L. infantum* were strikingly similar in terms of protein distribution and positioning. By using 18-cm IPG strips at the p*I*4--7 range in 2-DE comparative analysis, we successfully identified more than 700 protein spots for all the three species. These numbers represent about 9% of total proteins in *Leishmania* genome projects ([@B25]). The analysis and discussion about the detected proteins in this study could be categorized into three parts. First, we discuss comprehensively about the biological function of proteins whose expression abundances were common among each of the three species. In the second and third part, we represent the biological function of the proteins which are different and absent/present in the three different species. As it is described in [Table 1](#T1){ref-type="table"} and Venn diagram ([Fig. 2](#F2){ref-type="fig"}), in the first part we have studied a total number of 478 proteins among which we will discuss the common ones. Among the common proteins of three *Leishmania* species, which were surveyed, we would like to mention significant proteins such as beta tubulin and ADF/Coflin. These proteins could be categorized in the cell mobility and cytoskeleton group. ADF/cofilin is existent in all eukaryotic organisms and has been involved in cell motility and cytokinesis. *Leishmania* parasites express only one isoform of ADF/cofilin, which is essential for flagellar assembly and motility ([@B26]). Beta tubulin is known as one of the members of distinct microtubule networks in *Leishmania* and is implicated in locomotion, cell shape and division ([@B27]).

Within other common proteins we could name tryparedoxin, calpain-like cysteine peptidase, and calmodulin putative from the groups of oxidative stress defense, Intercellular survival/proteolysis, cell signalling, and vesicular trafficking, respectively. Tryparedoxins are special thiol disulfide oxidoreductases related to thioredoxins, which play a crucial role in hydroperoxide detoxification cascades of Kinetoplastida ([@B28]).

One of the cell signaling proteins was calmodulin-like preotein. Calmodulin is a kind of calcium binding protein. It participates in calcium signaling pathways that regulate multiple critical processes such as growth and proliferation ([@B29], [@B30]). The biological mechanism of elongation factor puts it in the protein synthesis group. Elongation factor1alpha plays a role in protein synthesis and assembly. It is a highly conserved GTP-binding protein involved in protein translation. In addition, it is an actin/microtubule-binding protein which interacts with the cytoskeleton ([@B31]--[@B33]). Another common protein could be calcium channel protein that is a member of the transporter group. A huge number of proteins exist under the group of hypothetical proteins whose biological mechanisms are still not well discovered.

The second group of proteins, which were studied, was the ones, which were different in all three *Leishmania* species. The differences in proteins mainly occur between *L. tropica* and *L. major* and the least differences occur between *L. major* and *L. infantum* ([Fig. 2](#F2){ref-type="fig"}). It is an interesting point that the proteins, which are different between *L. tropica* and *L. infantum* are less than those of *L. major* and *L. tropica*. *L. tropica* is the main cause of dry cutaneous leishmaniasis lesions; whereas, *L. major* is the main cause of wet CL lesion and *L. infantum* causes visceral leishmaniasis. Recently there have been reports of viscerotropic forms caused by *L. tropica* in either humans or dogs infected to visceral leishmaniasis in Iran and different parts of the world ([@B34], [@B2]). Among significant proteins which are different in *L. tropica* and *L. major* we could mention calcium channel proteins, tryparedoxin and Elongation factor from the groups of oxidative stress defense and protein synthesis respectively whose mechanism have been described previously. Recently, it was demonstrated that *Leishmania* EF-1alpha acts as a virulence factor ([@B35]). This protein could diffuse into the cytosol of infected macrophages, where it is able to activate tyrosine phosphatase-1 leading to macrophage deactivation ([@B31]). Among the proteins different in *L. tropica* and *L. infantum* we could mention proteins such as calmadolin and trypardoxin in which the former role is to transfer the material into the cells and the latter role is defending the host cell against oxidative factors and preventing its death. Other proteins are hypothetical and have an unknown mechanism. Nevertheless, some of these proteins have a specific domain and have a different mechanism such as hypothetical protein, conserved contains nucleoside 2-deoxyribosyltransferase domain, that have a role in nucleotide metabolism.

Among the most significant proteins between *L. major* and *L. infantum* we could mention Calmadolin. The remaining proteins are partially hypothetical. Another group of proteins, found in this study, was too rare and introduced as absent/present.

However, further studies are needed to define precise biological function for the mentioned proteins in the process of pathogenesis. Moreover, the examination of these species must be repeated with other strains in order to sanction the results. In addition, such differences must be evaluated and approved by other methods such as real time PCR and western blotting with monoclonal antibodies.

Conclusion
==========

The analysis of proteome mapping of 3 *Leishmania* species including, *L. tropica, L. major* and *L. infantum* by 2-DE and mass spectrometry demonstrated that the vast majority of *Leishmania* proteins are commonly expressed among 3 species. Therefore, differentiation, virulence and pathogenesis may be related not only to the immunity situation of the host but also to the differentiation expression of a number of proteins like Stress related proteins/Protein folding and protein ubiquitination/proteolysis. It must be pointed out that further studies must be undertaken using western blotting or real time PCR in order to support the results of the current study.

The authors would like to express their gratitude to Prof. Ghasem Hosseini Salekdeh from Agricultural Biotechnology Research Institute of Iran, (ABRI), for their kind cooperation in proteomics laboratory. We are grateful to Elham Sarhadi from ABRI for valuable technical assistant. Also special thanks from Prof. Burchmore and Ms Suzzane Eideman from Institute of Infection, Immunity and Inflammation, College of Medical, Veterinary and Life Science, University of Glasgow, United Kingdom for Mass spectrometry analysis. In addition, we are grateful to Dr S.R. Naddaf from Pasteur Institute of Iran for valuable assistant. The authors would like to express their gratitude to Dr. B. Akhoundi and Mrs. S. Charedar for their kind cooperation. This study was financially supported by the Vice Chancellor for Research, School of Public Health, and Tehran University of Medical Sciences, Iran (NO.10695). The authors declare that there is no conflict of interests.
