Mem Inst Oswaldo Cruz, Rio de Janeiro, VOLUME 115 | JANUARY 2020
Original Article

A new gene inventory of the ubiquitin and ubiquitin-like conjugation pathways in Giardia intestinalis

Isabel Cristina Castellanos1, Eliana Patricia Calvo2,+, Moisés Wasserman3

1Universidad EAN, Departamento de Ciencias Básicas, Bogotá, Colombia
2Universidad El Bosque, Laboratorio de Virología, Bogotá, Colombia
3Universidad Nacional de Colombia, Laboratorio de Investigaciones Básicas en Bioquímica, Bogotá, Colombia

DOI: 10.1590/0074-02760190242
405 views 334 downloads
ABSTRACT

BACKGROUND Ubiquitin (Ub) and Ub-like proteins (Ub-L) are critical regulators of complex cellular processes such as the cell cycle, DNA repair, transcription, chromatin remodeling, signal translation, and protein degradation. Giardia intestinalis possesses an experimentally proven Ub-conjugation system; however, a limited number of enzymes involved in this process were identified using basic local alignment search tool (BLAST). This is due to the limitations of BLAST’s ability to identify homologous functional regions when similarity between the sequences dips to < 30%. In addition Ub-Ls and their conjugating enzymes have not been fully elucidated in Giardia.

OBJETIVE To identify the enzymes involved in the Ub and Ub-Ls conjugation processes using intelligent systems based on the hidden Markov models (HMMs).

METHODS We performed an HMM search of functional Pfam domains found in the key enzymes of these pathways in Giardia’s proteome. Each open reading frame identified was analysed by sequence homology, domain architecture, and transcription levels.

FINDINGS We identified 118 genes, 106 of which corresponded to the ubiquitination process (Ub, E1, E2, E3, and DUB enzymes). The E3 ligase group was the largest group with 82 members; 71 of which harbored a characteristic RING domain. Four Ub-Ls were identified and the conjugation enzymes for NEDD8 and URM1 were described for first time. The 3D model for Ub-Ls displayed the β-grasp fold typical. Furthermore, our sequence analysis for the corresponding activating enzymes detected the essential motifs required for conjugation.

MAIN CONCLUSIONS Our findings highlight the complexity of Giardia’s Ub-conjugation system, which is drastically different from that previously reported, and provides evidence for the presence of NEDDylation and URMylation enzymes in the genome and transcriptome of G. intestinalis.

Ubiquitin (Ub) and ubiquitin-like modifiers (Ub-Ls) are small proteins that covalently attach to protein sub­strates and regulate various cellular processes such as the cell cycle, endocytosis, signaling pathways, intracellular trafficking, DNA repair and transcription, among others.(1)

Ub and Ub-Ls share two common features: a β-grasp fold composed of a five-stranded β-sheet and a C-termi­nal diglycine motif (GG) used for conjugation to target proteins. Currently, the Ub-L family includes 10 mem­bers: small ubiquitin modifier (SUMO), neural precur­sor cell expressed developmentally downregulated 8 (NEDD8) or Related to Ubiquitin 1 (RUB 1) in yeast, Ubiquitin-Related Modifier-1 (URM1), Ubiquitin-fold Modifier 1 (UFM1), autophagy-related proteins 8 and 12 (ATG8 and ATG12), interferon-stimulated gene 15 (ISG15), human leukocyte antigen (HLA)-F adja­cent transcript 10 (FAT10), fan ubiquitin-like protein 1 (FUB1), and histone mono-ubiquitination 1 (HUB1).(2)

The first step in the Ub-conjugation cascade is ac­tivation, which is mediated by the E1 protein (UBA1 in the budding yeast). Further, Ub is transferred to a Ubconjugating enzyme or E2 (UBC) through a transesteri­fication reaction. Finally, Ub ligase or E3 directly or in­directly transfers Ub to the substrate.(3)

