Comparative genomics and protein domain graph analyses link ubiquitination and RNA metabolism

Parkin inhibits proliferation and migration of bladder cancer via ubiquitinating Catalase

PRKN is a key gene involved in mitophagy in Parkinson's disease. However, recent studies have demonstrated that it also plays a role in the development and metastasis of several types of cancers, both in a mitophagy-dependent and mitophagy-independent manner. Despite this, the potential effects and underlying mechanisms of Parkin on bladder cancer (BLCA) remain unknown. Therefore, in this study, we investigated the expression of Parkin in various BLCA cohorts derived from human. Here we show that PRKN expression was low and that PRKN acts as a tumor suppressor by inhibiting the proliferation and migration of BLCA cells in a mitophagy-independent manner. We further identified Catalase as a binding partner and substrate of Parkin, which is an important antioxidant enzyme that regulates intracellular ROS levels during cancer progression. Our data showed that knockdown of CAT led to increased intracellular ROS levels, which suppressed cell proliferation and migration. Conversely, upregulation of Catalase decreased intracellular ROS levels, promoting cell growth and migration. Importantly, we found that Parkin upregulation partially restored these effects. Moreover, we discovered that USP30, a known Parkin substrate, could deubiquitinate and stabilize Catalase. Overall, our study reveals a novel function of Parkin and identifies a potential therapeutic target in BLCA.

Domain Architecture Based Methods for Comparative Functional Genomics Toward Therapeutic Drug Target Discovery

Genes duplicate, mutate, recombine, fuse or fission to produce new genes, or when genes are formed from de novo, novel functions arise during evolution. Researchers have tried to quantify the causes of these molecular diversification processes to know how these genes increase molecular complexity over a period of time, for instance protein domain organization. In contrast to global sequence similarity, protein domain architectures can capture key structural and functional characteristics, making them better proxies for describing functional equivalence. In Prokaryotes and eukaryotes it has proven that, domain designs are retained over significant evolutionary distances. Protein domain architectures are now being utilized to categorize and distinguish evolutionarily related proteins and find homologs among species that are evolutionarily distant from one another. Additionally, structural information stored in domain structures has accelerated homology identification and sequence search methods. Tools for functional protein annotation have been developed to discover, protein domain content, domain order, domain recurrence, and domain position as all these contribute to the prediction of protein functional accuracy. In this review, an attempt is made to summarise facts and speculations regarding the use of protein domain architecture and modularity to identify possible therapeutic targets among cellular activities based on the understanding their linked biological processes.

Linking the plant stress responses with RNA helicases

RNA helicases are omnipresent plant proteins across all kingdoms and have been demonstrated to play an essential role in all cellular processes involving nucleic acids. Currently, these proteins emerged as a new tool for plant molecular biologists to modulate plant stress responses. Here, we review the crucial role of RNA helicases triggered by biotic, abiotic, and multiple stress conditions. In this review, the emphasis has been given on the role of these proteins upon viral stress. Further, we have explored RNA helicase mediated regulation of RNA metabolism, starting from ribosome biogenesis to its decay upon stress induction. We also highlighted the cross-talk between RNA helicase, phytohormones, and ROS. Different overexpression and transgenic studies have been provided in the text to indicate the stress tolerance abilities of these proteins.

Origin and evolution of fungal HECT ubiquitin ligases

Ubiquitin ligases (E3s) are basic components of the eukaryotic ubiquitination system. In this work, the emergence and diversification of fungal HECT ubiquitin ligases is described. Phylogenetic and structural data indicate that six HECT subfamilies (RSP5, TOM1, UFD4, HUL4, HUL4A and HUL5) existed in the common ancestor of all fungi. These six subfamilies have evolved very conservatively, with only occasional losses and duplications in particular fungal lineages. However, an early, drastic reduction in the number of HECT genes occurred in microsporidians, in parallel to the reduction of their genomes. A significant correlation between the total number of genes and the number of HECT-encoding genes present in fungi has been observed. However, transitions from unicellularity to multicellularity or vice versa apparently had no effect on the evolution of this family. Likely orthologs or co-orthologs of all fungal HECT genes have been detected in animals. Four genes are deduced to be present in the common ancestor of fungi, animals and plants. Protein-protein interactions detected in both the yeast Saccharomyces cerevisiae and humans suggest that some ancient functions of HECT proteins have been conserved since the animals/fungi split.

