Ancient origin of animal U-box ubiquitin ligases

Background

The patterns of emergence and diversification of the families of ubiquitin ligases provide insights about the evolution of the eukaryotic ubiquitination system. U-box ubiquitin ligases (UULs) are proteins characterized by containing a peculiar protein domain known as U box. In this study, the origin of the animal UUL genes is described.

Results

Phylogenetic and structural data indicate that six of the seven main UUL-encoding genes found in humans (UBE4A, UBE4B, UIP5, PRP19, CHIP and CYC4) were already present in the ancestor of all current metazoans and the seventh (WDSUB1) is found in placozoans, cnidarians and bilaterians. The fact that only 4 - 5 genes orthologous to the human ones are present in the choanoflagellate Monosiga brevicollis suggests that several animal-specific cooptions of the U box to generate new genes occurred. Significantly, Monosiga contains five additional UUL genes that are not present in animals. One of them is also present in distantly-related protozoans. Along animal evolution, losses of UUL-encoding genes are rare, except in nematodes, which lack three of them. These general patterns are highly congruent with those found for other two families (RBR, HECT) of ubiquitin ligases.

Conclusions

Finding that the patterns of emergence, diversification and loss of three unrelated families of ubiquitin ligases (RBR, HECT and U-box) are parallel indicates that there are underlying, linage-specific evolutionary forces shaping the complexity of the animal ubiquitin system.

     

Redefining the catalytic HECT domain boundaries for the HECT E3 ubiquitin ligase family

There are 28 unique human members of the homologous to E6AP C-terminus (HECT) E3 ubiquitin ligase family. Each member of the HECT E3 ubiquitin ligases contains a conserved bilobal HECT domain of approximately 350 residues found near their C-termini that is responsible for their respective ubiquitylation activities. Recent studies have begun to elucidate specific roles that each HECT E3 ubiquitin ligase has in various cancers, age-induced neurodegeneration, and neurological disorders. New structural models have been recently released for some of the HECT E3 ubiquitin ligases, but many HECT domain structures have yet to be examined due to chronic insolubility and/or protein folding issues. Building on these recently published structural studies coupled with our in-house experiments discussed in the present study, we suggest that the addition of ∼50 conserved residues preceding the N-terminal to the current UniProt defined boundaries of the HECT domain are required for isolating soluble, stable, and active HECT domains. We show using in silico bioinformatic analyses coupled with secondary structural prediction software that this predicted N-terminal α-helix found in all 28 human HECT E3 ubiquitin ligases forms an obligate amphipathic α-helix that binds to a hydrophobic pocket found within the HECT N-terminal lobe. The present study brings forth the proposal to redefine the residue boundaries of the HECT domain to include this N-terminal extension that will likely be critical for future biochemical, structural, and therapeutic studies on the HECT E3 ubiquitin ligase family.     

Transgenic Arabidopsis expressing osmolyte glycine betaine synthesizing enzymes from halophilic methanogen promote tolerance to drought and salt stress

Gene duplication events exert key functions on gene innovations during the evolution of the eukaryotic genomes. A large portion of the total gene content in plants arose from tandem duplications events, which often result in paralog genes with high sequence identity. Ubiquitin ligases or E3 enzymes are components of the ubiquitin proteasome system that function during the transfer of the ubiquitin molecule to the substrate. In plants, several E3s have expanded in their genomes as multigene families. To gain insight into the consequences of gene duplications on the expansion and diversification of E3s, we examined the evolutionary basis of a cluster of six genes, duplC-ATLs, which arose from segmental and tandem duplication events in Brassicaceae. The assessment of the expression suggested two patterns that are supported by lineage. While retention of expression domains was observed, an apparent absence or reduction of expression was also inferred. We found that two duplC-ATL genes underwent pseudogenization and that, in one case, gene expression is probably regained. Our findings provide insights into the evolution of gene families in plants, defining key events on the expansion of the Arabidopsis Tóxicos en Levadura family of E3 ligases.     

RBR Ubiquitin Ligases: Diversification and Streamlining in Animal Lineages

The patterns of emergence and disappearance in animal species of genes encoding RBR ubiquitin ligases are described. RBR genes can be classified into subfamilies (Parkin, Ariadne, Dorfin, ARA54, etc.) according to sequence and structural data. Here, I show that most animal-specific RBR subfamilies emerged early in animal evolution, and that ancient animals, before the cnidarian/bilaterian split, had a set of RBR genes, which was as complex as the one currently found in mammals. However, some lineages (nematodes, dipteran insects) have recently suffered multiple losses, leading to a highly simplified set of RBR genes. Genes of a particular RBR subfamily, characterized by containing a helicase domain and so far found only in plants, are present also in some animal species. The meaning of these patterns of diversification and streamlining are discussed at the light of functional data. Extreme evolutionary conservation may be related to gene products having housekeeping functions. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Structure of the HHARI Catalytic Domain Shows Glimpses of a HECT E3 Ligase

