- Research article
The Schistosoma mansoni phylome: using evolutionary genomics to gain insight into a parasite’s biology
BMC Genomicsvolume 13, Article number: 617 (2012)
Schistosoma mansoni is one of the causative agents of schistosomiasis, a neglected tropical disease that affects about 237 million people worldwide. Despite recent efforts, we still lack a general understanding of the relevant host-parasite interactions, and the possible treatments are limited by the emergence of resistant strains and the absence of a vaccine. The S. mansoni genome was completely sequenced and still under continuous annotation. Nevertheless, more than 45% of the encoded proteins remain without experimental characterization or even functional prediction. To improve our knowledge regarding the biology of this parasite, we conducted a proteome-wide evolutionary analysis to provide a broad view of the S. mansoni’s proteome evolution and to improve its functional annotation.
Using a phylogenomic approach, we reconstructed the S. mansoni phylome, which comprises the evolutionary histories of all parasite proteins and their homologs across 12 other organisms. The analysis of a total of 7,964 phylogenies allowed a deeper understanding of genomic complexity and evolutionary adaptations to a parasitic lifestyle. In particular, the identification of lineage-specific gene duplications pointed to the diversification of several protein families that are relevant for host-parasite interaction, including proteases, tetraspanins, fucosyltransferases, venom allergen-like proteins, and tegumental-allergen-like proteins. In addition to the evolutionary knowledge, the phylome data enabled us to automatically re-annotate 3,451 proteins through a phylogenetic-based approach rather than solely sequence similarity searches. To allow further exploitation of this valuable data, all information has been made available at PhylomeDB (http://www.phylomedb.org).
In this study, we used an evolutionary approach to assess S. mansoni parasite biology, improve genome/proteome functional annotation, and provide insights into host-parasite interactions. Taking advantage of a proteome-wide perspective rather than focusing on individual proteins, we identified that this parasite has experienced specific gene duplication events, particularly affecting genes that are potentially related to the parasitic lifestyle. These innovations may be related to the mechanisms that protect S. mansoni against host immune responses being important adaptations for the parasite survival in a potentially hostile environment. Continuing this work, a comparative analysis involving genomic, transcriptomic, and proteomic data from other helminth parasites, other parasites, and vectors will supply more information regarding parasite’s biology as well as host-parasite interactions.
Schistosoma mansoni, S. haematobium, and S. japonicum (Platyhelminthes: Trematoda) are the main causative agents of human schistosomiasis, a neglected tropical disease that is endemic in 77 countries where more than 237 million people require preventive chemotherapy and other 779 million live in areas of risk of infection [1–4]. The genomes of these parasites have been recently published providing insights into parasite’s development, infection, and host-parasite interactions [5–7]. However, even with the progress made over the last years, schistosomiasis control depends primarily on the treatment of infected patients with Praziquantel®, the only drug available for mass treatment (e.g. [5, 8, 9]). Drawbacks of this drug are that it does not prevent against reinfection and its effectiveness varies depending on several factors such as the parasite’s gender, developmental stage, and the time of infection. Furthermore, Praziquantel®-resistant parasites have been found both in the laboratory and in the field, thus increasing the urgent need for new effective drugs and vaccines [10–13].
Schistosoma mansoni infects 7.1 million people in America, 95% of which in Brazil, and 54 million people in Sub-Saharan Africa causing intestinal and hepatosplenic schistosomiasis [14, 15]. The S. mansoni genome sequencing data was published in 2009 and a new version was recently released [5, 16]. The improved genome has 364.5 megabases (Mb) assembled in 885 scaffolds, half of which are represented in scaffolds greater than 2 kilobases . A total of 10,852 genes were identified, encoding over 11,000 proteins, 45% of which remain without known or predicted function [5, 16, 17]. 81% of the genome was assembled onto the parasite’s chromosomes, providing a partial genetic map [16, 18]. The availability of genomic data offers new opportunities for innovation in the control of schistosomiasis, by providing information that allows for the identification of novel drug targets and vaccine candidates through a system-wide perspective [5, 19, 20].
Making accurate functional predictions for genes or proteins is a key step in every genome sequencing project. However, on average, 30 to 50% of the predicted proteome remains uncharacterized while for the remaining set only general predictions are made. To deal with the gap between the rapid progress in genome sequencing and experimental characterization of genes and gene products, computational methods have been developed [21–23]. Two main approaches are generally used for functional prediction of genes and their products: one based on sequence similarity searches and another on phylogenetic analysis.
Owing to the computational cost and complexity of large scale phylogenetic analysis, the accurate identification of orthology relationships remains a challenge in comparative genomics and most of the orthology prediction methods rely on similarity-based search (e.g. BLAST , OrthoMCL , InParanoid ). In these cases, functional prediction is obtained based on the transfer of information from the most similar sequences in the database to the gene or protein of interest (e.g. ). However, several limitations are associated with this method, mainly the lack of a straightforward relationship between sequence similarity and protein function [21, 27–29]. Since this approach is fast, simple, and can be automated to analyze thousands of genes, it has been used frequently to predict functional products encoded by newly sequenced genomes. Over the last years this practice has generated systematic errors, the extent of which is not completely known [22, 27–32].
In an attempt to improve the accuracy of functional prediction at a large scale, phylogenetic methods may be applied [33, 34]. The advantage of such methods is that they focus on the evolutionary history of genes rather than merely on their sequence similarity [30, 35, 36]. Ideally, functional transfer in the genomic context or for specific genes/proteins should be performed only when there is any experimental evidence for those used as source of information. However, in databases as UniProt, only 3% of proteins have experimental support for their annotations . To deal with the absence of experimental support for most part of the available proteomes, transfer of functional annotation aiming to provide hints regarding the gene/protein function needs to follow strict requirements to avoid, as much as possible, misclassifications. In the last decade, the publication of a large number of genomic and proteomic data and the development of faster and powerful computers, new software, and automated pipelines have allowed for the reconstruction of phylogenetic trees of the complete set of proteins encoded in a genome – the so called phylome .
The phylome data may give a broad view of the evolution of an organism, since it comprises the phylogenies of all proteins encoded in its genome . Most notably, a phylome can be used to detect specific evolutionary scenarios, to quantify the fraction of individual phylogenies whose topologies are consistent with a given hypothesis, and to improve functional annotation of proteins and biological systems [38, 39]. Furthermore, comparing genomes or proteomes through an evolutionary perspective may provide insights to the understanding of the metabolism, physiology, pathogenicity, and the adaptation to a particular life style of organisms. In this context, the availability of S. mansoni genomic data provides the opportunity to study this parasite from a genome-wide perspective rather than from individual gene or protein analyses.
Taking advantage of the benefits provided by a genome-wide approach combined with an evolutionary perspective, we reconstructed the S. mansoni phylome with the goals of i) gaining insight into lineage-specific evolutionary events potentially related to the parasitic lifestyle, and ii) improving the functional annotation of the genome/proteome.
Phylogenetic techniques used in the present work included multiple sequence alignment [40–43] alignment trimming , neighbor-joining tree building , evolutionary model testing, and maximum likelihood analysis . The resulting phylome data contains 7,964 protein phylogenetic trees, covering the analysis of 11,763 S. mansoni proteins and their homologs in 12 other organisms, out of which we identified evolutionary events and homology relationships. The results provided useful information about the parasite’s genome evolution such as the identification of gene duplication events and expanded protein families such as proteases, tetraspanins, fucosyltransferases, venom allergen-like proteins (also called as SmVAL or SCP-like), tegumental-allergen-like proteins (SmTAL), among others. Altogether, the results obtained are likely to pave the way for a better understanding of the parasite’s biology including host-parasite interactions. This, in turn will accelerate the search for new drugs and vaccine directed toward the control and eradication of schistosomiasis.
Results and discussion
Reconstruction of the S. mansoni phylome
The S. mansoni phylome reconstructed in this work was derived from the comparative analysis of all proteins encoded in the parasite genome (predicted proteome) and their homologs in 12 other eukaryotic proteomes whose genomes were completely sequenced (Table 1). The set of selected species is particularly rich in metazoans (11 species), including ten invertebrates, one tunicate, and one vertebrate. One choanoflagellate, Monosiga brevicollis, was included as outgroup of the phylogenetic reconstruction. The metazoan species selected represent important evolutionary innovations, e.g. the origin of the third germ layer, the development of organs, systems, complex patterns of communication, and the emergence of the adaptive immune system, making this dataset set especially suitable for addressing the evolutionary innovations in S. mansoni in the context of metazoan evolution.
To perform the phylogenetic analyses, we applied an automated pipeline similar to the one used for the human phylome project . This pipeline is illustrated here (Figure 1). The resulting alignments, phylogenies, and orthology predictions can be accessed at PhylomeDB  (http://phylomedb.org).
Using this phylogenomic approach, we analyzed 11,763 S. mansoni proteins and obtained 7,964 phylogenetic trees covering 70% of the parasite’s proteome. This coverage is remarkably similar to that of other phylome data of newly sequenced genomes such as that of the pea aphid Acyrthosiphon pisum (67%) .
The absence of trees for the remaining 3,490 proteins is either due to a possible high degree of divergence between the S. mansoni proteins and their homologs in the other selected species, an indication of the uniqueness of the parasite’s proteome, or it reflects the presence of errors in gene models. Out of the 7,964 phylogenetic trees, 3.038 (38%) correspond to trees with “seed” proteins with a completely unknown function and without any GO  assignment in SchistoDB .
Phylogeny-based orthology prediction
In order to create a complete list of orthology and paralogy relationships among S. mansoni proteins and those encoded in the other eukaryotic proteomes included in this work, we analyzed the parasite’s phylome using a species-overlap algorithm as previously described . The comprehensive catalogue of phylogeny-based orthology and paralogy relationships among S. mansoni and other species was made publicly available at PhylomeDB .
Owing to the increasing rate at which new fully sequenced genomes are released, the accumulation of genomic and proteomic data has been much higher than the rates at which genes or proteins are experimentally characterized. Aiming at producing a high confidence set of functional predictions for S. mansoni proteins, we used the evolutionary relationships as inferred from phylogenetic trees to obtain subsets of one-to-one (single homolog in S. mansoni and in other species) homology relationships among S. mansoni proteins and the homologs from other species included in the present study (Figure 2).
By using such phylogeny-based approach, we transferred 10,175 functional annotations (GO terms ) to 3,451 S. mansoni proteins, from which 790 (7% of the parasite’s proteome) were previously annotated as “hypothetical protein”, corresponding to proteins whose function had not been predicted or experimentally tested before (Additional file 1 Table S1). The transfer was performed from each ortholog with known function in the selected taxa to the S. mansoni “seed” protein. For the other proteins that already had any functional prediction, the annotation was confirmed or improved. Consequently, a “seed” protein could receive more than one functional description. In these cases, all functional annotations were maintained allowing the user to choose the closest related transferred functional annotation, those that came from model organisms, or even to create a consensus based on all of them.
To validate the applied methodology, we retrieved reviewed S. mansoni proteins from UniProt , including experimentally confirmed ones, to evaluate the annotation transferred by the phylogenomic approach. The functional annotations performed by PhylomeDB correspond to known functions in the aforementioned database (Additional file 1 Table S2). Even though the BLAST search may detect distant homologs with additional domains, our subsequent phylogenetic reconstruction and our selection of orthologs will select those orthologs that are likely to have similar domain architecture. This is an additional reason why an orthology-based annotation is preferred over sequence similarity searches, since orthologs as compared to paralogs have a higher tendency to share a similar domain architecture .
Although less reliable than those based on one-to-one orthology relationships, annotation transfer based on more complex subsets (one-to-many, many-to-one, or many-to-many) may provide important hints to predict the biological function of S. mansoni proteins. However, in these cases, one or more genes are co-orthologous to a set of genes in another genome due to lineage-specific duplication(s) that can be associated with functional shifts, affecting the reliability of the functional transfer [38, 51]. An example of a one-to-one transfer from a Drosophila melanogaster protein to a S. mansoni protein comes from the phylogenetic reconstruction of the Phy000V14T_SCHMA (Smp_170950) protein, potentially related to the glycine cleavage system, and its homologs in the selected species (Figure 3). The analysis of this tree resulted in six transfers of functional annotation from homologous proteins to the S. mansoni “seed” protein. The GO terms in all six functional annotations are related to aminomethyltransferase activity and glycine catabolic process providing further support for the annotation transfer. In this example, to illustrate a case of a one-to-one transfer, we chose the functional annotation transferred from Drosophila melanogaster once, according to the information available in UniProt , it is one of the orthologs with known function and experimental validation. Tags for homologous sequences with experimental validation are not available in PhylomeDB . However, links to UniProt  and other databases are provided.
To explore the benefits offered by comparative genomics in order to improve functional annotation of genes and gene products, it is also necessary to consider the limitations involved in this approach. Although it is generally accepted that functional annotation through orthology, rather than just homology relationship, constitutes one of the most promising annotation approaches, these surveys are designed to provide predictions regarding the likely protein function, but it does not substitute experimental confirmation [36, 52]. Functional diversity is often associated with significant divergence at the sequence level, but high levels of identity do not ensure that two or more proteins perform the same function, since subtle changes in active sites are able to completely change the protein function .
As we previously mentioned, evolutionary analysis involving fully sequenced genomes/proteomes remains a challenge. Although the tools here applied were not originally designed for large scale phylogenetic analysis, we adapted them to work on a large scale, since we strongly believe that a system-wide perspective on evolutionary processes can greatly improve the understanding on how genomes came to be and what evolutionary process took them there. Functional prediction as described in the present work could be used as a starting point for future projects, prioritizing the selection of certain genes or proteins for new experimental studies.
Detection of gene duplications in S. mansoni
An additional advantage of the phylogeny-based approach is that it readily provides a collection of gene evolutionary histories that can be mined for particular events. Since gene duplication is considered one of the main mechanisms for functional innovation and diversification , we explored the S. mansoni phylome to identify protein families that have been specifically expanded in this lineage, since its diversification from the other sequenced metazoans. We used the above-mentioned species-overlap algorithm that identifies duplication nodes and also provides clues of the relative dating of the duplication event [39, 55].
Such analysis revealed that in 3,051 reconstructed phylogenetic trees there is at least one paralog connected to the “seed” protein through a duplication node (Additional file 1 Table S3). Among these, 211 phylogenies show lineage-specific duplications in the three Schistosoma species in comparison with the other taxa. These expansions are small-to-moderate in size, resulting in a total of two to ten paralogs, and include some of the most significant expansions as discussed below.
The inclusion of S. haematobium and S. japonicum proteomes gave us a high resolution within Schistosoma genus and allowed us to make comparisons across this taxon. In general, the expansions observed in S. mansoni can also be observed in the other two Schistosoma species, although with variable number of paralogs in each species. As previously observed by evolutionary relationships, cytogenetic data, and syntenic analyses, the present study shows that S. mansoni is more closely related to S. haematobium than to the S. japonicum[56–59]. Moreover, 170 evolutionary trees have only S. mansoni and S. haematobium proteins, while only six phylogenies have solely S. mansoni and S. japonicum proteins. Meanwhile, most of the homologous proteins shared by S. mansoni and S. haematobium are annotated as “hypothetical protein” and do not have any predicted function or significant hits with known proteins in public databases as UniProt , Pfam , or non-redundant (nr) NCBI database (ftp://ftp.ncbi.nih.gov/blast/db).
A small number of phylogenetic trees (1,45%) had only sequences of S. mansoni. These could be the result of very recent duplication events of proteins that are specific to this species. However, many of these genes were not found in the genetic map of S. mansoni[16, 18] and they do not contain protein domains traceable at Pfam . BLAST searches against the non-redundant (nr) NCBI database detected a few non-Schistosoma proteins as significant hits that were annotated as hypothetical in all cases. For these reasons we rather believe that these sequences correspond to spurious predictions. Further analyses will be conducted in the future in order to confirm or refute this hypothesis.
Among the most significant protein expansions in S. mansoni we identified tetraspanins, fucosyltransferases, venom allergen-like proteins (SmVAL), tegumental-allergen-like proteins (SmTAL), leishmanolysins, and elastases, which were previously proposed as drug targets, once they can be related to morphological or physiological specificities of this parasite [5, 20, 61–65]. In these cases, the protein family membership ranged from 6 to 23 paralogs encoded in the parasite’s genome.
Tetraspanins are small proteins with four transmembrane domains involved in the coordination of intra and intercellular processes, such as signal transduction, cell proliferation, adhesion, and migration, cell fusion and host-parasite interactions [66, 67]. The function of schistosome tetraspanins are not completely understood, but cell-cell interactions and maintenance of cell membrane integrity might be performed by these proteins as well as they can be receptors for host ligands, acting on immune evasion . The suppression of two tetraspanin genes (Sm-tsp-1 and Sm-tsp-2) by RNA interference in mice also suggests that these proteins play important structural roles in the parasite’s tegument, being a good target for anti-schistosomal vaccine . Figure 4 illustrates an example of tetraspanin lineage-specific duplications. In this case, the number of homologs in the three Schistosoma species varies from six to eight. Tree topology shows distinct well-supported clades suggesting that structural and/or functional variants might be present. Three proteins in this dataset have experimental evidence: Phy0048JNS_SCHHA (Q26499), Phy0048WJL_SCHMA (P19331), and Phy0005UU9_DROME (O46101) [49, 69, 70].
Venom allergen-like proteins (SmVAL), also called sperm-coating protein-like (SCP-like), are structurally related proteins members of the SCP/TAPS family. In Platyhelminthes, these proteins have been linked as potential modulators of immune function and components of sexual development . Although the specific function of each SmVAL family member is unknown, there is evidence suggesting potential roles in larval penetration, host immune response modulation, and adult worm development [63, 71]. Furthermore, analyses of SmVAL transcripts demonstrated that the corresponding genes are upregulated in infective stages of the parasite, highlighting SmVAL proteins as candidates for novel vaccine strategies [71, 72].
Fucosyltransferases are enzymes that catalyses the fucose transfer from the donor guanosine-diphosphate fucose to different acceptor molecules such as oligosaccharides, glycoproteins, and glycolipids . In schistosomes, fucosyltransferases are involved in producing immunomodulatory epitopes during infection, granuloma formation, egg/endothelium interactions, and were previously highlighted as anti-schistosomal candidates [63, 74].
Tegumental-Allergen-Like proteins (SmTALs) are members of a protein family present in parasitic Platyhelminthes [64, 75]. These proteins are located inside the tegument and have different life-cycle expression patterns . The tegumental protein Sm22.6 is considered the main target for human IgE in S. mansoni and human IgE response against this protein is associated with the development of age-dependent partial immunity to S. mansoni infections in endemic areas [64, 76].
Leishmanolysin, (also called invadolysin and SmPepM8), is a major surface protease member of the metallopeptidase M8 family. This protein can perform activities in schistosomes similar to those performed in Leishmania where these proteins are involved in different types of processes like degradation of the extracellular matrix and inhibition or perturbations of host cell interactions [63, 72]. In turn, elastases are serine proteases that in schistosomes play a pivotal role in the penetration by cercariae of host skin to initiate infection. Recent studies have also revealed that these proteases can be employed by schistosomes to overcome or evade the host immune response [77, 78]. Members of S. mansoni peptidase families such as leishmanolysins, cercarial elastases, and cathepsin D proteins were subjected to a detailed study in respect to their domain architectures, functional properties, and evolutionary relationships as described elsewhere .
Another specific feature of schistosomes is related to their tegument. Distinct from nematodes, which have a cuticle covering and protecting the organism body, schistosomes are covered by a living syncytium bounded by a complex multilaminate surface, which undergoes several adaptations soon after infection is initiated [79–81]. The external double membrane plays a crucial role in host-parasite interactions, being responsible for diverse mechanisms of survival [19, 82, 83]. The development of a tegument, highly specialized and resistant to immune damage, was accompanied by evolutionary adaptations, for example, the expansions of other protein families encoding annexins, cadherins, and innexins.
Annexins are widely distributed in eukaryotes performing a broad range of important biological processes related to tegument membrane [84–86]. In schistosomes, annexins appear to be involved in parasite’s stability protecting against immune attack by the host as well as against structural breakdown [85, 86]. Cadherins are adhesion molecules that mediate Ca2+-dependent cell-cell adhesion and whose duplication events happened probably in parallel to the advent of a third germ layer in flatworms [5, 87]. Innexins are components of gap-junction proteins, the intercellular channels that allow for the exchange of ions and other small signal molecules [88, 89]. In C. elegans, innexins have been implicated in different processes like electrical coupling between pharyngeal muscles, calcium propagation in the gut, gap junction-mediated oocyte, and sensory neuron identity .
In summary, we identified that approximately 45% of the S. mansoni predicted proteins that were covered by this phylogenomic analysis have, at least, one paralog encoded in the parasite genome that might have arisen by gene duplication events that occurred after its divergence from other selected taxa (Additional file 1 Table S3). In other eukaryotic genomes this value ranges from 30 and 65% , whereas in C. elegans this value is equal to 49% .
Altogether, the present results indicate that besides the exploitation of host endocrine and immune signals, the parasite genome exhibit multiple events of gene duplication which may be, at least partially, an adaptive response related to the parasitic lifestyle. These expansions probably reflect the intriguing complexity of evolutionary events that happened over time, resulting in important characteristics in schistosome’s biology with consequences to the disease it causes. Taking into account the host environment and the selective forces that it imposes to a parasite, the phylogeny of host(s) and parasite(s) are probably closely related, once this coevolution will be responsible for the continuity or elimination of such an interaction. Nonetheless, previous empirical experiments involving schistosomes and the intermediate host provide further support to suggest the potential for host-schistosome coevolution .
In this context, it is important to analyze the evolutionary history of protein families during screening for potential targets for drug and vaccine development. Incorporating the evolutionary perspective in drug development studies can improve our understanding regarding drug resistance and effectiveness, as well as to guide new strategies of drug discovery. Gene duplication events as well as adaptive evolution should be considered during this process, since an anti-parasitic drug could bind a single protein or in all proteins encoded by a multi-gene family . As a consequence, therapies which target a subset of genes that arose by duplication may not be effective at low doses. To solve this problem, the drug's effectiveness can be increased when a single-copy gene is targeted and its function is inactivated causing complete perturbation of a vital pathway [93, 94].
Through a systemic approach, we may accelerate the advance towards the understanding of schistosomiasis, its etiologic agents, and host-parasite interactions, optimizing the discovery of therapeutic targets to the development of new drugs and vaccines. Besides promoting a significant improvement in the functional annotation of the S. mansoni predicted proteome, our approach provided relevant information about the parasite’s genome evolution such as the identification of gene duplication events and expanded protein families, supplying important information regarding the mechanisms involved in Schistosoma’s genome evolution. Among the parasite paralog groups, we identified proteases, tetraspanins, fucosyltransferases, venom allergen-like proteins (also called as SmVAL or SCP-like), and tegumental-allergen-like proteins (SmTAL) that may be related to morphological or physiological specificities of this parasite. In addition, we strongly believe that the S. mansoni phylome data will pave the way for other, more detailed analysis, such as those that have been already performed on expanded peptidases families .
One of the remaining challenges is to understand which evasion strategies enable this parasite to survive for years in a potentially hostile environment, protected from the host immune system action and/or actively making the host response ineffective. Different mechanisms may be involved in these processes, including the generation of variant proteins by expression of micro-exon genes (MEG), which have been pointed as a potential strategy , and small non-coding RNAs which perform many essential regulatory functions .
Insights obtained through this phylogenomic approach will help us to guide forward genetic approaches to better understand the host-pathogen relationships toward to the elucidation of novel drug targets and vaccine candidates urgently needed to reduce the morbidity and mortality caused by schistosomiasis worldwide. Continuing this work, a comparative analysis involving genomic, transcriptomic, and proteomic data from other helminth species as Taenia solium, Echinococcus multiloculares, Echinococcus granulosus, Fasciola hepatica, other parasites, and vectors will provide valuable information from a system-wide perspective of a broad range of organisms, improving our understanding regarding the parasitic lifestyle.
Organisms and sequence data
Predicted proteomes from 13 fully sequenced eukaryotic genomes were downloaded from JGI Genome Projects, SchistoDB, Quest For Orthologs, WormBase, BeetleBASE, and FlyBase (Table 1). The taxon sampling was selected according to the availability of the predicted proteomes and based on the phylogenetic position of each species. The comprehensive taxa selected cover important evolutionary innovations making this dataset set especially suitable for addressing the evolutionary innovations in schistosomes in the context of metazoan evolution. Model organisms were also included to provide functional annotations that could be potentially transferred to S. mansoni homologous proteins.
To reconstruct the complete collection of phylogenetic trees for all S. mansoni proteins and their homologs in other 12 fully sequenced organisms (Table 1), we used a similar automated pipeline to that described earlier for the human proteome  (Figure 1). A local database was created containing data from the S. mansoni proteome and those of 12 other completely sequenced genomes/proteomes. Alternative splicing products and identical sequences from the S. mansoni proteome were filtered out. For each protein encoded in the S. mansoni genome (“seed”), a Smith-Waterman search  (E-value ≤ 10-5) was performed against the above mentioned database to retrieve proteins with significant sequence similarity. Sequences that aligned with a continuous region longer than 50% of the query sequence were selected and aligned using MUSCLE 3.6 , MAFFT , DIALIGN-TX , and M-Coffee  with default parameters. Positions in the alignment containing a high number of gaps were eliminated using trimAl , with a consistency cutoff of 0.1667 and a gap score cutoff of 0.1. Neighbor-joining trees were derived from the trimmed alignments using scoredist distances as implemented in BioNJ  and maximum likelihood trees were obtained as implemented in PhyML using the NJ tree as a starting point . For each “seed” protein phylogenetic reconstruction, we tested four different evolutionary models (JTT, WAG, BLOSUM62, VT, LG, CpREV, and DCMut). In all cases a discrete gamma-distribution model with four rate categories plus invariant positions was assumed, the gamma parameter and the fraction of invariant positions were estimated from the data. Tree support values were computed by approximate likelihood ratio test (aLTR) as implemented in PhyML [46, 97]. The evolutionary model best fitting the data was determined by comparing the likelihood of the used models according to the Akaike Information Criterion (AIC) .
Prediction of homology relationships
To derive orthology and paralogy relationships among S. mansoni proteins and those encoded in the other genomes included in this study we used a species-overlap algorithm as described in  and as implemented in ETE (Environment for Tree Exploration) . This algorithm uses the level of species overlap between the two daughter partitions of a given node to define it as duplication or a speciation event. The analysis starts at the protein used to generate the tree (“seed” protein) and runs through the internal nodes of the tree until it reaches the root. All the trees were rooted at the midpoint. If the two partitions share any species (if there is species overlap), the node is defined as a duplication node and the proteins are considered paralogous ones. Otherwise (if there is no overlap) the node is defined as a speciation node leading to orthologous proteins. Once all the nodes have been classified, the algorithm establishes the orthology and paralogy relationships between the “seed” protein and other proteins included in the tree according to the original definition of these terms [39, 100]. A previous study has shown that the species-overlap algorithm produces reliable orthology predictions with higher sensibility than a strict reconciliation method .
Orthology-based functional annotation
Based on the list of orthology and paralogy relationship we performed the transfer of functional annotation from each ortholog with known function to the S. mansoni “seed” proteins. To produce a confident set of functional predictions for S. mansoni proteins, we classified the list of orthologs in different subsets of orthology relationships (one-to-one, one-to-many, many-to-one, and many-to-many) between the S. mansoni proteins and the other proteins included in this phylome data. If no duplication has occurred since the speciation, the two genes form a one-to-one relationship. If subsequent duplications have occurred, other types of orthology relationships (one-to-many or many-to-many) were assigned . One example of this classification is provided (Figure 2).
To further analyze such large data set, we built the SchistoPhylomeSQL a Structured Query Language relational database using MySQL as a database management system. This local database integrates information from PhylomeDB (http://phylomedb.org) and SchistoDB (http://www.schistodb.net). Access to the database was obtained using DbVisualizer version 7.0.5 (http://www.dbvis.com), a graphical user interface that allows developing and accessing database management system (DBMSs) in different operating systems. The SchistoPhylomeSQL database was the main resource for data mining in this work. Perl scripts and SQL queries were implemented to parse the text files and load them to the database.
Detection of S. mansoni gene expansions
Using ETE , we analyzed the S. mansoni phylome data to identify protein families that were specifically expanded in the S. mansoni lineage since its diversification from the other metazoans (Additional file 1 Table S3). The duplication events defined by the species-overlap algorithm that only comprised paralogs from S. mansoni were considered lineage-specific duplications. In cases where the information extracted from more than one phylogenetic tree contained the same paralogous proteins, changing only the “seed” protein position, the data was filtered to obtain a non-redundant list of in-paralogs.
LLS: carried out the phylogenetic and functional annotation studies, and drafted the manuscript. MM: performed the phylome reconstruction and functional annotation transfer. LAN: participated in the coordination of this study, and co-wrote this manuscript. AZ: wrote the Perl scripts for data manipulation and provided computational support for this study. TG: participated in the coordination of this study, supervised the phylome reconstruction, and co-wrote this manuscript. GO: participated in the design and coordination of this study, and co-wrote this manuscript. All authors read and approved the final manuscript.
Engels D, Chitsulo L, Montresor A, Savioli L: The global epidemiological situation of schistosomiasis and new approaches to control and research. Acta Trop. 2002, 82 (2): 139-146.
Steinmann P, Keiser J, Bos R, Tanner M, Utzinger J: Schistosomiasis and water resources development: systematic review, meta-analysis, and estimates of people at risk. Lancet Infect Dis. 2006, 6 (7): 411-425.
Gryseels B: Schistosomiasis. Infect Dis Clin North Am. 2012, 26 (2): 383-397.
Organization WH: Schistosomiasis: population requiring preventive chemotherapy and number of people treated in 2010. Wkly Epidemiol Rec. 2012, 87 (4): 37-44.
Berriman M, Haas BJ, LoVerde PT, Wilson RA, Dillon GP, Cerqueira GC, Mashiyama ST, Al-Lazikani B, Andrade LF, Ashton PD, Aslett MA, Bartholomeu DC, Blandin G, Caffrey CR, Coghlan A, Coulson R, Day TA, Delcher A, DeMarco R, Djikeng A, Eyre T, Gamble JA, Ghedin E, Gu Y, Hertz-Fowler C, Hirai H, Hirai Y, Houston R, Ivens A, Johnston DA, Lacerda D, Macedo CD, McVeigh P, Ning Z, Oliveira G, Overington JP, Parkhill J, Pertea M, Pierce RJ, Protasio AV, Quail MA, Rajandream MA, Rogers J, Sajid M, Salzberg SL, Stanke M, Tivey AR, White O, Williams DL, Wortman J, Wu W, Zamanian M, Zerlotini A, Fraser-Liggett CM, Barrell BG, El-Sayed NM: The genome of the blood fluke Schistosoma mansoni. Nature. 2009, 460 (7253): 352-358.
SjGSaFA C: The Schistosoma japonicum genome reveals features of host-parasite interplay. Nature. 2009, 460 (7253): 345-351.
Young ND, Jex AR, Li B, Liu S, Yang L, Xiong Z, Li Y, Cantacessi C, Hall RS, Xu X, Chen F, Wu X, Zerlotini A, Oliveira G, Hofmann A, Zhang G, Fang X, Kang Y, Campbell BE, Loukas A, Ranganathan S, Rollinson D, Rinaldi G, Brindley PJ, Yang H, Wang J, Gasser RB: Whole-genome sequence of Schistosoma haematobium. Nat Genet. 2012, 44 (2): 221-225.
TDR: Tropical Disease Research: progress 2005–2006. 2007, Geneva: World Health Organization, Tropical Disease Research, 112-
Bruun B, Aagaard-Hansen J: The social context of schistosomiasis and its control: an introduction and annotated bibliography. 2008, Geneva: World Health Organization
Liang YS, Dai JR, Zhu YC, Coles GC, Doenhoff MJ: Genetic analysis of praziquantel resistance in Schistosoma mansoni. Southeast Asian J Trop Med Public Health. 2003, 34 (2): 274-280.
Pica-Mattoccia L, Cioli D: Sex- and stage-related sensitivity of Schistosoma mansoni to in vivo and in vitro praziquantel treatment. Int J Parasitol. 2004, 34 (4): 527-533.
Botros SS, Bennett JL: Praziquantel resistance. Expert Opin Drug Discov. 2007, 2: S35-S40.
Melman SD, Steinauer ML, Cunningham C, Kubatko LS, Mwangi IN, Wynn NB, Mutuku MW, Karanja DM, Colley DG, Black CL, Secor WE, Mkoji GM, Loker ES: Reduced susceptibility to praziquantel among naturally occurring Kenyan isolates of Schistosoma mansoni. PLoS Negl Trop Dis. 2009, 3 (8): e504-
Rokni MB: Schistosomiasis. 2002, InTech
van der Werf MJ, de Vlas SJ, Brooker S, Looman CW, Nagelkerke NJ, Habbema JD, Engels D: Quantification of clinical morbidity associated with schistosome infection in sub-Saharan Africa. Acta Trop. 2003, 86 (2–3): 125-139.
Protasio AV, Tsai IJ, Babbage A, Nichol S, Hunt M, Aslett MA, De Silva N, Velarde GS, Anderson TJ, Clark RC, Davidson C, Dillon GP, Holroyd NE, LoVerde PT, Lloyd C, McQuillan J, Oliveira G, Otto TD, Parker-Manuel SJ, Quail MA, Wilson RA, Zerlotini A, Dunne DW, Berriman M: A systematically improved high quality genome and transcriptome of the human blood fluke Schistosoma mansoni. PLoS Negl Trop Dis. 2012, 6 (1): e1455-
Zerlotini A, Heiges M, Wang H, Moraes RL, Dominitini AJ, Ruiz JC, Kissinger JC, Oliveira G: SchistoDB: a Schistosoma mansoni genome resource. Nucleic Acids Res. 2009, 37 (Database issue): 579-582.
Criscione CD, Valentim CL, Hirai H, LoVerde PT, Anderson TJ: Genomic linkage map of the human blood fluke Schistosoma mansoni. Genome Biol. 2009, 10 (6): R71-
Han ZG, Brindley PJ, Wang SY, Chen Z: Schistosoma genomics: new perspectives on schistosome biology and host-parasite interaction. Annu Rev Genomics Hum Genet. 2009, 10: 211-240.
DeMarco R, Verjovski-Almeida S: Schistosomes–proteomics studies for potential novel vaccines and drug targets. Drug Discov Today. 2009, 14 (9–10): 472-478.
Hawkins T, Kihara D: Function prediction of uncharacterized proteins. J Bioinform Comput Biol. 2007, 5 (1): 1-30.
Nahum LA, Pereira SL: Phylogenomics, Protein Family Evolution, and the Tree of Life: An Integrated Approach between Molecular Evolution and Computational Intelligence. Studies in Computational Intelligence (SCI). Edited by: Smolinski TG, Milanova MG, Hassanien AE. 2008, Berlin Heidelberg: Springer-Verlag, 259-279. vol. 122
Jiang Z: Protein function predictions based on the phylogenetic profile method. Crit Rev Biotechnol. 2008, 28 (4): 233-238.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215 (3): 403-410.
Li L, Stoeckert CJ, Roos DS: OrthoMCL: identification of ortholog groups for eukaryotic genomes. Genome Res. 2003, 13 (9): 2178-2189.
Ostlund G, Schmitt T, Forslund K, Köstler T, Messina DN, Roopra S, Frings O, Sonnhammer EL: InParanoid 7: new algorithms and tools for eukaryotic orthology analysis. Nucleic Acids Res. 2010, 38 (Database issue): 196-203.
Brenner SE: Errors in genome annotation. Trends Genet. 1999, 15 (4): 132-133.
Brown D, Sjölander K: Functional classification using phylogenomic inference. PLoS Comput Biol. 2006, 2 (6): e77-
Gabaldón T: Comparative genomics-based prediction of protein function. Methods Mol Biol. 2008, 439: 387-401.
Eisen JA, Wu M: Phylogenetic analysis and gene functional predictions: phylogenomics in action. Theor Popul Biol. 2002, 61 (4): 481-487.
Galperin MY, Koonin EV: Sources of systematic error in functional annotation of genomes: domain rearrangement, non-orthologous gene displacement and operon disruption. In Silico Biol. 1998, 1 (1): 55-67.
Koski LB, Golding GB: The closest BLAST hit is often not the nearest neighbor. J Mol Evol. 2001, 52 (6): 540-542.
Andrade LF, Nahum LA, Avelar LG, Silva LL, Zerlotini A, Ruiz JC, Oliveira G: Eukaryotic Protein Kinases (ePKs) of the Helminth Parasite Schistosoma mansoni. BMC Genomics. 2011, 12: 215-
Nahum LA, Goswami S, Serres MH: Protein families reflect the metabolic diversity of organisms and provide support for functional prediction. Physiol Genomics. 2009, 38 (3): 250-260.
Eisen JA: Phylogenomics: improving functional predictions for uncharacterized genes by evolutionary analysis. Genome Res. 1998, 8 (3): 163-167.
Gabaldón T, Dessimoz C, Huxley-Jones J, Vilella AJ, Sonnhammer EL, Lewis S: Joining forces in the quest for orthologs. Genome Biol. 2009, 10 (9): 403-
Sicheritz-Pontén T, Andersson SG: A phylogenomic approach to microbial evolution. Nucleic Acids Res. 2001, 29 (2): 545-552.
Huerta-Cepas J, Marcet-Houben M, Pignatelli M, Moya A, Gabaldón T: The pea aphid phylome: a complete catalogue of evolutionary histories and arthropod orthology and paralogy relationships for Acyrthosiphon pisum genes. Insect Mol Biol. 2010, 19 (Suppl 2): 13-21.
Huerta-Cepas J, Dopazo H, Dopazo J, Gabaldón T: The human phylome. Genome Biol. 2007, 8 (6): R109-
Edgar RC: MUSCLE: a multiple sequence alignment method with reduced time and space complexity. BMC Bioinforma. 2004, 5: 113-
Katoh K, Toh H: Recent developments in the MAFFT multiple sequence alignment program. Brief Bioinform. 2008, 9 (4): 286-298.
Subramanian AR, Kaufmann M, Morgenstern B: DIALIGN-TX: greedy and progressive approaches for segment-based multiple sequence alignment. Algorithms Mol Biol. 2008, 3: 6-
Wallace IM, O'Sullivan O, Higgins DG, Notredame C: M-Coffee: combining multiple sequence alignment methods with T-Coffee. Nucleic Acids Res. 2006, 34 (6): 1692-1699.
Capella-Gutiérrez S, Silla-Martínez JM, Gabaldón T: trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics. 2009, 25 (15): 1972-1973.
Gascuel O: BIONJ: an improved version of the NJ algorithm based on a simple model of sequence data. Mol Biol Evol. 1997, 14 (7): 685-695.
Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003, 52 (5): 696-704.
Huerta-Cepas J, Capella-Gutierrez S, Pryszcz LP, Denisov I, Kormes D, Marcet-Houben M, Gabaldón T: PhylomeDB v3.0: an expanding repository of genome-wide collections of trees, alignments and phylogeny-based orthology and paralogy predictions. Nucleic Acids Res. 2011, 39 (Database issue): 556-560.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM, Sherlock G: Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000, 25 (1): 25-29.
Apweiler R, Martin MJ, O'Donovan C, Magrane M, Alam-Faruque Y, Antunes R, Barrell D, Bely B, et al: Ongoing and future developments at the Universal Protein Resource. Nucleic Acids Res. 2011, 39 (Database issue): 214-219.
Forslund K, Pekkari I, Sonnhammer EL: Domain architecture conservation in orthologs. BMC Bioinforma. 2011, 12: 326-
Koonin EV: Orthologs, paralogs, and evolutionary genomics. Annu Rev Genet. 2005, 39: 309-338.
Kuzniar A, van Ham RC, Pongor S, Leunissen JA: The quest for orthologs: finding the corresponding gene across genomes. Trends Genet. 2008, 24 (11): 539-551.
Gerlt JA, Babbitt PC: Can sequence determine function?. Genome Biol. 2000, 1 (5): reviews0005.0001-reviews0005.0010.
Ohno S: Evolution by gene duplication. 1970, New York Heidelberg Berlin: Springer-Verlag, 1
Huerta-Cepas J, Gabaldón T: Assigning duplication events to relative temporal scales in genome-wide studies. Bioinformatics. 2011, 27 (1): 38-45.
Lawton SP, Hirai H, Ironside JE, Johnston DA, Rollinson D: Genomes and geography: genomic insights into the evolution and phylogeography of the genus Schistosoma. Parasit Vectors. 2011, 4: 131-
Webster BL, Southgate VR, Littlewood DT: A revision of the interrelationships of Schistosoma including the recently described Schistosoma guineensis. Int J Parasitol. 2006, 36 (8): 947-955.
Littlewood DT, Lockyer AE, Webster BL, Johnston DA, Le TH: The complete mitochondrial genomes of Schistosoma haematobium and Schistosoma spindale and the evolutionary history of mitochondrial genome changes among parasitic flatworms. Molecular phylogenetics and evolution. 2006, 39 (2): 452-467.
Swain MT, Larkin DM, Caffrey CR, Davies SJ, Loukas A, Skelly PJ, Hoffmann KF: Schistosoma comparative genomics: integrating genome structure, parasite biology and anthelmintic discovery. Trends Parasitol. 2011, 27 (12): 555-564.
Finn RD, Mistry J, Tate J, Coggill P, Heger A, Pollington JE, Gavin OL, Gunasekaran P, Ceric G, Forslund K, Holm L, Sonnhammer EL, Eddy SR, Bateman A: The Pfam protein families database. Nucleic Acids Res. 2010, 38 (Database issue): 211-222.
Tran MH, Pearson MS, Bethony JM, Smyth DJ, Jones MK, Duke M, Don TA, McManus DP, Correa-Oliveira R, Loukas A: Tetraspanins on the surface of Schistosoma mansoni are protective antigens against schistosomiasis. Nat Med. 2006, 12 (7): 835-840.
McManus DP, Loukas A: Current status of vaccines for schistosomiasis. Clin Microbiol Rev. 2008, 21 (1): 225-242.
Fitzpatrick JM, Peak E, Perally S, Chalmers IW, Barrett J, Yoshino TP, Ivens AC, Hoffmann KF: Anti-schistosomal intervention targets identified by lifecycle transcriptomic analyses. PLoS Negl Trop Dis. 2009, 3 (11): e543-
Fitzsimmons CM, Jones FM, Stearn A, Chalmers IW, Hoffmann KF, Wawrzyniak J, Wilson S, Kabatereine NB, Dunne DW: The Schistosoma mansoni tegumental-allergen-like (TAL) protein family: influence of developmental expression on human IgE responses. PLoS Negl Trop Dis. 2012, 6 (4): e1593-
Silva LL, Marcet-Houben M, Zerlotini A, Gabaldón T, Oliveira G, Nahum LA: Evolutionary Histories of Expanded Peptidase Families in Schistosoma mansoni. 2011, in press. Mem Inst Oswaldo Cruz
Levy S, Shoham T: Protein-protein interactions in the tetraspanin web. Physiology (Bethesda). 2005, 20: 218-224.
Yáñez-Mó M, Barreiro O, Gordon-Alonso M, Sala-Valdés M, Sánchez-Madrid F: Tetraspanin-enriched microdomains: a functional unit in cell plasma membranes. Trends Cell Biol. 2009, 19 (9): 434-446.
Tran MH, Freitas TC, Cooper L, Gaze S, Gatton ML, Jones MK, Lovas E, Pearce EJ, Loukas A: Suppression of mRNAs encoding tegument tetraspanins from Schistosoma mansoni results in impaired tegument turnover. PLoS Pathog. 2010, 6 (4): e1000840-
Lee KW, Shalaby KA, Medhat AM, Shi H, Yang Q, Karim AM, LoVerde PT: Schistosoma mansoni: characterization of the gene encoding Sm23, an integral membrane protein. Exp Parasitol. 1995, 80 (1): 155-158.
Inal J, Bickle Q: Sequence and immunogenicity of the 23-kDa transmembrane antigen of Schistosoma haematobium. Mol Biochem Parasitol. 1995, 74 (2): 217-221.
Chalmers IW, McArdle AJ, Coulson RM, Wagner MA, Schmid R, Hirai H, Hoffmann KF: Developmentally regulated expression, alternative splicing and distinct sub-groupings in members of the Schistosoma mansoni venom allergen-like (SmVAL) gene family. BMC Genomics. 2008, 9: 89-
Curwen RS, Ashton PD, Sundaralingam S, Wilson RA: Identification of novel proteases and immunomodulators in the secretions of schistosome cercariae that facilitate host entry. Mol Cell Proteomics. 2006, 5 (5): 835-844.
Ma B, Simala-Grant JL, Taylor DE: Fucosylation in prokaryotes and eukaryotes. Glycobiology. 2006, 16 (12): 158R-184R.
Marques ET, Ichikawa Y, Strand M, August JT, Hart GW, Schnaar RL: Fucosyltransferases in Schistosoma mansoni development. Glycobiology. 2001, 11 (3): 249-259.
Vichasri-Grams S, Subpipattana P, Sobhon P, Viyanant V, Grams R: An analysis of the calcium-binding protein 1 of Fasciola gigantica with a comparison to its homologs in the phylum Platyhelminthes. Mol Biochem Parasitol. 2006, 146 (1): 10-23.
Dunne DW, Butterworth AE, Fulford AJ, Kariuki HC, Langley JG, Ouma JH, Capron A, Pierce RJ, Sturrock RF: Immunity after treatment of human schistosomiasis: association between IgE antibodies to adult worm antigens and resistance to reinfection. Eur J Immunol. 1992, 22 (6): 1483-1494.
Salter JP, Choe Y, Albrecht H, Franklin C, Lim KC, Craik CS, McKerrow JH: Cercarial elastase is encoded by a functionally conserved gene family across multiple species of schistosomes. J Biol Chem. 2002, 277 (27): 24618-24624.
Aslam A, Quinn P, McIntosh RS, Shi J, Ghumra A, McKerrow JH, Bunting KA, Dunne DW, Doenhoff MJ, Morrison SL, Zhang K, Pleass RJ: Proteases from Schistosoma mansoni cercariae cleave IgE at solvent exposed interdomain regions. Mol Immunol. 2008, 45 (2): 567-574.
Hockley DJ, McLaren DJ: Schistosoma mansoni: changes in the outer membrane of the tegument during development from cercaria to adult worm. Int J Parasitol. 1973, 3 (1): 13-25.
Pearce EJ, Sher A: Mechanisms of immune evasion in schistosomiasis. Contrib Microbiol Immunol. 1987, 8: 219-232.
Braschi S, Borges WC, Wilson RA: Proteomic analysis of the schistosome tegument and its surface membranes. Mem Inst Oswaldo Cruz. 2006, 101 (Suppl 1): 205-212.
Mclaren DJ, Hockley DJ: Blood flukes have a double outer membrane. Nature. 1977, 269 (5624): 147-149.
Escobedo G, Roberts CW, Carrero JC, Morales-Montor J: Parasite regulation by host hormones: an old mechanism of host exploitation?. Trends Parasitol. 2005, 21 (12): 588-593.
Rescher U, Gerke V: Annexins–unique membrane binding proteins with diverse functions. J Cell Sci. 2004, 117 (Pt 13): 2631-2639.
Tararam CA, Farias LP, Wilson RA, Leite LC: Schistosoma mansoni Annexin 2: molecular characterization and immunolocalization. Exp Parasitol. 2010, 126 (2): 146-155.
Hofmann A, Osman A, Leow CY, Driguez P, McManus DP, Jones MK: Parasite annexins–new molecules with potential for drug and vaccine development. BioEssays. 2010, 32 (11): 967-976.
Takeichi M: Cadherins: a molecular family important in selective cell-cell adhesion. Annu Rev Biochem. 1990, 59: 237-252.
Phelan P, Bacon JP, Davies JA, Stebbings LA, Todman MG, Avery L, Baines RA, Barnes TM, Ford C, Hekimi S, Lee R, Shaw JE, Starich TA, Curtin KD, Sun YA, Wyman RJ: Innexins: a family of invertebrate gap-junction proteins. Trends Genet. 1998, 14 (9): 348-349.
Yeh E, Kawano T, Ng S, Fetter R, Hung W, Wang Y, Zhen M: Caenorhabditis elegans innexins regulate active zone differentiation. J Neurosci. 2009, 29 (16): 5207-5217.
Zhang J: Evolution by gene duplication: an update. Trends Ecol Evol. 2003, 18 (6): 292-298.
Rubin GM, Yandell MD, Wortman JR, Gabor Miklos GL, Nelson CR, Hariharan IK, Fortini ME, Li PW, Apweiler R, Fleischmann W, Cherry JM, Henikoff S, Skupski MP, Misra S, Ashburner M, Birney E, Boguski MS, Brody T, Brokstein P, Celniker SE, Chervitz SA, Coates D, Cravchik A, Gabrielian A, Galle RF, Gelbart WM, George RA, Goldstein LS, Gong F, Guan P, Harris NL, Hay BA, Hoskins RA, Li J, Li Z, Hynes RO, Jones SJ, Kuehl PM, Lemaitre B, Littleton JT, Morrison DK, Mungall C, O'Farrell PH, Pickeral OK, Shue C, Vosshall LB, Zhang J, Zhao Q, Zheng XH, Lewis S: Comparative genomics of the eukaryotes. Science. 2000, 287 (5461): 2204-2215.
Webster JP, Shrivastava J, Johnson PJ, Blair L: Is host-schistosome coevolution going anywhere?. BMC Evol Biol. 2007, 7: 91-
Emes RD, Yang Z: Duplicated paralogous genes subject to positive selection in the genome of Trypanosoma brucei. PLoS One. 2008, 3 (5): e2295-
DeMarco R, Mathieson W, Manuel SJ, Dillon GP, Curwen RS, Ashton PD, Ivens AC, Berriman M, Verjovski-Almeida S, Wilson RA: Protein variation in blood-dwelling schistosome worms generated by differential splicing of micro-exon gene transcripts. Genome Res. 2010, 20 (8): 1112-1121.
Oliveira KC, Carvalho ML, Maracaja-Coutinho V, Kitajima JP, Verjovski-Almeida S: Non-coding RNAs in schistosomes: an unexplored world. An Acad Bras Cienc. 2011, 83 (2): 673-694.
Smith TF, Waterman MS: Identification of common molecular subsequences. J Mol Biol. 1981, 147 (1): 195-197.
Anisimova M, Gascuel O: Approximate likelihood-ratio test for branches: A fast, accurate, and powerful alternative. Syst Biol. 2006, 55 (4): 539-552.
Akaike H: Information theory and an extension of the maximum likelihood principle. 2nd International Symposium on Information Theory 1973. Edited by: Petrov BN, Csàki F. 1973, Budapest: Akademiai kiado, 267-281.
Huerta-Cepas J, Dopazo J, Gabaldón T: ETE: a python Environment for Tree Exploration. BMC Bioinforma. 2010, 11: 24-
Fitch WM: Distinguishing homologous from analogous proteins. Syst Zool. 1970, 19 (2): 99-113.
Marcet-Houben M, Gabaldón T: The tree versus the forest: the fungal tree of life and the topological diversity within the yeast phylome. PLoS One. 2009, 4 (2): e4357-
The authors acknowledge the use of the computing resources of the Centre de Regulació Genòmica (CRG, Spain) and the Center for Excellence in Bioinformatics (CEBio/Fiocruz, Brazil). The authors are grateful to Jaime Huerta-Cepas (CRG), Eric Aguiar, and Francislon Silva (CEBio) for bioinformatics technical support and to Mariana de Oliveira (CEBio) for help with illustrations in this work. This work was funded by the National Institutes of Health - NIH/Fogarty International Center (TW007012 to GO), the National Council for Research and Development - CNPq (CNPq Research Fellowship 306879/2009-3 to GO, INCT-DT 573839/2008-5 to GO, and CNPq-Universal 476036/2010-0 to LAN), the Spanish Ministry of Science and Innovation (BFU2009-09168 to TG), and the Research Foundation of the State of Minas Gerais - FAPEMIG (CBB-1181/08 and PPM-00439-10 to GO).
The authors declare that they have no competing interests.
Electronic supplementary material
About this article
- Maximum likelihood analysis
- Homology prediction
- Functional annotation
- Paralogous families
- Parasite genomics