Open Access

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

  • Larissa Lopes Silva1, 2, 3,
  • Marina Marcet-Houben4, 5,
  • Laila Alves Nahum1, 2, 6,
  • Adhemar Zerlotini2, 7,
  • Toni Gabaldón4, 5 and
  • Guilherme Oliveira1, 2Email author
BMC Genomics201213:617

DOI: 10.1186/1471-2164-13-617

Received: 30 September 2012

Accepted: 22 October 2012

Published: 13 November 2012

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

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 [14]. The genomes of these parasites have been recently published providing insights into parasite’s development, infection, and host-parasite interactions [57]. 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 [1013].

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 [2123]. 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, 2729]. 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, 2732].

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 [4043] 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

Scientific Name

UniProt Species Code1

TaxID2

Proteins3

Source4

Download

Monosiga brevicollis

MONBE

81824

9,170

JGI

2011-06-01

Ciona Intestinalis

CIOIN

7719

14,048

UniProt Reference Proteomes

2011-07-09

Nematostella vectensis

NEMVE

45351

24,424

UniProt Reference Proteomes

2011-07-09

Schistosoma haematobium

SCHHA

6185

12,767

SchistoDB

2012-03-09

Schistosoma mansoni

SCHMA

6183

11,103

SchistoDB

2012-03-09

Schistosoma japonicum

SCHJA

6182

12,636

SchistoDB

2012-03-09

Caenorhabditis elegans

CAEEL

6239

19,758

UniProt Reference Proteomes

2011-07-09

Ascaris suum

ASCSU

6253

18,430

WormBase

2012-03-09

Brugia malayi

BRUMA

6279

19,916

WormBase

2012-03-09

Trichinella spiralis

TRISP

6334

15,878

WormBase

2012-03-09

Drosphila melanogaster

DROME

7227

11,794

FlyBase

2011-09-13

Tribolium castaneum

TRICA

7070

16,533

BeetleBASE - HGSC

2011-12-16

Homo sapiens

HUMAN

9606

20,965

UniProt Reference Proteomes

2011-07-09