E3 enzymes are a wide and diverse group of proteins that can be classified into three groups according to con­served structural domains and the transfer mechanism of Ub to the substrate. The family of Homologous to E6-associated protein carboxyl terminus (HECT) uses an indirect or two-step mechanism in which Ub is trans­ferred from E2 to E3 and then to the substrate. The fam­ily of Really Interesting New Gene (RING) and RING-related E3s have a domain containing short motifs of cysteines and histidines, which coordinate two zinc ions (Zf-C3HC4, Zf-UBR, Zf-B Box, PHD, and Zf-Mynd do­mains) and act as scaffolds for conjugation, thus promot­ing direct transfer of Ub from E2 to the substrate.(3)

Ubiquitination is a reversible process in which the deubiquitinating enzymes (DUBs) hydrolyze poly-Ub chains or remove Ub molecules. The human genome codifies approximately 90 DUBs that are classified into six families: ubiquitin C-terminal hydrolases (UCHs), ubiquitin-specific proteases (USPs), Machado-Joseph Disease (MJD), Permuted Papain fold Peptidases of DsRNA viruses and Eukaryotes (PPPDE), ovarian tu­mor (OTU), metalloproteases with a JAMM/MPN mo­tif, and the recently described motif interacting with Ub-containing protein (MINDY-4).(4)

Giardia intestinalis is a protozoan parasite that is considered to be an early divergent eukaryote; it lacks typical eukaryotic organelles such as mitochondria, peroxisomes, and Golgi apparatus. Giardia is an im­portant eukaryotic model because it could have only the key components of the principal regulation systems that characterize higher eukaryotes.(5) Our laboratory previously reported that a large number of ubiquitinated proteins exist during the motile, active metabolic, and replicative stage of Giardia (trophozoite); 151 proteins distributed over 14 functional categories were identified. However, in the infective stage (cyst), only 55 ubiquiti­nated substrates were observed. Despite this marked de­crease, ubiquitination of enzymes involved in cyst wall biogenesis suggested that Ub modification plays a cru­cial role in this stage of the cell cycle.(6) Therefore, Giar­dia might be a suitable biological model to define the fundamental elements of the Ub-conjugation pathway. The proteasome components have recently been anal­ysed using bioinformatics, confirming findings reported earlier where a remarkable conservation was observed.(7)

Previous studies have identified three genes for Ub, one E1 enzyme, 11 E2 enzymes, four E3 ligases, and 9 DUBs;(8, 9, 10, 11, 12, 13) however, the most divergent genes may have been overlooked. Herein, we performed an exhaustive search using an intelligent systems approach based on hidden Markov models (HMMs) with profiles from the Pfam and Superfamily databases. Approxi­mately 120 genes were identified, 88 of which corre­spond to new findings; among these genes, 76 were E3 ligases. Furthermore, we identified NEDD8 and URM1 conjugation pathways.

MATERIALS AND METHODS

