Plant Deubiquitinases and Their Role in the Control of Gene Expression Through Modification of Histones

Selective degradation of proteins in the cell occurs through ubiquitination, which consists of post-translational deposition of ubiquitin on proteins to target them for degradation by proteases. However, ubiquitination does not only impact on protein stability, but promotes changes in their functions. Whereas the deposition of ubiquitin has been amply studied and discussed, the antagonistic activity, deubiquitination, is just emerging and the full model and players involved in this mechanism are far from being completely understood. Nevertheless, it is the dynamic balance between ubiquitination and deubiquitination that is essential for the development and homeostasis of organisms. In this review, we present a detailed analysis of the members of the deubiquitinase (DUB) superfamily in plants and its division in different clades. We describe current knowledge in the molecular and functional characterisation of DUB proteins, focusing primarily on Arabidopsis thaliana. In addition, the striking function of the duality between ubiquitination and deubiquitination in the control of gene expression through the modification of chromatin is discussed and, using the available information of the activities of the DUB superfamily in yeast and animals as scaffold, we propose possible scenarios for the role of these proteins in plants.


INTRODUCTION
Ubiquitination is a post-translational modification that consists of the covalent binding of the small polypeptide ubiquitin to a target protein, either singly or sequentially (polyubiquitination; Akutsu et al., 2016).
Ubiquitination participates in several cell processes such as vesicular trafficking, signalling, subcellular location, regulation of transcription, chromatin structure, and the degradation of a tagged protein (Mukhopadhyay and Riezman, 2007).
Degradation either occurs via the 26S proteasome, a protein complex involved in the degradation of polyubiquitin substrates, while monoubiquitination (ub) or addition of short ub-chains generally targets proteins for degradation via the lysosome (Hicke, 2001;Vierstra, 2009). Ubiquitin-mediated 26S proteasome degradation is a sequential process that starts with the activation of inactive ubiquitin carried out by the E1 (ubiquitinactivating) enzyme in an ATP-dependent manner. The active ubiquitin is transferred from the E1 to the E2 (ubiquitin-conjugating) enzyme that acts as an intermediate and, finally, the E3 (ubiquitin ligase) enzyme mediates the deposition of the active ubiquitin to the target protein, mainly on a lysine residue. Protein ubiquitination can be direct or indirect depending on the E3 that mediates this process (Haas et al., 1982;Ishikura et al., 2010;Zheng and Shabek, 2017).
On the other hand, deubiquitination is carried out by an evolutionary conserved group of proteins known as ubiquitin deconjugating enzymes. Comprising one of the biggest superfamilies, the deubiquitinase superfamily (DUB) counter the action of E3 ligases. The DUBs have three molecular roles: (i) generation of ubiquitin monomers (Chung and Baek, 1999); (ii) regeneration of ubiquitin during the decomposition of ubiquitin-protein conjugates in the 26S proteasome (Amerik and Hochstrasser, 2004); and (iii) deubiquitination of conjugates by releasing intact both the ubiquitin and the target to prevent the degradation of the pre-targeted protein (Taya et al., 1999). This superfamily has five families: ubiquitin-specific proteases (USPs), also called ubiquitin-specific-processing proteases (UBPs) in organisms such as Arabidopsis thaliana (Arabidopsis), ubiquitin carboxy-terminal (UCH) proteases, the ovarian tumour proteases (OTUs), the Machado-Joseph disease protein domain proteases or joshephin (MJD) family and the JAB1/MPN + /MOV34 (JAMMs) proteases. The first four families are cysteine proteases while the JAMM family are zinc metalloisopeptidases. The DUB family in Arabidopsis contains an estimated 64 members (Yan et al., 2000). However, many of the putative members are still uncharacterised and their molecular activities are still poorly understood.
Ubiquitination/deubiquitination is a highly dynamic process that is ultimately essential for many processes including cell homeostasis, signal transduction, transcriptional gene regulation, protein degradation and endocytosis among others (Hershko and Ciechanover, 1998;Yan et al., 2000;Frappier and Verrijzer, 2011). Here we will summarise the latest plant-based research on this topic, and focus on the role of this dynamic mechanism in the regulation of gene expression.