1 - Code assigned to each species in the S. mansoni phylome. 2 - Taxonomic identifier at NCBI (TaxID). 3 - Number of proteins analyzed per species. 4 - Database from which the protein data were retrieved.

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).
https://static-content.springer.com/image/art%3A10.1186%2F1471-2164-13-617/MediaObjects/12864_2012_Article_4536_Fig1_HTML.jpg
Figure 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).
https://static-content.springer.com/image/art%3A10.1186%2F1471-2164-13-617/MediaObjects/12864_2012_Article_4536_Fig2_HTML.jpg
Figure 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.

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.
https://static-content.springer.com/image/art%3A10.1186%2F1471-2164-13-617/MediaObjects/12864_2012_Article_4536_Fig3_HTML.jpg
Figure 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[5659]. 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).

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, 6165]. 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].
https://static-content.springer.com/image/art%3A10.1186%2F1471-2164-13-617/MediaObjects/12864_2012_Article_4536_Fig4_HTML.jpg
Figure 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 [7981]. 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 [8486]. 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) 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 [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.

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.

Declarations

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

Authors’ 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
(2)
Centro de Excelência em Bioinformática, Fundação Oswaldo Cruz – FIOCRUZ
(3)
Instituto de Ciências Biológicas, Universidade Federal de Minas Gerais – UFMG
(4)
Bioinformatics and Genomics Programme, Centre for Genomic Regulation (CRG)
(5)
Universitat Pompeu Fabra (UPF)
(6)
Faculdade Infórium de Tecnologia
(7)
Laboratório Multiusuário de Bioinformática, Embrapa Informática Agropecuária

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.PubMed
  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.PubMed
  3. Gryseels B: Schistosomiasis. Infect Dis Clin North Am. 2012, 26 (2): 383-397.PubMed
  4. Organization WH: Schistosomiasis: population requiring preventive chemotherapy and number of people treated in 2010. Wkly Epidemiol Rec. 2012, 87 (4): 37-44.
  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.PubMed CentralPubMed
  6. SjGSaFA C: The Schistosoma japonicum genome reveals features of host-parasite interplay. Nature. 2009, 460 (7253): 345-351.
  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.PubMed
  8. TDR: Tropical Disease Research: progress 2005–2006. 2007, Geneva: World Health Organization, Tropical Disease Research, 112-
  9. Bruun B, Aagaard-Hansen J: The social context of schistosomiasis and its control: an introduction and annotated bibliography. 2008, Geneva: World Health Organization
  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.PubMed
  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.PubMed
  12. Botros SS, Bennett JL: Praziquantel resistance. Expert Opin Drug Discov. 2007, 2: S35-S40.PubMed
  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-PubMed CentralPubMed
  14. Rokni MB: Schistosomiasis. 2002, InTech
  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.PubMed
  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-PubMed CentralPubMed
  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.
  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-PubMed CentralPubMed
  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.PubMed
  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.PubMed
  21. Hawkins T, Kihara D: Function prediction of uncharacterized proteins. J Bioinform Comput Biol. 2007, 5 (1): 1-30.PubMed
  22. 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
  23. Jiang Z: Protein function predictions based on the phylogenetic profile method. Crit Rev Biotechnol. 2008, 28 (4): 233-238.PubMed
  24. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215 (3): 403-410.PubMed
  25. Li L, Stoeckert CJ, Roos DS: OrthoMCL: identification of ortholog groups for eukaryotic genomes. Genome Res. 2003, 13 (9): 2178-2189.PubMed CentralPubMed
  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.
  27. Brenner SE: Errors in genome annotation. Trends Genet. 1999, 15 (4): 132-133.PubMed
  28. Brown D, Sjölander K: Functional classification using phylogenomic inference. PLoS Comput Biol. 2006, 2 (6): e77-PubMed CentralPubMed
  29. Gabaldón T: Comparative genomics-based prediction of protein function. Methods Mol Biol. 2008, 439: 387-401.PubMed
  30. Eisen JA, Wu M: Phylogenetic analysis and gene functional predictions: phylogenomics in action. Theor Popul Biol. 2002, 61 (4): 481-487.PubMed
  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.PubMed
  32. Koski LB, Golding GB: The closest BLAST hit is often not the nearest neighbor. J Mol Evol. 2001, 52 (6): 540-542.PubMed
  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-PubMed CentralPubMed
  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.PubMed
  35. Eisen JA: Phylogenomics: improving functional predictions for uncharacterized genes by evolutionary analysis. Genome Res. 1998, 8 (3): 163-167.PubMed
  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-PubMed CentralPubMed
  37. Sicheritz-Pontén T, Andersson SG: A phylogenomic approach to microbial evolution. Nucleic Acids Res. 2001, 29 (2): 545-552.PubMed CentralPubMed
  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.PubMed
  39. Huerta-Cepas J, Dopazo H, Dopazo J, Gabaldón T: The human phylome. Genome Biol. 2007, 8 (6): R109-PubMed CentralPubMed
  40. Edgar RC: MUSCLE: a multiple sequence alignment method with reduced time and space complexity. BMC Bioinforma. 2004, 5: 113-
  41. Katoh K, Toh H: Recent developments in the MAFFT multiple sequence alignment program. Brief Bioinform. 2008, 9 (4): 286-298.PubMed
  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-PubMed CentralPubMed
  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.PubMed CentralPubMed
  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.PubMed CentralPubMed
  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.PubMed
  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.PubMed
  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.
  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.PubMed CentralPubMed
  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.
  50. Forslund K, Pekkari I, Sonnhammer EL: Domain architecture conservation in orthologs. BMC Bioinforma. 2011, 12: 326-
  51. Koonin EV: Orthologs, paralogs, and evolutionary genomics. Annu Rev Genet. 2005, 39: 309-338.PubMed
  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.PubMed
  53. Gerlt JA, Babbitt PC: Can sequence determine function?. Genome Biol. 2000, 1 (5): reviews0005.0001-reviews0005.0010.
  54. Ohno S: Evolution by gene duplication. 1970, New York Heidelberg Berlin: Springer-Verlag, 1
  55. Huerta-Cepas J, Gabaldón T: Assigning duplication events to relative temporal scales in genome-wide studies. Bioinformatics. 2011, 27 (1): 38-45.PubMed
  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-PubMed CentralPubMed
  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.PubMed
  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.PubMed
  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.PubMed CentralPubMed
  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.
  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.PubMed
  62. McManus DP, Loukas A: Current status of vaccines for schistosomiasis. Clin Microbiol Rev. 2008, 21 (1): 225-242.PubMed CentralPubMed
  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-PubMed CentralPubMed
  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-PubMed CentralPubMed
  65. 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
  66. Levy S, Shoham T: Protein-protein interactions in the tetraspanin web. Physiology (Bethesda). 2005, 20: 218-224.
  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.PubMed
  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-PubMed CentralPubMed
  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.PubMed
  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.PubMed
  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-PubMed CentralPubMed
  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.PubMed
  73. Ma B, Simala-Grant JL, Taylor DE: Fucosylation in prokaryotes and eukaryotes. Glycobiology. 2006, 16 (12): 158R-184R.PubMed
  74. Marques ET, Ichikawa Y, Strand M, August JT, Hart GW, Schnaar RL: Fucosyltransferases in Schistosoma mansoni development. Glycobiology. 2001, 11 (3): 249-259.
  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.PubMed
  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.PubMed
  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.PubMed
  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.PubMed
  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.PubMed
  80. Pearce EJ, Sher A: Mechanisms of immune evasion in schistosomiasis. Contrib Microbiol Immunol. 1987, 8: 219-232.PubMed
  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.PubMed
  82. Mclaren DJ, Hockley DJ: Blood flukes have a double outer membrane. Nature. 1977, 269 (5624): 147-149.PubMed
  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.PubMed
  84. Rescher U, Gerke V: Annexins–unique membrane binding proteins with diverse functions. J Cell Sci. 2004, 117 (Pt 13): 2631-2639.PubMed
  85. Tararam CA, Farias LP, Wilson RA, Leite LC: Schistosoma mansoni Annexin 2: molecular characterization and immunolocalization. Exp Parasitol. 2010, 126 (2): 146-155.PubMed
  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.PubMed
  87. Takeichi M: Cadherins: a molecular family important in selective cell-cell adhesion. Annu Rev Biochem. 1990, 59: 237-252.PubMed
  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.PubMed CentralPubMed
  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.PubMed
  90. Zhang J: Evolution by gene duplication: an update. Trends Ecol Evol. 2003, 18 (6): 292-298.
  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.PubMed CentralPubMed
  92. Webster JP, Shrivastava J, Johnson PJ, Blair L: Is host-schistosome coevolution going anywhere?. BMC Evol Biol. 2007, 7: 91-PubMed CentralPubMed
  93. Emes RD, Yang Z: Duplicated paralogous genes subject to positive selection in the genome of Trypanosoma brucei. PLoS One. 2008, 3 (5): e2295-PubMed CentralPubMed
  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.PubMed CentralPubMed
  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.PubMed
  96. Smith TF, Waterman MS: Identification of common molecular subsequences. J Mol Biol. 1981, 147 (1): 195-197.PubMed
  97. Anisimova M, Gascuel O: Approximate likelihood-ratio test for branches: A fast, accurate, and powerful alternative. Syst Biol. 2006, 55 (4): 539-552.PubMed
  98. 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.
  99. Huerta-Cepas J, Dopazo J, Gabaldón T: ETE: a python Environment for Tree Exploration. BMC Bioinforma. 2010, 11: 24-
  100. Fitch WM: Distinguishing homologous from analogous proteins. Syst Zool. 1970, 19 (2): 99-113.PubMed
  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-PubMed CentralPubMed

Copyright

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

This article is published under license to 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.