Inventory building - The full proteome database from G. intestinalis was downloaded from Eupath da­tabase version 5.0 (available at http://giardiadb.org). In addition, 66 Pfam HMM profiles associated with Ub and Ub-Ls conjugation systems were selected and down­loaded from the Pfam database version 31.0 (http://pfam. xfam.org/) (Table I). The HMMER package version 3.1 (http://hmmer.org) was used to search each Pfam profile against the entire proteome dataset using the hmmsearch tool and a threshold E-value ≤ 0.1. The repetitive tasks were automated using a perl script.

Each result was analysed for the respective domain, and other structural features were verified by basic local alignment search tool (BLAST) searches in the EMBL Pfam database (http://pfam.xfam.org/) and SMART analysis program (http://smart.embl-heidelberg.de/ smart/set_mode.cgi). Finally, each sequence identified in Giardia proteome was used as a query on BLASTp tool from www.giardiadb.org using the UniProtKB/ Swiss-Prot database to identify orthologs.

 

table01

 

Expression analysis - To analyse gene expression profiles, RNA-seq and microarray datasets were em­ployed. The data were downloaded from NCBI’s gene expression omnibus (GEO), with accession numbers GSE36490 and GSE25460 respectively, and parsed us­ing in-house perl scripts.

Protein model building and evaluation - 3D models for Ub and Ub-Ls were obtained using Phyre2 (www. sbg.bio.ic.ac.uk/phyre2/). 2LRW, 2QJL, 1YX5, and 1A5R PDB structures were used to model Ub, URM1, NEDD8, and SUMO, respectively. The predicted models were subjected to energy minimisation using YASARA (http://www.yasara.org/), and the stereochemical stabil­ity was verified using PROCHECK and ProSA analysis (http://www.ebi.ac.uk/thornton-srv/databases/pdbsum; https://prosa.services.came.sbg.ac.at/prosa.php). Ram­achandran plots were computed using Rampage to deter­mine the stereochemical quality and predicted accuracy of the structures.

Phylogenetic analysis - Protein sequences were re­trieved from Uniprot and aligned using CLUSTAL Omega (https://www.ebi.ac.uk/Tools/msa/clustalo/). A neighbor-joining phylogenetic tree was constructed us­ing the MEGA 7 program. Bootstrap values were ob­tained from 1000 replicates.

RESULTS

Our methodology for searching proteins that com­pose Ub and Ub-L conjugation systems in G. intesti­nalis identified 118 sequences that were classified into five groups: Ub and Ub-like, E1 and E1-like, E2, E3, and DUB enzymes (Table II).

The first group contains seven sequences: one for free Ub, two for fused Ub (Ub-L40 and Ub-S27), and four for Ub-like proteins: SUMO, URM1, RUB1, which is an ortholog of mammalian NEDD8, and UFM1. Among proteins identified, RUB1 is the closest to Ub (41% identical), whereas UFM1 is the least similar, with 16% identity. To characterise these sequences structur­ally, three-dimensional structure predictions for each protein were performed. The predicted structures for Ubiquitin, SUMO, NEDD8, and URM1 were similar to Ub-Ls as they possessed a β grasp fold; the characteris­tic diglycine motif at the C-terminus and the hydropho­bic core (Ile-66, Leu-65, and Leu-74) which is conserved in URM1 orthologs, was also identified.(2) Ramachan­dran plots were then built to assess the quality of the structures. The plots indicated that the distribution of residues in the allowed and disallowed regions as well as the plot analysis predicting the stability of the models are within reliable ranges (Fig. 1).

The second group includes proteins identified with the ThiF domain, which is characteristic of Ub activat­ing enzymes and members of the bacterial ThiF/MoeB/ HesA family. E1 enzymes for Ub (UBA1), SUMO (UBA2), NEDD8 (UBA3), URM1 (UBA4), and one ortholog of MoeB were identified (Table II). These en­zymes harbor two catalytic activities required for acti­vation; adenylation and thioester bond formation. The presence of the nucleotide binding motif, GXGXXGCE, and the catalytic cysteine motif, PZCTXXXXP, which are conserved among canonical E1s (UBA1, 2, and 3) in sequences GL50803_4083, GL50803_10661, and GL50803_6288 confirming our findings (Fig. 2). The sequences GL50803_12853 and GL50803_11436 exhib­ited some degree of similarity with the noncanonical E1 U BA4 a nd t he p rokaryotic p rotein M oeB, r espec­tively. UBA4 (MOCS3 in human) has dual functions in both protein urmylation and in sulfur transport within the tRNA thiolation pathways. MoeB activates MoaD through its C-terminal end during the first step to incor­porate sulfur during molybdenum cofactor biosynthesis. (14) The comparison of Giardia sequences with their or­thologs demonstrated that UBA4 is more closely related to MoeB than to canonical E1s (Fig. 2).

E2 enzymes are characterised by a highly conserved domain of approximately 150 amino acids or UBC, which contains a conserved catalytic cysteine and interacts with E1. E2 are classified into four groups based on the exis­tence of additional extensions flanking the UBC domain that confer functional differences. Class I has only the UBC domain, classes II and III have one extension (N- or C-terminal), and class IV has both. In Giardia, we found 11 class I conjugating enzymes and one class III enzyme. GL50803_5921 has a C-terminal extension of 170 amino acids and does not present any similarities within the pro­tein database. One of the 12 genes classified in this group is a new finding; GL50803_8638 shares 43% identity with UBC6 from Schizosaccharomyces pombe.

The fourth group includes 82 sequences corre­sponding to E3 ligases. Five are HECT enzymes: GL50803_137754 and GL50803_17386 showed 53% and 54% similarity to human E6-AP, which is the found­ing member of the HECT family.(15) GL50803_16321 is similar to Pub 1 from S. pombe, which ubiquitinates cdc25 (the mitotic phosphatase) in vivo. Furthermore, GL50803_32730 and GL50803_3117 are 31% identical to E3 ligases from Caenorhabditis elegans and Eimeria maxima, respectively (Table II).

The remaining sequences were classified as RING type, of which 71 displayed a typical C3HC4 domain. Within this class, we could distinguish three subgroups of metazoan E3 ligases. The first subgroup included 30 sequences containing ankyrin repeats and single C-ter­minal RING domain; this arrangement of domains has been reported in the XB3 family in plants.(16) Although 14 of these proteins show some degree of similarity with Mind bomb (MIB) proteins, they do not conserve the typ­ical modular architecture, including two substrate rec­ognition modules at the N-terminal, a series of ankyrin repeats, and multiple RING domains at the C-terminal. (17) Phylogenetic analysis was performed to determine whether Giardia ankyrin-RING proteins are members of the XB3 family. Our results revealed that the ankyrin- RING proteins were distributed into four distinct clus­ters, one of which appears to be related to XB3 (Fig. 3).

The second subgroup of RING ligases contains pro­teins with some degree of similarity to TRIM proteins, which are characterised by a tripartite motif composed of one RING domain, one or two B-Box domains, and a coiled-coil domain. TRIM proteins play important roles in various processes including cell growth, DNA damage signaling, senescence, tumor suppression, and innate anti­viral response;(15) the sequences identified share similarity with RING domains exclusively (Table II, Fig. 4).

In the last subgroup, there are several orthologs of ring finger and CHY zinc finger domain-containing pro­tein 1 (RCHY1) also known as p53-induced protein with a RING-H2 domain (Pirh2); these proteins regulate cell-cy­cle progression, cell proliferation, and cell death through the ubiquitination and degradation of diverse substrates such as p53, p27Kip1, p63, p73, c-Myc, and Chk2.(15)

Furthermore, there were E3 ligases involved in nu­clear functions, orthologs for Bre1, Pep5, and Hel 2 re­quired for the degradation of histones in yeast. Orthologs for RAD18, Pirh2, and the endonuclease SLX1 involved in DNA damage repair. RNF12, E3 ligase for c-Myc, and orthologs to E3s, which engage in apoptosis and cellular signaling (e.g., IAP and XIAP), were also identified.

 

table01

 

In contrast to the high number of sequences with a RING domain, only six harbor RING-related domains; one protein containing the Zf-UBR domain, which is involved in ubiquitination/degradation through the N-terminal rule; two proteins containing a B-Box domain; and two proteins with a PHD domain (C4HC3).

The fifth group corresponds to deubiquitinat­ing enzymes; three new findings were reported: One deNEDDylase (GL50803_10218); an OTU enzyme (GL50803_88556), and a member of the MINDY-4 fam­ily (GL50803_7349). Other OTU enzymes have been re­ported in protozoans such as Plasmodium falciparum, Cryptosporidium parvum, Toxoplasma gondii and Ei­meria acervulina,(18, 19) while orthologs of MINDY have not yet been reported in parasites.

Finally, to verify the expression of the genes identi­fied during the trophozoite stage, transcriptomics data from Franzen et al.(20) were used. These data were ob­tained from Illumina transcript analysis of WB strain assemblage A; genes with transcript levels formulated as fragments per kilobase per million fragments mapped (FPKM) < 0.5 were regarded as not expressed. From our analysis, transcription was detected in all of the genes but two (GL50803_16687 and GL50803_4430). FPKM ranged from 0.7 to 26.493, and a wide variation of ex­pression levels within the five groups of the ubiquitina­tion pathway genes was observed, with ubiquitin, sumo, nedd8, and urm1 with the highest levels of transcripts (Fig. 5, Supplementary data). Although we analysed the data available from encystation,(21) none of the genes were overexpressed during this differentiation process (data not shown).

DISCUSSION

Approximately 30 genes associated with the Ub-conjugation pathway have been identified in G. intes­tinalis, 12 of which were previously identified at our laboratory using BLAST.(10) This number is consider­ably lower than the hundreds of genes associated with this system reported in other eukaryotes. In humans, for example, there are two Ub activating enzymes (UBA1 and UBA6), approximately 40 conjugating proteins, and hundreds of E3 ligases.(1, 3) In this study, after an exhaus­tive search using the HMM, > 100 genes were identified.

In most organisms, Ub is codified in three forms: as a monomer, as polyUb, or as fusion proteins with ri­bosomal proteins. The Ub sequence is extensively con­served among all eukaryotes with similarities > 98% among humans, yeast, and apicomplexan parasites.(18) However, Giardia’s Ub gene is one of the most divergent sequences reported; it shares only 70% similarity with the human sequence.(8, 10)

Currently, the Ub-L family includes 10 members; some of them, such as FUB1, ISG15, and FAT10, are exclusive to metazoan and are involved in the immune response, T cell activation, and antiviral and antimicro­bial defense, respectively.(2) Here, we identified SUMO, NEDD8, URM1, and UFM1 homologs; however, we did not identify ATG8 and HUB1. Orthologs of HUB1 have been identified from yeasts to humans, in apicomplexan parasites, and in Entamoeba, (deep-branching eukary­ote, such as Giardia).(11, 18, 22) In S. cerevisiae, S. pombe, and humans, HUB1 is mainly involved in alternative splicing of pre-mRNA.(23) Notably, in Giardia, few in­trons have been reported, and splicing is not essential for the parasite’s viability.(24)

ATG8 plays a central role in the autophagy network, and its conjugation requires the activities of the ATG7 (activating) and ATG3 (conjugating) enzymes. In related organisms such as Entamoeba, Cryptosporidium, and Trichomonas vaginalis and in other protozoa parasites, the ATG8 system has been identified, and autophagy plays a principal role in the parasites’ survival.(11, 18, 22) In this study, ATG8, ATG7, and ATG3 enzymes were not identified, as previously reported by Bagchi et al.,(25) who suggested that in Giardia, autophagy does not oper­ate by the classical mechanism based on these proteins.

UFM1 is a highly conserved protein in metazoa and plants but not in yeasts. Initially, this modification was considered to be metazoa-specific and was associated with the endoplasmic reticulum stress response; how­ever, recently, Gannavaram et al.(26) demonstrated the existence of UFM1 and its conjugation enzymes: UBA5 (E1), UFC1 (E2), and UFL1 (E3) in Trypanosomatidae parasite proteomes. In Leishmania donovani, is associat­ed with the mitochondria, and plays an important role in pathogenesis.(26) UFM1 homologs have not been identi­fied in other unicellular parasites; although we identified a UFM1 domain-containing protein (GL50803_104982) with a high E-value (2.27e-30), we did not find the con­jugation enzymes and the predicted structure was not typical for Ub-L (Fig. 1).

URM1 acts as a protein modifier (urmylation) and be­longs to the superfamily MoaD and ThiS in prokaryotes; these are small sulfur carrier proteins involved in molyb­denum cofactor (MoaD) and thiamin (ThiS) biosynthe­sis, respectively. The URM1 sequence and structure are similar to those of bacterial proteins than those of Ub-Ls, and its conjugation process depends entirely on UBA4.(14) We identified URM1, UBA4, and bacterial MoeB homo­logs; however, we did not identify MoaD or ThiS protein (Table II). The sequences for UBA4 and MoeB share ap­proximately 24% sequence identity, and they are differ­entially expressed during encystation and temperature or redox stress response (data available at www.giardiadb. org), suggesting that they are involved in two different processes. However, whether URM1 can be activated by two enzymes, or if there is another member of the MoaD/ ThiS family to be identified, is yet unclear.

Among Ub-Ls, NEDD8 has the highest identity with Ub (approximately 60%) and is highly conserved from yeast to humans. NEDDylation is catalysed by the specif­ic enzymes: UBA3, UBE2F, UBC12, RING-box proteins (Rbx1 and Rbx2) and is reversed by specific proteases, such as DEN1/SENP8 and UCH L3.(2) This conjugation cascade has been identified in other protozoans, such as Plasmodium spp, T. gondii, C. parvum, Entamoeba spp, and Trypanosoma brucei.(18, 27) Although cullins are the most abundant substrates, proteins such as DNA damage binding protein 1 (DDB1), translation elongation factor α1, the chaperone DnaJ, in addition to the NEDDylation enzymes have been reported as targets.(27) In this study we found the NEDD8 ortholog and enzymes involved in the conjugation and de-NEDDylation processes (Table II).

The covalent modification of proteins by SUMO is the unique ubiquitination-like process that was described previously in Giardia; a single gene for SUMO, SUMO activating Enzyme subunit 2 (SAE2), SUMO conjugating enzyme (Ubc9), and one deSUMOylase were reported by Vranych et al.(13) The same sequences were identified here; however, this finding differs from that of higher eu­karyotes, in which multiple members of the SUMO fam­ily exist and SAE2 is a heterodimer of UBA1 and AOS1 subunits.(2) Although only two SUMOylated proteins, arginine deiminase and α-tubulin have been fully identi­fied in Giardia, SUMOylation of other proteins and the participation of SUMOylation in encystation, cell-cycle progression, cell growth, and morphology maintenance were recently demonstrated.(28)

Regarding the ubiquitination process, the extensive conservation exhibited by the E2 enzymes enabled the identification of most of Giardia’s enzymes using basic local alignment tools in the past;(10) however, a novel non­canonical ubiquitin-conjugating enzyme (NCUBE) was identified here (GL50803_8638); NCUBE enzymes are localised in the lumen of endoplasmic reticulum (ER) and participate in ER-associated degradation. Neverthe­less, Giardia’s protein and one ortholog of Entamoeba lack the hydrophobic C-terminal tail required for ER localisation; the function of this “truncated” Ubc6-like protein has not yet been established.(29)

The E3 Ub ligase family is the largest family of pro­teins involved in ubiquitination because it is required for specific substrate recognition. In Giardia, we identified approximately 80 enzymes, with five containing the HECT domain; our results agree with those reported for Entamoeba, apicomplexan, and yeast, which are organ­isms that codify for five or six HECT ligases.(22)

Considering the RING E3 ligases, our results largely differ from those previously reported by Gallego et al.;(10) they identified only one putative protein with the RING domain, whereas we identified approximately 70. Similar to other eukaryotes, in Giardia, ubiquitination is involved in numerous cellular processes such as protein quality control, metabolic pathways, endocytosis, cell signaling, DNA/RNA metabolism, and differentiation.(6) The large and diverse group of RING ligases reported here corre­sponds to multiple processes regulated by ubiquitination.

The over-representation of members of the XB3 family (Ankyrin repeat C3HC4 RING finger) (Figs 3, 4) suggests an important role in parasite growth; this type of protein has not been reported in other parasites or in yeast. Our data provide the first reference for the XB3 family in a single-cell eukaryote. In plants, these proteins play important roles in development, stress re­sponses, cell death induction, and pathogen response.(16)

Another interesting result is the absence of multisub­unit E3 enzymes, which are important cell-cycle regu­lators in higher eukaryotes. The anaphase promoter complex (APC) and SCF are E3 ligases that ubiquitinate cyclins involved in the M and G1/S phases, respectively. The APC/C holoenzyme comprises at least 14 different proteins distributed into three subcomplexes: a scaffolding subcomplex, a catalytic subcomplex (containing APC11 RING ligase, cullin-like subunit APC2, and APC10), and the substrate recognition subcomplex.(30) Our search iden­tified only one small protein (GL50803_8432; 78 amino acids) that is 44% similar to the APC11 from Mus muscu­lus (84 amino acids). The SCF complex consists of three subunits: Skp1/Cullin1 (scaffolding protein), Rbx1 (RING Ligase), and an interchangeable F-box-protein that deter­mines substrate specificity. This E3 complex ubiquiti­nates the S phase CDK inhibitor (Sic1p/p27kip1) involved in G1/S phase progression.(30) Our search failed to identify any sequence for Skp1 and Cullin1 homologs, although two Rbx proteins were identified (GL50803_8432 and GL50803_8241). These results agree with those reported by Eme et al.(31) and Gourguechon et al.;(32) none of them found any of the highly conserved components of these Ub ligases in Giardia. In addition, Gourguechon’s study showed that APC substrates as mitotic cyclin B, aurora and polo-like kinases were not ubiquitinated and that pro­teasome inhibition did not cause cell-cycle arrest, sug­gesting that no ubiquitination and protein degradation via proteasome are involved in Giardia cycle progression.

In 1994, Krebber suggested that the ubiquitination system in Giardia was a basal acquired system;(8) none­theless, the set of genes reported here demonstrates that the Giardia genome codifies members of the Ub and Ub-L conjugation system, similar to that described for higher eukaryotes.

AUTHORS’ CONTRIBUTION

ICC participated in the study design, accomplished gene inventory and 3D structure modeling; EPC participated in the study design, gene inventory validation, analysed the data, and wrote the manuscript; MW participated in the study design and reviewed the manuscript. All authors reviewed the results and approved the final version of the manuscript. The authors declare that there is no conflict of interest regarding the pub­lication of this paper.

REFERENCES
01. Ciechanover A. The unravelling of the ubiquitin system. Nat Rev Mol Cell Biol. 2015; 16(5): 322-24.
02. Cappadocia L, Lima CD. Ubiquitin-like protein conjugation: structures, chemistry, and mechanism. Chem Rev. 2018; 118(3): 889-918.
03. Metzger MB, Hristova VA, Weissman AM. HECT and RING finger families of E3 ubiquitin ligases at a glance. J Cell Sci. 2012; 125(3): 531-7.
04. Clague MJ, Urbé S, Komander D. Breaking the chains: deubiquitylating enzyme specificity begets function. Nat Rev Mol Cell Biol. 2019; 20: 338-52.
05. Carranza PG, Lujan HD. New insights regarding the biology of Giardia lamblia. Microbes Infect. 2010; 12(1): 71-80.
06. Niño CA, Chaparro J, Soffientini P, Polo S, Wasserman M. Ubiquitination dynamics in the early-branching eukaryote Giardia intestinalis. Microbiol Open. 2013; 2(3): 525-39.
07. Ray A, Sarkar S. The proteasome of the differently-diverged eukaryote Giardia lamblia and its role in remodeling of the microtubule-based cytoskeleton. Crit Rev Microbiol. 2017; 43(4): 481-92.
08. Krebber H, Wöstmann C, Bakker-Grunwald T. Evidence for the existence of a single ubiquitin gene in Giardia lamblia. FEBS Lett. 1994; 343(3): 234-6.
09. Catic A, Sun Z, Ratner D, Misaghi S, Spooner E, Samuelson J, et al. Sequence and structure evolved separately in a ribosomal ubiquitin variant. EMBO J. 2007; 26(14): 3474-83.
10. Gallego E, Alvarado M, Wasserman M. Identification and expression of the protein ubiquitination system in Giardia intestinalis. Parasitol Res. 2007; 101(1): 1-7.
11. Ponder E, Bogyo M. Ubiquitin-like modifiers and their deconjugating enzymes in medically important parasitic protozoa. Eukaryot Cell. 2007; 6(11): 1943-52.
12. Iyer LM, Burroughs AM, Aravind L. The prokaryotic antecedents of the ubiquitin-signaling system and the early evolution of ubiquitin-like β-grasp domains. Genome Biol. 2006; 7(7): R60.
13. Vranych CV, Merino MC, Zamponi N, Touz MC, Rópolo AS. SUMOylation in Giardia lamblia: a conserved post-translational modification in one of the earliest divergent eukaryotes. Biomolecules. 2012; 2(3): 312-30.
14. Jüdes A, Bruch A, Klassen R, Helm M, Schaffrath R. Sulfur transfer and activation by ubiquitin-like modifier system Uba4•URM1 link protein urmylation and tRNA thiolation in yeast. Microb Cell. 2016; 3(11): 554-64.
15. Uchida C, Kitagawa M. RING-, HECT-, and RBR-type E3 ubiquitin ligases: involvement in human cancer. Current Cancer Drug Targets. 2016; 16(2): 157-74.
16. Yuan X, Zhang S, Liu S, Yu M, Su H, Shu H, et al. Global analysis of ankyrin repeat domain C3HC4-type RING finger gene family in plants. PLoS One. 2013; 8(3): e58003.
17. Guo B, McMillan BJ, Blacklow SC. Structure and function of the Mind bomb E3 ligase in the context of Notch signal transduction. Curr Opin Struct Biol. 2016; 41: 38-45.
18. Ponts N, Yang J, Chung DW, Prudhomme J, Girke T, Horrocks P, et al. Deciphering the ubiquitin-mediated pathway in apicomplexan parasites: a potential strategy to interfere with parasite virulence. PLoS One. 2008; 3(6): e2386.
19. Wang P, Gong P, Wang W, Li J, Ai Y, Zhang X. An Eimeria acervulina OTU protease exhibits linkage-specific deubiquitinase activity. Parasitol Res. 2019; 118(1): 47-55.
20. Franzen O, Jerlström-Hultqvist J, Einarsson E, Ankarklev J, Ferella M, Andersson B, et al. Transcriptome profiling of Giardia intestinalis using strand-specific RNASeq. PLoS Comput Biol. 2013; 93: e1003000.
21. Morf L, Spycher C, Rehrauer H, Fournier CA, Morrison HG, Hehl AB. The transcriptional response to encystation stimuli in Giardia lamblia is restricted to a small set of genes. Eukaryot Cell. 2010; 9(10): 1566-76.
22. Arya S, Sharma G, Gupta P, Tiwari S. In silico analysis of ubiquitin / ubiquitin-like modifiers and their conjugating enzymes in Entamoeba species. Parasitol Res. 2012; 111(1): 37-51.
23. Mishra SK, Ammon T, Popowicz GM, Krajewski M, Nagel RJ, Ares M, et al. Role of the ubiquitin-like protein Hub1 in splice-site usage and alternative splicing. Nature. 2011; 474(7350): 173-8.
24. Gómez V, Wasserman M. Interactions between Giardia duodenalis Sm proteins and their association with spliceosomal snRNAs. Parasitol Res. 2017; 116(2): 617-26.
25. Bagchi S, Oniku AE, Topping K, Mamhoud ZN, Paget TA. Programmed cell death in Giardia. Parasitology. 2012; 139(7): 894-903.
26. Gannavaram S, Connelly PS, Daniels MP, Duncan R, Salotra P, Nakhasi HL. Deletion of mitochondrial associated ubiquitin fold modifier protein UFM1 in Leishmania donovani results in loss of β-oxidation of fatty acids and blocks cell division in the amastigote stage. Mol Microbiol. 2012; 86(1): 187-98.
27. Liao S, Hu H, Wang T, Tu X, Li Z. The protein neddylation pathway in Trypanosoma brucei functional characterization and substrate identification. J Biol Chem. 2017; 292(3): 1081-91.
28. Di Genova BM, da Silva RC, da Cunha JPC, Gargantini PR, Mortara RA, Tonelli RR. Protein SUMOylation is involved in cell-cycle progression and cell morphology in Giardia lamblia. J Eukaryot Microbiol. 2017; 64(4): 491-503.
29. Kumari R, Gupta P, Tiwari S. Ubc7/Ube2g2 ortholog in Entamoeba histolytica: connection with the plasma membrane and phagocytosis. Parasitol Res. 2018; 117(5): 1599-611.
30. Nakayama KI, Nakayama K. Ubiquitin ligases: cell-cycle control and cancer. Nat Rev Cancer. 2006; 6(5): 369-81.
31. Eme L, Trilles A, Moreira D, Brochier-Armanet C. The phylogenomic analysis of the anaphase promoting complex and its targets points to complex and modern-like control of the cell cycle in the last common ancestor of eukaryotes. BMC Evol Biol. 2011; 11(1): 265.
32. Gourguechon S, Holt LJ, Cande WZ. The Giardia cell cycle progresses independently of the anaphase-promoting complex. J Cell Sci. 2013; 126(10): 2246-55.

+ Corresponding author: calvoeliana@unbosque.edu.co
https://orcid.org/0000-0002-9135-0748
Received 10 July 2019
Accepted 02 January 2020

Our Location

Memórias do Instituto Oswaldo Cruz

Av. Brasil 4365, Castelo Mourisco 
sala 201, Manguinhos, 21040-900 
Rio de Janeiro, RJ, Brazil

Tel.: +55-21-2562-1222

This email address is being protected from spambots. You need JavaScript enabled to view it.

Support Program

ioc

fiocruz governo
faperj cnpq capes