Email updates

Keep up to date with the latest news and content from BMC Genomics and BioMed Central.

Open Access Highly Accessed Research article

The Schistosoma mansoni phylome: using evolutionary genomics to gain insight into a parasite’s biology

Larissa Lopes Silva123, Marina Marcet-Houben45, Laila Alves Nahum126, Adhemar Zerlotini27, Toni Gabaldón45 and Guilherme Oliveira12*

Author Affiliations

1 Grupo de Genômica e Biologia Computacional, Centro de Pesquisas René Rachou. Instituto Nacional de Ciência e Tecnologia em Doenças Tropicais. Fundação Oswaldo Cruz - FIOCRUZ, Belo Horizonte, MG, 30190-002, Brazil

2 Centro de Excelência em Bioinformática, Fundação Oswaldo Cruz – FIOCRUZ, Belo Horizonte, MG, Brazil

3 Instituto de Ciências Biológicas, Universidade Federal de Minas Gerais – UFMG, Belo Horizonte, MG, Brazil

4 Bioinformatics and Genomics Programme, Centre for Genomic Regulation (CRG), Dr. Aiguader, 88, 08003, Barcelona, Spain

5 Universitat Pompeu Fabra (UPF), 08003, Barcelona, Spain

6 Faculdade Infórium de Tecnologia, Belo Horizonte, MG, 30130-180, Brazil

7 Laboratório Multiusuário de Bioinformática, Embrapa Informática Agropecuária, Campinas, São Paulo, Brazil

For all author emails, please log on.

BMC Genomics 2012, 13:617  doi:10.1186/1471-2164-13-617

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2164/13/617


Received:30 September 2012
Accepted:22 October 2012
Published:13 November 2012

© 2012 Silva et al.; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

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.

Results

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 webcite).

Conclusions

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.

Keywords:
Phylogenomics; Maximum likelihood analysis; Homology prediction; Functional annotation; Paralogous families; Parasite genomics; Schistosomiasis

Background

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 [16]. 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 [24], OrthoMCL [25], InParanoid [26]). 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. [24]). 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 [28]. 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 [37].

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 [37]. 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 [44], neighbor-joining tree building [45], evolutionary model testing, and maximum likelihood analysis [46]. 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.

Table 1. Proteomes selected for the S. mansoni phylome reconstruction

To perform the phylogenetic analyses, we applied an automated pipeline similar to the one used for the human phylome project [39]. This pipeline is illustrated here (Figure 1). The resulting alignments, phylogenies, and orthology predictions can be accessed at PhylomeDB [47] (http://phylomedb.org webcite).

thumbnailFigure 1. Pipeline used to reconstruct and analyze the S. mansoni phylome. Each protein sequence encoded in the parasite genome was compared against a database of proteins from other 12 fully sequenced eukaryotic proteomes (Table 1) to select putative homologous proteins. Groups of potential homologs were aligned and subsequently trimmed to remove gap-rich regions. The refined alignment was used to build a NJ tree, which was then used as a “seed” tree to perform a ML likelihood analysis as implemented in PhyML. In the ML analysis, up to five different evolutionary models were tested and the model best fitting to the data was determined by the Akaike Information Criterion (AIC). Different algorithms were used to identify homology relationships and lineage-specific duplications. To extract and interpret the large data set obtained a Structured Query Language (SQL) relational database was built. This database was the main resource for data mining in this work. Adapted from [39].

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%) [38].

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 [48] assignment in SchistoDB [17].

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 [39]. The comprehensive catalogue of phylogeny-based orthology and paralogy relationships among S. mansoni and other species was made publicly available at PhylomeDB [47].

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).

thumbnailFigure 2. Homology relationships and evolutionary events inferred from the analysis of a S. mansoni protein.A) Phylogenetic tree reconstructed for the parasite “seed” protein Phy000V0I5_SCHMA (Smp_175750). B) Homology relationships identified between the “seed” protein and its homologs in the other species.

By using such phylogeny-based approach, we transferred 10,175 functional annotations (GO terms [48]) 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.

Additional file 1. Table S1: Functional prediction for theS. mansoni "seed" proteins. Table S2: Comparison between reviewed proteins with functional annotation at UniProt and GO terms transferred by phylogenomic approach. Table S3: Homologous proteins to the S. mansoni "seed" proteins.

