Improved orthologous databases to ease protozoan targets inference
Parasites & Vectors volume 8, Article number: 494 (2015)
Homology inference helps on identifying similarities, as well as differences among organisms, which provides a better insight on how closely related one might be to another. In addition, comparative genomics pipelines are widely adopted tools designed using different bioinformatics applications and algorithms. In this article, we propose a methodology to build improved orthologous databases with the potential to aid on protozoan target identification, one of the many tasks which benefit from comparative genomics tools.
Our analyses are based on OrthoSearch, a comparative genomics pipeline originally designed to infer orthologs through protein-profile comparison, supported by an HMM, reciprocal best hits based approach. Our methodology allows OrthoSearch to confront two orthologous databases and to generate an improved new one. Such can be later used to infer potential protozoan targets through a similarity analysis against the human genome.
The protein sequences of Cryptosporidium hominis, Entamoeba histolytica and Leishmania infantum genomes were comparatively analyzed against three orthologous databases: (i) EggNOG KOG, (ii) ProtozoaDB and (iii) Kegg Orthology (KO). That allowed us to create two new orthologous databases, “KO + EggNOG KOG” and “KO + EggNOG KOG + ProtozoaDB”, with 16,938 and 27,701 orthologous groups, respectively.
Such new orthologous databases were used for a regular OrthoSearch run. By confronting “KO + EggNOG KOG” and “KO + EggNOG KOG + ProtozoaDB” databases and protozoan species we were able to detect the following total of orthologous groups and coverage (relation between the inferred orthologous groups and the species total number of proteins): Cryptosporidium hominis: 1,821 (11 %) and 3,254 (12 %); Entamoeba histolytica: 2,245 (13 %) and 5,305 (19 %); Leishmania infantum: 2,702 (16 %) and 4,760 (17 %).
Using our HMM-based methodology and the largest created orthologous database, it was possible to infer 13 orthologous groups which represent potential protozoan targets; these were found because of our distant homology approach.
We also provide the number of species-specific, pair-to-pair and core groups from such analyses, depicted in Venn diagrams.
The orthologous databases generated by our HMM-based methodology provide a broader dataset, with larger amounts of orthologous groups when compared to the original databases used as input. Those may be used for several homology inference analyses, annotation tasks and protozoan targets identification.
Historically, the very definition of a Protozoa represents an open debate. Despite many classifications and changes provided over history [1–4], in this article we will refer to Protozoa as eukaryotic organisms, apart from those who do not have a primitive mitochondria, peroxisomes (Archezoa) and the shared characteristics which define the Animalia, Fungi, Plantae and Chromista kingdoms .
There are over 200,000 described protozoan organisms and among them over 10,000 parasites of invertebrate organisms and nearly all vertebrate ones . There are several protozoan related diseases, which affect more than 25 % of the world population, such as Chagas’ disease, Human African Trypanosomosis, Leishmaniosis, Amoebiosis, Giardiosis, Toxoplasmosis, Cryptosporidiosis, Theileriosis, Babesiosis among many others [5–10].
Neglected Tropical Diseases (NTDs) are diseases caused by a variety of organisms and are usually associated to developing countries, which suffer from poor sanitation, hygiene, social and financial conditions. Over 1 billion people are affected by such diseases, in 149 countries worldwide . Among the 17 NTDs listed by WHO, three are caused by protozoan organisms: Chagas’ disease (Trypanosoma cruzi), Human African Trypanosomosis (Trypanosoma brucei) and Leishmaniosis (Leishmania spp.) .
According to the 3rd WHO report on NTDs, even though several advances have been achieved in the recent years, there is a permanent need for research and innovation in improved diagnosis, next-generation treatments and interventions for such NTDs [12, 13].
Leishmaniosis is a neglected disease caused by the Leishmania spp. and transmitted by phlebotomine sandflies . More than 1.3 million people are infected worldwide, especially those who live in poor sanitation, hygiene and social conditions and WHO estimates about 20,000 to 30,000 deaths occur yearly . Such disease has three distinct presentations: cutaneous, visceral and mucocutaneous, each of them related to different Leishmania spp. and world regions.
Leishmaniosis is hard to diagnose and treat. So far, the available drugs and vaccines are either toxic or present poor efficiency . Also, its elevated treatment cost (up to US$252/patient, depending on the applied treatment and drugs) eventually becomes prohibitive to the most affected and poor countries .
Besides allowing for a better comprehension of such Protozoa organism, many molecular studies have been done over the last few years using DNA/RNA sequencing methodologies [16–18], which have been used in order to infer new drug targets. These data are available in several public databases and allow for comparative genomics studies among either closely or distant related organisms. Also, that might increase the odds of discovering relevant information applied to drug manufacturing or reuse, which could be later applied to disease treatments.
Comparative genomics mainly refers to homology and evolutionary dynamics between organisms, genes and proteins, which provides better understanding on how species evolved through comparing either their complete genomes or specific genes . Homologous genes share a common ancestry, either intra- or inter-species. Several scenarios relate to homology, such as orthology, paralogy, horizontal gene transfer, gene loss, orphan genes and others; for this study, we will focus on orthology aspects only .
Orthology might be inferred when the same genes or proteins are present in distinct species, and this was due to a speciation event .
Homology inference has become an important issue when inferring function to recently sequenced genes because orthologs tend to preserve their ancestor function. Besides that, such studies provide a better insight on genes evolutionary history and consequently, to the species evolution [20, 21].
Inferring putative function is one of the particular benefits in orthologous group (OG) assignment, especially when dealing with recently sequenced genome data . In addition, OGs may provide us a better comprehension on species evolutionary relationships , since it is through such data that one might provide information that could help on both evolutionary and functional analysis .
Moreover, several tasks could benefit from OGs, such as genome annotation, gene conservation, protein family identification, phylogenetic tree reconstruction, pharmacology and many others [22, 24–27]. Topics as positional orthology and synteny conservation among orthologs are also appealing to those who aggregate genomic context in their homology inference methods .
There are several available methodologies to aid on homology detection. Besides a simple categorization effort , we will follow Dalquen’s proposition . Briefly, three distinct approaches are available: (i) the one which use multiple sequence alignment (MSA) scores along with reciprocal best hits, such as OrthoSearch , OrthoMCL  and InParanoid ; (ii) that which rely on evolutionary distance calculus, as RSD [32, 33]; (iii) and that based on phylogenetic trees reconstruction, as SPIMAP .
OrthoSearch  is a scientific workflow  for homology inference among species. Initially conceived as a Perl-based routine, it uses a reciprocal best hits, HMM-based approach. OrthoSearch has already proven to be effective inferring orthology among five protozoan genomes, using COG and KOG ODs .
In this work, we propose an update and a new functionality for OrthoSearch, showing it as an effective tool in providing means to create new ODs (n-ODs). So far, we tested our methodology in a controlled, three steps scenario: (i) Protozoa orthology inference and (ii) n-ODs creation, both supported by publicly available ODs used as input; and (iii) improved Protozoa orthology inference, supported by such recently created n-ODs.
With our methodology and generated n-ODs, we expect to be able to provide ODs with broader data sets, which in turn can be applied in target identification for protozoan organisms, such as stated by Timmers et al.  review on research efforts related to genomic database development for protozoan parasites.
Moreover, previous initiatives, such as the study performed by Tschoeke et al.  regarding the Leishmania amazonensis parasite, as well as the Leishmania donovani comparative genomics analysis performed by Satheesh et al.  corroborate the benefits provided by the use of broader orthologous data sets.
OrthoSearch improvements and analyses scenarios
In order to reach our main methodological goal, which is to provide OrthoSearch with means to create n-ODs, we revisited its original pipeline. Notably (i) we adopted HMMER version 3 and (ii) changed from a Perl-based routine to C++ 4.43 and Ruby 1.8.7 modules. A dedicated Ubuntu 12.04 single-server machine with 64 cores and 32GB RAM was used for all assembled scenarios.
OrthoSearch for protozoa orthology inference
OrthoSearch needs as input data an (i) OD and (ii) an organism multifasta protein data. We used Kegg Orthology (KO) [44, 45], EggNOG KOG and ProtozoaDB as input ODs. KO, downloaded via FTP, contains data from all life domains – Archaea, Bacteria and Eukarya. EggNOG KOG is a eukaryotic-only groups Eggnog subset , downloaded directly from its website. ProtozoaDB OGs  (which contain only protozoan species) were also used for our analyses. Details about each OD are available in Additional file 1.
We randomly selected three protozoan species as OrthoSearch organisms input data: Cryptosporidium hominis, Entamoeba histolytica and Leishmania infantum, each with 3,885, 7,973 and 7,872 proteins and downloaded such data from ProtozoaDB .
OrthoSearch for orthologous database building
Figure 1 depicts OrthoSearch pipeline with its two possibilities, which could be (i) a standard orthology inference or (ii) a n-OD creation. In order to build the n-ODs, we used as input data an OD in its original composition; and another OD data subset, which enacts as an organism multifasta protein data.
That was called an impersonated proteome and was generated by choosing a representative protein for each OG at the confronted OD. The selection and extraction of such proteins was performed with the support of a Python script kindly developed by Salvador Capella (personal communication, URL: https://github.com/scapella/trimal/blob/dev/scripts/get_sequence_representative_from_alignment.py) and an internally developed Ruby script. Each representative protein identifier and its amino acid sequence were stored in a single multifasta file (see Additional file 2).
We started with KO as our fixed, complete OD, against an impersonated EggNOG KOG multifasta protein data. Reciprocal best hits between both of them were processed using internal scripts developed in both Ruby and Unix/POSIX shell script languages, so that, new OGs were created and arranged, generating the n-OD conveniently named “KO + EggNOG KOG”.
During such n-OD creation, there were three possible scenarios related to the fixed database OGs – KO – and the impersonated proteome. Briefly, (i) those OGs from KO database which did not have any reciprocal best hit with the impersonated OGs from EggNOG KOG; and (ii) those OGs from EggNOG KOG that did not have any reciprocal best hit against KO database were identified, selected and incorporated to our n-OD (“KO + EggNOG KOG”), without any changes; (iii) those OGs from KO database which presented a reciprocal best hit with the impersonated OGs from EggNOG KOG were expanded, by adding up every respective EggNOG KOG OG proteins into such KO group.
Therefore, at the end of this run, we had a n-OD called “KO + EggNOG KOG” which comprised original, unaltered KO and EggNOG KOG OGs, as well as expanded KO OGs, which now also contain EggNOG KOG protein data. Once such n-OD was created, we performed the same steps described above, confronting “KO + EggNOG KOG” against ProtozoaDB impersonated OGs, which generated our second n-OD called “KO + EggNOG KOG + ProtozoaDB”.
Inferring protozoan orthologs with OrthoSearch and the n-ODs
Having built the two n-ODs: “KO + EggNOG KOG” and “KO + EggNOG KOG + ProtozoaDB”, we analyzed OrthoSearch in a standard orthology inference against the three above mentioned protozoan species.
Comparison between the n-ODs and OrthoMCLDB
The three protozoan species were confronted against OrthoMCLDB through online phyletic pattern search queries, in order to infer its orthologous proteins, in the same way as we did with the n-ODs created by our proposed methodology.
Quantitative results obtained while executing both OrthoSearch (with the two n-ODs) and OrthoMCLDB against the three protozoan species were compared in order to offer a better understanding of the proposed methodology behavior and to analyze if we were able to provide better results or not.
Potential Leishmania spp. targets against the human proteome
In order to identify potential protozoan targets that are not available at the human genome, a BlastP  was performed between the largest n-OD generated by our methodology - “KO + Eggnog KOG + ProtozoaDB” (details on the orthologous groups proteins are available at Additional file 3) and the human proteome, downloaded via RefSeq . We used BlastP 2.2.28+ with 0.1 as e-value, extracted and analyzed the orthologous groups which did not perform any hit against the human proteome but provided results against Leishmania spp. and therefore could represent potential targets. A BlastP was also performed against KO, Eggnog KOG and ProtozoaDB orthologous databases separately.
OrthoSearch for protozoa orthology inference
The protein data of the three protozoan species were confronted against (i) KO, (ii) EggNOG KOG and (iii) ProtozoaDB ODs (Fig. 2). ProtozoaDB performed best, with: 2,830 OGs against Cryptosporidium hominis, 4,952 for Entamoeba histolytica and 3,821 for Leishmania infantum.
With such data, we extracted coverage percentage information, which shows the total number of OGs inferred by OrthoSearch versus how many OGs are contained within each OD. For Cryptosporidium hominis, which has the smallest number of proteins of the three protozoan species studied, EggNOG KOG performed best, with 33 % coverage. Entamoeba histolytica also performed well with EggNOG KOG (37 %), but showed very similar results with ProtozoaDB (36 %), while showing a poor coverage with KO (13 %). Finally, Leishmania infantum had the best coverage (47 %), with EggNOG KOG.
Internal scripts, developed with the R language and its Venn Diagram library, processed reciprocal best hits for such protozoan species. We identified species-specific, pair-to-pair and core OGs, depicted at Fig. 3.
EggNOG KOG presented the best core (ratio between OGs shared by the three protozoan species and all the inferred OGs), corresponding to 26.72 % of the total inferred best hits (828/3,099 OGs), followed by KO - 20.94 % (719/3,434 OGs) and ProtozoaDB - 4.65 % (464/9,979 OGs).
In addition, ProtozoaDB presented the best species-specific results, with Entamoeba histolytica performing 40.77 % of the total OGs (4,086/9,979); Leishmania infantum with 27.93 % (2,787/9,979); and Cryptosporidium hominis with 19.68 % (1,964/9,979).
OrthoSearch for Orthologous Database building
Table 1 shows details on how many OGs remained intact and directly migrated to the n-ODs created by our methodology as well as those that were expanded.
After “KO + EggNOG KOG” building, we had a 14.02 % increase in the total number of OGs when compared to KO (16,938/14,856 OGs). In addition, 18.63 % KO OGs were expanded (2,769/14,856 OGs). When this n-OD was confronted against the impersonated ProtozoaDB OGs, we had a 63.54 % increase in the total number of OGs (27,701/16,938) and 17.88 % “KO + Eggnog KOG” OGs were expanded. Table 2 summarizes how many OGs, proteins, the average proteins per group and each n-OD size. “KO + EggNOG KOG + ProtozoaDB”, when compared to KO, provides 86.46 % more OGs (27,701/14,856 OGs) and 22.45 % more proteins (2,582,631/2,109,027 proteins). At last, Table 3 shows protozoan species representation at each created n-OD.
Inferring protozoan orthologs with OrthoSearch and the n-ODs
With these recently created n-ODs, on our second scenario we executed OrthoSearch using as input such n-ODs and the same three protozoan species then compared the obtained results against previous KO analysis. Figure 4 depicts coverage percentage data for each of the OG databases created by the methodology itself, for each organism adopted.
Our methodology provided an 86.47 % increase on the total number of OGs and a 22.45 % on the total amount of proteins when comparing “KO + EggNOG KOG + ProtozoaDB” against KO. Although there was a relevant increase in the number of inferred OGs for Cryptosporidium hominis (from 1,499 up to 3,254 groups), coverage increase was very subtle (10 %-12 %). Entamoeba histolytica, on the other hand, shows a relevant increase in coverage (19 %), especially when confronted against “KO + EggNOG KOG + ProtozoaDB” n-OD. Leishmania infantum had a very similar behavior to Cryptosporidium hominis, with a total of up to 4,760 OGs and from 15 % up to 17 % coverage.
We also consolidated the reciprocal best hits obtained in Venn diagrams, so that we might have a glimpse on species-specific, pair-to-pair and core OGs, as shown in Fig. 5.
Additional file 4 shows how many OGs were inferred as best hits with OrthoSearch when protozoan species were confronted against each of the generated n-ODs. “KO + EggNOG KOG” provided the best OGs core coverage ratio – 18.76 % (784/4,178) and “KO + EggNOG KOG + ProtozoaDB”, 5.65 % (627/11,089). “KO + EggNOG KOG + ProtozoaDB” had the best results in species-specific OGs, with Entamoeba histolytica at a 36.80 % ratio (4,081/11,089); Leishmania infantum with 29.94 % (3,320/11,089); and Cryptosporidium hominis at 18.81 % (2,086/11,089).
Comparison between the n-ODs and OrthoMCLDB
After submitting the same protozoan species to OrthoMCLDB online phyletic pattern search, 3,516 (Cryptosporidium hominis), 6,107 (Entamoeba histolytica) and 7,538 (Leishmania infantum) OGs were inferred. OrthoMCLDB inferred a 760 OGs core, which represents 5.12 % of the total best hits (760/14,846).
Concerning species-specific OGs, OrthoMCLDB detected 2,340 – 15.76 % (2,340/14,846) OGs for Cryptosporidium hominis; 4,816 – 32.44 % (4,816/14,846) for Entamoeba histolytica; and 6,145 – 41.39 % for Leishmania infantum; at last, pairwise shared OGs corresponded to 162 (Cryptosporidium hominis and Entamoeba histolytica), 254 (Cryptosporidium hominis and Leishmania infantum) and 369 (Entamoeba histolytica and Leishmania infantum) OGs respectively. Figure 6 shows a Venn diagram with obtained results.
Potential Leishmania spp. targets against the human genome
A BlastP against our largest created n-OD, “KO + Eggnog KOG + ProtozoaDB” (27,701 orthologous groups) allowed us to infer 7,622 (27.5 %) orthologous groups which did not perform any hit against the human proteome. Among such, 6.5 % (1,805/27,701) groups belong to KO or Eggnog KOG, but are not available in ProtozoaDB, which contains only protozoan organisms (Leishmania spp. included). Furthermore, 13 orthologous groups (0.05 %) contain at least one Leishmania spp. (Table 4), that should be considered as potential targets for further analysis.
The same BlastP query against each of the original ODs provided us the results listed in Table 5. These groups have no similarity with the human proteome and have at least one Leishmania spp. sequence.
In this analysis, we adopted new programming languages and updated the OrthoSearch pipeline with several bioinformatics tools, rewriting it to be later used in homology inference analyses and n-ODs creation.
OrthoSearch uses an algorithm based on reciprocal best hits calculation via HMM profiles, with Mafft being responsible for providing MSA’s and HMMER3 tools for generating HMM models, profiles and statistics. It benefits from multithreading provided by both Mafft and HMMER3 and accepts any of both tools available parameters. In a regular run, OrthoSearch uses an e-value cut-off flexible enough in order to aid on later profiles calculation.
Each of the studied ODs has its particular characteristics (e.g., while KO contains OGs from all life domains, EggNOG KOG contains only eukaryotes OGs and ProtozoaDB only protozoan OGs). That might influence the obtained results when running OrthoSearch with such ODs and organisms.
OrthoSearch execution with KO OD provides a significantly small core compared to KO size and the total number of best hits. That could be explained as KO contains proteins from many evolutionarily distant organisms, what could pose a challenge in the identification of closely related OGs.
Later, EggNOG KOG OD provided a discrete increase in the obtained protein core, most likely due to EggNOG KOG having only eukaryotic organisms’ data. While ProtozoaDB OD provided the smallest core among the three ODs, the total number of species-specific protein is extremely higher. This could be due to the reduced number of species in ProtozoaDB OD, along with the fact that all of those are protozoan organisms. Basically, the odds of obtaining a hit with a protein belonging to the own species being analyzed within OrthoSearch could increase.
We opted to choose a representative protein for each OG at the confronted OD and impersonate an organism multifasta protein data because that could minimize the required computational power and time needed to run OrthoSearch analyses.
Since an OD contains several OGs, which also contains several proteins, that would easily escalate the required time to confront such ODs. In addition, as each OG contains two or more proteins usually from closely related organisms, that could imply the possibility of two (or more) distinct proteins from the same OG obtaining a hit with distinct OGs at the confronted OD.
Our scenarios for n-OD creation were based on KO, EggNOG KOG and ProtozoaDB ODs. According to the literature, each of these ODs were created through particular methodologies: the use of metabolic pathways (KO), heuristic approaches and Gene Ontology  support (EggNOG KOG) and OrthoMCL algorithm (ProtozoaDB).
Our methodology allowed us to create n-ODs that either contain intact OGs which originated from the source or the confronted ODs or expanded OGs from the obtained reciprocal best hits inferred by OrthoSearch. The intact OGs contribution relates to offering more OGs for further analyses, while expanded ones provide more variability than those OGs from the original databases.
Besides providing a means to improve ProtozoaDB orthology inference, we opted to begin our n-OD creation tasks with KO and Eggnog due to both database variability – proteins from organisms from all life domains – and size. We also decided to maintain the original orthologous groups database identifiers, as well as their functional annotation. That might ease further steps related to information provenance.
Our proposed methodology works as a non-intrusive approach to a HMM-based pipeline – OrthoSearch – without changing its core functions. It uses ODs as input data and is capable to create n-ODs without requiring extensive computational power.
When looking at how protozoan species data fit to our proposed n-OD creation methodology, we observe a 61.98 % increase from KO to “KO + EggNOG KOG” OD, scaling from 3,612 up to 5,851 OGs with at least one Protozoa protein. In addition, the total number of protozoan proteins had a 162 % increase (from 46,027 up to 74,630) in such n-OD.
Furthermore, when KO and “KO + EggNOG KOG + ProtozoaDB” ODs were compared, we identified a 379 % increase in the number of OGs that contain at least one Protozoa protein (from 3,612 up to 17,305) and a 300 % increase in the total number of protozoan proteins (from 46,027 up to 138,814).
A broader dataset is usually desirable, as it may increase the odds of obtaining hits while inferring homology. As more organisms contribute to a n-OD, one might be able to obtain more hits with regular OrthoSearch runs confronting n-ODs and organisms multifasta protein data.
OrthoSearch uses a Markov chain based approach in order to create the n-OD OGs, which tend to comprise more evolutionary distant orthologous proteins than BLAST-based methodologies, such as OrthoMCLDB.
With the n-ODs created within our methodology, our initiative is another step to reinforce possibilities to build a gold, reference dataset  for orthology inference.
As more OGs (and respective proteins) from distinct species are added up to the n-ODs created, the methodology offers a broader dataset, with more data variability. Such data may be used for further homology inference analyses, which is a very desirable aspect in several comparative genomics applications. For example, phylogenomic studies which try to address gene conflicts and allow for optimal tree construction  or even review species’ definition .
The obtained results point towards the success of our proposed methodology, which encourage us in refining and creating more n-ODs. Such n-ODs might also be used in order to improve future functional annotation, re-annotation and also potential targets identification.
OrthoMCLDB was able to provide a broader species-specific OGs dataset, varying from 12 % (Cryptosporidium hominis - 2,340/2,086) up to 85 % (Leishmania infantum – 6,145/3,320) more OGs than our methodology. That may be due to several OrthoMCLDB groups containing only in-paralogs (20,853/124,740) from a single species lineage .
When looking at pairwise groups, there is a change in this scenario. Our methodology either provides the same quantitative results as OrthoMCLDB (Cryptosporidium hominis and Entamoeba histolytica – 162 OGs) or better, with 17.88 % more OGs (Entamoeba histolytica and Leishmania infantum – 435/369) and up to 48.81 % more OGs (Cryptosporidium hominis and Leishmania infantum – 378/254).
OrthoMCLDB inferred a larger absolute number of OGs in the three protozoan species core (760) than our methodology (627), which may be related to a broader seed of OGs (124,740/27,701). On the other hand, OrthoMCLDB performed poorly in coverage aspect, with only 0.06 % OGs in its core (760/124,740), while our methodology covered 2.26 % OGs (627/27,701), with our largest n-OD (“KO + EggNOG KOG + ProtozoaDB”). This may be due to the fact that the studied species are not so closely related. In addition, OrthoMCLDB uses a Blast-based algorithm, in a less sensitive approach than our methodology (OrthoSearch uses a protein-profile comparison) [9, 32].
Our methodology provides means for improved orthologous database creation using a HMM-based approach. Those new databases may contain a greater set of evolutionary distant homologous proteins, which could further extend the odds of inferring knowledge regarding the target organisms.
Specifically, our analyses allowed for a better comprehension on three protozoan species, as well as a deeper analysis on potential targets. For example, the obtained protozoan core orthologous proteins may allow us to evaluate which of these are housekeeping proteins and how they relate to the organism fitness.
Also, the species specific proteins – those which do not belong to the core, or those shared between two of the three studied protozoan organisms might be explored either as species-specific or group-specific targets, respectively.
The obtained BlastP results allowed us to infer orthologous groups which contain protozoan proteins – specifically Leishmania spp. - that could be used as potential targets for further analysis, as they posed no hit against the human proteome.
Among the Leishmania spp. inferred orthologous groups without hits against the human proteome (Table 4) are proteins already described in the literature as possible drug targets, briefly: trypanothione  (K01833.cdhit) – which relates to defense against oxidative stress ; and alpha-1,3-mannosyltransferase [55, 56] (K13690.cdhit) – enzyme essential to add mannose on the glycosylphosphatidyl, relates to the growing resistance to miltefosine. However, there are also other proteins, not yet described as drug targets, which should be further studied, briefly: the energy-converting hydrogenase B subunit I  (K14118.cdhit), found in the Archaea organism Methanothermobacter thermautotrophicus, which belongs to a domain related to MnhB subunit of Na+/H+ antiporter and is predicted as an integral membrane protein [57, 58]; and galactofuranosyltransferase  (K13672.cdhit), related to the LPG1 gene, which acts as a major ligand for macrophage adhesion [60, 61].
Our methodology also provided means to allocate new and evolutionary distant proteins to the original orthologous groups’ databases, identifying orthology relationships which have not been previously described.
Even though this is a preliminary analysis, it allowed us to evaluate the applied methodology and to forecast how its results may be used for protozoan target identification, either in a species-specific or shared point-of-view. This methodology will be later applied to all of 22 ProtozoaDB  protozoan organisms.
Our analyses used initially KO and EggNOG KOG databases as a starting point, adding later ProtozoaDB. We were able to create two n-ODS, “KO + EggNOG KOG” and “KO + EggNOG KOG + ProtozoaDB”, with each providing a larger amount of OGs when compared to the original. Those represent a broader dataset that can be used in future homology inference, annotation transfer and protozoan targets identification.
So far, our methodology allowed the identification of 13 potential targets in protozoan that would not have been identified without the distant homology approach provided by OrthoSearch and the n-ODs created by our methodology.
new orthologous database
Imam T. The complexities in the classification of protozoa: a challenge to parasitologists. Bayero J Pure Appl Sci. 2009;2:159–64.
Cavalier-Smith T. Predation and eukaryote cell origins: A coevolutionary perspective. Int J Biochem Cell Biol. 2009;41:307–22 [Molecular and Cellular Evolution: A Celebration of the 200th Anniversary of the Birth of Charles Darwin].
Cavalier-Smith T. Kingdoms Protozoa and Chromista and the eozoan root of the eukaryotic tree. Biol Lett. 2009;6(3):rsbl20090948.
Cavalier-Smith T. Kingdom Protozoa and Its 18 Phyla. Microbiol Rev. 1993;57:953–94.
Widmer G, London E, Zhang L, Ge G, Tzipori S, Carlton J, et al. Preliminary Analysis of the Cryptosporidium muris Genome. In: Ortega-Pierres G, Cacciò Simone M, Fayer R, Mank TG, Smith HV, Thompson RCA, editors. GIARDIA AND CRYPTOSPORIDIUM From Molecules to Disease. Wallingford, UK: CAB International; 2009. p. 320–7.
Thompson RCA. The Impact of Giardia on Science and Society. In: Ortega-Pierres, Guadalupe; Cacciò, Simone M; Fayer, Ronald; Mank, Theo G; Smith, Huw V; Thompson R, editors. GIARDIA AND CRYPTOSPORIDIUM From Molecules to Disease. Wallingford, UK: CAB International; 2009. p. 1–11.
Pain A, Renauld H, Berriman M, Murphy L, Yeats C, Weir W, et al. Genome of the host-cell transforming parasite Theileria annulata compared with T. parva. Science. 2005;309:131–3.
Elmore S, Jones JL, Conrad P, Patton S, Lindsay DS, Dubey JP. Toxoplasma gondii: epidemiology, feline clinical aspects, and prevention. Trends Parasitol. 2010;26:190–6.
Carlton JM, Adams JH, Silva JC, Bidwell SL, Lorenzi H, Caler E, et al. Comparative genomics of the neglected human malaria parasite Plasmodium vivax. Nature. 2008;455:757–63.
Brayton K, Lau AOT, Herndon DR, Hannick L, Kappmeyer LS, Berens SJ, et al. Genome sequence of Babesia bovis and comparative analysis of apicomplexan hemoprotozoa. PLoS Pathog. 2007;3:1401–13.
WHO | Neglected Diseases [http://www.who.int/neglected_diseases/diseases/en/].
WHO/Department of control of neglected tropical diseases. Investing to Overcome the Global Impact of Neglected Tropical Diseases. Geneva, Switzerland: World Health Organization; 2015.
The London declaration on Neglected Tropical Diseases [http://unitingtocombatntds.org/sites/default/files/resource_file/london_declaration_on_ntds.pdf].
WHO | Leishmaniasis [http://www.who.int/mediacentre/factsheets/fs375/en/].
WHO Technical Report Series [http://whqlibdoc.who.int/trs/WHO_TRS_949_eng.pdf].
Singh N, Chikara S, Sundar S. SOLiDTM sequencing of genomes of clinical isolates of leishmania donovani from india confirm leptomonas co-infection and raise some key questions. PLoS One. 2013;8, e55738.
Brotherton M-C, Bourassa S, Leprohon P, Légaré D, Poirier GG, Droit A, et al. Proteomic and genomic analyses of antimony resistant leishmania infantum mutant. PLoS One. 2013;8, e81899.
Tschoeke DA, Nunes GL, Jardim R, Lima J, Dumaresq AS, Gomes MR, et al. The comparative genomics and phylogenomics of leishmania amazonensis parasite. Evol Bioinforma Online. 2014;10:131–53.
Hardison RC. Comparative genomics. PLoS Biol. 2003;1:156–60.
Koonin EV. Orthologs, paralogs, and evolutionary genomics. Annu Rev Genet. 2005;39:309–38.
Sequence - Evolution - Function - NCBI Bookshelf [http://www.ncbi.nlm.nih.gov/books/NBK20260/].
Dessimoz C, Gabaldón T, Roos DS, Sonnhammer ELL, Herrero J, Quest for Orthologs Consortium. Toward community standards in the quest for orthologs. Bioinforma Oxf Engl. 2012;28:900–4.
Delsuc F, Brinkmann H, Philippe H. Phylogenomics and the reconstruction of the tree of life. Nat Rev Genet. 2005;6:361–75.
Li L, Stoeckert CJ, Roos DS. OrthoMCL: Identification of ortholog groups for eukaryotic genomes. Genome Res. 2003;13:2178–89.
Dessimoz C. Editorial: Orthology and applications. Brief Bioinform. 2011;12:375–6.
Kristensen DM, Wolf YI, Mushegian AR, Koonin EV. Computational methods for Gene Orthology inference. Brief Bioinform. 2011;12:379–91.
Cuadrat RRC, Cruz SM da S, Tschoeke DA, Silva E, Tosta F, Jucá H, et al. An Orthology-Based Analysis of Pathogenic Protozoa Impacting Global Health: An Improved Comparative Genomics Approach with Prokaryotes and Model Eukaryote Orthologs. OMICS J Integr Biol. 2014;140624130015005.
Dewey CN. Positional orthology: putting genomic evolutionary relationships into context. Brief Bioinform. 2011;12:401–12.
Dalquen DA, Altenhoff AM, Gonnet GH, Dessimoz C. The Impact of Gene Duplication, Insertion, Deletion, Lateral Gene Transfer and Sequencing Error on Orthology Inference: A Simulation Study. PLoS One. 2013;8, e56925.
Da Cruz SMS, Batista V, Silva E, Tosta F, Vilela C, Cuadrat R, et al. Detecting distant homologies on protozoans metabolic pathways using scientific workflows. Int J Data Min Bioinforma. 2010;4:256–80.
Remm M, Storm CEV, Sonnhammer ELL. Automatic clustering of orthologs and in-paralogs from pairwise species comparisons. J Mol Biol. 2001;314:1041–52.
DeLuca TF, Cui J, Jung J-Y, Gabriel KCS, Wall DP. Roundup 2.0: Enabling comparative genomics for over 1800 genomes. Bioinformatics. 2012;28:715–6.
Wall DP, Fraser HB, Hirsh AE. Detecting putative orthologs. Bioinformatics. 2003;19:1710–1.
Rasmussen MD, Kellis M. A Bayesian Approach for Fast and Accurate Gene Tree Reconstruction. Mol Biol Evol. 2011;28:273–90.
Chen F, Mackey AJ, Stoeckert CJ, Roos DS. OrthoMCL-DB: querying a comprehensive multi-species collection of ortholog groups. Nucleic Acids Res. 2006;34 suppl 1:D363–8.
O’Brien KP, Remm M, Sonnhammer ELL. Inparanoid: a comprehensive database of eukaryotic orthologs. Nucleic Acids Res. 2005;33(Database issue):D476–80.
Tatusov RL, Fedorova ND, Jackson JD, Jacobs AR, Kiryutin B, Koonin EV, et al. The COG database: an updated version includes eukaryotes. BMC Bioinformatics. 2003;4:41.
Powell S, Forslund K, Szklarczyk D, Trachana K, Roth A, Huerta-Cepas J, et al. eggNOG v4.0: nested orthology inference across 3686 organisms. Nucleic Acids Res 2013;42(D1):gkt1253.
Powell S, Szklarczyk D, Trachana K, Roth A, Kuhn M, Muller J, et al. eggNOG v3.0: orthologous groups covering 1133 organisms at 41 different taxonomic ranges. Nucleic Acids Res. 2011;40:D284–9.
Da Cruz SMS, Batista V, Dávila AMR, Silva E, Tosta F, Vilela C, et al. OrthoSearch: a scientific workflow approach to detect distant homologies on protozoans. New York, NY, USA: ACM; 2008. p. 1282–6 [SAC ’08].
Taylor IJ, Deelman E, Gannon DB, Shields M. Workflows for E-Science. 2007.
Timmers LFSM, Pauli I, Barcellos GB, Rocha KB, Caceres RA, de Azevedo WF, et al. Genomic databases and the search of protein targets for protozoan parasites. Curr Drug Targets. 2009;10:240–5.
S SK, R.K. G, Ghosh M. Comparative in-silico genome analysis of Leishmania (Leishmania) donovani: A step towards its species specificity. Meta Gene. 2014;2:782–98.
KEGG Orthology [http://www.genome.jp/kegg/ko.html].
KEGG. Kyoto Encyclopedia of Genes and Genomes [http://nar.oxfordjournals.org/content/28/1/27.abstract?ijkey=8863b3d5385d8722cdd87f19eccd8a62b567a0b4&keytype2=tf_ipsecsha].
Dávila AMR, Mendes PN, Wagner G, Tschoeke DA, Cuadrat RRC, Liberman F, et al. ProtozoaDB: dynamic visualization and exploration of protozoan genomes. Nucleic Acids Res. 2008;36(Database issue):D547–52.
Camacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K, et al. BLAST+: architecture and applications. BMC Bioinformatics. 2009;10:421.
RefSeq: NCBI Reference Sequence Database [http://www.ncbi.nlm.nih.gov/refseq/].
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, et al. Gene Ontology: tool for the unification of biology. Nat Genet. 2000;25:25–9.
Chen F, Mackey AJ, Vermunt JK, Roos DS. Assessing Performance of Orthology Detection Strategies Applied to Eukaryotic Genomes. PLoS One. 2007;2.
McCormack JE, Faircloth BC, Crawford NG, Gowaty PA, Brumfield RT, Glenn TC. Ultraconserved elements are novel phylogenomic markers that resolve placental mammal phylogeny when combined with species-tree analysis. Genome Res. 2012;22:746–54.
Konstantinidis KT, Tiedje JM. Genomic insights that advance the species definition for prokaryotes. Proc Natl Acad Sci. 2005;102:2567–72.
Colotti G, Baiocco P, Fiorillo A, Boffi A, Poser E, Chiaro FD, et al. Structural insights into the enzymes of the trypanothione pathway: targets for antileishmaniasis drugs. Future Med Chem. 2013;5:1861–75.
Krauth-Siegel RL, Meiering SK, Schmidt H. The Parasite-Specific Trypanothione Metabolism of Trypanosoma and Leishmania. Biol Chem. 2003;384.
Shinde S, Mol M, Jamdar V, Singh S. Molecular modeling and molecular dynamics simulations of GPI 14 in Leishmania major: insight into the catalytic site for active site directed drug design. J Theor Biol. 2014;351:37–46.
Garami A, Ilg T. The role of phosphomannose isomerase in Leishmania mexicana glycoconjugate synthesis and virulence. J Biol Chem. 2001;276:6566–75.
Ito M, Guffanti AA, Krulwich TA. Mrp-dependent Na+/H+ antiporters of Bacillus exhibit characteristics that are unanticipated for completely secondary active transporters. FEBS Lett. 2001;496:117–20.
Hiramatsu T, Kodama K, Kuroda T, Mizushima T, Tsuchiya T. A putative multisubunit Na+/H+ antiporter from Staphylococcus aureus. J Bacteriol. 1998;180:6642–8.
Huang C, Turco SJ. Defective galactofuranose addition in lipophosphoglycan biosynthesis in a mutant of Leishmania donovani. J Biol Chem. 1993;268:24060–6.
Zhang K, Barron T, Turco SJ, Beverley SM. The LPG1 gene family of Leishmania major. Mol Biochem Parasitol. 2004;136:11–23.
Spath GF, Epstein L, Leader B, Singer SM, Avila HA, Turco SJ, et al. Lipophosphoglycan is a virulence factor distinct from related glycoconjugates in the protozoan parasite Leishmania major. Proc Natl Acad Sci. 2000;97:9258–63.
Our sincere thanks to everyone involved, not only from our laboratory, but to everyone who inspired, stimulated and provided us with feedback. To Salvador Capella-Gutierrez (CRG) for the script to obtain representative sequences from multiple alignments. A special thank you to our fellow colleagues: Diogo Tschoeke, Rafael Cuadrat and Sérgio Serra.
The authors declare that there is no conflict of interest regarding the publication of this paper.
AMRD designed and coordinated the analyses. NPKF and RJ were responsible for programming tasks, experiment design and writing the manuscript. NPKF performed and collected experiment data. All authors revised the final version of the manuscript and AMRD approved it.
Orthologous databases main characteristics. This file provides a detailed view on the orthologous databases characteristics, regarding: how many organisms (total and protozoan only), orthologous groups, total proteins and the database total size. (XLS 6 kb)
Sample impersonated multifasta protein sample; This file provides a quick view on how the organisms’ representative proteins were stored in a single multifasta file, which impersonates an organism proteome and actually contains proteins from several orthologous groups from a single orthologous database. (FASTA 1 kb)
“KO + Eggnog KOG + ProtozoaDB” orthologous groups (protein identifiers only). For each ortholog group, we list all protein identifiers (corresponding organism and accession number). (BZ2 9971 kb)
Inferred RBH between protozoan organisms and orthologous databases. This file provides information on inferred RBH between protozoan species against KO database and the two n-ODs created by the methodology as well. It is detailed on species-specific, pairwise and core RBH (shared between all organisms) absolute RBH and their percentage when compared to the total amount of obtained RBH. (XLS 27 kb)
About this article
Cite this article
Kotowski, N., Jardim, R. & Dávila, A.M.R. Improved orthologous databases to ease protozoan targets inference. Parasites Vectors 8, 494 (2015). https://doi.org/10.1186/s13071-015-1090-0
- Comparative genomics
- Homology inference
- Target identification
- Orthologous database
- Distant homology