THE UBP FAMILY: COMPONENTS AND MOLECULAR ACTIVITIES
Arabidopsis UBP members have redundant functions, but also specific roles in plant development (Liu et al., 2008). This family of UBPs, which possess a highly similar sequence to human USPs proteins, has 27 members in Arabidopsis divided in 14 subfamilies based on specific protein domains (Yan et al., 2000;Zhou et al., 2017). All UBPs in Arabidopsis contain a UBP domain (although these vary in length depending on the protein) and one or more domains that are speculated to be involved in protein-protein interactions (Komander et al., 2009). Among these subfamilies, the cluster composed of UBP10/9/11/5 and 8 all contain DUSP domains, which in humans seem to participate either in protein-protein interactions or in substrate recognition (De Jong et al., 2006). The second, larger, subfamily comprises UBP15/16/17/18 and 19. All these proteins have a Myeloid-type zinc finger domain, Nervy and DEAF-1 (ZnF-MYND) domain, which is found in Eukaryotes such as yeast, Drosophila, mammals as well as plants and might have a role in protein-protein interaction (Gross and McGinnis, 1996). The subfamily formed by UBP12 and 13 have another evolutionary conserved domain called a meprin and TRAF homology (MATH) domain that can act as a folding unit and is sufficient for interaction with receptors (Ye et al., 1999;Sunnerhagen et al., 2002). Another characteristic domain is a zinc finger UBP type (ZnF-UBP) that is shared between the UBP1 and UBP2 subfamily and UBP14, the only member of another subfamily that additionally has an ubiquitin-associated (UBA) domain related with (poly)ubiquitin binding. The ZnF-UBP evolutionary conserved domain has only been reported in some members of the DUB superfamily and maintains functional similarities with other zinc finger domains. In addition, the ubiquitin-like (UBL) domain, which is found in proteins related to the ubiquitin pathway as well as other, non-related, proteins, is found in two subfamilies, the first one composed of UBP6 and UBP7 and the second one represented only by UBP26 (Larsen et al., 1998;Liu et al., 2008). A schematic representation of the members of this family has been previously shown by Liu et al. (2008).
UBP1 and UBP2 share 62% sequence similarity and also show high sequence conservation with putative members of this subfamily in other plant species (Figure 1). UBP2 breaks ubiquitin chains in vivo and both UBP1 and UBP2 can also break multi-ubiquitin chains in vitro. Analysis of ubp1 and ubp2 single mutants producing truncated proteins without functional catalytic domains did not show and altered phenotype when compared to wild-type plants, even under stress conditions. Nevertheless, in the presence of canavanine, an argininelike compound, the mutant plants display several phenotypic alterations such as shorter roots, chlorotic leaves and reduced growth rate, indicating that UBP1 and UBP2 are involved in the metabolism of this anti-herbivory compound (Yan et al., 2000).
UBP3 and UBP4 localise in the nucleus and are necessary for proper pollen development and transmission (Chandler et al., 1997). Both proteins share redundant function since the single mutants present wild type phenotype, but ubp3;ubp4 double mutant is not able to mature. Indeed, UBP3 and UBP4 are required for pollen germination because defective sperm production was reported in ubp3,ubp4 double mutants (Doelling et al., 2007).
UBP6 has a non-canonical Calmodulin motif-binding domain (CaMBD). As calmodulin is involved in the Ca 2+ signalling pathway as a sensor, this motif in UBP6 suggests the role of Ca 2+ in protein degradation and/or plant homeostasis. Arabidopsis UBP6 can complement the yeast orthologue ubp6 mutant background in vivo. In addition, the Arabidopsis UBP6 protein allows the recovery of canavanine resistance in the complemented ubp6 yeast, suggesting a role of the Arabidopsis UBP6 in external stimulus response (Byeong et al., 2005).
UBP12 and UBP13 are involved in plant immunity. Under infection of virulent and avirulent Pseudomonas syringae pv. tomato DC3000 (PstDC3000) UBP12 and 13 genes were found to be upregulated. The study of the single mutants indicates that UBP12 and UBP13 share a redundant function in plant defence response; the ubp12; ubp13 double mutant shows FIGURE 1 | Ubiquitin-specific-processing proteases (UBPs) are highly conserved in plants. Arabidopsis lyrata (AL), Capsella rubella (CRU), Brassica rapa (BRO), Popolus trichocarpa (PT), Physcomitrella patens (PP). The analysis was conducted in Clustal Omega (Sievers et al., 2014). The tree was edited for easier visualization in MEGA 7. Sequences obtained from PLAZA (Proost et al., 2015).
lethality (Brazma et al., 2003;Ewan et al., 2011). UBP12 and UBP13 have deubiquitinating activity in vitro, localise in the cytoplasm as well as in the nucleus and are encoded by ubiquitously expressed genes, which indicates a more general role of these proteins in the plant. This hypothesis was confirmed by analyses of a ubp12-mild;ubp13 double mutant, which displays an altered flowering time and changes in the circadian rhythms due to misregulation of CONSTANS (CO), a photoperiod pathway regulator, as well as clock genes such as LATE ELONGATED HYPOCOTYL (LHY), CIRCADIAN CLOCK ASSOCIATED 1 (CCA1), and TIMING OF CAB EXPRESSION 1 (TOC1) (Cui et al., 2013). UBP12 and UBP13 also regulate the ubiquitination of the transcription factor MYC2, a regulator of jasmonate (JA)-mediated responses, which places them as components of the JA signalling pathway (Jeong et al., 2017). Finally, UBP12 and 13 are involved in the process of gene repression that we will discuss in depth in the next chapter. Taken together, the data are consistent with UBP12 and 13 playing key roles in the regulation of plant development and plant defence.
UBP14, which has a highly similar protein sequence to UBP14 in the yeast model organism Saccharomyces cerevisiae and USP5 in humans (see Figure 2), also seems to have an important role in plant development. This ubiquitin protease has in vitro activity against multi-ubiquitin chains and in vivo activity against hexameric polyubiquitin. In Arabidopsis, UBP14 is necessary for embryo development, as ubp14 mutant plants are embryo lethal, displaying an increase in multi-ubiquitin chains as well as ubiquitinated proteins in the seeds resulting in the arrest of embryo development at the globular stage (Doelling et al., 2001). Decrease in UBP14 expression (ubp14-mild/per1) affects root development in the absence of phosphate, specifically reducing root hair elongation, suggesting that UBP14 has a role in plant adaptation to different environmental conditions (Li et al., 2010). GFP-UBP14 fusion protein is localised to the FIGURE 2 | Consensus phylogenetic tree of UBPs/USPs. The tree shows phylogenetic conservation of proteins from Arabidopsis, Drosophila, and humans. The analysis was conducted in MEGA7 (Kumar et al., 2016). The tree was constructed by the Neighbor-Joining method with bootstrap resampling (5000 replicates). The numbers in the nodes indicate bootstrap values. The distances were computed using the Poisson method. Scale bars indicate units of the number of amino acid substitutions per site.
nucleus where it is thought to regulate different molecular and developmental process such as endoreduplication and cell and organ growth. Specifically, the UBP domain of UBP14 genetically and physically interacts with UV-B-INSENSITIVE4 (UVA4) to repress anaphase-promoting complex/cyclosome (APC/C) and acts antagonistically with the APC/C activator CELL CYCLE SWITCH52 A1 (CCS52A1), which ultimately represses endoreduplication and so influences cell and organ growth. In addition, ubp14/da3 mutant presents lower amount of cyclin A2;3 (CDKB1;1) and cyclin-dependent kinase B1;1 (CYCA2;3) that together act as a endoreduplication repressor complex downstream of APC/C (Xu et al., 2016). These results indicate that UBP14 is a determinant for Arabidopsis growth and development through the regulation of the cell cycle.
GFP-UBP15 is localised to the cytosol and nucleus in transient expression experiment in onion cells and, although ubiquitous, its expression is stronger in leaves and flowers where it regulates leaf shape, cell number and flowering time. UBP15 de-ubiquitination activity in vitro has also been reported (Liu et al., 2008). Because of its role in controlling cell number, the ubp15 mutant displays several morphological defects such as shorter roots, smaller flowers and shorter fruits. These results suggests a role of UBP15 in plant size as well as leaf development. The activity of UBP15 seems to be partially redundant to UBP16 despite having higher sequence similarity to UBP17 (see Liu et al., 2008). In addition to cell number, UBP15 also affects seed size promoting cell proliferation in the integuments of ovules. UBP15 overexpression lines develop bigger seeds compared to wild type, while the ubp15 mutants form smaller seeds, again highlighting its role in plant development (Du et al., 2014).
On the other hand, UBP16 has a role in tolerance to saline abiotic stress conditions by regulating the plasma membrane Na + /H + antiport activity and promoting the accumulation of K + instead of Na + . ubp16 mutant displays hypersensitivity to salt stress and accumulate more sodium and less potassium. Besides, UBP16 stabilises SERINE HYDROXYMETHYLTRANSFERASE1 (SHM1), reported as a salt tolerance regulator, and also directly affects the Na + /H + antiport/K + activity. In addition, UBP16 regulates cell death via the SHM1-mediated reactive oxygen species inhibitor pathway (Zhou et al., 2012). Another member of this subfamily, UBP19, may have an essential role in plant development as ubp19 mutants are embryo lethal (Liu et al., 2008).
Abscisic acid signalling and salt and drought tolerance is mediated by the USP domain of UBP24, which affects stomata closure. In addition, UBP24 can form homodimers and cleave UBQ1 into ubiquitin monomers in vitro (Zhao et al., 2016).
The activities of UBP26 have been related to transcriptional control through deubiquitination of specific histone tails that will be discussed below.
Finally, UBP27 seems to be involved in mitochondria morphogenesis. UBP27 localises to the mitochondria and cleaves polyubiquitin in vitro. Although the ubp27 mutant does not show any phenotypic variation even in the mitochondria, UBP27 overexpression lines display changes in mitochondria shape (Pan et al., 2014).

