- Research article
- Open Access
Evolutionary comparison of prenylation pathway in kinetoplastid Leishmania and its sister Leptomonas
BMC Evolutionary Biology volume 15, Article number: 261 (2015)
Leptomonas is monogenetic kinetoplastid parasite of insects and is primitive in comparison to Leishmania. Comparative studies of these two kinetoplastid may share light on the evolutionary transition to dixenous parasitism in Leishmania. In order to adapt and survive within two hosts, Leishmania species must have acquired virulence factors in addition to mechanisms that mediate susceptibility/resistance to infection in the pathology associated with disease. Rab proteins are key mediators of vesicle transport and contribute greatly to the evolution of complexity of membrane transport system. In this study we used our whole genome sequence data of these two divergent kinetoplastids to analyze the orthologues/paralogues of Rab proteins.
During change of lifestyle from monogenetic (Leptomonas) to digenetic (Leishmania), we found that the prenyl machinery remained unchanged. Geranylgeranyl transferase-I (GGTase-I) was absent in both Leishmania and its sister Leptomonas. Farnesyltransferase (FTase) and geranylgeranyl transferase-II (GGTase-II) were identified for protein prenylation. We predict that activity of the missing alpha-subunit (α-subunit) of GGTase-II in Leptomonas was probably contributed by the α-subunit of FTase, while beta-subunit (β-subunit) of GGTase-II was conserved and indicated functional conservation in the evolution of these two kinetoplastids. Therefore the β-subunit emerges as an excellent target for compounds inhibiting parasite activity in clinical cases of co-infections. We also confirmed that during the evolution to digenetic life style in Leishmania, the parasite acquired capabilities to evade drug action and maintain parasite virulence in the host with the incorporation of short-chain dehydrogenase/reductase (SDR/MDR) superfamily in Rab genes.
Our study based on whole genome sequences is the first to build comparative evolutionary analysis and identification of prenylation proteins in Leishmania and its sister Leptomonas. The information presented in our present work has importance for drug design targeted to kill L. donovani in humans but not affect the human form of the prenylation enzymes.
Visceral leishmaniasis (VL) is endemic in 98 countries with 350 million people at risk around the world and 300,000 are infected/year (http://www.dndi.org/diseases/leishmaniasis.php). More than 90 % of visceral leishmaniasis (VL) occurs in seven countries: India, Bangladesh, Nepal, Sudan, Ethiopia, Kenya and Brazil (http://www.dndi.org/diseases/leishmaniasis.php). Leishmaniasis is caused by protozoan parasites of Leishmania genus. VL also known as kala-azar in India, is the most severe form of leishmaniasis (http://www.dndi.org/diseases-projects/diseases/vl.html) and anthroponotic transmission (human to vector to human) occurs. Whole genome sequencing of clinical isolates of Leishmania donovani and Leptomonas has previously been completed in our lab . An important insight gained through this effort was the confirmation of co-infection of Leptomonas with the visceralizing Leishmania species. This disturbing trend has been reported so far only in India [1–3] hinting towards zoonotic spread of the disease and indicates that Leptomonas may be a new pathogen. In such a situation the question remains that a changed paradigm in chemotherapy should be adopted targeting the co-infections.
Leptomonas is monogenetic kinetoplastid parasite of insects and primitive in comparison to Leishmania  and these two parasites arose from a common ancestor as shown in Fig. 1 [4, 5]. At some time during the evolution, Leishmania appears to have lost the ability to be transmitted in nature from invertebrate host to another and adapted to a life cycle alternating between invertebrate and vertebrate host [6, 7]. With the availability of complete genome sequence with us of these two kinetoplastids  representing important evolutionary branch points viz Leishmania and Leptomonas, in this study we have analyzed the enzymes prenyltransferases (farnesyltransferase and geranylgeranyltransferase) of these kinetoplastids. Cellular organization and signaling in both unicellular and multicellular organisms is heavily influenced by the Ras superfamily of small GTP-binding proteins. These proteins have a structurally and mechanistically preserved GTP-binding core despite considerable divergence in sequence and function . Protein trafficking pathways are frequently exploited by human pathogens to gain entry and survive within host cells . Studying the evolution of protein trafficking is essential to understand the origins of eukaryotes. Among trafficking-associated proteins, the Rab family expanded most in evolution, suggesting that it provided the primary diversification element in the evolution of trafficking . The prenylation pathway remains largely unexplored in Kinetoplastida. Therefore, it was of interest for us to carry out the phylogenetic profiling of Rab geranylgeranyltransferase (RabGGTase) and Rab proteins in both Leishmania and its sister Leptomonas. In this study we used a classification and phylogeny based approach to define subfamilies of RabGGTase and Rab proteins, we built evolutionary trees using maximum likelihood and distance based methods. We provide a comprehensive view of RabGGTase and Rab evolution for sister Leptomonas and with various Leishmania species. Since, drugs that inhibit enzymes involved in protein prenylation could be interesting antiparasitic agents; we used one specific known inhibitor of RabGGTase to predict the structure activity relationship with the parasite enzyme.
Results and discussion
Kinetoplastid genome phylogeny
We compared the genomes of Leishmania species to Leptomonas in order to obtained information on kinetoplastid phylogeny (Fig. 2). Crithidia was taken as an outgroup. We observed that the L. donovani isolate from India was placed at the crown in comparison to Leptomonas.
Leishmania donovani and Leptomonas farnesyltransferase
In Leishmania donovani, α-subunit of farnesyltransferase (LdFTase alpha) was located on chromosome 29 (LdBPK_292070) with accession number XP_003862625.1. In Leptomonas α-subunit of FTase (Leptomonas FTase alpha) was identified from the contig (contig_2652) of Leptomonas genome. Using ClustalW, percent identity of the α-subunit of LdFTase and Leptomonas FTase with L. major (XP_003722277.1) and L. infantum (XP_001466722.1) was summarized in Fig. 3a and it can be seen that the region of sequence identity was found extensively at the C terminus (Additional file 1: Figure S1).
β-subunit of FTase in Leishmania donovani (LdFTase beta) was located on chromosome 26 (LdBPK_261450) with accession number XP_003861732.1. In Leptomonas, β-subunit of FTase (Leptomonas FTase beta) was identified from the contig (contig_1135) of Leptomonas genome. Using ClustalW, percent identity of the β-subunit of LdFTase and Leptomonas FTase with L. major (XP_001684151.1) and L. infantum (XP_001470492.1) was summarized in Fig. 3a. In comparison with the α-subunit of FTase among L. major, L. infantum and L. donovani the β-subunit showed more sequence identity (Additional file 1: Figure S2). It is also reported that the β-subunit of prenyltransferases are more conserved than α-subunit . The β-subunit appears to have invariable domain architecture throughout evolution and its invariable nature is due to functional rather than structural constraints.
Leishmania donovani and Leptomonas geranylgeranyltransferase
Using ClustalW, percent identity of the α-subunit of LdGGTase-II (XM_001468149) with L. major (XM_001685808) and L. infantum (XP_001468186) was summarized in Fig. 3b. α-subunit of this enzyme was totally conserved between the two visceral species L. donovani and L. infantum but quite variable between the cutaneous L. major and visceral species (Additional file 1: Figure S3). Non conservative amino acids alanine and serine at positions 82, 88, 214, 440 and 235, of the visceral species were replaced by valine and threonine in the cutaneous species respectively (Additional file 1: Figure S3). It has been reported by Peterson et al. that parasites harboring a pair of point mutations from Ala-16 to Val-16 and from Ser-108 to Thr-108 are resistant to cycloguanil . A similar substitution in α-subunit of GGTase-II could be responsible for its survival and multiplication within the host cell. No sequence identity was observed to α-subunit of GGTase-II in Leptomonas genome with L. major (XM_001685808), L. infantum (XP_001468186) and LdGGTase-II (XM_001468149).
In our study, the β-subunit of GGTase-II of L. donovani was totally conserved in both the Leishmania species responsible for the visceral manifestation of the disease viz. L. infantum (XP_001468743) and L. donovani (XP_003864545). The carboxy terminal of this enzyme was completely conserved between the visceral (L. infantum, XP_001468743) and cutaneous species (L. major, XP_001686510) of Leishmania (Additional file 1: Figure S4) however in the N-terminal, differences in the cutaneous species were observed. Non conservative amino acids aspartate and glycine at positions 57 and 61 of the visceral species were substituted in the cutaneous species (Additional file 1: Figure S4). The susceptibility to infection by Leishmania NRAMP1 gene has been similarly reported  to be substitution in non conservative amino acid aspartate and glycine. A similar substitution in β-subunit of GGTase-II maybe due to the different requirement for mode of protein trafficking.
Percent identity of the β-subunit of GGTase-II in Leptomonas to Old World L. major (XP_001686510) and L. infantum (XP_001468743) is shown (Figs. 3b and Fig. 4). New World, L. panamensis (XP_010698581) and L. braziliensis (XP_001564646), showed 81 % identity with Leptomonas whereas in L. donovani identity with only Old World was observed.
The evolutionary relationship between FTase and GGTase-II of Leishmania and Leptomonas was reconstructed by using Neighbor-Joining method (Fig. 5). In our study, α-subunit of FTase between Leishmania and Leptomonas were in close phylogeny. But as regards to the holoenzyme, α-subunit was placed very differently to the β-subunit of these enzymes. We concluded that α-subunit of GGTase-II is almost similar to the corresponding α-subunit of FTase in both kinetoplastids. β-subunit of GGTase-II is very similar to the corresponding β-subunit of FTase in both kinetoplastids.
Rab sequence mining in L. major and L. infantum
Among the prenylation proteins, Rabs are small GTP-binding proteins which are the largest family within the Ras superfamily. These act as membrane associated molecular switches and regulate vesicular trafficking [14–16]. The number of Rab proteins and gene complexity in an organism has been correlated with multicellularity of organism. A multicellular organism will require/have higher number of Rabs encoded by its genome . Plasmodium falciparum has been shown to encode 11 Rabs, Toxoplasma gondii encodes 15 Rabs, Leishmania encodes 12 Rabs, Drosophila melanogaster encodes 29 Rabs, and Homo sapiens encodes > 60 Rab proteins, which correlates with their multicelluarity and organism structural complexity.
Until now, despite the importance of the intracellular vesicle transport in Leishmania pathogenicity, very few publications are available regarding number, diversity, complexity and functions of Leishmania Rab gene family [18–20]. Therefore, it was of interest for us to carry out the phylogenetic profiling of the Rab repertoire in the lower eukaryote Leishmania. The exhaustive mining of Leishmania Rab proteins in GeneDB (www.genedb.org) and NCBI (http://www.ncbi.nlm.nih.gov/) database revealed many full length proteins, some were annotated as Rabs of various families, but most were annotated as putative Rabs or Rab like proteins, in this aspect the L. major Rabs are better annotated compared to L. infantum Rabs. Our search retrieved 14 completely annotated Rabs in L. major whereas 18 Rabs are reported in L. major by Berriman et al. . Similarly in L. infantum, we found 13 proteins annotated as Rabs and few other proteins annotated as putative Rab or Rab like proteins. The details of Rabs of both L. infantum and L. major, along with name, accession number, gene ID and chromosomal location in Additional file 2: Table S1.
Rab sequence mining in Leishmania and Leptomonas
From our whole genome sequence of L. donovani and Leptomonas some Rabs were annotated as Rabs of various families, but most were annotated as putative Rabs or Rab like proteins. We retrieved 13 completely annotated Rabs in L. donovani (Additional file 1: Figure S5). Similarly in Leptomonas, we found 11 proteins annotated as Rabs (Additional file 1: Figure S6). The details of Rabs of both L. donovani and Leptomonas are presented in Additional file 2: Table S2 and S3 respectively. A pair wise comparison of the Rab proteins retrieved from the L. major (Fig. 6a) and L. infantum (Fig. 6b), L. donovani (Additional file 2: Table S4A) and Leptomonas (Additional file 2: Table S4B) were calculated and phylogenetic tree was constructed using the Rabs orthologues sequences of L. major, L. infantum, L. donovani and Leptomonas (Additional file 1: Figure S7).
To visualize the synteny in Rab genes between Leishmania donovani and Leptomonas, we constructed a Circos plot linking the homologous Rab genes (Additional file 1: Figure S8) The Rab genes of Leptomonas are represented on the left side and Rab genes of Leishmania donovani are represented on the right side. Rab4, Rab7 GTP and Rab18 showed extensive synteny in both the species. Rab1 and Rab7 of Leishmania donovani have no similarity with any genes in Leptomonas. Whereas few genes in Leishmania donovani are partially similar such as Rab28, Rab14, Rab11B, Rab11, Rab6, Rab5, Rab2A and GTP binding Rab1 to Leptomonas.
Thirteen annotated Rabs of L. major and L. infantum were exactly the same as in L. donovani, only Rab 21 was present in L. major but was not identified in either of the visceral Leishmania. What does Rab 21 do? The explanation for this has been addressed in the Phylogeny section. In Leptomonas 11 Rabs could be annotated by us, all of them same as found in Leishmania.
Phylogenetic analysis of Rab GTPase
An important role in evolution is played by small GTPases belonging to the Ras superfamily. Rab proteins form the largest branch of the Ras superfamily. This family is much diversified . Each Rab protein has a distinct subcellular location and is responsible for a specific transport step [22, 23]. Rab proteins are key regulators of intracellular vesicular transport and membrane trafficking in exocytic and endocytic pathways [14–16]. The evolution of Rab proteins has been analyzed extensively in Trypanosomes  but not in Leishmania. We believe that molecular phylogenetic analysis of Rab proteins would provide a unique opportunity for finding the relationship between evolution of Leishmania and Leptomonas. Leishmania is believed to have evolved from Leptomonas where life cycle changed from monogenetic to digenetic . Our study is the first extensive data mining and phylogenetic construction of clinical isolates of Leishmania and Leptomonas which is based on annotation of Rabs obtained from our whole genome sequencing .
The orthologues which correspond to genes separated by divergence of the Leishmania Rab superfamily in comparison with Leptomonas is shown in Fig. 7a. We found that the core Rab repertoire associated with basic endocytosis and exocytosis is fully conserved , the endomembrane system appears to be very ancient and maintained same for Leishmania and Leptomonas, whilst gene duplication has facilitated the building of additional, lineage-specific complexity .
Out of 20 Rabs reported to be present in the last eukaryotic common ancestor (LECA) , five Rabs are said to be present in any well characterized genome. The core set includes a conserved mechanism for ER to Golgi transport (Rab1 and Rab2), trafficking through the early endosome (Rab5), recycling (Rab4 and Rab11), delivery to the late endosome/lysosome (Rab7) and retrograde transport through the Golgi complex (Rab6) . Inspection of the tree in Fig. 7a reveals that major groups of Rab proteins in the Leishmania genome can be broadly related to Leptomonas Rabs by homology thus indicating similar function and/or subcellular localization. This co-segregation of Rab GTPases according to common functions rather than to taxonomic relationships indicates a conserved mechanism of Rab interaction with regulators/effectors across evolution . Orthologues of Rabs in Leishmania and Leptomonas e.g. Rab2, Rab4, Rab5, Rab6, Rab7, Rab18 and Rab28 co-segregated with a strong bootstrap support. This observation suggests that strict phylogeny of function in the Rab superfamily has been maintained and steps required for the basic vesicular transport were almost same in both these kinetoplastids. These core set of Rabs also appear to be orthologues of mammalian Rabs . Remarkably high level sequence conservation is retained between these two kinetoplastids and higher eukaryote orthologues. Parasites (T. bruci, T. curci and L. major) have ~ 20 Rabs gene and most of which are clearly orthologues with metazoan Rabs . Rab28 is ancient Rab protein that arose before the speciation event separating trypanosomes from the mammalian lineage  and is maintained in the kinetoplastid Leishmania and its sister species Leptomonas. Rab8 is lost in kinetoplastids .
Non-essentiality of the different Rab isoforms was observed in both the kinetoplastid genomes but which is observed to be a major contribution in the Rab repertoire of other organisms. We found isoform of Rab2 (2A) and its orthologues with Leptomonas in the crown position of phylogenetic tree. Isoform of Rab11 (11B) in Leishmania and Leptomonas takes basal position in the tree. PrePS analysis (PrePS- Prenylation Prediction Suite, http://mendel.imp.ac.at/sat/PrePS/index.html) of these protein showed that Rab11B lacks the cysteine residue at C terminus of CAAX signal sequences which is essential for prenylation (Fig. 7b). On blastp analysis with the well annotated human Rab family, this Rab11B of L. donovani, showed 58 % identity with Rab like protein 2A isoform X7 (XP_006712281.1). Therefore we determined that the Rab11B as annotated in Leishmania and Leptomonas genome is actually not Rab11B but an ancient isoform of Rab2A. Maybe this protein has conserved GTPase activity and as it evolved to Rab2A (crown position) it acquired specific prenylation activity (Fig. 7c). Alternatively spliced transcript variants encoding different isoforms have been reported. Evolutionary distance of the lineage requires experimental evidence to further support these sequence data. It has been established by Diniz et al. : that Rab2 and Rab11B (which we now identify to be isoform of Rab2A) are actin interacting proteins (AIP) in L. infantum and L. major.
Rab5 is an early endosome marker as established by Real et al. [28, 29]: in Leishmania major and Leishmania amazonensis. From our study we found only Rab5 annotation in the genomes of both the kinetoplatids. In L. donovani Rab5 has been shown to be is involved in hemoglobin endocytosis . Two isoforms of Rab5 (Rab5A and Rab5B) were reported by them. The strain used by Singh et al. : was UR6 which is a very long term laboratory cultivated avirulent strain and has been established to be L. tropica  whereas our study of Rab annotations is based on whole genome sequences of recently cultivated clinical isolates of L. donovani and Leptomonas. On blastp analysis of the Rab5 of both L. donovani and Leptomonas with human (taxid 9606, accession number NP_002859.1) showed 45 % and 50 % identity to isoform 5B respectively and Rab5b may act as a marker for early endosomes in L. donovani . No isoform A in the genome sequence of these two kinetoplastids was present thereby indicating only Rab5 orthologues.
Rab1 GTPase and Rab7 GTPase are actually true Rab1 and Rab7 respectively on the basis of bastp analysis of human Rab family and these were identified as orthologues in both L. donovani and Leptomonas. Rab7 and Rab9 are two related subfamilies, which show overlapping localization to late endosomes . Rab 9 was absent in the kinetoplastid genomes. Rab9 protein has not yet been reported in unicellular eukaryotes, plants and fungi . Rab7 subfamily belongs to the ancestral set of Rab and arose before the radiation of eukaryotes. Rab7b isoforms are found only in representatives of amphibians, birds and mammals . The Rab1 and Rab7 identified only in Leishmania were branched together which suggests paralogous gene duplications. This is a major force in evolution which provides new gene function. Retroposition or lack of introns, as observed in kinetoplastids, is a major contributor for these phenomena. However on blastp analysis with human Rab family, these showed identity 43 and 37 % to Rab 13 and Rab 43 respectively.
RAB11 was found to be a useful marker to track flagellar pocket division and to follow mitochondrial DNA (kinetoplast) segregation. Leishmania major antibody raised against T. brucei RAB11 which also cross-reacts with the L. major homologous . Rab14, appears to regulate endo- lysosomal fusion in Dictyostelium  primarily found on the contractile vacuole network of membranes, an organelle important in osmotic regulation and, in lesser amounts, in the endo-lysosomal pathway .
Rabs in intracellular parasites like Leishmania are associated with pathogenesis, as vesicle transport is essential for extracellular nutrient acquisition, release of virulence factors, microbicidal resistance and evasion of host immune responses , these requirements being different in the monogenetic ancestor Leptomonas hence the difference in phylogeny between Leishmania and Leptomonas for some Rabs like Rab 11 and Rab 14.
In our study Rab 4 and 14 in Leishmania are paralogues and possibly associated with endosomal recycling. Rabs sharing similar functions and/or locations are more closely related than Rab proteins with distinct functions . Rab3 and Rab27, two examples of tissue-restricted Rabs  with specialized functions are rightly not detected in our genomes. Rabs are involved in resistance to chemotherapeutic drugs . Rab4 interacts with P-gp, (P-glycoprotein), a large transmembrane protein localized at the plasma membrane that extrudes anticancer drugs and thus decreases their concentration, and toxicity inside the cells . Rab8 has also been reported to be involved in resistance .
Flow cytometry experiments showed that sensitivity to doxorubicin was associated with increased drug accumulation in cells expressing the RAB6A variant. WTH3 gene’s product is homologous to the Rab6 and Rab6c genes . WTH3 is a house-keeping gene and its product is capable of binding to GTP molecules. WTH3  and Rab5  both are important genes involved in the cellular multi drug resistance (MDR) phenotype development. It is a well established fact that the H locus of Leishmania codes for a short chain dehydrogenase (SDR) gene that is involved in antilfolate resistance . We found SDR to be present in Rab4, Rab6 and Rab18 of Leishmania but not in the corresponding Rab genes of Leptomonas. A study in our lab has provided evidence supporting the important role played in drug resistance by the Rab6 gene of Leishmania having SDR at C terminus of gene (manuscript in communication).
Structure activity relationship of L. donovani GGTase-II with specific inhibitor
Since most prenylated proteins in cells are modified by geranylgeranyl moiety rather than farnesyl  it offers Rab prenylation as a lucrative target for chemotherapeutic intervention in Leishmania cell. Impairment of Rabs (a GGTase-II substrate) geranylgeranylation by isoprenoid pathway inhibitors (lovastatin, zoledronate and DGBP) have been demonstrated . NE10790 is a phosphonocarboxylate analogue of the potent bisphosphonate risedronate. NE10790 inhibited incorporation of [(14)C] mevalonic acid into Rab but not into H-Ras or Rap1, proteins that are modified by FTase and GGTase I, respectively . NE10790 therefore appears to be the first specific inhibitor of Rab prenylation.
The active site architecture of GGTase-II (β-subunit) of L. donovani was compared with other organisms by alignment. The protein sequences of GGTase-II (β-subunit) of Leishmania donovani were aligned with the GGTase-II (β-subunit) sequences of Homo sapiens (NP_004573.2), Plasmodium falciparum (ACT90629.1), Trypanosoma brucei (AAX80948.1) and Drosophila melanogaster (AAF51183.1). Sequence similarity was seen within the kinetoplastid Trypanosoma brucei 65 %, Mycobacterium tuberculosis which is also an intracellular parasite like Leishmania, however, did not show any presence of GGTase-II enzyme on data mining. Available literature indicated that Mycobacterium is predicted to express secreted effectors containing CAAX motifs that may be prenylated by host protein prenyltransferases . While Drosophila melanogaster, Homo sapiens and Plasmodium falciparum were showed 52, 48 and 29 % sequence similarity to GGTase-II (β-subunit) of Leishmania donovani respectively (Fig. 8).
Template for GGTase-II (β-subunit) of Leishmania donovani was selected on the basis of sequence similarity and lowest E-value. Blastp analysis against PDB database, GGTase-II (β-subunit) of Rattus norvegicus shows 48 % identity and minimum E-value with GGTase-II (β-subunit) of Leishmania donovani and this template was further used for molecular modeling studies.
We have ascertained that the β-subunit of GGTase-II in both Leishmania and Leptomonas lie in close proximity done through Neighbor-Joining method (Fig. 5). Therefore the β-subunit emerges as the excellent target for compounds inhibiting parasite activity in clinical cases of co-infections. The crystal structure of Rab escort protein-1 in complex with geranylgeranyltransferase-II (Fig. 9a) and isoprenoid from Rattus showed 48 % of sequence similarity (Fig. 9b) with β-subunit of Leishmania GGTase-II thus selected as the template for homology modeling . Ten models were generated with Modeller. The validation of the resulting models was performed using the SAVS server and the best validated model was selected . In the selected model, the majority of the residues (89.4 %) occupy the most favored region of Ramachandran Plot generated by PROCHECK and 9.9, 0.7 and 0.0 % residues lie in additional allowed region (yellow), generously allowed (light yellow) and disallowed (white) region respectively (Fig. 10a). The wiring diagram of β-subunit of LdGGTase-II was performed using the PDBsum server (http://www.ebi.ac.uk/thornton-srv/databases/pdbsum/); a schematic diagram showed the protein’s secondary structure elements (alpha-helices and beta-sheets) together with various structural motifs such as beta- and gamma-turns, and beta-hairpins. Helices are labeled as H1, H2, etc. (Fig. 10b). The superimposition of the modeled protein (L. donovani RabGGTase) with template 1LTX shows an RMSD of 0.205 Å (Fig. 10c). The farnesyl binding site was selected as the active site for docking of the inhibitor. The NE10790 has a docking score of −6.10 kcal/mol. The inhibitor is involved in formation of three hydrogen bonds with Asn37 and one hydrogen bond with Trp41 as shown by binding mode (Fig. 10d).
By data mining of the prenylation pathway from the genomes of L. donovani and Leptomonas, we have tried to analyse the evolutionary and functional similarities between both these kinetoplastids (Fig. 11). When compared to Rab proteins defined and elucidated in Homo sapiens, the presence of Rabs in Leishmania suggests that Leishmania had a sophisticated endo membrane system, which provides its ability to survive in various intracellular and extracellular stages. This detailed analysis of Rab proteins has provided information of degree of complexity of exo and endocytic trafficking in this organism. The presence of SDR motif known to be associated with drug resistance in several Rab genes of Leishmania indicates that as the parasite evolved from monogenetic to digenetic lifestyle it acquired capabilities to evade drug action for its survival in host. The disruption of prenylation pathway of the parasite using highly specific inhibitors will disable the parasite to reside and survive intracellularly within host phagosome microenvironment and prevent it from modulating host defense pathways for its survival.
Whole genome sequence
L. major and L. infantum genomes are available at (www.ncbi.nlm.nih.gov) with accession number GCA_000002725.2 and GCA_000002875.2 respectively.
Whole genome sequence (WGS) of kinetoplastid species has been deposited by us for Leishmania donovani (BHU 1220, AVPQ01000001) http://www.ncbi.nlm.nih.gov/Traces/wgs/?val=AVPQ01000001and Leptomonas (BHU 1095, ANAF00000000.1) (http://www.ncbi.nlm.nih.gov/Traces/wgs/?val=ANAF01#contigs) .
Kinetoplastid genome phylogeny
Following genomes sequences were retrieved from TritrypDB database.
L. sp. MAR LEM2494
For whole-genome multiple sequence alignments of all these 15 genomes were carried out using Mugsy (http://mugsy.sf.net). Core segments of the alignment that are shared among all isolates included in the analysis were identified and concatenated. The phylogenetic analysis was carried out by considering BHU1220 as reference to the other 14 genome using Phylomark and dendrogram.
Database extraction of geranylgeranyltransferase from whole genome sequence
Rab geranylgeranyltransferase is a heterodimer, consisting of α and β subunits. When we searched for GGTase-I and FTase –I or its homologous sequences in the genomes of L. major, L. infantum, L. donovani and Leptomonas, we found no GGTase-I sequences (EC. 126.96.36.199) in any of the Leishmania species but putative farnesyltransferase sequences (EC. 188.8.131.52) of Leishmania major was found on chromosome 31 (Lmj F31.2940) with accession no. XP_001685284, and in Leishmania infantum also it is located on chromosome 31 (LinJ_31_3050) with accession no. XP_00146718.
When we searched for GGTase-II and its homologous sequences in Leishmania, we found GGTasae-II equences but it was putative and only β-subunit was annotated. In Leishmania major it was located on chromosome 34 (Lmj F34.4030) with accession no. XP_001686510, in Leishmania infantum it was located on chromosome 34 (LinJ_34_3860) with accession no. XP_001468743 and in Leishmania donovani it was also located on chromosome 34 (LdBPK_343860) with accession no. XP_003864545. In Leishmania donovani there is only one protein of 330 amino acids (Accession no. XP_003864545), annotated as putative Rab-GGTase.
For annotation of α-subunit of GGTase-II of L. donovani, it was located on chromosome 33 with accession no XM_001468149.
When we searched for GGTase-II and its homologous sequences in Leptomonas genome we found no α-subunit of GGTase-II but β-subunit of GGTase-II was obtained. To annotate the β-subunit of GGTase-II in Leptomonas, the DNA sequence was retrieved from the contig (genes_5507) of its genome and was put in the translate tools of Expasy. Open reading frame of the β-subunit of GGTase-II from Leptomonas was thus obtained. Blastp analysis of β-subunit of GGTase-II ORF was then done by NCBI and we obtained the sequences producing significant alignment with L. panamensis (XP_010698581), L. braziliensis (XP_001564646) and L. major (XP_001686510).
Multiple sequence alignments of farnesyltransferase and geranylgeranyltransferase were performed with ClustalW. The evolutionary history was inferred using the Neighbor-Joining method . Phylogenetic reconstruction was done using 1000 bootstraps and evolutionary distances were computed using the JTT matrix-based method [48, 49], are in the units of the number of amino acid substitutions per site. The analysis involved 8 amino acid sequences. All positions containing gaps and missing data were eliminated. Evolutionary analysis were conducted in MEGA5 . Trees were drawn using in FigTree v. 1.0.
Rab sequence mining in L. major and L. infantum
Initially in first set we extracted all the annotated H. sapiens Rab sequences reported previously from NCBI . Based upon their annotation they were grouped in families, subfamilies and isoforms. Then each Rab protein sequence was used as query to blastp against the NCBI protein database of L. major and L. infantum, a cutoff for inclusion of protein in dataset was set to be having > 30 % similarity with the query sequence. For generation of second dataset all the sequences retrieved in first dataset being already well annotated as a Rab protein of a family were grouped and Clustal was done to find the conserved domains in these Rab families .
Further analysis was done by NCBI conserved domain search (http://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi) and motif scan (http://hits.isb-sib.ch/cgi-bin/PFSCAN) server to check the conserved domains in these proteins. Based on these results a consensus profile sequence for each family of the Rab was generated which included conserved domains, these consensus sequences were used as query sequence in PSI-BLAST against L. major and L. infantum database at NCBI to retrieve Rab family sequences. Similarly a cutoff of >30 % similarity was set for including sequence in the dataset. Sequences in all the datasets were having similarity more than cut off value were retrieved in fasta format, they were checked for redundancy and duplication or incorrect annotation with rejection of sequence from dataset if it showed 100 % redundancy with any other sequence. The protein datasets were treated separately using the jalview software . Then all the Rab groups in the datasets were compared with each other for redundancy at sequence level. The sequences thus obtained were then fed to the PrePS server (http://mendel.imp.ac.at/PrePS/index.html)  to check that proteins can be predicted to be of Rab family and gets prenylated by GGTase-II enzyme.
Sequence similarity between Rab proteins of L. major and L. infantum was calculated by MatGAT (Matrix Global Alignment Tool. 2.02) software. It is a simple, easy to use computer application that generates similarity/identity matrices for protein sequences without needing pre-alignment of the data .
Rab sequence mining in L. donovani and Leptomonas
The Rabs in L. donovani were retrieved from its genome based on homology search for the corresponding protein sequences from the reference genome of L. infantum. 17 Rabs were identified from its genome as Rabs and Rabs like sequences; whereas 13 out of 17 Rabs were fully annotated (Additional file 2: Table S2).
The Rabs in Leptomonas were retrieved from the contigs of its genome and were based on homology search for the corresponding protein sequences from the reference genome of L. infantum. To annotate the Rabs in these contigs, the DNA sequence of these contigs was put in the translate tools of Expasy. Open reading frames of the Rab proteins of Leptomonas were thus obtained. Blastp analysis of these Rab ORFs was then done by NCBI and we obtained the Rab sequences producing significant alignment with L. infantum. We were able to retrieve fully annotated 11 Rabs out of 53 contigs (Additional file 2: Table S3).
Multiple sequence alignment and phylogenetic reconstruction
The protein sequences were aligned by using ClustalW algorithm in MEGALIGN, DNASTAR package (DNASTAR, Inc. Madison WI, USA). Phylogenetic trees were generated with MEGALIGN which implements a Maximum Likelihood probabilistic model, using standard parameters (gap penalty of 15 and gap length penalty of 6.66) and 1000 bootstraps. The graphically enhanced alignment was obtained using ESPript 2.2 server. The branch lengths are drawn proportional to the evolutionary distances. The numbers at the nodes show the bootstrap support for that node.
Molecular modeling studies of GGTase-II with specific drug target molecule
Since Rab proteins play critical roles in all aspects of intracellular membrane trafficking, this makes them particularly attractive therapeutic targets in diseases. Proper function of Rab proteins depends upon correct membrane localization which is achieved through geranylgeranylation. The quest to determine agents which disrupt Rab geranylgeranylation therefore represents an intriguing therapeutic strategy by which to induce cellular stress and apoptosis. Two lines of attack may achieve diminished cellular levels of Rab geranylgeranylation. One approach would be to deplete cells of the isoprenoid substrate geranylgeranyl diphosphate. This might be achieved by inhibition of any of the individual steps in the isoprenoid biosynthetic pathway that lead to geranylgeranyldiphosphate (GGDP), from inhibition of HMGCoA reductase by a statin, to inhibition of GGPP synthase. This suggests that a long, lipophilic tail can enhance the potency of potential GGTase-II inhibitors, presumably through interaction with the enzyme site that holds the tail of the natural substrate geranylgeranyl diphosphate. The second direct strategy involves inhibition of GGTase-II itself. While inhibitors of GGTase-II are rare, a few have been reported . In silico analysis we used NE10790 for its ability to inhibit Leishmania GGTase-II using GGPP, Rab6 substrate and REP. Our results show that NE10790 is an active compound which bears a geranylgeranyl chain.
Our search for Leishmania GGTase-II crystal structure in databases showed that there are no crystal structures of Leishmania GGTase-II for structural studies, so we carried out the exercise of molecular modeling of Leishmania donovani GGTase-II. Three-dimensional (3D) homology model of Leishmania donovani geranylgeranyltransferase-II (LdGGTase-II) was generated by comparative modeling using MODELLER 9v8  with crystal structure of Rab escort protein-1 in complex with Rab geranylgeranyltransferase and isoprenoid from Rattus norvegicus (PDB ID-1LTX) was selected as a template for structure modeling. ClustalW server was used to align the query sequence with the template sequence . The graphically enhanced alignment was generated with the help of ESPript 2.2 server . All models were assessed stereo-chemically by PROCHECK . All the graphical visualization and image production were performed using Chimera . The docking studies were performed using Autodock4.2 software . The structure of the inhibitor NE10790 was drawn using sketch module of Sybyl7.1 and geometry was optimized by performing energy minimization using MMFF94 force field with 1000 iterations (Sybyl, version 7.1; Tripos, Inc.: St. Louis, MO, 2005).
Availability of supporting data
Whole genome sequence (WGS) of Leishmania donovani (BHU 1220, AVPQ01000001) and Leptomonas (BHU 1095, ANAF00000000.1) supporting the results of this article is available at http://www.ncbi.nlm.nih.gov/Traces/wgs/?val=AVPQ01000001 http://www.ncbi.nlm.nih.gov/Traces/wgs/?val=ANAF01#contigs.
All the supporting data are included as additional files.
Whole genome sequencing
Rab escort protein
Protein Data Bank
Last eukaryotic common ancestor
Actin interacting protein
Short-chain dehydrogenase/ reductase
Open reading frame
Singh N, Chikara S, Sundar S. SOLiD™ sequencing of genomes of clinical isolates of Leishmania donovani from India confirm leptomonas co-infection and raise some key questions. PLoS One. 2013;8(2), e55738.
Srivastava P, Prajapati VK, Vanaerschot M, Van der Auwera G, Dujardin JC, Sundar S. Detection of Leptomonas sp. parasites in clinical isolates of Kala-azar patients from India. Infect Genet Evol. 2010;10(7):1145–50.
Ghosh S, Banerjee P, Sarkar A, Datta S, Chatterjee M. Coinfection of Leptomonas seymouri and Leishmania donovani in Indian Leishmaniasis. J Clin Microbiol. 2012;50(8):2774–8.
Garnham P. The genus leishmania. Bull World Health Organ. 1971;44(4):477.
Baker J. The evolution of parasitic protozoa. Evo Parasites. 1965;1965:1–28.
Killick-Kendrick R, Ward R. Leishmania in phlebotomid sandflies. II. The insusceptibility of sandfly larvae to leishmania. Proc R Soc Lond Ser B Biol Sci. 1975;188(1091):229–31.
Lynn MA, McMaster WR. Leishmania: conserved evolution–diverse diseases. Trends Parasitol. 2008;24(3):103–5.
Rojas AM, Fuentes G, Rausell A, Valencia A. The Ras protein superfamily: evolutionary tree and role of conserved amino acids. J Cell Biol. 2012;196(2):189–201.
Gruenberg J, Van der Goot FG. Mechanisms of pathogen entry through the endosomal compartments. Nat Rev Mol Cell Biol. 2006;7(7):495–504.
Gurkan C, Koulov AV, Balch WE. An evolutionary perspective on eukaryotic membrane trafficking, Eukaryotic Membranes and Cytoskeleton. New York: Springer; 2007. p. 73–83.
Rasteiro R, Pereira-Leal JB. Multiple domain insertions and losses in the evolution of the Rab prenylation complex. BMC Evol Biol. 2007;7(1):140.
Peterson DS, Milhous WK, Wellems TE. Molecular basis of differential resistance to cycloguanil and pyrimethamine in Plasmodium falciparum malaria. Proc Natl Acad Sci. 1990;87(8):3018–22.
Malo D, Vogan K, Vidal S, Hu J, Cellier M, Schurr E, et al. Haplotype mapping and sequence analysis of the mouse Nramp gene predict susceptibility to infection with intracellular parasites. Genomics. 1994;23(1):51–61.
Barbacid M. Ras genes. Annu Rev Biochem. 1987;56(1):779–827.
Bourne HR, Sanders DA, McCormick F. The GTPase superfamily: a conserved switch for diverse cell functions. Nature. 1990;348(6297):125–32.
Bourne HR, Sanders DA, McCormick F. The GTPase superfamily: conserved structure and molecular mechanism. Nature. 1991;349(6305):117–27.
Pereira-Leal JB, Seabra MC. The mammalian Rab family of small GTPases: definition of family and subfamily sequence motifs suggests a mechanism for functional specificity in the Ras superfamily. J Mol Biol. 2000;301(4):1077–87.
Patel N, Singh SB, Basu SK, Mukhopadhyay A. Leishmania requires Rab7-mediated degradation of endocytosed hemoglobin for their growth. Proc Natl Acad Sci. 2008;105(10):3980–5.
Chenik M, Chaabouni N, Achour-Chenik YB, Ouakad M, Lakhal-Naouar I, Louzir H, et al. Identification of a new developmentally regulated Leishmania major large RAB GTPase. Biochem Biophys Res Commun. 2006;341(2):541–8.
Singh SB, Tandon R, Krishnamurthy G, Vikram R, Sharma N, Basu SK, et al. Rab5‐mediated endosome–endosome fusion regulates hemoglobin endocytosis in Leishmania donovani. EMBO J. 2003;22(21):5712–22.
Berriman M, Ghedin E, Hertz-Fowler C, Blandin G, Renauld H, Bartholomeu DC, et al. The genome of the African trypanosome Trypanosoma brucei. Science. 2005;309(5733):416–22.
Pereira-Leal JB, Seabra MC. Evolution of the Rab family of small GTP-binding proteins. J Mol Biol. 2001;313(4):889–901.
Mackiewicz P, Wyroba E. Phylogeny and evolution of Rab7 and Rab9 proteins. BMC Evol Biol. 2009;9(1):101.
Ackers JP, Dhir V, Field MC. A bioinformatic analysis of the RAB genes of Trypanosoma brucei. Mol Biochem Parasitol. 2005;141(1):89–97.
Field MC, Natesan SKA, Gabernet‐Castello C, Lila Koumandou V. Intracellular trafficking in the trypanosomatids. Traffic. 2007;8(6):629–39.
Klöpper TH, Kienle N, Fasshauer D, Munro S. Untangling the evolution of Rab G proteins: implications of a comprehensive genomic analysis. BMC Biol. 2012;10(1):71.
Diniz MC, Costa MP, Pacheco AC, Kamimura MT, Silva SC, Carneiro LD, et al. Actin-interacting and flagellar proteins in Leishmania spp.: Bioinformatics predictions to functional assignments in phagosome formation. Genet Mol Biol. 2009;32(3):652–65.
Real F, Mortara RA. The diverse and dynamic nature of Leishmania parasitophorous vacuoles studied by multidimensional imaging. PLoS Negl Trop Dis. 2012;6(2), e1518.
Bucci C, Parton RG, Mather IH, Stunnenberg H, Simons K, Hoflack B, et al. The small GTPase rab5 functions as a regulatory factor in the early endocytic pathway. Cell. 1992;70(5):715–28.
Ghosh SS, Mukerjee S, Adhya S. Chromosome profile of Leishmania donovani: Interstrain and interspecific variations. J Biosci. 1998;23(3):247–54.
Marotta DE, Gerald N, Dwyer DM. Rab5b localization to early endosomes in the protozoan human pathogen Leishmania donovani. Mol Cell Biochem. 2006;292(1–2):107–17.
Ambit A, Woods KL, Cull B, Coombs GH, Mottram JC. Morphological events during the cell cycle of Leishmania major. Eukaryot Cell. 2011;10(11):1429–38.
Bush J, Temesvari L, Rodriguez-Paris J, Buczynski G, Cardelli J. A role for a Rab4-like GTPase in endocytosis and in regulation of contractile vacuole structure and function in Dictyostelium discoideum. Mol Biol Cell. 1996;7(10):1623–38.
Harris E, Cardelli J. RabD, a Dictyostelium Rab14-related GTPase, regulates phagocytosis and homotypic phagosome and lysosome fusion. J Cell Sci. 2002;115(18):3703–13.
Frei SB, Rahl PB, Nussbaum M, Briggs BJ, Calero M, Janeczko S, et al. Bioinformatic and comparative localization of Rab proteins reveals functional insights into the uncharacterized GTPases Ypt10p and Ypt11p. Mol Cell Biol. 2006;26(19):7299–317.
Recchi C, Seabra MC. Rab GTPases and their interacting proteins in health and disease: novel functions for Rab GTPases in multiple aspects of tumour progression. Biochem Soc Trans. 2012;40(Pt 6):1398.
Shan J, Mason JM, Yuan L, Barcia M, Porti D, Calabro A, et al. Rab6c, a new member of the rab gene family, is involved in drug resistance in MCF7/AdrR cells. Gene. 2000;257(1):67–75.
Tian K, Wang Y, Huang Y, Sun B, Li Y, Xu H. Methylation of WTH3, a possible drug resistant gene, inhibits p53 regulated expression. BMC Cancer. 2008;8(1):327.
Fu D, van Dam EM, Brymora A, Duggin IG, Robinson PJ, Roufogalis BD. The small GTPases Rab5 and RalA regulate intracellular traffic of P-glycoprotein. Biochim Biophys Acta, Mol Cell Res. 2007;1773(7):1062–72.
Papadopoulou B, Roy G, Ouellette M. Frequent amplification of a short chain dehydrogenase gene as part of circular and linear amplicons in methortexate resistant Lesihmania. Nucleic Acids Res. 1993;21(18):4305–12.
Farnsworth CC, Gelb MH, Glomset JA. Identification of geranylgeranyl-modified proteins in HeLa cells. Science. 1990;247(4940):320–2.
Dudakovic A, Wiemer AJ, Lamb KM, Vonnahme LA, Dietz SE, Hohl RJ. Inhibition of geranylgeranyl diphosphate synthase induces apoptosis through multiple mechanisms and displays synergy with inhibition of other isoprenoid biosynthetic enzymes. J Pharmacol Exp Ther. 2008;324(3):1028–36.
Coxon FP, Helfrich MH, Larijani B, Muzylak M, Dunford JE, Marshall D, et al. Identification of a novel phosphonocarboxylate inhibitor of Rab geranylgeranyl transferase that specifically prevents Rab prenylation in osteoclasts and macrophages. J Biol Chem. 2001;276(51):48213–22.
Zverina EA, Lamphear CL, Wright EN, Fierke CA. Recent advances in protein prenyltransferases: Substrate identification, regulation, and disease interventions. Curr Opin Chem Biol. 2012;16(5):544–52.
Pylypenko O, Rak A, Reents R, Niculae A, Sidorovitch V, Cioaca M-D, et al. Structure of Rab escort protein-1 in complex with Rab geranylgeranyltransferase. Mol Cell. 2003;11(2):483–94.
Laskowski RA, MacArthur MW, Moss DS, Thornton JM. PROCHECK: a program to check the stereochemical quality of protein structures. J Appl Crystallogr. 1993;26(2):283–91.
Saitou N, Nei M. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987;4(4):406–25.
Felsenstein J. Confidence limits on phylogenies: an approach using the bootstrap. Evolution. 1985;783–791.
Jones DT, Taylor WR, Thornton JM. The rapid generation of mutation data matrices from protein sequences. CABIOS. 1992;8(3):275–82.
Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011;28(10):2731–9.
Larkin MA, Blackshields G, Brown N, Chenna R, McGettigan PA, McWilliam H, et al. Clustal W and Clustal X version 2.0. Bioinformatics. 2007;23(21):2947–8.
Waterhouse AM, Procter JB, Martin DM, Clamp M, Barton GJ. Jalview Version 2—a multiple sequence alignment editor and analysis workbench. Bioinformatics. 2009;25(9):1189–91.
Maurer-Stroh S, Eisenhaber F. Refinement and prediction of protein prenylation motifs. Genome Biol. 2005;6(6):R55.
Campanella JJ, Bitincka L, Smalley J. MatGAT: an application that generates similarity/identity matrices using protein or DNA sequences. BMC Bioinf. 2003;4(1):29.
Zhou X, Hartman SV, Born EJ, Smits JP, Holstein SA, Wiemer DF. Triazole-based inhibitors of geranylgeranyltransferase II. Bioorg Med Chem Lett. 2013;23(3):764–6.
Šali A, Blundell TL. Comparative protein modelling by satisfaction of spatial restraints. J Mol Biol. 1993;234(3):779–815.
Chenna R, Sugawara H, Koike T, Lopez R, Gibson TJ, Higgins DG, et al. Multiple sequence alignment with the Clustal series of programs. Nucleic Acids Res. 2003;31(13):3497–500.
Gouet P, Courcelle E, Stuart DI. ESPript: analysis of multiple sequence alignments in PostScript. Bioinformatics. 1999;15(4):305–8.
Pettersen EF, Goddard TD, Huang CC, Couch GS, Greenblatt DM, Meng EC, et al. UCSF Chimera—a visualization system for exploratory research and analysis. J Comput Chem. 2004;25(13):1605–12.
Morris GM, Huey R, Lindstrom W, Sanner MF, Belew RK, Goodsell DS, et al. AutoDock4 and AutoDockTools4: Automated docking with selective receptor flexibility. J Comput Chem. 2009;30(16):2785–91.
Real F, Vidal RO, Carazzolle MF, Mondego JMC, Costa GGL, Herai RH, et al. The genome sequence of Leishmania (Leishmania) amazonensis: functional annotation and extended analysis of gene models. DNA Res. 2013;20(6):567–81.
We are grateful to Mr. Arun Kumar for bioinformatic analysis. This work was supported by Council of Scientific and Industrial Research (CSIR) funded network project “Host Interactome analysis: Understanding the Role of Host molecules in Parasitic Infection (HOPE)”. This is CDRI manuscript number −159/2015/NS.
The authors declare that they have no competing interests.
NS conceived the study. NS and ISC designed the study. ISC performed the experiments. SK performed the molecular modeling studies of GGTase-II. MIS analyzed the molecular modeling data. AG performed and analyzed the genome phylogeny. PS analyzed the chemical entities. NS, ISC and JK analyzed the data. NS and ISC wrote the manuscript. All authors read and approved the final manuscript.
Additional file 1: Figure S1.
Multiple sequence alignment of FTase (α-subunit) of Leishmania donovani (accession number XP_003862625.1, putative protein) with FTase (α-subunit) of Leishmania major (accession number XP_003722277.1, putative protein), Leishmania infantum (accession number XP_001466722.1, putative protein) and its sister Leptomonas (contig_2654, putative protein) showing conserved regions (red colored) based on percentage identity. Figure S2. Multiple sequence alignment of FTase (β-subunit) of Leishmania donovani (accession number XP_003861732.1) with FTase (β-subunit) of Leishmania major (accession number XP_001684151), Leishmania infantum (accession number XP_001470492) and its sister Leptomonas (contig_1135) showing conserved regions (red colored) based on percentage identity. Figure S3. Multiple sequence alignment of GGTase-II (α-subunit) of Leishmania donovani (accession number XM_001468149) with GGTase-II (α-subunit) of Leishmania major (accession number XM_001685808) and Leishmania infantum (accession number XP_001468186) showing conserved regions (red colored) based on percentage identity. Its sister Leptomonas has no α -subunit of GGTase-II. Figure S4. Multiple sequence alignment of GGTase-II (β-subunit) of Leishmania donovani (accession number XP_003864545) with GGTase-II (β-subunit) of Leishmania major (accession number XP_001686510) and Leishmania infantum (accession number XP_001468743) showing conserved regions (red colored) based on percentage identity. Figure S5. Alignment of various Leishmania donovani Rab protein sequences showing sequence similarities and domain conservation among them colored red based in their percentage similarity. Figure S6. Alignment of various Leptomonas Rab protein sequences showing sequence similarities and domain conservation among them colored red based in their percentage similarity. Figure S7. Phylogenetic relationship among Rab sequences of L. major, L. infantum, L. donovani and Leptomonas which is based on multiple sequence alignment (ClustalW algorithm in MEGALIGN with 1000 bootstraps). Figure S8. Circos plot was generated using Circos version 0.60, which is used to determine the similarity between Rab genes of Leishmania donovani and its sister Leptomonas. Where cyan color for all Rabs protein obtained from Leptomonas genome and green color for all Rabs protein obtained from Leishmania donovani genome. Other colors for each different Rabs. (ZIP 81526 kb)
Additional file 2: Table S1.
L. major and L. infantum Rabs with their accession numbers, gene IDs and chromosome location. Table S2. Rab GTPase sequences identified from Leishmania donovani. Table S3. Rab GTPase sequences identified from Leptomonas. Table S4. (A) Heat-map depicting percent identity between L. donovani and (B) Leptomonas. (DOC 102 kb)
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Chauhan, I.S., Kaur, J., Krishna, S. et al. Evolutionary comparison of prenylation pathway in kinetoplastid Leishmania and its sister Leptomonas . BMC Evol Biol 15, 261 (2015). https://doi.org/10.1186/s12862-015-0538-3
- Leishmania donovani
- Prenylation pathway