Genome-wide identification and characterization of RBR ubiquitin ligase genes in soybean

RBR (RING1-IBR-RING2) proteins play an important role in protein ubiquitination and are involved in many cellular processes. Recent studies showed plant RBR genes were induced by abiotic and biotic stresses. However, detailed studies on RBR genes in the important oil crop, soybean (Glycine max (L.) Merr.), is still lacking. Here we performed a genome-wide search and identified 24 RBR domain-containing genes from the soybean genome sequence and cloned 11 of them. Most soybean RBR proteins contain a highly conserved RBR supra-domain. Phylogenetic analyses indicated all 24 soybean RBR proteins are most related to the RBR proteins from Phaseolus vulgaris, and could be classified into seven groups including Ariadne A, Ariadne B, ARA54, Plant IIA, Plant IIB, Plant IIC, and Helicase. Tandem duplication and block duplication were found among the Ariadne B and Plant IIC group of soybean RBR genes. Despite the conserved RBR supra-domain, there are extensive variations in the additional protein motifs and exon-intron structures between different groups, which indicate they might have diverse functions. Most soybean RBR proteins are predicted to localize in nucleus, and four of them were experimentally confirmed by GFP fusion proteins. Soybean RBR genes are broadly expressed in many tissue types with a little more abundant in the roots and flowers than leaves, stems, and seeds. The expression of GmRTRTP3 (Plant IIB) and GmRTRTP5 (Plant IIC) are induced by NaCl treatment, which suggests these RBR genes might be involved in soybean response to abiotic stresses.

A molecular explanation for the recessive nature of parkin-linked Parkinson's disease

Mutations in the park2 gene, encoding the RING-inBetweenRING-RING E3 ubiquitin ligase parkin, cause 50% of autosomal recessive juvenile Parkinsonism cases. More than 70 known pathogenic mutations occur throughout parkin, many of which cluster in the inhibitory amino-terminal ubiquitin-like domain, and the carboxy-terminal RING2 domain that is indispensable for ubiquitin transfer. A structural rationale showing how autosomal recessive juvenile Parkinsonism mutations alter parkin function is still lacking. Here we show that the structure of parkin RING2 is distinct from canonical RING E3 ligases and lacks key elements required for E2-conjugating enzyme recruitment. Several pathogenic mutations in RING2 alter the environment of a single surface-exposed catalytic cysteine to inhibit ubiquitination. Native parkin adopts a globular inhibited conformation in solution facilitated by the association of the ubiquitin-like domain with the RING-inBetweenRING-RING C-terminus. Autosomal recessive juvenile Parkinsonism mutations disrupt this conformation. Finally, parkin autoubiquitinates only in cis, providing a molecular explanation for the recessive nature of autosomal recessive juvenile Parkinsonism.

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.

Evolution of plant HECT ubiquitin ligases

HECT ubiquitin ligases are key components of the ubiquitin-proteasome system, which is present in all eukaryotes. In this study, the patterns of emergence of HECT genes in plants are described. Phylogenetic and structural data indicate that viridiplantae have six main HECT subfamilies, which arose before the split that separated green algae from the rest of plants. It is estimated that the common ancestor of all plants contained seven HECT genes. Contrary to what happened in animals, the number of HECT genes has been kept quite constant in all lineages, both in chlorophyta and streptophyta, although evolutionary recent duplications are found in some species. Several of the genes found in plants may have originated very early in eukaryotic evolution, given that they have clear similarities, both in sequence and structure, to animal genes. Finally, in Arabidopsis thaliana, we found significant correlations in the expression patterns of HECT genes and some ancient, broadly expressed genes that belong to a different ubiquitin ligase family, called RBR. These results are discussed in the context of the evolution of the gene families required for ubiquitination in plants.