Format: XLS Size: 10.3MB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

To validate the applied methodology, we retrieved reviewed S. mansoni proteins from UniProt [49], 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 [50].

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 [49], it is one of the orthologs with known function and experimental validation. Tags for homologous sequences with experimental validation are not available in PhylomeDB [47]. However, links to UniProt [49] and other databases are provided.

thumbnailFigure 3. Example of functional prediction based on phylogenetic analysis. The protein sequences are represented by the internal identifier in PhylomeDB. Relationships among the parasite Phy000V14T_SCHMA “seed” protein (Smp_170950) and its homologs in other species (Table 1) as inferred by maximum likelihood method implemented in PhyML. Support values were computed by approximate likelihood ratio test (aLTR). Curly brackets hold Gene Ontology (GO) terms for proteins in this dataset.

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 [53].

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 [54], 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 [49], Pfam [60], or non-redundant (nr) NCBI database (ftp://ftp.ncbi.nih.gov/blast/db webcite).

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 [60]. 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 [61]. 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 [68]. 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].

thumbnailFigure 4. Phylogenetic relationships of schistosome lineage-specific duplicated tetraspanins. Analysis was performed with trimmed sequence alignment by using the maximum likelihood method as implemented in PhyML. Best fit model (WAG) and support values for each node were estimated by the Akaike Likelihood Ratio Test (aLRT). Sequence labels follow the PhylomeDB internal identifier. For details, see supplementary data (Additional file 1 Table S3).

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 [71]. 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 [73]. 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 [64]. 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 [65].

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 [89].

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% [90], whereas in C. elegans this value is equal to 49% [91].

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 [92].

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 [93]. 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].

Conclusions

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 [65].

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 [94], and small non-coding RNAs which perform many essential regulatory functions [95].

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.

Methods

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.

Phylome reconstruction

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 [39] (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 [96] (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 [40], MAFFT [41], DIALIGN-TX [42], and M-Coffee [43] with default parameters. Positions in the alignment containing a high number of gaps were eliminated using trimAl [44], 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 [45] and maximum likelihood trees were obtained as implemented in PhyML using the NJ tree as a starting point [46]. 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) [98].

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 [39] and as implemented in ETE (Environment for Tree Exploration) [99]. 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 [101].

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 [51]. 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 webcite) and SchistoDB (http://www.schistodb.net webcite). Access to the database was obtained using DbVisualizer version 7.0.5 (http://www.dbvis.com webcite), 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 [99], 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.

Competing interests

The authors declare that they have no competing interests.

Author’s contributions

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.

Acknowledgments

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).