The ubiquitin-signaling pathway utilizes E1 activating, E2 conjugating, and E3 ligase enzymes to sequentially transfer the small modifier protein ubiquitin to a substrate protein. During the last step of this cascade different types of E3 ligases either act as scaffolds to recruit an E2 enzyme and substrate (RING), or form an ubiquitin-thioester intermediate prior to transferring ubiquitin to a substrate (HECT). The RING-inBetweenRING-RING (RBR) proteins constitute a unique group of E3 ubiquitin ligases that includes the Human Homologue of Drosophila Ariadne (HHARI). These E3 ligases are proposed to use a hybrid RING/HECT mechanism whereby the enzyme uses facets of both the RING and HECT enzymes to transfer ubiquitin to a substrate. We now present the solution structure of the HHARI RING2 domain, the key portion of this E3 ligase required for the RING/HECT hybrid mechanism. The structure shows the domain possesses two Zn²⁺-binding sites and a single exposed cysteine used for ubiquitin catalysis. A structural comparison of the RING2 domain with the HECT E3 ligase NEDD4 reveals a near mirror image of the cysteine and histidine residues in the catalytic site. Further, a tandem pair of aromatic residues exists near the C-terminus of the HHARI RING2 domain that is conserved in other RBR E3 ligases. One of these aromatic residues is remotely located from the catalytic site that is reminiscent of the location found in HECT E3 enzymes where it is used for ubiquitin catalysis. These observations provide an initial structural rationale for the RING/HECT hybrid mechanism for ubiquitination used by the RBR E3 ligases.     

UPL1 and 2, two 405 kDa ubiquitin-protein ligases from Arabidopsis thaliana related to the HECT-domain protein family

The ubiquitin/26S proteasome pathway is a major route for degrading abnormal and important short-lived regulatory proteins in eukaryotes. Covalent attachment of ubiquitin, which triggers entry of target proteins into the pathway, is accomplished by an ATP-dependent reaction cascade involving the sequential action of three enzymes, E1s, E2s and E3s. Although much of the substrate specificity of the pathway is determined by E3s (or ubiquitin-protein ligases, UPLs), little is known about these enzymes in plants and how they choose appropriate targets for ubiquitination. Here, we describe two 405 kDa E3s (UPL1 and 2) from Arabidopsis thaliana related to the HECT-E3 family that is essential in yeast and animals. UPL1 and 2 are encoded by 13 kbp genes 26 cM apart on chromosome I, that are over 95% identical within both the introns and exons, suggesting that the two loci arose from a recent gene duplication. The C-terminal HECT domain of UPL1 is necessary and sufficient to conjugate ubiquitin in vitro in a reaction that requires the positionally conserved cysteine within the HECT domain, E1, and an E2 of the UBC8 family. Given that HECT E3s help define target specificity of the ubiquitin conjugation, a continued characterization of UPL1 and 2 should be instrumental in understanding the functions of ubiquitin-dependent protein turnover in plants and for identifying pathway substrates.     

Structure of the HECT:ubiquitin complex and its role in ubiquitin chain elongation

Several mechanisms have been proposed for the synthesis of substrate-linked ubiquitin chains. HECT ligases directly catalyse protein ubiquitination and have been found to non-covalently interact with ubiquitin. We report crystal structures of the Nedd4 HECT domain, alone and in complex with ubiquitin, which show a new binding mode involving two surfaces on ubiquitin and both subdomains of the HECT N-lobe. The structures suggest a model for HECT-to-substrate ubiquitin transfer, in which the growing chain on the substrate is kept close to the catalytic cysteine to promote processivity. Mutational analysis highlights differences between the processes of substrate polyubiquitination and self-ubiquitination.     

Diversification and Specialization of Plant RBR Ubiquitin Ligases

Background

RBR ubiquitin ligases are components of the ubiquitin-proteasome system present in all eukaryotes. They are characterized by having the RBR (RING – IBR – RING) supradomain. In this study, the patterns of emergence of RBR genes in plants are described.

Methodology/Principal Findings

