- Short report
- Open Access
Overview of the organization of protease genes in the genome of Leishmania spp
Parasites & Vectorsvolume 7, Article number: 387 (2014)
The genus Leishmania includes protozoan parasites that are able to infect an array of phlebotomine and vertebrate species. Proteases are related to the capacity of these parasites to infect and survive in their hosts and are therefore classified as virulence factors.
By analyzing protease genes annotated in the genomes of four Leishmania spp [Leishmania (Leishmania) infantum, L. (L.) major, L. (L.) mexicana and L. (Viannia) braziliensis], these genes were found on every chromosome of these protozoa. Four protease classes were studied: metallo-, serine, cysteine and aspartic proteases. Metalloprotease genes predominate in the L. (V.) braziliensis genome, while in the other three species studied, cysteine protease genes prevail. Notably, cysteine and serine protease genes were found to be very abundant, as they were found on all chromosomes of the four studied species. In contrast, only three aspartic protease genes could be detected in these four species. Regarding gene conservation, a higher number of conserved alleles was observed for cysteine proteases (42 alleles), followed by metalloproteases (35 alleles) and serine proteases (15 alleles).
The present study highlights substantial differences in the organization of protease genes among L. (L.) infantum, L. (L.) major, L. (L.) mexicana and L. (V.) braziliensis. We observed significant distinctions in many protease features, such as occurrence, quantity and conservation. These data indicate a great diversity of protease genes among Leishmania species, an aspect that may be related to their adaptations to the peculiarities of each microenvironment they inhabit, such as the gut of phlebotomines and the immune cells of vertebrate hosts.
The World Health Organization classifies the leishmaniases, infections caused by parasites of the genus Leishmania, among emerging diseases that lack effective control. Annually, an estimated 1.3 million new cases occur and 20,000 to 30,000 deaths are attributed to these diseases. The clinical forms range in severity and are classified as follows: punctuate skin lesions to oronasal disfigurement are classified as cutaneous leishmaniasis (CL), whereas fatal systemic infections are classified as visceral leishmaniasis (VL). Leishmania spp are distributed worldwide and are organized into subgenera and species complexes. Their transmission to mammalian hosts occurs during the blood meal of infected sandflies, which in turn acquire the parasites when feeding on an infected host, thus maintaining the cycle of the disease. The species grouped into the Leishmania (Leishmania) donovani complex, including L. (L.) infantum, are the agents of VL. As for the species commonly associated with CL, L. (L.) major is reported in the Old World, whereas L. (L.) mexicana and L. (Viannia) braziliensis are the main species reported in the New World. This latter species is also associated with the mucocutaneous form of the disease.
In a recent review study, we have highlighted the pivotal roles of proteases as virulence factors for Leishmania spp. Such enzymes have been implicated in many parasitic activities, such as tissue invasion, survival in macrophages and host immune response modulation.
Proteases are classified according to their physicochemical features as: optimal pH for activity, kind of catalytic activity, nature of catalytic site and homology with reference structure. According to the enzymatic nomenclature committee, the Joint Commission on Biochemical Nomenclature (JCBN), peptidases are allocated into the Enzyme Class (EC) 3 (hydrolases) and subclass 3.4 (peptidases). They can be subdivided into exopeptidases (EC 4.11 - EC.4.19) and endopeptidases (EC 3.4.21 - 3.4.25), and the latter are organized according to the amino acids related to catalysis and the nature of catalytic site. In addition, endopeptidases are further divided into classes according to the main catalytic mechanism involved in their hydrolytic activities, e.g., serine, threonine, aspartate, metallo- and cysteine proteases.
The aim of the present study is to analyse the genomic organization of proteases in four Leishmania species known to cause disease in humans: L. (V.) braziliensis, L. (L.) major, L. (L.) mexicana and L. (L.) infantum, and, concomitantly, to evaluate their diversity among these species. Due to the importance of these enzymes in the life cycle of these parasites, the genomic data gathered here would be very useful as a basis for further studies correlating infection characteristics of each of the studied species with their protease richness. Understanding how these enzymes are organized and conserved (or diverged) in the different Leishmania subgenera and species is very useful in helping to identify new targets with the most potential for chemotherapy or vaccination strategies.
Findings and discussion
We performed a comparative genomic analysis on the organization of protease genes in four species, a methodology we applied to identify species-specific features that may account for phenotypic or virulence differences among the studied species. Gene divergence, acquisition, loss, and rearrangement within and between syntenic regions have shaped the genomes of the trypanosomatids and can explain the organization and diversity of the degradome (the complete set of protease genes encoded by the genome of a certain organism) of Leishmania spp. Initially, we performed a survey of the predicted protease sequences present in the annotated genomes of L. (V.) braziliensis, L. (L.) major, L. (L.) mexicana and L. (L.) infantum in the GeneDB genome database. This survey was conducted using the following keywords: protease, peptidase, proteinase, aspartic protease, cysteine protease, serine protease and metalloprotease.
In an initial analysis of the data retrieved by the methodology above, the abundance of protease genes in the genomes of each of the studied species was defined. While protease genes account for 2.18% of the total genes in L. (V.) braziliensis, these genes account for smaller percentages in the other species: 1.61% in L. (L.) infantum, 1.52% in L. (L.) mexicana and 1.41% in L. (L.) major.
Metalloprotease genes predominate in L. (V.) braziliensis, while in the other species the cysteine protease genes prevail. Our analysis showed that 52% of the protease genes in L. (V.) braziliensis are metalloproteases and this same class accounts for 40% of protease genes in L. (L.) infantum and 35% in L. (L.) major and L. (L.) mexicana. The percentages of cysteine and serine protease gene are close among the studied species: cysteine protease genes represent 36 to 47% of the total protease genes, whereas serine protease genes represent 10 to 16%. Very few aspartic protease genes were identified, amounting to only three in each of the four species (Figure 1).
A very interesting finding we observed is that protease genes are present in every chromosome of the studied Leishmania spp, but occur in different frequencies (Figure 1). This discovery is consistent with the previously reported importance of proteases for these parasitic organisms, as it reveals that genes encoding these enzymes are abundantly scattered among the Leishmania spp genomes, and is, complementarily, an indication of the distinct pattern evolution has impinged over the different species.
Other studies regarding gene organization in Leishmania spp have been conducted before and related the structural configuration of the genes with important functional features. The organization of genes in tandem repeats allows parasites to quickly generate a high number of transcripts that may be needed in large amounts. Other authors hypothesise that Leishmania spp. might have a strategy to increase mRNA levels by duplicating genes on disomic chromosomes or by forming supernumerary chromosomes.
Of the chromosomes that we identified as containing metalloprotease genes, 18 are common for all studied species. Notably, the presence of metalloprotease genes on chromosomes 8 and 30 is exclusive to L. (L.) mexicana. Similar exclusiveness for the presence of metalloprotease genes was observed for chromosome 22 in L. (L.) infantum and chromosome 20 in L. (V.) braziliensis. Cysteine protease genes are present in 22 chromosomes common to all four species studied. Cysteine protease genes are also present on chromosome 7 exclusively in L. (L.) mexicana, on chromosome 28 in L. (V.) braziliensis and on chromosome 35 in L. (L.) infantum and L. (L.) major. Serine protease genes are present in 9 chromosomes common to all four species studied and the number of these genes does not exceed three per chromosome. The presence of genes for this protease class on chromosome 29 is exclusive to L. (L.) major and on chromosome 20 to L. (V.) braziliensis. The protease class found to have the fewest coding genes was aspartic proteases: only three genes for this class were observed, but the chromosomes on which they are present are common to all studied species. These genes are located on chromosomes 1, 15 and 29 (Figure 1).
Regarding genes for different protease classes that occur on the same chromosome, most of the studied chromosomes were found to contain genes for multiple protease classes. The exceptions were chromosomes 3 and 6, which were found to contain only serine protease genes and chromosomes 5, 11 and 22, which were found to contain only metalloprotease genes.
Due to fusion events that occurred in Leishmania chromosomes, we observed an interesting pattern of organization of protease genes where the same arrangement of alleles is maintained across different species but is located on different chromosomes. Graphical representations of such fusion events were developed using the Artemis and ACT software (Additional file1: Figure S1 to S8).
Nevertheless, there is a trend of conservation of some alleles in the same chromosomes across the studied Leishmania species. We observed 42 conserved alleles of cysteine proteases, 35 of metalloproteases and 15 of serine proteases (Figure 2). The conserved alleles are predominantly grouped on chromosome 10 for cysteine proteases, chromosome 30 for metalloproteases and chromosome 28 for serine proteases.
Among all the analysed protease genes, only two alleles were found to be conserved on the same chromosome for all four species: alleles of cysteine protease genes coding for ubiquitin carboxyl-terminal hydrolases (Clan CA, family C12) located on chromosomes 24 and 25 (alleles 0420 and 0190, respectively) of all species.
Notably, L. (L.) major and L. (L.) mexicana were found to show more synteny than the other species, containing 23, 15 and 13 conserved alleles for cysteine, metallo- and serine proteases, respectively. Conversely, L. (V.) braziliensis was not found to show synteny for serine protease genes of any other species. Although this absence of synteny was observed in the only species in our analysis classified into a different subgenus, it has been proposed by Peacock et al. that such absence would not necessarily indicate a lineage-specific diversity in Leishmania spp.
One of the first comparative genomic studies of Leishmania showed that despite phenotypic variations among species, only a few genes are truly species-specific. In agreement with such reports, we also observed few genes that do not show similarity to any others. They show sequence identities lower than 80% to other genes (Additional file2: Table S1). This is an important finding, as these exclusive genes can help explain why these species cause different forms of diseases and are present in specific vectors and hosts. Previously, it was reported that more than 99% of genes are conserved between L. (V.) braziliensis, L. (L.) infantum and L. (L.) major, revealing a high degree of synteny for genomes of different Leishmania species. Our data indicates that, when analysing strictly protease genes, this same scenario holds up, as we also observed high synteny between the studied species.
When contemplating the usefulness of parasite proteases as new targets for chemotherapies, it is very important to consider the hypothesis that these enzymes are unique to the Leishmania species and quite different from corresponding enzymes in their mammalian hosts, such as humans and dogs. Thus, to verify this hypothesis, we conducted a BLAST (Basic Local Alignment Search Tool) analysis to compare the genes that show synteny among the greater number of the four species (represented in the intersection of the Venn diagram – Figure 2) with mammalian protease genes (taxid:40674). The genes 05.0960 and 11.0630 of L. (L.) major, L. (L.) mexicana, L. (L.) major show the highest degree of relational similarity with mammalian genes, with approximately 69% sequence identity and a query coverage of up to 39%. However, in general, the query coverage was very low, with a mean value of 2%. In addition, to perform a similar study with other proteases that did not show synteny among all the studied species, we used a different approach.
Initially, a multiple alignment analysis was carried out on the sequences of protease genes of the four species (software CD-HIT), using a cutoff of 80% sequence identity to cluster them. As result, we were able to establish 28 clusters of metalloprotease genes, 27 of cysteine protease genes, 11 of serine protease genes and 1 of aspartic protease genes.
The consensus sequences (Additional file2: Table S2) of each cluster were then used in the BLAST analysis to find similarity with mammalian genes. We identified sequences of O-sialoglycoprotein endopeptidase genes of hamster, dog, wolf and mice with 69% sequence identity to a consensus sequence of Leishmania metalloprotease genes LbrM.31.0100, LinJ.31.0110, LmjF.31.0100 and LmxM.31.0100. Sequences of 26S subunit ATPase genes of a lagomorph Ochotona princeps and of mice showed 65% sequence identity to a consensus sequence of serine protease genes LbrM.03.0450, LinJ.03.0520, LmjF.03.0540 and LmxM.03.0540. Additionally, we could not find any similarity among sequences of cysteine and aspartic protease genes of mammals and Leishmania spp.
As proteases can be grouped into different families and clans depending on intrinsic evolutionary relationships, we classified and organized the protease genes surveyed in this study applying criteria from MEROPS (up to December 2013) (Figure 3). This classification is based on structural and functional similarities between these proteolytic enzymes. Clans contain enzymes with related structures and families contain enzymes with related sequences. This classification is highly relevant to understanding the organization of these parasites’ degradomes.
Cysteine proteases and metalloproteases are the major representative classes of proteases in this study, corresponding to 43% and 42%, respectively of the protease genes in the studied Leishmania spp. In this survey, three clans of cysteine proteases were observed in the studied species: clan CA, CD and CF. These cysteine proteases are distributed among 11 families from which C1, C2 and C19 have more members. The MPs observed in the study belong to the clans MA, MC, ME, MF, MG, MH and MP and are further distributed among 14 families (Figure 3). The diversity of protease genes observed in the analysis reinforces the idea that this class of enzyme is crucial to the parasite lifecycle, although until now the role of most of these proteases can only be predicted based on current knowledge of homologous enzymes, therefore pointing to the necessity of more studies characterising proteases.
The high number of metalloprotease genes in L. (V.) braziliensis relates to the 36 distinct genes of the zinc metalloprotease gp63. This metalloprotease is a very well-characterised virulence factor for L. (L.) braziliensis and has several reported functions in the interactions of this parasite with its hosts. In L. (L.) major, L. (L.) mexicana and L. (L.) infantum, the diversity of gp63 genes is much lower: only 6, 7 and 8 genes, respectively, of this protease could be found (Figure 3). The organization of metalloprotease genes in species of the subgenus Viannia is rather different than that of species of the subgenus Leishmania. The predominance of metalloprotease genes in L. (V.) braziliensis, a peculiarity also observed in L. (V.) guyanensis, has a biological significance not completely understood[8, 18]. Amplification of genes is a common phenomenon in Leishmania[19–21] and is a likely source of the differences between the two subgenera. Such interesting variation might have fundamental implications for the way each species interacts with its hosts.
Our study highlights the informative potential of analysing genome databases for understanding the gene organization of parasites. However, one should be aware that not all annotated proteases have described roles in the Leishmania life cycle. Thus, the picture observed here is not yet complete.
It is still unclear how the current organization of Leishmania spp genomes evolved, but the set of results gathered here emphasises the capacity of Leishmania species to use the plasticity of their genomes to modulate their phenotypes and increase their odds of survival within hosts, among other biological processes. The diversity of protease genes described by our present study points to their potential importance as survival and adaptation tools and, consequently, as important targets in vaccination and therapy strategies.
WHO: Fact sheet n° 375. 2014, updated January 2014 [http://www.who.int/mediacentre/factsheets/fs375/en/]
Silva-Almeida M, Pereira BA, Ribeiro-Guimarães ML, Alves CR: Proteinases as virulence factors in Leishmania spp. infection in mammals. Parasit Vectors. 2012, 7 (5): 160-
Barrett AJ: Classification of peptidases. Methods Enzymol. 1994, 244: 1-15.
Rawlings ND, Waller M, Barrett AJ, Bateman A: MEROPS: the database of proteolytic enzymes, their substrates and inhibitors. Nucleic Acids Res. 2014, 42 (Database issue): D503-509.
El-Sayed NM, Myler PJ, Blandin G, Berriman M, Crabtree J, Aggarwal G, Caler E, Renauld H, Worthey EA, Hertz-Fowler C, Ghedin E, Peacock C, Bartholomeu DC, Haas BJ, Tran AN, Wortman JR, Alsmark UC, Angiuoli S, Anupama A, Badger J, Bringaud F, Cadag E, Carlton JM, Cerqueira GC, Creasy T, Delcher AL, Djikeng A, Embley TM, Hauser C, Ivens AC: Comparative genomics of trypanosomatid parasitic protozoa. Science. 2005, 309: 404-409. 10.1126/science.1112181.
López-Otín C, Overall CM: Protease degradomics: a new challenge for proteomics. Nat Rev Mol Cell Biol. 2002, 3 (7): 509-519. 10.1038/nrm858.
Victoir K, Dujardin JC: How to succeed in parasitic life without sex? Asking Leishmania. Trends Parasitol. 2002, 18 (2): 81-85. 10.1016/S1471-4922(01)02199-7.
Rogers MB, Hilley JD, Dickens NJ, Wilkes J, Bates PA, Depledge DP, Harris D, Her Y, Herzyk P, Imamura H, Otto TD, Sanders M, Seeger K, Dujardin JC, Berriman M, Smith DF, Hertz-Fowler C, Mottram JC: Chromosome and gene copy number variation allow major structural change between species and strains of Leishmania. Genome Res. 2011, 21 (12): 2129-42. 10.1101/gr.122945.111.
Carver T, Berriman M, Tivey A, Patel C, Bohme U, Barrell BG, Parkhill J, Rajandream MA: Artemis and ACT: viewing, annotating and comparing sequences stored in a relational database. Bioinformatics. 2008, 24: 2672-2676. 10.1093/bioinformatics/btn529.
Peacock CS, Seeger K, Harris D, Murphy L, Ruiz JC, Quail MA, Peters N, Adlem E, Tivey A, Aslett M, Kerhornou A, Ivens A, Fraser A, Rajandream MA, Carver T, Norbertczak H, Chillingworth T, Hance Z, Jagels K, Moule S, Ormond D, Rutter S, Squares R, Whitehead S, Rabbinowitsch E, Arrowsmith C, White B, Thurston S, Bringaud F, Baldauf SL: Comparative genomic analysis of three Leishmania species that cause diverse human disease. Nat Genet. 2007, 7: 839-847. 10.1038/nri2207.
Li W, Jaroszewski L, Godzik A: Clustering of highly homologous sequences to reduce the size of large protein databases. Bioinformatics. 2001, 17 (3): 282-283. 10.1093/bioinformatics/17.3.282.
Mottram JC, Coombs GH, Alexander J: Cysteine peptidases as virulence factors of Leishmania. Curr Opin Microbiol. 2004, 7: 375-381. 10.1016/j.mib.2004.06.010.
Olivier M, Atayde VD, Isnard A, Hassani K, Shio MT: Leishmania virulence factors: focus on the metalloprotease GP63. Microbes Infect. 2012, 14 (15): 1377-1389. 10.1016/j.micinf.2012.05.014.
Voth BR, Kelly BL, Joshi PB, Ivens AC, McMaster WR: Differentially expressed Leishmania major gp63 genes encode cell surface leishmanolysin with distinct signals for glycosylphosphatidylinositol attachment. Mol Biochem Parasitol. 1998, 93 (1): 31-41. 10.1016/S0166-6851(98)00013-9.
Steinkraus HB, Greer JM, Stephenson DC, Langer PJ: Sequence heterogeneity and polymorphic gene arrangements of the Leishmania guyanensis gp63 genes. Mol Bioch Parasitol. 1993, 62: 173-185. 10.1016/0166-6851(93)90107-9.
Victoir K, Arevalo J, De Doncker S, Barker DC, Laurent T, Godfroid E, Bollen A, Le Ray D, Dujardin JC: Complexity of the major surface protease (msp) gene organization in Leishmania (Viannia) braziliensis: evolutionary and functional implications. Parasitology. 2005, 131: 207-214. 10.1017/S0031182005007535.
Iovannisci DM, Beverley SM: Structural alterations of chromosome 2 in Leishmania major as evidence for diploidy, including spontaneous amplification of the mini-exon array. Mol Bioch Parasitol. 1989, 34: 177-188. 10.1016/0166-6851(89)90009-1.
Inga R, De Doncker S, Gomez J, Lopez M, Garcia R, Le Ray D, Arevalo J, Dujardin JC: Relation between variation in copy number of ribosomal RNA encoding genes and size of harbouring chromosomes in Leishmania of subgenus Viannia. Mol Bioch Parasitol. 1998, 92: 219-228. 10.1016/S0166-6851(98)00009-7.
Kebede A, De Doncker S, Arevalo J, Le Ray D, Dujardin JC: Size-polymorphism of mini-exon gene-bearing chromosomes among natural populations of Leishmania, subgenus Viannia. Int J Parasitol. 1999, 29: 549-557. 10.1016/S0020-7519(99)00010-7.
We thank FAPERJ (E-26/102.413/2010, E-26/110.592/2012) and CAPES (23038.007057-98) for partial financial support of this research. Mariana Silva-Almeida and Franklin Souza-Silva are doctoral students of Fiocruz institution; Dr. Bernardo A. S. Pereira and Dr Michelle Lopes Ribeiro-Guimarães are postdoctoral researcher’s fellow of CAPES/FAPERJ; and Dr. Carlos R. Alves is fellow researcher of CNPq institution.
The authors declare that they have no competing interests.
MSA, FSS and CRA formulated the idea and wrote the manuscript; MSA, FSS, and MFM performed the analysis processes. MSA, FSS, CRA, BASP, provided critical comments to the methods and the discussion. All authors approved the final version of this manuscript.
Mariana Silva-Almeida, Franklin Souza-Silva contributed equally to this work.
Electronic supplementary material
Additional file 1: Figure S1: Representation of fusion events between chromosomes 29 and 8 of L. (L.) major (LmjF) and L. (L.) mexicana (LmxM), respectively. Figure S2. Representation of allelic transpositions between chromosomes 30 and 29 of L. (L.) major (LmjF) and L. (L.) mexicana (LmxM), respectively. Figure S3. Representation of allelic transpositions between chromosomes 31 and 30 of L. (L.) major (LmjF) and L. (L.) mexicana (LmxM), respectively. Figure S4. Representation of allelic transpositions between chromosomes 32 and 31 of L. (L.) major (LmjF) and L. (L.) mexicana (LmxM), respectively. Figure S5. Representation of allelic transpositions between chromosomes 33 and 32 of L. (L.) major (LmjF) and L. (L.) mexicana (LmxM), respectively. Figure S6. Representation of allelic transpositions between chromosomes 34 and 33 of L. (L.) major (LmjF) and L. (L.) mexicana (LmxM), respectively. Figure S7. Representation of allelic transpositions between chromosomes 35 and 34 of L. (L.) major (LmjF) and L. (L.) mexicana (LmxM), respectively. Figure S8. Representation of fusion events between chromosomes 36 and 20 of L. (L.) major (LmjF) and L. (L.) mexicana (LmxM), respectively. (PDF 732 KB)
About this article
- Leishmania (Viannia) braziliensis
- Leishmania (Leishmania) infantum
- Leishmania (Leishmania) major
- Leishmania (Leishmania) mexicana