UCH FAMILY: COMPONENTS AND MOLECULAR ACTIVITIES
The Arabidopsis UCH family has only three members (UCH1, 2, 3). All three share an N-terminal as well as the catalytic domain; nevertheless, only UCH1 and UCH2 show a high similarity of the C-terminal sequence (see Figure 3). In addition, the three members have a C-terminal extension, which in humans was reported to be an active-site crossover loop that is necessary for substrate size specificity, blocking the catalysis of multimeric ubiquitin (Popp et al., 2009).
UCH2 can release ubiquitin in vitro from the targets linked as peptides as well as isopeptides. UCH1/2 have a strong impact on shoot development, producing changes in the architecture of the inflorescence. The analysis of the uch1-1;uch2-1 double mutant showed an enhancement of AUX/IAA protein turn-over, affecting the shoot architecture due to changes in auxin signalling patterns. The opposite effect was observed using the UCH1 overexpression line. Both UCH1/2 are expressed ubiquitously, although UCH2 expression is higher than UCH1. UCH1 and 2 localise in the cytoplasm as well as in the nucleus . Strikingly, in contrast to the role of these proteins in yeast and humans, UCH1/2 does not seem to form part of the Arabidopsis 26S proteasome, although a truncated version of the protein was able to interact with subunits of the 26S proteasome in vivo. In addition, UCH1/2 interact with each other in vivo forming heterodimers in addition to subunits of the transcription-export complex 2 (TREX-2), which mediates the export of mature mRNA and tagged genes that are being transcribed (Tian et al., 2012). Future research will be required to elucidate the activities of the third member of the family in Arabidopsis and other plant species.