Phylogenetic and structural data confirm that just four RBR subfamilies (Ariadne, ARA54, Plant I/Helicase and Plant II) exist in viridiplantae. All of them originated before the split that separated green algae from the rest of plants. Multiple genes of two of these subfamilies (Ariadne and Plant II) appeared in early plant evolution. It is deduced that the common ancestor of all plants contained at least five RBR genes and the available data suggest that this number has been increasing slowly along streptophyta evolution, although losses, especially of Helicase RBR genes, have also occurred in several lineages. Some higher plants (e. g. Arabidopsis thaliana, Oryza sativa) contain a very large number of RBR genes and many of them were recently generated by tandem duplications. Microarray data indicate that most of these new genes have low-level and sometimes specific expression patterns. On the contrary, and as occurs in animals, a small set of older genes are broadly expressed at higher levels.

Conclusions/Significance

The available data suggests that the dynamics of appearance and conservation of RBR genes is quite different in plants from what has been described in animals. In animals, an abrupt emergence of many structurally diverse RBR subfamilies in early animal history, followed by losses of multiple genes in particular lineages, occurred. These patterns are not observed in plants. It is also shown that while both plants and animals contain a small, similar set of essential RBR genes, the rest evolves differently. The functional implications of these results are discussed.     

Evolutionary analysis of CBL-interacting protein kinase gene family in plants

Calcineurin B-like proteins (CBLs) and CBL-interacting protein kinases (CIPKs) form the CBL-CIPK complexes, perceiving calcium signals and relaying the signals to downstream responses in plants. To further understand the CBL-CIPK signaling system, here we focused on the evolutionary analysis of CIPKs. We re-evaluated eight plant genomes and identified 146 CIPKs, providing several new CIPKs in rice and poplar. A phylogenetic tree was constructed, showing that these 146 CIPKs are grouped into intron-rich and intron-less clades. Furthermore, all the CIPKs from the non-angiosperm species were found in intron-rich clade. We identified 30 conserved protein motifs among these 146 CIPKs. Analysis of gene duplication showed that the expansion of CIPKs in both clades is partly contributed by segmental duplications, however, tandem duplicates were found only in intron-less clade. Ka/Ks ratios showed that CIPK genes have experienced purifying selective pressure. Additionally, clustering of gene expression revealed that some CIPK genes in two clades share similar expression patterns under abiotic stresses and four CIPKs in intron-less clade form a distinct cluster (i.e., different expression patterns), suggesting the complexity of CIPK gene expression under abiotic stresses. Taken together, our results provided some new insights into the evolution of CIPKs and the hint that the expansion of CIPKs in intron-less clade is adaptive to environmental stresses.     

The Ubiquitin Binding Region of the Smurf HECT Domain Facilitates Polyubiquitylation and Binding of Ubiquitylated Substrates

Mono- and polyubiquitylation of proteins are key steps in a wide range of biological processes. However, the molecular mechanisms that mediate these different events are poorly understood. Here, we employed NMR spectroscopy to map a non-covalent ubiquitin binding surface (UBS) on the Smurf ubiquitin ligase HECT domain. Analysis of mutants of the HECT UBS reveal that interfering with the UBS surface blocked Smurf-dependent degradation of its substrate RhoA in cells. In vitro analysis revealed that the UBS was not required for UbcH7-dependent charging of the HECT catalytic cysteine. Surprisingly, although the UBS was required for polyubiquitylation of both Smurf itself and the Smurf substrate RhoA, it was not required for monoubiquitylation. Furthermore, we show that mutating the UBS interfered with efficient binding of a monoubiquitylated form of RhoA to the Smurf HECT domain. Our findings suggest the UBS promotes polyubiquitylation by stabilizing ubiquitylated substrate binding to the HECT domain.     

Knots in the family tree: evolutionary relationships and functions of knox homeobox genes

Knotted-like homeobox (knox) genes constitute a gene family in plants. Class I knox genes are expressed in shoot apical meristems, and (with notable exceptions) not in lateral organ primordia. Class II genes have more diverse expression patterns. Loss and gain of function mutations indicate that knox genes are important regulators of meristem function. Gene duplication has contributed to the evolution of families of homeodomain proteins in metazoans. We believe that similar mechanisms have contributed to the diversity of knox gene function in plants. Knox genes may have contributed to the evolution of compound leaves in tomato and could be involved in the evolution of morphological traits in other species. Alterations in cis-regulatory regions in some knox genes correlate with novel patterns of gene expression and distinctive morphologies. Preliminary data from the analysis of class I knox gene expression illustrates the evolution of complex patterns of knox expression is likely to have occurred through loss and gain of domains of gene expression.     

Progressive Purkinje Cell Degeneration in tambaleante Mutant Mice Is a Consequence of a Missense Mutation in HERC1 E3 Ubiquitin Ligase