Origin and diversification of TRIM ubiquitin ligases

Most proteins of the TRIM family (also known as RBCC family) are ubiquitin ligases that share a peculiar protein structure, characterized by including an N-terminal RING finger domain closely followed by one or two B-boxes. Additional protein domains found at their C termini have been used to classify TRIM proteins into classes. TRIMs are involved in multiple cellular processes and many of them are essential components of the innate immunity system of animal species. In humans, it has been shown that mutations in several TRIM-encoding genes lead to diverse genetic diseases and contribute to several types of cancer. They had been hitherto detected only in animals. In this work, by comprehensively analyzing the available diversity of TRIM and TRIM-like protein sequences and evaluating their evolutionary patterns, an improved classification of the TRIM family is obtained. Members of one of the TRIM subfamilies defined, called Subfamily A, turn to be present not only in animals, but also in many other eukaryotes, such as fungi, apusozoans, alveolates, excavates and plants. The rest of subfamilies are animal-specific and several of them originated only recently. Subfamily A proteins are characterized by containing a MATH domain, suggesting a potential evolutionary connection between TRIM proteins and a different type of ubiquitin ligases, known as TRAFs, which contain quite similar MATH domains. These results indicate that the TRIM family emerged much earlier than so far thought and contribute to our understanding of its origin and diversification. The structural and evolutionary links with the TRAF family of ubiquitin ligases can be experimentally explored to determine whether functional connections also exist.

Deciphering network community structure by surprise

The analysis of complex networks permeates all sciences, from biology to sociology. A fundamental, unsolved problem is how to characterize the community structure of a network. Here, using both standard and novel benchmarks, we show that maximization of a simple global parameter, which we call Surprise (S), leads to a very efficient characterization of the community structure of complex synthetic networks. Particularly, S qualitatively outperforms the most commonly used criterion to define communities, Newman and Girvan's modularity (Q). Applying S maximization to real networks often provides natural, well-supported partitions, but also sometimes counterintuitive solutions that expose the limitations of our previous knowledge. These results indicate that it is possible to define an effective global criterion for community structure and open new routes for the understanding of complex networks.

The plexus model for the inference of ancestral multidomain proteins

Interactions of protein domains control essential cellular processes. Thus, inferring the evolutionary histories of multidomain proteins in the context of their families can provide rewarding insights into protein function. However, methods to infer these histories are challenged by the complexity of macroevolutionary events. Here, we address this challenge by describing an algorithm that computes a novel network-like structure, called plexus, which represents the evolution of domains and their combinations. Finally, we demonstrate the performance of this algorithm with empirical data sets.

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.

Jerarca: efficient analysis of complex networks using hierarchical clustering

Background

How to extract useful information from complex biological networks is a major goal in many fields, especially in genomics and proteomics. We have shown in several works that iterative hierarchical clustering, as implemented in the UVCluster program, is a powerful tool to analyze many of those networks. However, the amount of computation time required to perform UVCluster analyses imposed significant limitations to its use.

Methodology/Principal Findings

We describe the suite Jerarca, designed to efficiently convert networks of interacting units into dendrograms by means of iterative hierarchical clustering. Jerarca is divided into three main sections. First, weighted distances among units are computed using up to three different approaches: a more efficient version of UVCluster and two new, related algorithms called RCluster and SCluster. Second, Jerarca builds dendrograms based on those distances, using well-known phylogenetic algorithms, such as UPGMA or Neighbor-Joining. Finally, Jerarca provides optimal partitions of the trees using statistical criteria based on the distribution of intra- and intercluster connections. Outputs compatible with the phylogenetic software MEGA and the Cytoscape package are generated, allowing the results to be easily visualized.

Conclusions/Significance