OTU FAMILY: COMPONENTS AND MOLECULAR ACTIVITIES
The OTU domain, that gives its name to this family of proteases, was first described in Drosophila through studies focused on the ovarian tumour (otu) gene (Steinhauer et al., 1989) and later on the domain was recognised in proteins of eukaryotes, viruses and bacteria (Makarova et al., 2000). In Arabidopsis 12 genes have been found to encode OTU domain proteins, which share similarity with other OTU proteases in metazoans (Figure 4). Some Arabidopsis OTU proteins have isoforms with a disrupted OTU domain due to alternative splicing. In addition to the OTU domain, with its characteristic catalytic Cys-His-Asp/Asn triad, Arabidopsis OTUs proteases can also have other domains, such as a UBA-like domain in OTUBAIN-LIKE DEUBIQUITINASE 1 (OTLD1) isoforms, a coiled-coil domain in OTU5 isoforms and a putative nucleic acid binding SEC-C motif in OTU7a. In addition the OTU domain displayed differences both with substrate and specificity. OTU1/2/3/4c/7a/9/10 cleave lysine-48 (K-48) linked ub chains in vitro, suggesting a role in the proteasome degradation pathway, and can also hydrolyse lysine-63 (K-63) linked ub chains, involved in other biological processes such as endocytosis, DNA repair and complex assembly. OTU1 also cleaves ub linear dimers, trimers and tetramers; whereas, OTU3/4c has only been shown to hydrolyze trimers and tetramers (Radjacommare et al., 2014).
The only OTU protease described at the molecular level in Arabidopsis so far is OTLD1. Using OTLD1 overexpression lines, different defects in Arabidopsis were reported such as increased number of flowers per inflorescence and enhanced plant growth due to the increase of cell size. OTLD1 represses genes involved in plant growth, cell expansion, hormone signalling and transition from cell division to cell expansion such us WUS, OSR2 or ABI5 among others through H2Bub1 deubiquitination at these loci. In addition, a decrease of H3K4me3 was also reported in these genes (Keren and Citovsky, 2016). These results suggest that OTU domain proteases may have different roles and further research will be needed in the characterisation of these proteins.

MJD FAMILY: COMPONENTS AND MOLECULAR ACTIVITIES
The Machado Joseph disease protein domain (MJD) proteases have not been described in Arabidopsis. This family has four members in humans and evolutionary analysis demonstrates the presence of at least two genes with high sequence similarity to the human MJD protease Ataxin-3 (ATXN3) in Arabidopsis, which acts as a histone-binding protein that mediates gene repression interacting with key regulators of transcription (Figure 5) (Li et al., 2002). Therefore, considering this data from humans, Arabidopsis MJD proteins may be involved in the epigenetic regulation of gene expression, an exciting hypothesis that will need further characterisation.