The HERC gene family encodes proteins with two characteristic domains: HECT and RCC1-like. Proteins with HECT domains have been described to function as ubiquitin ligases, and those that contain RCC1-like domains have been reported to function as GTPases regulators. These two activities are essential in a number of important cellular processes such as cell cycle, cell signaling, and membrane trafficking. Mutations affecting these domains have been found associated with retinitis pigmentosa, amyotrophic lateral sclerosis, and cancer. In humans, six HERC genes have been reported which encode two subgroups of HERC proteins: large (HERC1-2) and small (HERC3-6). The giant HERC1 protein was the first to be identified. It has been involved in membrane trafficking and cell proliferation/growth through its interactions with clathrin, M2-pyruvate kinase, and TSC2 proteins. Mutations affecting other members of the HERC family have been found to be associated with sterility and growth retardation. Here, we report the characterization of a recessive mutation named tambaleante, which causes progressive Purkinje cell degeneration leading to severe ataxia with reduced growth and lifespan in homozygous mice aged over two months. We mapped this mutation in mouse chromosome 9 and then performed positional cloning. We found a G⇔A transition at position 1448, causing a Gly to Glu substitution (Gly483Glu) in the highly conserved N-terminal RCC1-like domain of the HERC1 protein. Successful transgenic rescue, with either a mouse BAC containing the normal copy of Herc1 or with the human HERC1 cDNA, validated our findings. Histological and biochemical studies revealed extensive autophagy associated with an increase of the mutant protein level and a decrease of mTOR activity. Our observations concerning this first mutation in the Herc1 gene contribute to the functional annotation of the encoded E3 ubiquitin ligase and underline the crucial and unexpected role of this protein in Purkinje cell physiology.     

Identification of the highly accumulated microRNA*s in Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa)

Plant microRNAs (miRNAs) are crucial for the regulation of gene expression, which is involved in almost all the important biological processes. In the cytoplasm, the miRNA strand is selectively incorporated into a specific Argonaute (AGO)-associated gene silencing complex, while the miRNA* is degraded rapidly. Thus, most miRNA*s were thought to be biologically meaningless. Interestingly, several recent reports in both plants and animals have shaken this notion. Many miRNA*s were demonstrated to possess regulatory roles in gene expression. However, the low accumulation levels of most miRNA*s raise the question whether the activities of this small RNA (sRNA) species are widespread in plants. Here, by using publicly available sRNA high-throughput sequencing data, we found that the accumulation levels of several miRNA*s could be much higher than those of their miRNA partners in certain organs, mutants and/or AGO-associated silencing complexes of both Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa). Based on target prediction and degradome sequencing data-based validation, some of these highly accumulated miRNA*s were indicated to possess cleavage-based potential regulatory role on certain targets. Besides, some interesting biological interpretations were obtained based on the accumulation patterns of the miRNA*s, the annotations of the target genes, and literature mining. Taken together, the expanded list of the highly accumulated miRNA*s along with their potential target genes discovered in this study further strengthened the current notion that certain members of the miRNA* species are biologically relevant, which needs further inspection.     

Plant ubiquitylation and proteolytic systems, the key elements in hormonal signaling pathways

The general function of the ubiquitylation systems is to conjugate ubiquitin to lysine residues within substrate proteins, thus targeting them for degradation by the proteasome. In Arabidopsis thaliana more than 1300 genes (approximately 5% of the proteome) encode components of the ubiquitin/26S proteasome pathway. Approximately 90% of these genes encode subunits of the E3 ubiquitin ligases, which confer substrate specificity to the ubiquitin/26S proteasome pathway. The plant E3 ubiquitin ligases comprise a large and diverse family of proteins or protein complexes containing either a HECT domain, a RING-finger or U-box domain. The SCF class of E3 ligases is the most thoroughly studied in plants because some of them participate in regulation of hormone signaling pathways. The role of the SCF is to ubiquitylate repressors of hormone response (auxin, gibberellins), whereas in response to ethylene, abscisic acid and brassinosteroids the SCF participate in degradation of positive regulators in the absence of the hormone.     

HECT ubiquitin ligases link viral and cellular PPXY motifs to the vacuolar protein-sorting pathway