The four main advantages of Jerarca in respect to UVCluster are: 1) Improved speed of a novel UVCluster algorithm; 2) Additional, alternative strategies to perform iterative hierarchical clustering; 3) Automatic evaluation of the hierarchical trees to obtain optimal partitions; and, 4) Outputs compatible with popular software such as MEGA and Cytoscape.

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.

Animal HECT ubiquitin ligases: evolution and functional implications

Background

HECT ubiquitin ligases (HECT E3s) are key components of the eukaryotic ubiquitin-proteasome system and are involved in the genesis of several human diseases. In this study, I analyze the patterns of diversification of HECT E3s since animals emerged in order to provide the right framework to understand the functional data available for proteins of this family.

Results

I show that the current classification of HECT E3s into three groups (NEDD4-like E3s, HERCs and single-HECT E3s) is fundamentally incorrect. First, the existence of a "Single-HECT E3s" group is not supported by phylogenetic analyses. Second, the HERC proteins must be divided into two subfamilies (Large HERCs, Small HERCs) that are evolutionarily very distant, their structural similarity being due to convergence and not to a common origin. Sequence and structural analyses show that animal HECT E3s can be naturally classified into 16 subfamilies. Almost all of them appeared either before animals originated or in early animal evolution. More recently, multiple gene losses have occurred independently in some lineages (nematodes, insects, urochordates), the same groups that have also lost genes of another type of E3s (RBR family). Interestingly, the emergence of some animal HECT E3s precedes the origin of key cellular systems that they regulate (TGF-β and EGF signal transduction pathways; p53 family of transcription factors) and it can be deduced that distantly related HECT proteins have been independently co-opted to perform similar roles. This may contribute to explain why distantly related HECT E3s are involved in the genesis of multiple types of cancer.

Conclusions

The complex evolutionary history of HECT ubiquitin ligases in animals has been deciphered. The most appropriate model animals to study them and new theoretical and experimental lines of research are suggested by these results.

Interactome and Gene Ontology provide congruent yet subtly different views of a eukaryotic cell

BACKGROUND

The characterization of the global functional structure of a cell is a major goal in bioinformatics and systems biology. Gene Ontology (GO) and the protein-protein interaction network offer alternative views of that structure.

RESULTS

This study presents a comparison of the global structures of the Gene Ontology and the interactome of Saccharomyces cerevisiae. Sensitive, unsupervised methods of clustering applied to a large fraction of the proteome led to establish a GO-interactome correlation value of +0.47 for a general dataset that contains both high and low-confidence interactions and +0.58 for a smaller, high-confidence dataset.

CONCLUSION

The structures of the yeast cell deduced from GO and interactome are substantially congruent. However, some significant differences were also detected, which may contribute to a better understanding of cell function and also to a refinement of the current ontologies.

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.

Plant RNA helicases: linking aberrant and silencing RNA

RNA helicases are ATPases that are capable of rearranging RNA and ribonucleoprotein (RNP) structure, and they can potentially function in any aspect of RNA metabolism. The RNA helicase gene family of plant genomes is larger and more diverse than genome families observed in other systems and provides an ideal model for investigation of the physiological importance of RNA secondary structure rearrangement in plant development. Numerous plant RNA helicases are associated with a variety of physiological functions, but this review will focus on the thirteen RNA helicases associated with the metabolism of aberrant and silencing RNAs. The results emphasize the crucial role RNA helicase activity has in the regulation of mRNA quality control and gene expression in plant development.

Identification of a novel Zn2+-binding domain in the autosomal recessive juvenile Parkinson-related E3 ligase parkin