References

  1. 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. OpenURL

  2. 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. OpenURL

  3. Gryseels B: Schistosomiasis.

    Infect Dis Clin North Am 2012, 26(2):383-397. OpenURL

  4. Organization WH: Schistosomiasis: population requiring preventive chemotherapy and number of people treated in 2010.

    Wkly Epidemiol Rec 2012, 87(4):37-44. OpenURL

  5. 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. OpenURL

  6. SjGSaFA C: The Schistosoma japonicum genome reveals features of host-parasite interplay.

    Nature 2009, 460(7253):345-351. OpenURL

  7. 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. OpenURL

  8. TDR: Tropical Disease Research: progress 2005–2006. Geneva: World Health Organization, Tropical Disease Research; 2007:112. OpenURL

  9. Bruun B, Aagaard-Hansen J: The social context of schistosomiasis and its control: an introduction and annotated bibliography. Geneva: World Health Organization; 2008. OpenURL

  10. 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. OpenURL

  11. 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. OpenURL

  12. Botros SS, Bennett JL: Praziquantel resistance.

    Expert Opin Drug Discov 2007, 2:S35-S40. OpenURL

  13. 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. OpenURL

  14. Rokni MB: Schistosomiasis. InTech; 2002. OpenURL

  15. 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. OpenURL

  16. 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. OpenURL

  17. 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. OpenURL

  18. 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. OpenURL

  19. 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. OpenURL

  20. DeMarco R, Verjovski-Almeida S: Schistosomes–proteomics studies for potential novel vaccines and drug targets.

    Drug Discov Today 2009, 14(9–10):472-478. OpenURL

  21. Hawkins T, Kihara D: Function prediction of uncharacterized proteins.

    J Bioinform Comput Biol 2007, 5(1):1-30. OpenURL

  22. Nahum LA, Pereira SL: Phylogenomics, Protein Family Evolution, and the Tree of Life: An Integrated Approach between Molecular Evolution and Computational Intelligence. In Studies in Computational Intelligence (SCI). vol. 122 edition. Edited by Smolinski TG, Milanova MG, Hassanien AE. Berlin Heidelberg: Springer-Verlag; 2008:259-279. OpenURL

  23. Jiang Z: Protein function predictions based on the phylogenetic profile method.

    Crit Rev Biotechnol 2008, 28(4):233-238. OpenURL

  24. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool.

    J Mol Biol 1990, 215(3):403-410. OpenURL

  25. Li L, Stoeckert CJ, Roos DS: OrthoMCL: identification of ortholog groups for eukaryotic genomes.

    Genome Res 2003, 13(9):2178-2189. OpenURL

  26. 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. OpenURL

  27. Brenner SE: Errors in genome annotation.

    Trends Genet 1999, 15(4):132-133. OpenURL

  28. Brown D, Sjölander K: Functional classification using phylogenomic inference.

    PLoS Comput Biol 2006, 2(6):e77. OpenURL

  29. Gabaldón T: Comparative genomics-based prediction of protein function.

    Methods Mol Biol 2008, 439:387-401. OpenURL

  30. Eisen JA, Wu M: Phylogenetic analysis and gene functional predictions: phylogenomics in action.

    Theor Popul Biol 2002, 61(4):481-487. OpenURL

  31. 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. OpenURL

  32. Koski LB, Golding GB: The closest BLAST hit is often not the nearest neighbor.

    J Mol Evol 2001, 52(6):540-542. OpenURL

  33. 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. OpenURL

  34. 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. OpenURL

  35. Eisen JA: Phylogenomics: improving functional predictions for uncharacterized genes by evolutionary analysis.

    Genome Res 1998, 8(3):163-167. OpenURL

  36. 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. OpenURL

  37. Sicheritz-Pontén T, Andersson SG: A phylogenomic approach to microbial evolution.

    Nucleic Acids Res 2001, 29(2):545-552. OpenURL

  38. 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. OpenURL

  39. Huerta-Cepas J, Dopazo H, Dopazo J, Gabaldón T: The human phylome.

    Genome Biol 2007, 8(6):R109. OpenURL

  40. Edgar RC: MUSCLE: a multiple sequence alignment method with reduced time and space complexity.

    BMC Bioinforma 2004, 5:113. OpenURL

  41. Katoh K, Toh H: Recent developments in the MAFFT multiple sequence alignment program.

    Brief Bioinform 2008, 9(4):286-298. OpenURL

  42. Subramanian AR, Kaufmann M, Morgenstern B: DIALIGN-TX: greedy and progressive approaches for segment-based multiple sequence alignment.

    Algorithms Mol Biol 2008, 3:6. OpenURL

  43. 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. OpenURL

  44. 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. OpenURL

  45. 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. OpenURL

  46. Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood.

    Syst Biol 2003, 52(5):696-704. OpenURL

  47. 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. OpenURL

  48. 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. OpenURL

  49. 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. OpenURL

  50. Forslund K, Pekkari I, Sonnhammer EL: Domain architecture conservation in orthologs.

    BMC Bioinforma 2011, 12:326. OpenURL

  51. Koonin EV: Orthologs, paralogs, and evolutionary genomics.

    Annu Rev Genet 2005, 39:309-338. OpenURL

  52. 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. OpenURL

  53. Gerlt JA, Babbitt PC: Can sequence determine function?

    Genome Biol 2000, 1(5):reviews0005.0001-reviews0005.0010. OpenURL

  54. Ohno S: Evolution by gene duplication. 1st edition. New York Heidelberg Berlin: Springer-Verlag; 1970. OpenURL

  55. Huerta-Cepas J, Gabaldón T: Assigning duplication events to relative temporal scales in genome-wide studies.

    Bioinformatics 2011, 27(1):38-45. OpenURL

  56. 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. OpenURL

  57. 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. OpenURL

  58. 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. OpenURL

  59. 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. OpenURL

  60. 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. OpenURL

  61. 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. OpenURL

  62. McManus DP, Loukas A: Current status of vaccines for schistosomiasis.

    Clin Microbiol Rev 2008, 21(1):225-242. OpenURL

  63. 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. OpenURL

  64. 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. OpenURL

  65. Silva LL, Marcet-Houben M, Zerlotini A, Gabaldón T, Oliveira G, Nahum LA: Evolutionary Histories of Expanded Peptidase Families in Schistosoma mansoni. in press. Mem Inst Oswaldo Cruz; 2011. OpenURL

  66. Levy S, Shoham T: Protein-protein interactions in the tetraspanin web.

    Physiology (Bethesda) 2005, 20:218-224. OpenURL

  67. 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. OpenURL

  68. 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. OpenURL

  69. 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. OpenURL

  70. Inal J, Bickle Q: Sequence and immunogenicity of the 23-kDa transmembrane antigen of Schistosoma haematobium.

    Mol Biochem Parasitol 1995, 74(2):217-221. OpenURL

  71. 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. OpenURL

  72. 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. OpenURL

  73. Ma B, Simala-Grant JL, Taylor DE: Fucosylation in prokaryotes and eukaryotes.

    Glycobiology 2006, 16(12):158R-184R. OpenURL

  74. Marques ET Jr, Ichikawa Y, Strand M, August JT, Hart GW, Schnaar RL: Fucosyltransferases in Schistosoma mansoni development.

    Glycobiology 2001, 11(3):249-259. OpenURL

  75. 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. OpenURL

  76. 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. OpenURL

  77. 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. OpenURL

  78. 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. OpenURL

  79. 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. OpenURL

  80. Pearce EJ, Sher A: Mechanisms of immune evasion in schistosomiasis.

    Contrib Microbiol Immunol 1987, 8:219-232. OpenURL

  81. 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. OpenURL

  82. Mclaren DJ, Hockley DJ: Blood flukes have a double outer membrane.

    Nature 1977, 269(5624):147-149. OpenURL

  83. 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. OpenURL

  84. Rescher U, Gerke V: Annexins–unique membrane binding proteins with diverse functions.

    J Cell Sci 2004, 117(Pt 13):2631-2639. OpenURL

  85. Tararam CA, Farias LP, Wilson RA, Leite LC: Schistosoma mansoni Annexin 2: molecular characterization and immunolocalization.

    Exp Parasitol 2010, 126(2):146-155. OpenURL

  86. 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. OpenURL

  87. Takeichi M: Cadherins: a molecular family important in selective cell-cell adhesion.

    Annu Rev Biochem 1990, 59:237-252. OpenURL

  88. 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. OpenURL

  89. 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. OpenURL

  90. Zhang J: Evolution by gene duplication: an update.

    Trends Ecol Evol 2003, 18(6):292-298. OpenURL

  91. 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. OpenURL

  92. Webster JP, Shrivastava J, Johnson PJ, Blair L: Is host-schistosome coevolution going anywhere?

    BMC Evol Biol 2007, 7:91. OpenURL

  93. Emes RD, Yang Z: Duplicated paralogous genes subject to positive selection in the genome of Trypanosoma brucei.

    PLoS One 2008, 3(5):e2295. OpenURL

  94. 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. OpenURL

  95. 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. OpenURL

  96. Smith TF, Waterman MS: Identification of common molecular subsequences.

    J Mol Biol 1981, 147(1):195-197. OpenURL

  97. Anisimova M, Gascuel O: Approximate likelihood-ratio test for branches: A fast, accurate, and powerful alternative.

    Syst Biol 2006, 55(4):539-552. OpenURL

  98. Akaike H: Information theory and an extension of the maximum likelihood principle. In 2nd International Symposium on Information Theory 1973. Edited by Petrov BN, Csàki F. Budapest: Akademiai kiado; 1973:267-281. OpenURL

  99. Huerta-Cepas J, Dopazo J, Gabaldón T: ETE: a python Environment for Tree Exploration.

    BMC Bioinforma 2010, 11:24. OpenURL

  100. Fitch WM: Distinguishing homologous from analogous proteins.

    Syst Zool 1970, 19(2):99-113. OpenURL

  101. 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. OpenURL