Many enveloped viruses exploit the class E vacuolar protein-sorting (VPS) pathway to bud from cells, and use peptide motifs to recruit specific class E VPS factors. Homologous to E6AP COOH terminus (HECT) ubiquitin ligases have been implicated as cofactors for PPXY motif-dependent budding, but precisely which members of this family are responsible, and how they access the VPS pathway is unclear. Here, we show that PPXY-dependent viral budding is unusually sensitive to inhibitory fragments derived from specific HECT ubiquitin ligases, namely WWP1 and WWP2. We also show that WWP1, WWP2, or Itch ubiquitin ligase recruitment promotes PPXY-dependent virion release, and that this function requires that the HECT ubiquitin ligase domain be catalytically active. Finally, we show that several mammalian HECT ubiquitin ligases, including WWP1, WWP2, and Itch are recruited to class E compartments induced by dominant negative forms of the class E VPS ATPase, VPS4. These data indicate that specific HECT ubiquitin ligases can link PPXY motifs to the VPS pathway to induce viral budding.     

Evolution of the F-Box Gene Family in Euarchontoglires: Gene Number Variation and Selection Patterns

F-box proteins are substrate adaptors used by the SKP1–CUL1–F-box protein (SCF) complex, a type of E3 ubiquitin ligase complex in the ubiquitin proteasome system (UPS). SCF-mediated ubiquitylation regulates proteolysis of hundreds of cellular proteins involved in key signaling and disease systems. However, our knowledge of the evolution of the F-box gene family in Euarchontoglires is limited. In the present study, 559 F-box genes and nine related pseudogenes were identified in eight genomes. Lineage-specific gene gain and loss events occurred during the evolution of Euarchontoglires, resulting in varying F-box gene numbers ranging from 66 to 81 among the eight species. Both tandem duplication and retrotransposition were found to have contributed to the increase of F-box gene number, whereas mutation in the F-box domain was the main mechanism responsible for reduction in the number of F-box genes, resulting in a balance of expansion and contraction in the F-box gene family. Thus, the Euarchontoglire F-box gene family evolved under a birth-and-death model. Signatures of positive selection were detected in substrate-recognizing domains of multiple F-box proteins, and adaptive changes played a role in evolution of the Euarchontoglire F-box gene family. In addition, single nucleotide polymorphism (SNP) distributions were found to be highly non-random among different regions of F-box genes in 1092 human individuals, with domain regions having a significantly lower number of non-synonymous SNPs.     

The life-cycle of operons

Operons are a major feature of all prokaryotic genomes, but how and why operon structures vary is not well understood. To elucidate the life-cycle of operons, we compared gene order between Escherichia coli K12 and its relatives and identified the recently formed and destroyed operons in E. coli. This allowed us to determine how operons form, how they become closely spaced, and how they die. Our findings suggest that operon evolution may be driven by selection on gene expression patterns. First, both operon creation and operon destruction lead to large changes in gene expression patterns. For example, the removal of lysA and ruvA from ancestral operons that contained essential genes allowed their expression to respond to lysine levels and DNA damage, respectively. Second, some operons have undergone accelerated evolution, with multiple new genes being added during a brief period. Third, although genes within operons are usually closely spaced because of a neutral bias toward deletion and because of selection against large overlaps, genes in highly expressed operons tend to be widely spaced because of regulatory fine-tuning by intervening sequences. Although operon evolution may be adaptive, it need not be optimal: new operons often comprise functionally unrelated genes that were already in proximity before the operon formed.     

Evolution of Stress-Regulated Gene Expression in Duplicate Genes of Arabidopsis thaliana

Due to the selection pressure imposed by highly variable environmental conditions, stress sensing and regulatory response mechanisms in plants are expected to evolve rapidly. One potential source of innovation in plant stress response mechanisms is gene duplication. In this study, we examined the evolution of stress-regulated gene expression among duplicated genes in the model plant Arabidopsis thaliana. Key to this analysis was reconstructing the putative ancestral stress regulation pattern. By comparing the expression patterns of duplicated genes with the patterns of their ancestors, duplicated genes likely lost and gained stress responses at a rapid rate initially, but the rate is close to zero when the synonymous substitution rate (a proxy for time) is >~0.8. When considering duplicated gene pairs, we found that partitioning of putative ancestral stress responses occurred more frequently compared to cases of parallel retention and loss. Furthermore, the pattern of stress response partitioning was extremely asymmetric. An analysis of putative cis-acting DNA regulatory elements in the promoters of the duplicated stress-regulated genes indicated that the asymmetric partitioning of ancestral stress responses are likely due, at least in part, to differential loss of DNA regulatory elements; the duplicated genes losing most of their stress responses were those that had lost more of the putative cis-acting elements. Finally, duplicate genes that lost most or all of the ancestral responses are more likely to have gained responses to other stresses. Therefore, the retention of duplicates that inherit few or no functions seems to be coupled to neofunctionalization. Taken together, our findings provide new insight into the patterns of evolutionary changes in gene stress responses after duplication and lay the foundation for testing the adaptive significance of stress regulatory changes under highly variable biotic and abiotic environments.