Missense mutations in park2, encoding the parkin protein, account for approximately 50% of autosomal recessive juvenile Parkinson disease (ARJP) cases. Parkin belongs to the family of RBR (RING-between-RING) E3 ligases involved in the ubiquitin-mediated degradation and trafficking of proteins such as Pael-R and synphillin-1. The proposed architecture of parkin, based largely on sequence similarity studies, consists of N-terminal ubiquitin-like and C-terminal RBR domains. These domains are separated by a approximately 160-residue unique parkin sequence having no recognizable domain structure. We used limited proteolysis experiments on bacterially expressed and purified parkin to identify a new domain (RING0) within the unique parkin domain sequence. RING0 comprises two distinct, conserved cysteine-rich clusters between Cys(150)-Cys(169) and Cys(196)-His(215) consisting of CX(2)-(3)CX(11)CX(2)C and CX(4-6)CX(10-16)-CX(2)(H/C) motifs. The positions of the cysteine/histidine residues in this region bear similarity to parkin RING1 and RING2 domains, as well as other E3 ligase RING domains. However, in parkin a 26-residue linker region separates the motifs, which is not typical of other RING domain structures. Further, the RING0 domain includes all but one of the known ARJP mutation sites between the ubiquitin-like and RBR regions of parkin. Using electrospray ionization mass spectrometry and inductively coupled plasma-atomic emission spectrometry analysis, we determined that the RING0, RING1, IBR, and RING2 domains each bind two Zn(2+) ions, the first observation of an E3 ligase with the ability to bind eight metal ions. Removal of the zinc from parkin causes near complete unfolding of the protein, an observation that rationalizes cysteine-based ARJP mutations found throughout parkin, including RING0 (C212Y) that form cellular inclusions and/or are defective for ubiquitination likely because of poor zinc binding and misfolding. The identification of the RING0 domain in parkin provides a new overall domain structure for the protein that will be important in assessing the roles of ARJP mutations and designing experiments aimed at understanding the disease.

Basic networks: definition and applications

We define basic networks as the undirected subgraphs with minimal number of units in which the distances (geodesics, minimal path lengths) among a set of selected nodes, which we call seeds, in the original graph are conserved. The additional nodes required to draw the basic network are called connectors. We describe a heuristic strategy to find the basic networks of complex graphs. We also show how the characterization of these networks may help to obtain relevant biological information from highly complex protein-protein interaction data.

Just how versatile are domains?

Background

Creating new protein domain arrangements is a frequent mechanism of evolutionary innovation. While some domains always form the same combinations, others form many different arrangements. This ability, which is often referred to as versatility or promiscuity of domains, its a random evolutionary model in which a domain's promiscuity is based on its relative frequency of domains.

Results

We show that there is a clear relationship across genomes between the promiscuity of a given domain and its frequency. However, the strength of this relationship differs for different domains. We thus redefine domain promiscuity by defining a new index, DV I ("domain versatility index"), which eliminates the effect of domain frequency. We explore links between a domain's versatility, when unlinked from abundance, and its biological properties.

Conclusion

Our results indicate that domains occurring as single domain proteins and domains appearing frequently at protein termini have a higher DV I. This is consistent with previous observations that the evolution of domain re-arrangements is primarily driven by fusion of pre-existing arrangements and single domains as well as loss of domains at protein termini. Furthermore, we studied the link between domain age, defined as the first appearance of a domain in the species tree, and the DV I. Contrary to previous studies based on domain promiscuity, it seems as if the DV I is age independent. Finally, we find that contrary to previously reported findings, versatility is lower in Eukaryotes. In summary, our measure of domain versatility indicates that a random attachment process is sufficient to explain the observed distribution of domain arrangements and that several views on domain promiscuity need to be revised.

Arrangements in the modular evolution of proteins

It has been known for the last couple of decades that proteins evolve partly through rearrangements of larger fragments, typically domains. These units are considered the basic modules of protein structure, evolution and function. In the last few years, the analysis of protein-domain rearrangements has provided us with functional and evolutionary insights and has aided improved functional predictions and domain assignments to previously uncharacterised genes and proteins. Although some mechanisms that govern modular rearrangements of protein domains have been uncovered, such as the addition or deletion of a single N- or C-terminal domain, much is still unknown about the genetics behind these arrangements.

Ancient origin of the Parkinson disease gene LRRK2

