- Open Access
Ubiquitin-specific proteases are differentially expressed throughout the Schistosoma mansoni life cycle
Parasites & Vectorsvolume 8, Article number: 349 (2015)
The ubiquitination process can be reversed by deubiquitinating enzymes (DUBs). These proteases are involved in ubiquitin processing, in the recovery of modified ubiquitin trapped in inactive forms, and in the recycling of ubiquitin monomers from polyubiquitinated chains. The diversity of DUB functions is illustrated by their number and variety of their catalytic domains with specific 3D architectures. DUBs can be divided into five subclasses: ubiquitin C-terminal hydrolases (UCHs), ubiquitin-specific proteases (USPs or UBPs), ovarian tumour proteases (OTUs), Machado-Joseph disease proteases (MJDs) and JAB1/MPN/Mov34 metalloenzymes (JAMMs).
Considering the role that the ubiquitin-proteasome system has been shown to play during the development of Schistosoma mansoni, our main goal was to identify and characterize SmUSPs. Here, we showed the identification of putative ubiquitin-specific proteases using bioinformatic approaches. We also evaluated the gene expression profile of representative USP family members using qRT-PCR.
We reported 17 USP family members in S. mansoni that present a conservation of UCH domains. Furthermore, the putative SmUSP transcripts analysed were detected in all investigated stages, showing distinct expression during S. mansoni development. The SmUSPs exhibiting high expression profiles were SmUSP7, SmUSP8, SmUSP9x and SmUSP24.
S. mansoni USPs showed changes in expression levels for different life cycle stages indicating their involvement in cellular processes required for S. mansoni development. These data will serve as a basis for future functional studies of USPs in this parasite.
The covalent modification of proteins by the addition or removal of ubiquitin, a highly conserved protein comprising of 76 amino acids, changes the molecular function of the target protein and can therefore influence its interactions with other proteins. In turn, these interactions can regulate many biological processes, including DNA repair, cell-cycle control, endocytosis, transcription and protein degradation by the proteasome [1–4].
Deubiquitinating enzymes (DUBs) catalyse the removal of ubiquitin from ubiquitin-conjugated proteins and from its precursor proteins [5, 6]. DUBs can be divided into five subclasses, of which four are cysteine proteases and one comprises a group of metalloproteases. They are named as ubiquitin C-terminal hydrolases (UCHs), ubiquitin-specific proteases (USPs or UBPs), ovarian tumour proteases (OTUs), Machado-Joseph disease proteases (MJDs) and JAB1/MPN/Mov34 metalloenzymes (JAMMs) . Of these, USPs represent the largest subclass with approximately 56 members in humans. The USP catalytic domain is highly divergent in size (295–850 residues). High sequence homology is mainly observed in two regions that surround the catalytic Cys and His residues: the so-called Cys Box domains, containing 19 amino acids, and the His Box domains, containing 60–90 amino acids [8–10].
As is the case for most cellular enzymes, the activity of DUBs can be controlled through multiple mechanisms. Several DUBs require assembly into large multimolecular complexes for full activation; this is exemplified by proteasomal DUBs. Other DUBs, including USP1, USP7, USP12 and USP46, are allosterically regulated by co-activator complexes or proteins, such as the WD40-repeat proteins, RAD50 and cullin-3 [7, 11]. Cross-talk between phosphorylation and ubiquitination is a significant aspect of intracellular signalling networks [12, 13]. Differential phosphorylation and dephosphorylation of some USPs can result in enhanced or reduced activity of these enzymes. Other post-translational modifications are emerging as modifiers of their activity, including ubiquitination and SUMOylation . Furthermore, DUBs are prone to reversible inactivation caused by reactive oxygen species (ROS) .
Schistosomes are parasitic worms that require several coordinated morphological and biochemical changes that guarantee adaptation to various environments such as water and the internal milieu of their vertebrate and invertebrate hosts [15, 16]. Our group was the first to observe that the ubiquitin–proteasome system plays a crucial role in regulating cercariae to schistosomula transition in Schistosoma mansoni [17, 18]. Subsequent studies also revealed both differential expression of 20S proteasome subunits and specific patterns of ubiquitinated proteins during S. mansoni egg maturation, highlighting the importance of controlled protein turnover during embryo development .
Considering the aforementioned results, the main objective of this work was to identify the putative and non-annotated genes encoding ubiquitin-specific enzymes using bioinformatic approaches. For this we took advantage of the available S. mansoni sequences as 81 % of the parasite’s genome has been assigned to specific chromosomes . Additionally, we evaluated the gene expression profile of 17 identified members of the USP family. These may be involved in important cellular processes during the life cycle of this parasite.
All experiments involving animals were authorized by the Ethics Committee for Animal Care of the Federal University of Ouro Preto (CEUA/UFOP protocol no. 2011/55). The experiments were performed in accordance with national and international regulations accepted for laboratory animal use and care. Mice (Balb/c strain, age 6 weeks, weight ~ 16–18 g) were kept under environmentally controlled conditions (temperature ~ 25 °C; humidity ~70 %) with free access to water and rodent diet.
The LE strain was maintained by routine passage through Biomphalaria glabrata snails and BALB/c mice. The infected snails were induced to shed cercariae under light exposure for 2 h, and the larvae were recovered by sedimentation on ice. Adult parasites were obtained by liver perfusion of mice after 50 days of infection. Livers of infected mice were macerated in phosphate buffer (Na2HPO4 63 mM, KH2PO4 330 mM, pH 8.2) and trypsinised, and the homogenate was then incubated for 2.5 h at 37 °C in a water bath. Eggs were recovered in saline solution after sequential sieving of the liver homogenate through 360- and 180-μm meshes. Mechanically transformed schistosomula (MTS) were prepared as described by Harrop and Wilson (1993) , using a protocol that mimics skin-transformed S. mansoni schistosomula . Briefly, cercariae were recovered and washed in RPMI 1640 medium (Invitrogen, Sao Paulo, Brazil) before vortexing at maximum speed for 90 s. The cercariae were immediately cultured in 169 medium for 3.5 h at 37 °C with 5 % CO2. Then, the recovered schistosomula were washed in RPMI 1640 until no tails were detected. For subsequent incubations, the parasites were maintained in M169 medium supplemented with 10 % FBS, penicillin/streptomycin at 100 μg/mL and 5 % Schneider's medium  at 37 °C in a 5 % CO2 incubator for 3.5, 24, 48 and 72 hours.
Computational analysis of USPs
The putative SmUSPs were identified and selected by mining S. mansoni sequences in GeneDB database (version 5.0, available at http://www.genedb.org/genedb/smansoni/) using the basic local alignment search tool (BLAST) algorithm BLASTp and Homo sapiens reference USP proteins as queries . Best Blastp hits showing cut-off values < e−12 were selected. Reference proteins from Mus musculus, Rattus norvegicus, Drosophila melanogaster and Caenorhabditis elegans were searched in the NCBI (National Center for Biotechnology Information) using BLASTp tool and non-redundant database to obtain a full set of putative orthologue USP proteins to compare with the S. mansoni putative proteins. Analyses of protein families, domains and active sites were performed using the PFAM (version 27.0, available at http://pfam.sanger.ac.uk) and Conserved Domains Database (CDD) (http://www.ncbi.nlm.nih.gov/cdd/) . The entire protein sequences were used to perform multiple sequence alignments using CLUSTALX2 with the default settings (available at http://www.clustal.org/clustal2/) [26, 27]. A phylogenetic tree was inferred using the neighbour-joining method and the Jones-Taylor-Thornton model . A bootstrap consensus tree inferred from 1,000 replicates was used to represent the evolutionary history of the taxa analysed. The molecular phylogenetic analyses were conducted using MEGA 5 software . All positions containing gaps and missing data were eliminated from the dataset. Furthermore, the Kyoto Encyclopedia of Genes and Genomes (KEGG) database (available at http://www.genome.jp/kegg/) was used to search for orthologue proteins in S. mansoni compared with Eukaryote protein-coding genes generated from the GENES database in KEGG . The protein domain logos were generated using WebLogo 2.8.2 at http://weblogo.berkeley.edu/logo.cgi .
Expression analysis of identified USPs
Total RNA from cercariae, schistosomula, adult worms and eggs was obtained using a combination of TRIzol (GIBCO, Sao Paulo, Brazil) and chloroform extraction; it was then column-purified using the SV Total RNA Isolation System (Promega, Belo Horizonte, Brazil). The preparation was treated three times with RNase-free DNase I (1 unit each treatment), as described by the manufacturer (RQ1 DNase; Promega). The RNA was quantified using a spectrophotometer, and an aliquot containing 1 μg of total RNA was reverse transcribed using an oligodT primer from the ThermoScript RT-PCR System (Invitrogen) as described by the manufacturer. The efficiency of DNAse I treatment was evaluated by PCR amplification of a cDNA reaction mix lacking the ThermoScript enzyme. S. mansoni-specific primers were designed using the GeneRunner® program. Despite the fact that USPs have well-conserved sequences, primers were designed to less conserved regions. The sequence accession numbers and their primer pairs are shown in Additional file 1: Table S1. Reverse-transcribed cDNA samples were used as templates for PCR amplification using SYBR Green Master Mix UDG-ROX® (Invitrogen) and a 7300 Real-Time PCR System (Applied Biosystems, Rio de Janeiro, Brazil). S. mansoni EIF4E was used as an endogenous control (GeneDB ID: Smp_001500) (forward primer: 5’-TGTTCCAACCACGGTCTCG-3’, reverse primer: 5’-TCGCCTTCCAATGCTTAGG-3’) . The efficiency of each pair of primers was evaluated according to the protocol developed by the Applied Biosystems application (the cDNA dilutions used were 1:4, 1:16, 1:64, 1:256 and 1:1024). The absence of non-specific products was confirmed by the presence of a single peak in the dissociation curves. For all investigated transcripts, three biological and technical replicates were performed, and their gene expression levels normalized using the EIF4E transcript as a reference according to the 2−ΔCt method  using the Applied Biosystems 7300 software.
Statistical analysis was performed using the GraphPad Prism software package, version 5.0 (Irvine, CA, USA). The normality of the data was established using one-way analysis of variance (ANOVA). Tukey post-tests were used to investigate significant differences in the expression of transcripts throughout the investigated stages. In all cases, the differences were considered significant when the p values were < 0.05.
Results and Discussion
Here, we report 17 USP family members in S. mansoni and a particular emphasis was given to their structures and conserved domains. Phylogenetic analysis was conducted to understand how closely related these proteases are when compared to their orthologs. In addition, the S. mansoni orthologue proteins were compared throughout the parasite’s life cycle and their expression profile evaluated by qRT-PCR.
Conservation of UCH domains in S. mansoni USPs
Putative members of the USP family in S. mansoni were retrieved by mining the parasite databases. Our analyses revealed that SmUSPs are conserved at the amino acid level compared to their orthologues from diverse organisms, such as D. melanogaster, C. elegans, H. sapiens, M. musculus and R. norvegicus. Furthermore, by comparing the S. mansoni predicted genes and their related ESTs (Expressed Sequence Tag) from S. japonicum in NCBI, up to 54 % sequence similarity could be found (Additional file 1: Table S2). A total of 18 USPs were identified in S. mansoni. The number of USPs annotated for H. sapiens (~50) is much higher possibly due to the differences in complexity between these organisms (Table I). In addition, the S. mansoni orthologue proteins were compared to eukaryotic protein-coding genes, generated from the GENES database in KEGG, confirming their identification. Structures of the UCH domains differed comparing fifteen subclades, and their length varied from 300 to 700 amino acids among orthologs. We report that this parasite possesses 12 out of 18 USP genes with complete UCH domains. This difference in UCH domains may indicate diversity in substrate specificity as described for other organisms , the presence of non-functional USP or incomplete sequences deposited on parasite databases. The S. mansoni USP proteins exhibiting a complete UCH domain displayed the well-conserved Cys and His boxes, which include the catalytic triad formed by Cys, His and Asp residues (Fig. 1). It is reported that USPs constitute the largest DUB sub-family with approximately 56 members in humans and 16 in yeast. Such diversity highlights their involvement in diverse cellular processes .
In addition, conserved domains present in each USP were investigated (Fig. 2). USPs contain a diverse range of ancillary domains whose roles are poorly characterized for the majority of catalogued and annotated USPs . This characteristic has been conserved for S. mansoni USPs. SmUSP2, SmUSP9x, SmUSP10, SmUSP14, SmUSP16, SmUSP22, SmUSP30, SmUSP36-42 and SmUSP46 contain only the UCH domain (PF00443). SmUSP5, SmUSP39 and SmUSP49-44 exhibit the zf-UBP domain (PF02148), which is a relatively small motif responsible for docking interactions with their target molecules . Furthermore, SmUSP5 and SmUSP24 possess the UBA (ubiquitin associated) domain (PF00627) found in diverse proteins involved in the ubiquitin/proteasome pathway, DNA excision-repair, and cell signalling via protein kinases . SmUSP7 contains three unique domains: MATH (Meprin and TRAF-Homology) (PF00917), USP7_ICPO (PF12436) and USP7_C2 (PF14533). The MATH domain, which participates in protein-protein interactions, is also found in intracellular TRAFs (tumour necrosis receptor-associated factors) and extracellular meprins . USP7_ICP0, found only in SmUSP7, is known to interact with the herpesvirus 1 trans-acting transcriptional protein ICP0/VMW110 . The third domain, USP7_C2, is found at the C-terminus of USP7, and its function is unclear. Only SmUSP8 contains the IPPT domain (PF01715) related to ATP binding and the USP8_dimer (PF08969) found at its N-terminus and likely involved in homodimer formation . SmUSP15 and SmUSP20 also contain the DUSP (domain present in ubiquitin-specific protease) (PF06337), located at the N- and C-terminal sides of their UCH domains, respectively. The function of DUSP is unknown; however, it may play a role in protein/protein interaction or substrate recognition . SmUSP48 is unique among the investigated molecules as it contains the ubiquitin domain (PF00240) at its C-terminus. Understanding the interaction of these protein domains provides insights into the diverse biological activities exhibited by SmUSP enzymes. Furthermore, the presence of the DUSP domain is not indicative of catalytic activity. A recent report showed that Usp39 lacks catalytic activity in vitro and is unable to cleave ubiquitin from ubiquitylated Aurora B in vivo .
A phylogenetic tree generated with the neighbour-joining method was used to separate the putative SmUSPs from their respective orthologues in H. sapiens, M. musculus, R. norvegicus, C. elegans and D. melanogaster (Fig. 3). Most of the S. mansoni USP proteins clustered with their orthologues into 15 subclades: USPs 2, 5, 7, 8, 10, 14, 15, 22, 30, 36–42, 39, 46, 48 and 49–44. The putative USP proteins 9x, 16 and 24 did not cluster in the same clade as their orthologues due to the incompleteness of their UCH domains. One exception was SmUSP20s, which exhibits the complete domain but did not cluster with their orthologs.
Distinct expression of SmUSPs is observed during S. mansoni development
SmUSP transcript levels were analysed by qRT-PCR at different developmental stages of S. mansoni: cercariae, MTS-3.5 h, MTS-24 h, MTS-48 h, MTS-72 h, paired adult worms and eggs. Three biological and technical replicates were performed. The putative SmUSP transcripts analysed were detected in all investigated stages (Fig. 4). A 2-fold change in the expression of USP genes, relative to EIF4E transcript levels, was considered of biological relevance [42, 43]. According to their expression profile, SmUSPs were divided into three groups, to account for their high (0.5-5), medium (0.2-1.5) and low expression (0.02-0.5) levels. The SmUSPs exhibiting high expression profiles were SmUSP7, SmUSP8, SmUSP9x and SmUSP24 (Fig. 4a). Those showing medium expression were SmUSP2, SmUSP15, SmUSP20, SmUSP30, SmUSP36-42, SmUSP39, SmUSP46 and SmUSP49-44 (Fig. 4b). The SmUSPs displaying low expression profiles were SmUSP10, SmUSP14, SmUSP16, SmUSP22 and SmUSP48 (Fig. 4c).
The expression profile of Smusp5 has been previously described . Of the 17 analysed transcripts, none of the USPs was up-regulated in the cercariae and MTS-72 h stages (USP2, USP5, USP7, USP8, USP9x, USP10, USP14, USP15, USP16, USP20, USP22, USP24, USP30, USP36-42, USP39, USP46, USP48, USP49-44). In contrast, Smusp2, Smusp39, Smusp46 and Smusp49-44 were up-regulated in early schistosomula (p < 0.05). Extending the culture period to 24 h, Smusp8 transcripts were observed at increased levels (p < 0.05). Concerning MTS-48 h, Smusp24 was up-regulated in this stage (p < 0.05). Smusp7, Smusp9x and Smusp15 transcripts presented at high levels in adult worms (p < 0.05). Smusp20 was up-regulated in three stages: MTS-3.5 h, MTS-48 h and adult worms (p < 0.05). Smusp30 was also up-regulated in three stages: MTS-24 h, MTS-48 h and adult worms (p < 0.05). Smusp36-42 was up-regulated in early schistosomula (MTS-3.5 h MTS-24 h and MTS-48 h) (p < 0.05). Furthermore, Smusp7, Smusp8, Smusp9x and Smusp24 were more highly expressed USPs in S. mansoni. Overall, upon examining the expression levels of all USP genes, it was found that the SmUSPs were more abundant in schistosomula and adult worms when compared with cercariae and eggs. This finding corroborates with previous data from our group, which demonstrated low levels of free, non-conjugated ubiquitin in the cercariae stage followed by its increase during schistosomula development up to adult worms . In parallel, western blot analysis revealed the accumulation of ubiquitinated proteins in cercariae and schistosomula cultured up to 5 days relative to other intra-mammalian stages, suggesting that SmUSPs are less active at this stage . This also coincides with the reported lower proteasomal proteolysis observed during early schistosomula compared with that found for the adult worm stage . Whether the accumulation of ubiquitinated proteins are due to low levels of USPs transcripts, decreased proteasomal activity or their combined effects remains to be elucidated. These data reinforce the hypothesis that a major function of USP members in S. mansoni is the regulation of protein stability during the cercariae to adult worm development. The high levels of ubiquitinated conjugates found in early schistosomula likely indicates ubiquitination linked to alternative protein fates, such as sub-cellular localization and lysosome mediated degradation [45–47]. The importance of the ubiquitin cycle to protein stability is now recognized in all eukaryotic cells as a key mechanism for maintenance of cell viability .
The cellular roles of DUBs are as wide as that of the ubiquitin-proteasome system itself given the involvement of the Ub system in intracellular signalling . Here, we selected some examples to illustrate their broader functional categorization. The proteasome has both ubiquitin ligases and DUBs that associate with it and several DUB-ligase pairs interact directly, including BRCC36-BRCA1, BAP1-BRCA1, USP4-Ro52, USP7-MDM2, USP8-GRAIL, USP20-pVHL, USP33-pVHL and USP44-APC . Although the E3 repertoire in S. mansoni is not known, previous analyses from our group suggest the conservation of the MDM2, GRAIL and APC in the parasite genome (data not shown), indicating a role for USP49-44 in cell cycle progression. USP7 removes ubiquitin not only from p53 itself but also from the p53 E3 ubiquitin-ligase MDM2. These combined effects determine functional p53 levels, creating an important role for USP7 in p53-dependent stress responses. Our group showed that the p53 orthologues, p63 and p73 , are up-regulated in MTS-3.5 h and MTS-5 days . We also observed a differential expression profile for Smusp7. Smusp8 and Smusp9x were highly expressed throughout the parasite cycle. USP8 and USP9 appear to be involved in the control of mammalian cell proliferation, inducing apoptosis , indicating a role of these enzymes in parasite remodelling. Smusp14 transcript level is higher in cercariae when compared with adult worms, suggesting that SmUSP14 can inhibit proteasome function noncatalytically, as previously observed for its yeast orthologue Ubp6 . These data could explain, at least in part, the low rate of proteasome-dependent proteolysis detected in cercariae when compared with adult worms . Smusp20 expression is up-regulated in intra-mammalian stages, particularly in MTS-3.5 h. Orthologues of USP20 seem to be involved in hypoxia signalling, and the possibility of this protease to regulate the schistosomula development will be investigated . Additionally, Smusp24 and Smusp30 were three times more abundant in MTS-48 h when compared with adult worms. Considering that USP24 was recently identified as a novel regulator of DDB2 (damage-specific DNA-binding protein 2) stability  and the implication of USP30 in the maintenance of mitochondrial morphology , SmUSP24 and SmUSP30 could be necessary during parasite remodelling. Furthermore, the transcription levels of Smusp36-42 are similar when comparing MTS-3.5 h, 24 and 48 h, reinforcing the hypothesis that specific SmUSPs can be important for early schistosomula development.
In conclusion, we describe the USP enzyme repertoire in S. mansoni and their regulated expression in the parasite life cycle. This differential profile of transcripts can reflect the stage-specific subset of their target ubiquitinated substrates during the parasite’s life cycle. Our results raise a number of interesting questions concerning the regulation of SmUSP activities and their role during schistosome development. USPs are temporally and spatially controlled and most often act as part of multi-protein complexes [56, 55]. Moreover, a single USP can act upon various substrates, and its activity can be regulated by several DUBs. Further experiments shall clarify the functions of putative SmUSPs in S. mansoni.
Hershko A, Ciechanover A. The ubiquitin system. Annu Rev Biochem. 1998;67:425–79.
Hicke L, Dunn R. Regulation of membrane protein transport by ubiquitin and ubiquitin-binding proteins. Annu Rev Cell DevBiol. 2003;19:141–72.
Muratani M, Tansey WP. How the ubiquitin-proteasome system controls transcription. Nat Rev Mol Cell Biol. 2003;4(3):192–201.
Sun L, Chen ZJ. The novel functions of ubiquitination in signalling. Curr Opin Cell Biol. 2004;16(2):119–26.
Fraile JM, Quesada V, Rodríguez D, Freije JM, López-Otín C. Deubiquitinases in cancer: new functions and therapeutic options. Oncogene. 2012;31:2373–88.
Hutchins AP, Liu S, Diez D, Miranda-Saavedra D. The Repertoires of Ubiquitinating and Deubiquitinating Enzymes in Eukaryotic Genomes. Mol Biol Evol. 2013;30(5):1172–87.
Clague MJ, Barsukov I, Coulson JM, Liu H, Rigden DJ, Urbé S. Deubiquitylases from genes to organism. Physiol Rev. 2013;93(3):1289–315.
Papa FR, Hochstrasser M. The yeast DOA4 gene encodes a deubiquitinating enzyme related to a product of the human tre-2 oncogene. Nature. 1993;366:313–19.
D’Andrea A, Pellman D. Deubiquitinating enzymes. A new class of biological regulators. Crit Rev Biochem Mol Biol. 1998;33:337–52.
Eletr ZM, Wilkinson KD. Regulation of proteolysis by human deubiquitinating enzymes. Biochim Biophys Acta. 2014;1843(1):114–28.
Sowa ME, Bennett EJ, Gygi SP, Harper JW. Defining the human deubiquitinating enzyme interaction landscape. Cell. 2009;138:389–403.
Hunter T. The age of crosstalk: phosphorylation, ubiquitination, and beyond. Mol Cell. 2007;28:730–38.
Kessler BM, Edelmann MJ. PTMs in Conversation: Activity and Function of Deubiquitinating Enzymes Regulated via Post-Translational Modifications. Cell Biochem Biophys. 2011;60(1–2):21–38.
Cotto-Rios XM, Bekes M, Chapman J, Ueberheide B. Huang TT. Cell Reports: Deubiquitinases as a signalling target of oxidative stress; 2012.
Stirewalt MA. Schistosoma mansoni: Cercaria toschistosomule. Adv Parasitol. 1974;12:115–82.
Fishelson Z, Amiri P, Friend DS, Marikovsky M, Petitt M, Newport G, et al. Schistosoma mansoni: Cell-specific expression and secretion of a serine protease during development of cercariae. Exp Parasitol. 1992;75:87–98.
Guerra-Sá R, Castro-Borges W, Evangelista EA, Kettelhut IC, Rodrigues V. Schistosoma mansoni: Functional proteasomes are required for development in the vertebrate host. Exp Parasitol. 2005;109:228–36.
Castro-Borges W, Cartwright J, Ashton PD, Braschi S, Guerra-Sá R, Rodrigues V, et al. The 20S proteasome of Schistosoma mansoni: A proteomic analysis. Proteomics. 2007;7:1065–75.
Mathieson W, Castro-Borges W, Wilson RA. The proteasome-ubiquitin pathway in the Schistosoma mansoni egg has development- and morphology-specific characteristics. Mol Biochem Parasitol (Print). 2011;175:118–25.
Protasio AV, Tsai IJ, Babbage A, Nichol S, Hunt M, Aslett MA, et al. A Systematically Improved High Quality Genome and Transcriptome of the Human Blood Fluke Schistosoma mansoni. PLoS Negl Trop Dis. 2012;6(1), e1455.
Harrop R, Wilson RA. Protein synthesis and release by cultured schistosomula of Schistosoma mansoni. Parasitology. 1993;107:265–74.
Protasio AV, Dunne DW, Berriman M. Comparative Study of Transcriptome Profiles of Mechanical- and Skin-Transformed Schistosoma mansoni Schistosomula. PLoS Negl Trop Dis. 2013;7(3), e2091.
Basch P, DiConza J. In vitro development of Schistosoma mansoni cercariae. J Parasitol. 1977;63:245–49.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215(3):403–10.
Sonnhammer EL, Eddy SR, Durbin R. Pfam: a comprehensive database of protein domain families based on seed alignments. Proteins. 1997;28(3):405–20.
Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG. The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 1997;25(24):4876–82.
Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, et al. Clustal W and Clustal X version 2.0. Bioinformatics. 2007;23:2947–8.
Saitou N, Nei M. The neighbour-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987;4:406–25.
Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: Molecular Evolutionary Genetics Analysis using Maximum Likelihood, Evolutionary Distance, and Maximum Parsimony Methods. Mol Biol Evol. 2011;28(10):2731–9.
Ogata H, Goto S, Fujibuchi W, Kanehisa M. Computation with the KEGG pathway database. Biosystems. 1998;47(1–2):119–28.
Crooks GE, Hon G, Chandonia JM, Brenner SE. WebLogo: a sequence logo generator. Genome Res. 2004;14(6):1188–90.
Liu S, Cai P, Hou N, Piao X, Wang H, Hung T, et al. Genome-wide identification and characterization of a panel of house-keeping genes in Schistosoma japonicum. Mol Biochem Parasitol. 2012;182(1–2):75–82.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta DeltaC(T)) Method”. Methods. 2001;25:402–8.
Nijman SM, Luna-Vargas MP, Velds A, Brummelkamp TR, Dirac AM, Sixma TK, et al. A genomic and functional inventory of deubiquitinating enzymes. Cell. 2005;123(5):773–86.
Ye Y, Scheel H, Hofmann K, Komander D. Dissection of USP catalytic domains reveals five common insertion points. Mol Biosyst. 2009;5(12):1797–808.
Gamsjaeger R, Liew CK, Loughlin FE, Crossley M, Mackay JP. Sticky fingers: zinc-fingers as protein-recognition motifs. Trends Biochem Sci. 2007;32(2):63–70.
Hofmann K, Bucher P. The UBA domain: a sequence motif present in multiple enzyme classes of the ubiquitination pathway. Trends Biochem Sci. 1996;21(5):172–3.
Sunnerhagen M, Pursglove S, Fladvad M. The new MATH: homology suggests shared binding surfaces in meprin tetramers and TRAF trimers. FEBS Lett. 2002;530(1–3):1–3.
Holowaty MN, Sheng Y, Nguyen T, Arrowsmith C, Frappier L. Protein interaction domains of the ubiquitin-specific protease, USP7/HAUSP. J Biol Chem. 2003;278(48):47753–61.
Avvakumov GV, Walker JR, Xue S, Finerty Jr PJ, Mackenzie F, Newman EM, et al. Amino-terminal dimerization, NRDP1-rhodanese interaction, and inhibited catalytic domain conformation of the ubiquitin-specific protease 8 (USP8). J Biol Chem. 2006;281(49):38061–70.
van Leuken RJ, Luna-Vargas MP, Sixma TK, Wolthuis RM, Medema RH. Usp39 is essential for mitotic spindle checkpoint integrity and controls mRNA-levels of aurora B. Cell Cycle. 2008;7(17):2710–9.
McCarthy DJ, Smyth GK. Testing significance relative to a fold-change threshold is a TREAT. Bioinformatics. 2009;25(6):765–71.
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–40.
Pereira RV, Vieira HG, Oliveira VF, Gomes MD, Passos LK, Borges WD, et al. Conservation and developmental expression of ubiquitin isopeptidases in Schistosoma mansoni. Mem Inst Oswaldo Cruz. 2013;109(1):1–8.
Komander D, Clague MJ, Urbé S. Breaking the chains: structure and function of the deubiquitinases. Nat Rev Mol Cell Biol. 2009;10(8):550–63.
Mizuno T, Hayashi H, Naoi S, Sugiyama Y. Ubiquitination is associated with lysosomal degradation of cell surface-resident ATP-binding cassette transporter A1 (ABCA1) through the endosomal sorting complex required for transport (ESCRT) pathway. Hepatology. 2011;54(2):631–43.
Zhang L, Xu M, Scotti E, Chen ZJ, Tontonoz P. Both K63 and K48 ubiquitin linkages signal lysosomal degradation of the LDL receptor. J Lipid Res. 2013;54(5):1410–20.
Su V, Lau AF. Ubiquitination, intracellular trafficking, and degradation of connexins. Arch Biochem Biophys. 2012;524(1):16–22.
Wilkinson KD. DUBs at a glance. J Cell Sci. 2009;122(14):2325–29.
Rutkowski R, Hofmann K, Gartner A. Phylogeny and function of the invertebrate p53 superfamily. Cold Spring Harb Perspect Biol. 2010;2(7):a001131.
Pereira RV, Gomes MS, Olmo RP, Souza DM, Jannotti-Passos LK, Baba EH, et al. NEDD8 conjugation in Schistosoma mansoni: genome analysis and expression profiles. Parasitol Int. 2013;62(2):199–207.
Ramakrishna S, Suresh B, Baek KH. The role of deubiquitinating enzymes in apoptosis. Cell Mol Life Sci. 2011;68:15–26.
Hanna J, Hathaway NA, Tone Y, Crosas B, Elsasser S, Kirkpatrick DS, et al. Deubiquitinating enzyme Ubp6 functions noncatalytically to delay proteasomal degradation. Cell. 2006;127(1):99–111.
Li Z, Wang D, Messing EM, Wu G. VHL protein-interacting deubiquitinating enzyme 2 deubiquitinates and stabilizes HIF-1alpha. EMBO Rep. 2005;6:373–78.
Zhang L, Lubin A, Chen H, Sun Z, Gong F. The deubiquitinating protein USP24 interacts with DDB2 and regulates DDB2 stability. Cell Cycle. 2012;11(23):4378–84.
Nakamura N, Hirose S. Regulation of mitochondrial morphology by USP30, a deubiquitinating enzyme present in the mitochondrial outer membrane. Mol Biol Cell. 2008;19(5):1903–11.
Felsenstein J. Confidence limits on phylogenies: An approach using the bootstrap. Evolution. 1985;39:783–91.
The authors thank the following transcriptome initiatives: the São Paulo Transcriptome Consortium, the Minas Gerais Genome Network and the Wellcome Trust Genome Initiative (UK). This work was supported by the following Brazilian research agencies: FAPEMIG (Fundação de Amparo à Pesquisa do Estado de Minas Gerais, CBB - APQ-02101-11), NuBio UFOP (Núcleo de Bioinformática da Universidade Federal de Ouro Preto) and CNPq (Conselho Nacional de Desenvolvimento Científico e Tecnológico).
The authors declare that they have no competing interests.
RVP performed transcriptome analyses, RNA isolation and qRT-PCR experiments. MSG performed transcriptome and bioinformatics analyses. RPO and DMS helped with qRT-PCR experiments. LKJP participated in S. mansoni life cycle maintenance and schistosomula culture. EHB, LKJP, WCB and RGS conceived the study and participated in design and coordination. RVP, MSG and RGS wrote the manuscript. All authors read and approved the manuscript.
Sequence accession numbers and primer pairs. Table S2. Putative ubiquitin-specific proteases in S. mansoni and S. japonicum.