JAMM FAMILY: COMPONENTS AND MOLECULAR ACTIVITIES
Whereas the previous three families are cysteine proteases, the JAMMs are zinc metalloproteases in which their characteristic MPN+ domain have two zinc ions. One of the two zinc ions activates a water molecule to promote the cleavage of the isopeptide bond on a K-63 linked ub chains (Maytal-Kivity et al., 2002;Sato et al., 2008). The phylogenetic tree in Figure 6 shows the eight JAMM proteases of Arabidopsis.
ASSOCIATED MOLECULE WITH THE SH3 DOMAIN OF STAM3 (AMSH3), one of the three Arabidopsis proteins that show high similarity to human JAMM protease AMSH, cleaves K-48 and K-63 linked ub chain in vitro. The amsh3 mutant shows lethality at seedling stage and this was related to the proper vacuolar biogenesis at late embryogenesis stage. The lack of AMSH3 impairs the transport from the plasma membrane to the vacuole that is going to ultimately affect the ubiquitin-dependent endocytosis (Isono et al., 2010). AMHS3 also interacts with ENDOSOMAL COMPLEX REQUIRED FOR TRANSPORT-III (ESCRT-III) subunits VPS2.1 and VPS24.1, both involved in intracellular trafficking, through its Microtubule Interacting and Transport (MIT) domain, which is essential for AMSH3 function in vivo (Katsiarimpa et al., 2011(Katsiarimpa et al., , 2014. AMSH1 also interacts with the ESCRT-III subunit VPS2.1 which is important for the proper autophagic degradation process (Katsiarimpa et al., 2013).
BRCA2-CONTAINING COMPLEX 36 HOMOLOG A (BRCC36A) and BRCC36B are two other JAMM proteases in Arabidopsis, which display high sequence similarity to the mammalian BRCC36 that is involved in the DNA damage response pathway as subunit of the BRCA1-A complex, one of the most characterised complex involved in tumour suppression in humans (Her et al., 2016). BRCC36A/B are expressed ubiquitously and neither the single mutants nor the brcc36a/b double mutant shows any developmental problems. On the other hand, brcc36a presents defects in the intra-and interchromosomal homologous recombination, BRRC36A localises in the nucleus after genotoxic stress and is involved in the DNA crosslink repair in somatic cells showing an epistatic effect to BRCA1 in this process (Block-Schmidt et al., 2011). These results suggest that BRCC36A/B could affect different complexes, as it has been shown in other organisms, to regulate different molecular processes.
Two paralogous genes, CSN5A and CSN5B, encode CSN5 JAMM protease isoforms in Arabidopsis (Kwok et al., 1998). CSN5 is the core subunit of the COP9 signalosome (CSN), a conserved complex involved in signalling as well as development, and promotes protein degradation through removal of the Related to Ubiquitin (RUB) [also known as neddylin (NEDD8)] protein (Cope et al., 2002). csn5b single mutants do not show phenotypic defects, csn5a mutant has growth defects and cns5a/b double mutants are lethal (Dohmann et al., 2005). CSN5A expression is higher compared to CSN5B in Arabidopsis, both CSN5A and CSN5B are the core subunit of the CSN complex. Furthermore, CSN5A/B may have different functions since VTC1, an enzyme involved in the synthesis of vitamin C, interacts with CSN5B but not with CSN5A. (Gusmaroli et al., 2004;Wang et al., 2013). These results suggest that CSN5A/B have only partial redundant activities and, therefore, further investigation will be required to fully decipher their combined and exclusive roles. RPN11 was found in the purification of the proteasome complex in yeast (Glickman et al., 1998). RPN11 is a subunit of the 26S yeast proteasome lid subcomplex involved in the deubiquitination of proteasome substrates when the proteins arrive for their degradation (Verma, 2002). RPN11 activity is supported by RPN8. These two proteins form a heterodimer which alters the RPN11 MNP + domain to allow RPN11mediated cleavage of ubiquitin from different tagged substrates in a non-specific manner (Worden et al., 2014). In Arabidopsis, a RPN11 homolog was co-purified with the 26S proteasome, specifically as part of the 19S Regulatory Particle (Book et al., 2010).
Further research is necessary in the characterisation of this family and, considering that several JAMM proteases work as subunits in complexes that are not fully described in plants, a deeper knowledge of their molecular activities will be essential to elucidate the role of these complexes.

DUBs INFLUENCE ON GENE EXPRESSION REGULATION: H2A and H2B
A crosstalk between DUBs and the regulation of gene expression by epigenetic changes has been described in humans, Drosophila and yeast and, based on the evolutionary conservation of the superfamily (Figures 1, 2), future advances of this exciting topic in plants may be expected. Indeed, novel information on the role of different Arabidopsis DUBs in the epigenetic regulation of gene expression is emerging.
Epigenetics studies the mechanisms that produce stable and heritable changes in gene expression patterns without changes in the DNA sequence (Wolffe and Matzke, 1999). Different epigenetic processes have been observed acting at several levels, such as chromatin remodelling, modifications in non-coding RNAs, histones variants, DNA methylation and post-translational histone tail modifications (Allis and Jenuwein, 2016).
Regarding histone tail modifications, more than 30 different epigenetic marks have an impact on the structure of the chromatin, by increasing or decreasing its level of compaction and, subsequently, preventing or promoting gene transcription (Berger, 2007;Pfluger and Wagner, 2007). The deposition of these marks is mediated by chromatin remodelling complexes (Rando and Ahmad, 2007). One of these post-translational histone tail modifications is ubiquitination, which can be deposited on the tails of H2A and/or H2B causing opposite transcriptional results (Figure 7).
H2A monoubiquitination (H2Aub1) is related to transcriptional repression. This mark is deposited on the H2AK121 in Arabidopsis, H2AK119 in humans and H2AK118 in Drosophila but is absent in yeast. In Arabidopsis, PRC1 RING-finger orthologues were identified, specifically three Bmi1 proteins (AtBMI1A, AtBMI1B and AtBMI1C) and two Ring1 proteins (AtRING1a and AtRING1b) (Sanchez-Pulido et al., 2008). These five PRC1 RING-fingers can monoubiquitinate the H2A.1 isoform in vitro and AtBMI1A/1B are necessary for H2Aub1 in vivo due to the decrease of H2Aub1 in AtBMI1A/1B double mutant background. However, in vivo monoubiquitination of H2A by AtRING1a/1b has yet to be demonstrated (Bratzel et al., 2010).
H2A deubiquitination is usually associated with transcriptional activation (Nakagawa et al., 2008). In mammals, MYSM1, a member of the JAMM DUB family was the first enzyme reported to affect H2Aub1 levels (Figures 6, 7B). MYSM1 overexpression shows a decreased level of H2Aub1 and the mysm1 mutant in human embryonic kidney cell lines (HEK293T) displays an accumulation of H2Aub1. Changes in the levels of H2Aub1 also alter the enrichment of other epigenetic marks. For instance, H2A deubiquitination correlates with an increase of H1 phosphorylation, which is related with gene activation . The phylogenetic analyses of the plant JAMM family did not allow the identification of a MYSM1-like protein (Figure 6) and, thus, molecular data will be required to demonstrate whether a plant JAMM protein might develop a similar activity.
In HeLa cells, USP16 deubiquitinates H2A in vitro and in vivo and USP16 knock-down RNA line shows accumulation of H2Aub1, which affects the cell growth ratio and regulates the expression of a HOX gene, HOXD10 (Joo et al., 2007). The closer proteins in Arabidopsis based on high sequence similarity are UBP1 and UBP2 (Figure 2). Nevertheless, the role of these two proteins in H2A deubiquitination have not been characterised.
USP3 also affects the cell cycle progression and at a molecular level mediates the deubiquitination of H2Aub1 and γ-H2AX under DNA damage response (Nicassio et al., 2007;Sharma et al., 2014). Based on our phylogenetic analysis, we do not detect a clear high similarity candidate in Arabidopsis's UBPs.
The PR-DUB complex counteracts PRC1-mediated repression, as it removes H2Aub1 (Scheuermann et al., 2012) (Figure 7B). In contrast to PRC1 and PRC2, a functional PR-DUB has not been described in Arabidopsis or any other plant.
The UCH family has an important role in the levels of the H2Aub1 mark. The UCH protein CALYPSO is the catalytic subunit of the PR-DUB of Drosophila ( Figure 7B). PR-DUB mutants display an accumulation of H2Aub1 and misregulation of homeotic genes, indicating the essential role of this complex in the dynamics of H2Aub1 (Alonso et al., 2007;Scheuermann et al., 2010). The CALYPSO orthologue in humans is the tumour suppressor BRCA-1-associated protein 1 (BAP1) (Sahtoe et al., 2016). Strikingly, the ubiquitination and deubiquitination at H2A is a highly dynamic process in crosstalk with other PRCs as well as with other epigenetic marks. In Arabidopsis there are no clear orthologues of CALYPSO (see Figure 3) and complexes involved in H2A deubiquitination have indeed not been reported. However, the conservation of the H2Aub1 mark suggests a similar dynamic control by still unidentified DUB proteins or DUB protein complexes. Considering the UCH domain is conserved in the three Arabidopsis's UCH members, the regulatory role of BAP1/CALYPSO might have been acquired by one of these proteases in Arabidopsis.

H2B
While H2Aub1 is a repressive transcriptional mark, H2B monoubiquitination (H2Bub1) generally plays a role in transcriptional activation. While in Arabidopsis this mark is deposited on the K143 of the H2B tail, in other organisms H2Bub1 occurs on different lysine residues of the H2B. Similar to the situation in mammals and yeasts, genome-wide distribution of H2Bub1 in Arabidopsis correlates with other active epigenetic marks, such as H3K4me3 and H3K36me3 (Roudier et al., 2011).
Radiation sensitivity protein 6 (Rad6) was the first yeast protein reported to have E2 activity in the deposition of the H2Bub1 in vitro and in vivo (Robzyk et al., 2000) (Figure 7C). Rad6 co-operates with the E3 enzyme Bre1, which is essential for H2Bub1 in vivo. In the monoubiquitination of H2B several support complexes are needed. The Paf1 complex is necessary for proper H2Bub1 since defective single mutants of components of this complex showed a loss of H2Bub1 (Wood et al., 2003) ( Figure 7C). In humans these factors and their functions are conserved (Kim et al., 2005;Zhu et al., 2005). The activity of the BUR complex is also required since the defective mutant of Bur2, one of the complex components, shows a decrease in H2Bub1 (Wood et al., 2005) (Figure 7C). In Arabidopsis, both, Rad6 (UBC1, UBC2 and UBC3) and Bre1 (HUB1 and HUB2) play a repressive role in the control of flowering time through the activation of FLOWERING LOCUS C (FLC), the main repressor of flowering in Arabidopsis, by the deposition of ubiquitin from UBC1/2 to H2B guided by the E3 HUB1/2 (Xu et al., 2009). Subsequently, H2Bub1 enrichment in FLC promotes the deposition of H3K4me3 and H3K36me3 (Cao et al., 2008).
Considering the role of H2Bub1, it is obvious that removal of this mark entails a reduction of transcription. In yeast, Ubp8 and Ubp10 deubiquitinate H2Bub1. Ubp10 acts independently, while Ubp8 is a subunit of the Spt-Ada-Gcn5-acetyltransferase (SAGA) complex (Gardner et al., 2005) (Figure 7D). Orthologues of the yeast Ubp8 exist in Drosophila (Non-stop) and humans (USP22), as well as other SAGA complex subunits Weake et al., 2008). The Arabidopsis SAGA complex has a highly similar UBP protease (UBP22), but no H2B deubiquitination activity has yet been reported for either UBP22 or the SAGA complex. On the other hand, the role of SAGA in the control of gene expression through histone acetylation is conserved (Moraga and Aquea, 2015;Kim et al., 2015) and, hence, conservation of its involvement in H2B deubiquitination is a plausible hypothesis. In Drosophila, USP7 contributes to homeotic gene silencing guided by Polycomb (Pc).
In humans, USP7 and USP11 physically interact with members of PRC1 in vivo, such as Mel18, Bmi1 and Ring1. USP7 deubiquitinates H2A and H2B in vitro and changes in Bmi1 and Ring1 ubiquitin levels were reported in USP7 and USP11 overexpression lines. usp7 and usp11 mutants in human fibroblasts result in de-repression and loss of PRC1 binding to the tumour suppressor INK4α locus (Maertens et al., 2010). These results suggest that USP7 and USP11 have a double role in PRC1 functions: (i) regulating the activity of PRC1 through modification of the ubiquitin levels of the catalytic subunits: (ii) acting on PcG targets as direct partners of PRC1. In Arabidopsis, UBP12 and UBP13, which share a similar protein sequence to USP7 (see Figure 2), interact in vivo with LHP1, a PcG protein. It has been shown that UBP12 binds to PcG targeted chromatin and is needed for the deposition of H3K27me3 in PcG target genes (Derkacheva et al., 2016). Finally, Derkacheva and colleagues showed how UBP12 and UBP13 also contribute to gene silencing in heterochromatin, although its role is smaller compared to UBP26 in this process, sharing this function with the Drosophila USP7.
UBP26 and OTLD1 are the only DUBs described with H2Bub1 deubiquitination activity in Arabidopsis ( Figure 7D). toc1 was identified as a suppressor of mutations affecting REPRESSOR OF SILENCING1 (ROS1), a DNA demethylase involved in suppressing gene silencing (Sridhar et al., 2007). ubp26 shows higher levels of H2B monoubiquitination (H2Bub1) as well as decreased non-CpG DNA methylation. These results indicate that UBP26 does indeed deubiquitinate H2B ( Figure 7D) and furthermore that this post-translational modification is required for the deposition of the repressive mark H3K9me2, which in turn is needed for gene silencing through DNA methylation in heterochromatin. Mutations in UBP26 arrest embryo development, similar to some PcG members mutants, upregulating the expression of PHERES1 (PHE1) due to low enrichment of H3K27me3 at the PHE1 locus (Luo et al., 2008). UBP26 affects flowering time, mediating in the silencing of FLC expression due to H2Bub1 deubiquitination on FLC chromatin. As a consequence, ubp26 mutant displays an early flowering phenotype as well as higher global level of H2Bub1 (Schmitz et al., 2009). Methylation levels of H3K36 at FLC also decreases in ubp26; whereas, H3K27me3 levels increase. Thus, these results suggest that UBP26 might regulate FLC expression by decreasing the repressive mark H3K27me3 through H2B deubiquitination. Finally, it was shown that the PcG target gene AT1G80160 is also upregulated in ubp26 mutant (Derkacheva et al., 2016). Taken together, these data show that UBP26 plays an important role in the regulation of the expression of loci located in both heterochromatin and euchromatin.
OTLD1 was found to be interacting with the histone lysine demethylase KDM1C in planta. Indeed, OTLD1 has H2B deubiquitination activity in vitro and the KDM1C-OTLD1 complex represses gene expression due to H2Bub1 deubiquitination (Krichevsky et al., 2011).
Other interactions and regulation between DUB members and PcG components have not been reported yet and will require further investigation.

REMARKS AND FUTURE PERSPECTIVES
The characterisation of DUB proteins in different model organisms has accelerated but our understanding of their cellular functions in Arabidopsis remains limited.
In other organisms, DUBs commonly interact physically with partner proteins which determine their protein conformation and target specificity (Larsen et al., 1998;Amerik and Hochstrasser, 2004). AtUBP12 and UBP13 need to interact with the transcription factor MYC2, a regulator of the jasmonate-response pathway, to stabilise MYC2 as a substrate and ultimately affect the JA responses (Jeong et al., 2017), which suggests that protein-protein interactions will also be important in plant DUB function. Similarly, BAP-1 needs to be activated by the DEUBAD domain of ASXL1-3 in order to acquire its H2A deubiquitination activity (Sahtoe et al., 2016). In this review, we have highlighted that in yeast and animals DUBs usually function as part of complexes and these interactions are strictly necessary for correct biological function. However, in plants these complexes remain to be described. Therefore, identification of functional Arabidopsis DUB complexes and the proteins of which they are composed will be paramount to fully understanding DUB-mediated mechanisms.
To study Arabidopsis DUBs a promising strategy can be the use of knock-down and overexpression lines instead of knockout lines as it has been exemplified in this review. This strategy is especially interesting if we consider that several Arabidopsis DUBs mutants show an embryo-lethal phenotype. In addition, phylogenetic analyses demonstrate that many DUBs might share functional redundancy indicating the difficulty of studying them by forward genetic strategies (Figures 3-6).
In addition, more research on the DUB family will need to be performed on plant species other than Arabidopsis. So far, the second plant species with the most information on this family is rice, another angiosperm, where the preliminary characterisation of DUBs is emerging. OsUBP6 is involved in rice plant growth and is a member of the putative OsUBP family with another 20 members (Moon et al., 2009). WTG1, which shows high sequence similarity to the human OTUB1 and controls kernel size and shape (Huang et al., 2017). Finally, OsOTUB1, another high sequence similarity to the human OTUB1 is involved in grain yield affecting meristematic activity, seeing that OsOTUB1 mutants increased grain number and grain weight as well as reducing tiller number per plant . All these phenotypes increase plant architecture for farming purposes. Considering the results from rice and Arabidopsis and sequence conservation in other plant species (Figure 1), we can hypothesise that DUB proteins will have a key impact on the regulation of agronomic traits.
The second part of this review was focused on the role of the ubiquitin levels in the control of gene expression via histone tail modifications. As mentioned before, in Drosophila H2A deubiquitination by CALYPSO plays an unexpected repressive role in transcription. In Arabidopsis, OTLD1, involved in H2B deubiquitination, can promote gene activation as well as gene repression (Keren and Citovsky, 2017). These examples indicate that this mark can share functional similarities to other epigenetic mechanisms and the study of multi-epigenetic modifications in certain loci in a specific cell type during development will be the next challenge in the field.
Finally, we believe that there is still much to do regarding the role of DUBs in the modification of different histone variants. So far only the role of the humans USP3 in H2A and γH2AX deubiquitination, BRCC36 in γH2AX deubiquitination as well and the yeast Ubp8 in Cse4 deubiquitination were reported (Shao et al., 2009;Sharma et al., 2014;Canzonetta et al., 2015) and no information comes from plants despite the large number of H2A and H2B variants (Bönisch and Hake, 2012;Molden et al., 2015). Taking into account the essential role of histone variants in control of cell cycle, development and diseases in different organisms, including Arabidopsis (Talbert and Henikoff, 2016), deepening the study of this regulatory link between DUBs and histone variants will add another missing piece to the function of DUBs in the epigenetic regulation of gene expression.

AUTHOR CONTRIBUTIONS
EM and SF contributed equally to the design and writing of this article.

FUNDING
This work was supported by funding from the School of Natural Science at the National University of Ireland Galway (NUI Galway) and EM is funded by a scholarship from the College of Science (NUI Galway).