Dominant mutations in the LRRK2 gene, a member of the Roco family, cause both familial and sporadic Parkinson disease. LRRK genes had so far been detected only in bilaterian animals. In deuterostomes, including humans, two LRRK genes (LRRK1 and LRRK2) exist, while in protostomes a single LRRK gene has been found. In this study, I combine structural and phylogenetic analyses to show that the cnidarian Nematostella vectensis has four LRRK genes. One of them is a bona fide orthologue of the human LRRK2 gene, demonstrating that this gene has an ancient origin. Two others are, respectively, orthologues of the deuterostome LRRK1 and the protostome LRRK genes. The fourth gene is probably cnidarian-specific. This precise characterization of the early evolution of LRRK genes in animals has important implications, because it indicates that the Drosophila and Caenorhabditis LRRK genes, which are studied to gain an understanding of LRRK2 function, are not true orthologues of the human Parkinson disease gene. Novel functional insights are also gained by comparison of the structures of LRRK2 genes in distantly related species.

SPIN1, a K homology domain protein negatively regulated and ubiquitinated by the E3 ubiquitin ligase SPL11, is involved in flowering time control in rice

The rice (Oryza sativa) E3 ligase SPOTTED LEAF11 (SPL11) negatively regulates programmed cell death and disease resistance. We demonstrate here that SPL11 also regulates flowering via interaction with SPIN1 (for SPL11-interacting protein1), a Signal Transduction and Activation of RNA family member. SPIN1 binds RNA and DNA in vitro and interacts with SPL11 in the nucleus. Spl11 mutants have delayed flowering under long-day conditions. Spin1 overexpression causes late flowering independently of daylength; expression analyses of flowering marker genes in these lines suggested that SPIN1 represses flowering by downregulating the flowering promoter gene Heading date3a (Hd3a) via Hd1-dependent mechanisms in short days and by targeting Hd1-independent factors in long days. Both Spin1 and Spl11 are regulated diurnally in opposing phases. SPL11 negatively regulates Spin1 transcript levels, while SPIN1 also affects Spl11 expression. Moreover, we show that coincidence of high accumulation of Spin1 mRNA with the light in the morning and early evening is needed to repress flowering. SPIN1 is monoubiquitinated by SPL11, suggesting that it is not targeted for degradation. Our data are consistent with a model in which SPIN1 acts as a negative regulator of flowering that itself is negatively regulated by SPL11, possibly via ubiquitination.

A general strategy to determine the congruence between a hierarchical and a non-hierarchical classification

Background

Classification procedures are widely used in phylogenetic inference, the analysis of expression profiles, the study of biological networks, etc. Many algorithms have been proposed to establish the similarity between two different classifications of the same elements. However, methods to determine significant coincidences between hierarchical and non-hierarchical partitions are still poorly developed, in spite of the fact that the search for such coincidences is implicit in many analyses of massive data.

Results

We describe a novel strategy to compare a hierarchical and a dichotomic non-hierarchical classification of elements, in order to find clusters in a hierarchical tree in which elements of a given "flat" partition are overrepresented. The key improvement of our strategy respect to previous methods is using permutation analyses of ranked clusters to determine whether regions of the dendrograms present a significant enrichment. We show that this method is more sensitive than previously developed strategies and how it can be applied to several real cases, including microarray and interactome data. Particularly, we use it to compare a hierarchical representation of the yeast mitochondrial interactome and a catalogue of known mitochondrial protein complexes, demonstrating a high level of congruence between those two classifications. We also discuss extensions of this method to other cases which are conceptually related.

Conclusion

Our method is highly sensitive and outperforms previously described strategies. A PERL script that implements it is available at .

The ring between ring fingers (RBR) protein family

An overview of the large and functionally diverse RBR protein family that mediates protein-protein interactions of various kinds in development and disease.

Summary

Proteins of the ring between ring fingers (RBR)-domain family are characterized by three groups of specifically clustered (typically eight) cysteine and histidine residues. Whereas the amino-terminal ring domain (N-RING) binds two zinc ions and folds into a classical cross-brace ring finger, the carboxy-terminal ring domain (C-RING) involves only one zinc ion. The three-dimensional structure of the central ring domain, the IBR domain, is still unsolved. About 400 genes coding for RBR proteins have been identified in the genomes of uni- and multicellular eukaryotes and some of their viruses, but the family has not been found in archaea or bacteria. The RBR proteins are classified into 15 major subfamilies (besides some orphan cases) by the phylogenetic relationships of the RBR segments and the conservation of their sequence architecture. The RBR domain mediates protein-protein interactions and a subset of RBR proteins has been shown to function as E3 ubiquitin ligases. RBR proteins have attracted interest because of their involvement in diseases such as parkinsonism, dementia with Lewy bodies, and Alzheimer's disease, and in susceptibility to some intracellular bacterial pathogens. Here, we present an overview of the RBR-domain containing proteins and their subcellular localization, additional domains, function, specificity, and regulation.

The Parkinson Disease Gene LRRK2: Evolutionary and Structural Insights

Mutations in the human leucine-rich repeat kinase 2 (LRRK2) gene are associated with both familial and sporadic Parkinson disease (PD). LRRK2 belongs to a gene family known as Roco. Roco genes encode for large proteins with several protein domains. Particularly, all Roco proteins have a characteristic GTPase domain, named Roc, plus a domain of unknown function called COR. In addition, LRRK2 and several other Roco proteins also contain a protein kinase domain. In this study, I use a combination of phylogenetic and structural analyses of the COR, Roc, and kinase domains present in Roco proteins to describe the origin and evolutionary history of LRRK2. Phylogenetic analyses using these domains demonstrate that LRRK2 emerged from a duplication that occurred after the protostome-deuterostome split. The duplication was followed by the acquisition by LRRK2 proteins of a specific type of N-terminal repeat, described here for the first time. This repeat is absent in the proteins encoded by the paralogs of LRRK2, called LRRK1 or in protostome LRRK proteins. These results suggest that Drosophila or Caenorhabditis LRRK genes may not be good models to understand human LRRK2 function. Genes in the slime mold Dictyostelium discoideum with structures very similar to those found in animal LRRK genes, including the protein kinase domain, have been described. However, phylogenetic analyses suggest that this structural similarity is due to independent acquisitions of distantly related protein kinase domains. Finally, I confirm in an extensive sequence analysis that the Roc GTPase domain is related but still substantially different from small GTPases, such as Rab, Ras, or Rho. Modeling based on known kinase structures suggests that mutations in LRRK2 that cause familiar PD may alter the local 3-dimensional folding of the LRRK2 protein without affecting its overall structure.

UVPAR: fast detection of functional shifts in duplicate genes

BACKGROUND

The imprint of natural selection on gene sequences is often difficult to detect. A plethora of methods have been devised to detect genetic changes due to selective processes. However, many of those methods depend heavily on underlying assumptions regarding the mode of change of DNA sequences and often require sophisticated mathematical treatments that made them computationally slow. The development of fast and effective methods to detect modifications in the selective constraints of genes is therefore of great interest.

RESULTS

We describe UVPAR, a program designed to quickly test for changes in the functional constraints of duplicate genes. Starting with alignments of the proteins encoded by couples of duplicate genes in two different species, UVPAR detects the regions in which modifications of the functional constraints in the paralogs occurred since both species diverged. Sequences can be analyzed with UVPAR in just a few minutes on a standard PC computer. To demonstrate the power of the program, we first show how the results obtained with UVPAR compare to those based on other approaches, using data for vertebrate Hox genes. We then describe a comprehensive study of the RBR family of ubiquitin ligases in which we have performed 529 analyses involving 14 duplicate genes in seven model species. A significant increase in the number of functional shifts was observed for the species Danio rerio and for the gene Ariadne-2.

CONCLUSION

These results show that UVPAR can be used to generate sensitive analyses to detect changes in the selection constraints acting on paralogs. The high speed of the program allows its application to genome-scale analyses.