The ubiquilin gene family: evolutionary patterns and functional insights
© Marín; licensee BioMed Central Ltd. 2014
Received: 26 September 2013
Accepted: 17 March 2014
Published: 28 March 2014
Ubiquilins are proteins that function as ubiquitin receptors in eukaryotes. Mutations in two ubiquilin-encoding genes have been linked to the genesis of neurodegenerative diseases. However, ubiquilin functions are still poorly understood.
In this study, evolutionary and functional data are combined to determine the origin and diversification of the ubiquilin gene family and to characterize novel potential roles of ubiquilins in mammalian species, including humans. The analysis of more than six hundred sequences allowed characterizing ubiquilin diversity in all the main eukaryotic groups. Many organisms (e. g. fungi, many animals) have single ubiquilin genes, but duplications in animal, plant, alveolate and excavate species are described. Seven different ubiquilins have been detected in vertebrates. Two of them, here called UBQLN5 and UBQLN6, had not been hitherto described. Significantly, marsupial and eutherian mammals have the most complex ubiquilin gene families, composed of up to 6 genes. This exceptional mammalian-specific expansion is the result of the recent emergence of four new genes, three of them (UBQLN3, UBQLN5 and UBQLNL) with precise testis-specific expression patterns that indicate roles in the postmeiotic stages of spermatogenesis. A gene with related features has independently arisen in species of the Drosophila genus. Positive selection acting on some mammalian ubiquilins has been detected.
The ubiquilin gene family is highly conserved in eukaryotes. The infrequent lineage-specific amplifications observed may be linked to the emergence of novel functions in particular tissues.
KeywordsUbiquitination Spermiogenesis Ubiquitin receptors Gene duplication Neurodegeneration
Mutations in the Saccharomyces cerevisiae gene DSK2 were found to be suppressors of temperature-sensitive mutations in KAR1, a gene involved in duplication of the yeast spindle pole body . Soon it became clear that DSK2 was a member of a family of evolutionary conserved genes, present not only in yeasts but also in animals and plants, whose protein products are characterized by having a N-terminal ubiquitin-like (UBL) domain, an C-terminal ubiquitin-associated (UBA) domain and a variable number of internal Sti1 repeats [2–6] These proteins are today commonly known as ubiquilins. Mammals have several ubiquilin genes [4, 7–9] Early studies demonstrated that three of them, called UBQLN1 (formerly known also as PLIC-1), UBQLN2 (a. k. a. PLIC-2) and UBQLN4 (a. k. a. A1Up, UBIN, CIP75), are widely expressed in human, mouse or rat, while a fourth one, UBQLN3, is testis-specific in both human and mouse [4, 6, 8–12]. A fifth ubiquilin gene, called UBQLNL, was later detected in humans (first mentioned in ).
Ubiquilins are functionallly linked to the ubiquitin-proteasome system [7, 14]. The UBL domain interacts with the proteasome and also with proteins containing ubiquitin-interacting motifs (UIMs), while the UBA domain serves to interact with polyubiquitinated proteins, at the same time protecting ubiquilins from proteasomal degradation [14–20]. UBL and UBA domains can also mutually interact . These results suggested that ubiquilins might function as ubiquitin receptors , i. e. they would contact with ubiquitinated proteins either to deliver them to the proteasome for degradation or to make them enter other destruction pathways (e. g. autophagy). This has been shown to be true not only for ubiquilins, but for proteins with related structures, also containing UBL and UBA domains, such as yeast Rad23 and its ortholog proteins in other eukaryotes (reviewed in [23–26]).
Multiple results linked ubiquilins to several neurodegenerative diseases. One of them is Alzheimer’s disease (reviewed in [27, 28]). UBQLN1 interacts with presenilins  and overexpression of either UBQLN1 or UBQLN2 protects presenilins from degradation [6, 29]. Also, a particular polymorphism in a UBQLN1 intron may increase the risk of suffering Alzheimer’s disease [30–37]. Additional results linking UBQLN1 with the quality control of Alzheimer’s disease-related proteins have been found recently [38–40]. Finally, reduced UBQLN1 levels were found in the brain cortex of Alzheimer’s disease patients . Drosophila ubiquilin, encoded by the Ubqn (a. k. a. dUbqln) gene, may have related functions [41, 42]. A second disease is amyotrophic lateral sclerosis (ALS). Missense mutations in the proline residues of PXX repeats present in UBQLN2 were found to cause sex-linked, dominant ALS, often associated to frontotemporal dementia . Later, additional missense mutations outside those repeats were also linked to ALS [44–46]. In spinal motor neurons, UBQLN2 appears in ubiquitin-rich protein aggregates typical of ALS, both in patients with UBQLN2 mutations and in patients lacking those mutations . UBQLN1 also interacts with TDP-43, a protein involved in ALS-specific protein aggregates, and TDP-43 aggregates with either UBQLN1 or UBQLN2 in cell systems [44–47]. Similar results have been found in Drosophila. Characteristic ubiquilin-containing aggregates are also found in ALS patients with hexanucleotide expansions in the non-coding region of the C9orf72 gene, which is a common mutation found in both familial and sporadic ALS . The fact that ubiquilins interact also with proteins involved in spinocerebellar ataxia type 1 (UBQLN4 with ataxin-1;[9, 50]] and Huntington’s disease (UBLQN1 with huntingtin ), that UBLQN1 and UBQLN2 proteins are found in protein aggregates in Huntington disease models and human brains affected by Huntington and other neurodegenerative diseases [52–54] and the finding of UBQLN1 mutations in Brown-Vialetto-Van Laere syndrome, a rare neurological disease , further suggested important roles of these proteins in neural tissues. However, it must be emphasized that, no matter how interesting all these results are, they probably reflect just a small fraction of the range of ubiquilin functions in humans and other mammals. For example, the broad expression patterns of UBQLN1, UBQLN2 and UBQLN4 suggest that it is likely that other tissues or organs, besides the brain, may be affected by mutations in those genes. It is also significant that the testis-specific roles of UBLQN3 or the functions of UBQLNL are totally unknown. For these reasons, to determine additional roles for members of this family of proteins is a significant goal.
In this study, I analyze the patterns of diversification of ubiquilin-encoding genes in all eukaryotes, with emphasis in mammals, in which it has been found a unique expansion of this gene family. UBQLN4 turns to be the oldest among the ubiquilin genes present in humans and other vertebrates, corresponding to the ancestral gene present in many other eukaryotes. Vertebrate species often have quite different numbers of ubiquilin-encoding genes, due to duplications and losses of these genes in different lineages. Several of these duplications, occurred in mammals, have generated a group of genes that are expressed only in testis. Functional data suggest a postmeiotic role, in spermiogenesis. These results are the first comprehensive analysis of the ubiquilin gene family available and solve the main questions regarding the origin and diversification of this family in eukaryotes. In addition, they suggest significant new views of ubiquilin functional roles.
Global patterns of ubiquilin family evolution
Comprehensive searches, summarized in the Methods section, determined the presence in the databases of 643 full-length or almost complete ubiquilin sequences, all of them derived from eukaryotic species. So far already described in animals, plants and fungi [2, 4, 6], the ancient origin of this family of proteins is confirmed by the fact that they can be detected in most eukaryotic groups. They are present in both unikonts such as animals, fungi, amoebozoans, choanoflagellates or ichthyosporeans and bikonts such as plants, alveolates, stramenopiles or excavates. The presence of multiple ubiquilin genes was detected both in some plant and in some animal species. Given that organisms belonging to the sister groups of plants (green algae) and the sister group of animals (choanoflagellates) have single ubiquilin genes, it seemed likely that the expansion of the ubiquilin family in the lineages that gave rise to those plants or animals with multiple genes occurred relatively recently. This interesting question will be examined in more detail in the next Sections. A single ubiquilin gene was also found in 127 fungal species. The only fungus for which two sequences were detected was Batrachochytrium dendrobatidis (Chytridiomycota). However, while one of the sequences resembles the rest of fungal ubiquilins, the other one (accession number GL882891.1) encodes for a ubiquilin protein which is extremely different from the other fungal ubiquilins, moreover not resembling any other sequence in our database. Assuming it is not just a sequencing artifact, how this gene originated, i. e. whether simply emerged in a Batrachochytrium-specific duplication followed by drastic sequence changes or, perhaps, by horizontal transmission from another, unknown organism, is uncertain. Finally, two very different ubiquilins were detected in some excavates and in some alveolates, another result that will be described in more detail below.
An alignment of the 643 sequences was obtained (available as Additional file 1) and I performed phylogenetic analyses, either based on full-length sequences or only in the highly conserved UBL and UBA sequences. In general, the second type of analysis must be preferred when the goal is to compare very different ubiquilins. The reason is that ubiquilin sequences have a variable number of Sti1 repeats, in a way that is often unrelated to the phylogenetic proximity among species. This is an important problem if the full-length sequences are used for phylogenetic reconstruction, because the presence of exactly the same number of repeats may cause a spurious, convergent similarity among very distant sequences. However, whenever all sequences have the same structure (e. g. plants, see below), the full-length sequences can be safely used, thus increasing the amount of useful information.
Diversification of ubiquilin genes in animal species
The peculiar Drosophila results can be explained quite simply. One of the Drosophila species ubiquilin genes (called Ubqn) is typical, i. e. very similar to those found in all other insects (Figure 2). This gene has been the one examined in previous functional papers using D. melanogaster as a model [41, 42, 56]. The second Drosophila gene (named CG31528) clearly corresponds to a highly divergent but recently emerged, genus-specific duplicate. According to data compiled in FlyBase (http://www.flybase.org), Ubqn has high levels of expression in multiple tissues, while CG31528 is highly expressed only in male testis.
A more complex situation is detected in vertebrates. Results indicate that most vertebrate sequences fit into five main classes with good bootstrap support. This is shown in Figure 3, which is simply the section of the tree shown in Figure 2 that corresponds to the vertebrate sequences, expanded. From top to bottom, the first class, indicated as “UBQLN1/2 genes” in Figure 3, includes two human genes, UBQLN1 and UBQLN2, and their orthologs, which are present in all other vertebrates except actinopterygians. The second corresponds to the set of UBQLN4 orthologs, which appear in the tree as close relatives of the UBQLN1/2 genes. The third one surprisingly corresponds to a hitherto undescribed ubiquilin gene. Probably, the reason for not having been detected as such before is that it is present in many mammals, but not in humans. I have called that gene UBQLN5. Finally, the fourth and fifth correspond to the two remnant known ubiquilin genes, UBQLNL and UBQLN3. Only a few reptilian and bird sequences are so highly divergent that do not fit well in any of those five main classes. They may correspond to a sauropsid-specific duplicate, conserved in just a handful of the species for which data are currently available (bottom of Figure 3). That kind of genes can be named UBQLN6.
As indicated in Figure 3, not all genes are detected in all vertebrates. On the contrary, only a single gene, corresponding to the UBQLN4 class, was detected in all main types of vertebrates, including several actinopterygian fish species, such as Danio rerio or Salmo salar, which have this single ubiquilin gene. When we examine a closer relative of humans, the sarcopterygian Latimeria chalumnae (coelacanth), two genes can already be detected; one of them a typical UBQLN4 and the second one belonging to the UBQLN1/2 class. This situation, with two genes, is also found in amphibians, (e. g. Xenopus), and birds (e. g. Gallus). In mammals, additional gene amplifications are observed, with one exception, namely the monotreme Ornithorhynchus anatinus. Only two bona fide ubiquilin genes were detected in O. anatinus, namely a UBQLN1/2 gene and a typical UBQLNL gene. These results indicate when UBQLNL may have originated after the split which separated the mammalian ancestors from the rest of vertebrates. They also suggest that UBQLN4 (which is present in species of all vertebrate groups) have been lost in O. anatinus. In marsupials, such as Sarcophilus harrisii or Monodelphis domestica, five ubiquilin genes (UBQLN1/2, UBQLN3, UBQLN4, UBQLNL and UBQLN5) were detected. On the other hand, most eutherians have six: UBLQN1, UBQLN2, UBQLN3, UBQLN4, UBQLNL and UBQLN5. Given that the presence of the two different genes, UBQLN1 and UBQLN2, is restricted to this lineage, it means that they derive from a recent, eutherian-specific duplication of the precursor UBQLN1/2 gene present in other vertebrates. Finally, as I have already indicated, UBQLN5 -- which most likely emerged after the split that separated monotremes from the rest of mammals, given its presence in both marsupials and eutherians -- has been lost in some primates. More specifically, UBQLN5 is found in the genomes of prosimians, such as Otolemur garnettii or Microcebus murinus. However, it has not been detected either in platyrrhines or in catarrhines, including our own species. All these results, put together, indicate that vertebrates have increased their number of ubiquilin genes from a single original one (which would correspond to UBQLN4) to up to 6 genes, as found today in many mammals.
Analyses of the genomic locations of these genes in multiple organisms were performed at the Ensembl web page (see Methods) and provided significant complementary information to understand their diversification in vertebrates. I first analyzed the mouse genome, finding the significant result, confirmed later in other species, that UBQLNL, UBQLN3 and UBQLN5 are located in tandem. This indicates that these three genes are evolutionary closely related, being the most likely that UBQLN3 and UBQLN5, exclusive of marsupials and eutherian mammals, emerged by tandem duplications of UBQLNL, which is the only one detected also in monotremes. I will call from now on these three genes as the “UBQLNL group”. The second important result obtained is that the ortholog of one of the genes adjacent to UBQLN4 in the mouse genome, called Lamtor2, is also adjacent to the putative UBQLN4 gene of Danio rerio. This is additional evidence supporting the conclusions obtained from the phylogenetic reconstructions, indicating that all the genes that I have been hitherto calling UBQLN4 are true orthologs and that the first UBQLN4 gene originated before the split of actinopterygians and the rest of vertebrates. The third interesting result indicates that UBQLN2 originated from UBQLN1. This derives from the study of the Latimeria genome. It turns out that the two coelacanth genes, which I defined above as UBQLN4 and UBQLN1/2 are adjacent in the genome to, respectively, Lamtor2 (as expected, again confirming the ancient origin of UBQLN4) and Idnk. Given that this Idnk gene is in other mammals just adjacent to UBQLN1 (e. g. in mouse, they are both together in chromosome 13), but not adjacent to UBQLN2 (which is X-linked), we can conclude that the Latimeria gene named so far UBQLN1/2 most likely corresponds to UBQLN1, with UBQLN2 being thus an eutherian-specific duplicate. Additional confirmation is obtained from the fact that, in other species in which a single UBQLN1/2 gene is present (e. g. Gallus gallus), that gene is also adjacent to the Idnk ortholog.
A final type of information which is relevant here concerns the protein domain structure of animal ubiquilins. As indicated in the Introduction, in addition to the terminal UBL and UBA domains, ubiquilin typically have one to a few Sti1 domains. I have explored the structures of all these proteins using the integrated tool InterProScan (see Methods). The conclusion is that UBQLN1, UBQLN2 and UBQLN4 proteins are very similar, typically having 4 Sti1 domains (although less than four are detected in some cases), while UBQLNL proteins generally have 2 Sti1 domains and UBQLN3 and UBQLN5 usually a single one, being the Sti1 domains in UBQLN3 proteins particularly divergent. Examining then invertebrate animal proteins, it was observed that almost all have 4 Sti1 domains, being thus structurally more similar to UBQLN1, UBQLN2 and UBQLN4 than to the proteins of the UBQLNL group. Given that we have deduced that UBQLN4 is the oldest vertebrate gene, this coincidence is not surprising, and provides additional evidence supporting the monophyly of all animal ubiquilins.
Evolution of ubiquilins in other organisms
In this section, I will first describe the results for the Viridiplantae, excluding the chlorophytes, which have highly divergent ubiquilins (Figure 1). The only relevant information obtained from chlorophyte sequences is that a single gene was found in seven different species and that the proteins encoded by those genes have 3 or 4 Sti1 repeats.
The simplest explanation for the results obtained is that a single ubiquilin gene was present when viridiplantae originated. After that, a few independent duplications have occurred in many lineages. Considering that the evolutionary history of many of these plants has involved multiple rounds of genome duplication, some of them ancient, it is significant that only two old duplications can be deduced from Figure 5. One of them is, which is highly supported by bootstrap analyses, is observed in the poaceae (see “Poaceae I” and “Poaceae II” in Figure 5). The second one is a duplication in dicots, produced before the splits that separated asterids (see Lactuca species), rosids (e. g. Arabidopsis) and saxifragales (see the two Paeonia genes). This second duplication, indicated as “Dicot I” and “Dicot II” in Figure 5, has a more limited bootstrap support, but it is the simplest way to reconcile the observed tree with the known phylogenetic relationships among all these dicot species. All the other increases in the number of ubiquilin genes that have been observed in plants can be explained by independent, very recent duplications. Notice for example the two almost identical ubiquilins found in Selaginella or the three similar ubiquilins found in Physcomitrella. The same occurs in many other species, such as Arabidopsis thaliana, in which there are two very similar ubiquilin genes (called Dsk2A/At2g17190 and Dsk2B/At2g17200), which are located in tandem. The species with four ubiquilin genes (Glycine max, Brassica rapa, Zea mays) obtained that number also due to recent duplications, not observed in close relatives (as can be easily deduced from Figure 5).
Compared with the relatively complex evolutionary patterns described in vertebrate animals and in plants, the rest is much simpler, given that multiplications of the genes of this family are not detected in any other organism for which data are available. First, a single ubiquilin gene has been detected in all the remnant opisthokonts analyzed. This includes the fungal species, a choanoflagellate (Monosiga brevicollis) and an icthyosporean (Capsaspora owczarzaki). Single genes were also found in 5 amoebozoan species (from the genera Entamoeba, Dictyostelium and Polysphondylium). Also, a single gene has been detected in all stramenopile species characterized so far. Alveolates for which data are available have 1 or 2 ubiquilins. In species with two genes (e. g. which belong to the apicomplexan genera Plasmodium, Cryptosporidium or Neospora), it is clear that they are highly divergent, appearing in two distant groups in the general trees (see details in Additional file 2). This result suggests that they may derive from ancient duplications. Also, 1–2 genes are detected in excavate species, duplicates being detected in Trypanosoma and Leishmania species.
Functional data for mammalian ubiquilins
Expression patterns of mouse and human ubiquilin genes, in arbitrary units
Bronchial Epithelial Cells
CD56+ NK Cells
Mast cells IgE + antigen 1hr
Granulocytes mac1 + gr1+
Mast cells IgE
Average ± s. e. m.
2651.98 ± 123.15
Average ± s. e. m.
57.58 ± 1.24
Cerebral cortex prefrontal
Thymocyte SP CD4+
B-cells marginal zone
Dorsal Root Ganglion
Granulocytes mac1 + gr1+
Testis Leydig Cell
Testis Seminiferous Tubule
Average ± s. e. m.
3071.72 ± 336.37
Average ± s. e. m.
78.22 ± 10.35
Embryonic stem line Bruce4 p13
Cerebral cortex prefrontal
Superior Cervical Ganglion
Embryonic stem line V26 2 p16
Leukemia lymphoblastic (MOLT-4)
Granulocytes mac1 + gr1+
C3H 10T1 2
Macrophage bone marrow 6hr LPS
Average ± s. e. m.
309.12 ± 19.63
Average ± s. e. m.
7.39 ± 0.18
Cerebral cortex prefrontal
Testis Germ Cell
Testis Seminiferous Tubule
Testis Leydig Cell
Average ± s. e. m.
87.92 ± 83.09
Average ± s. e. m.
23.80 ± 9.76
UBQLN5 Not present in humans
Embryonic stem line Bruce4 p13
Average ± s. e. m.
11.93 ± 7.24
Granulocytes mac1 + gr1+
Average ± s. e. m.
26.14 ± 21.10
Average ± s. e. m.
7.56 ± 0.15
The results shown in Table 1 and Figure 6 complement what was known about ubiquilin expression in mouse, adding some interesting new information. First, these results confirm Northern blot expression data for mouse UBQLN1-4[4, 8, 12] which indicated that the genes of the UBQLN4 group, UBQLN1, UBQLN2 and UBQLN4, are expressed in multiple tissues and UBQLN3 is testis-specific. Second, it was found that UBQLNL and UBLQN5 are also testis-specific (Table 1). Thus, the three UBQLNL group genes are not only evolutionarily but also functionally related. Third, it turned out that the lowest level of expression of both UBQLN1 and UBQLN2 among all samples analyzed actually corresponds also to the male testis (Table 1). This is not the case for UBQLN4: “testis” does not appear in the section of Table 1 corresponding to that gene because the level of expression in that organ (263.24) was very similar to the average level for the whole set of samples (309.12 ± 19.63; Table 1). Actually, the low standard error of the mean of UBQLN4 values indicates a similar expression in all tissues, including testis.
Additional useful information can be obtained from Figure 6. By adding together in a single column the values of expression for all ubiquilin genes in a given tissue, some patterns become evident. Dominant in Figure 6 are the orange and brown segments, which respectively correspond to UBQLN1 and UBQLN2. These two have by far the highest levels of expression among all ubiquilin genes in most tissues. The consistent but quite low expression of UBQLN4 is, by comparison, dwarfed. It is also easily noticeable in Figure 6 how radically different from the rest of samples is the one corresponding to the testis, in which the UBQLNL group genes account for most of the expression detected. Among the genes of the UBQLN4 group, only UBQLN1 has a relatively high level of expression in testis. Finally, it can be also appreciated in that figure that many among the samples with the highest total levels of expression, obtained adding together all ubiquilin genes, come from the nervous system (see e. g. hypothalamus, cerebral cortex, cerebellum, etc.) in good agreement with previous data [6, 11].
Considering now human ubiquilins, it is important first to notice that the available information is a priori somewhat less convincing than the data for mouse genes, given that the levels of expression observed for all genes are much lower and therefore are closer to background levels (Table 1, right panels). Even with that caveat, the fact is that results very similar to those found in mouse are observed for UBQLN2 (i. e. a broad pattern of expression, with high level in nervous system samples and low levels in testis) and UBLQN3 (testis-specific expression, which is confirmed also by independent results ). The highest values for UBQLNL are also found in testis, although here the specificity is not as high as in mouse. Finally, some of the results for UBQLN1 are UBQLN4 are somewhat incongruent with those found in mouse. On one hand, UBQLN1 is broadly expressed, but no particularly low expression in testis in detected (this is clearly seen in Figure 7). This is probably a real result, given that a relatively high level of expression in testis was observed before . On the other hand, expression in whole brain samples for UBQLN4 appears as one of the lowest. However, the values of UBLQN4 are, as those for UBLQNL just mentioned, so low in all samples that is unclear to what point they are reliable. Actually, other experiments showed a relatively high level of expression of UBLQN4 in brain . In any case, even with these differences, both the obvious uniqueness of the patterns of expression observed in the testis and the fact that several nervous system samples are among the ones with the highest levels of expression, are two general results detected in both mouse and human (Figures 6 and 7).
Test for positive selection acting on human ubiquilin genes
Significant positive selection acting on particular codons in the whole set of sequences was not detected. Analyses implementing the six main site models implemented in PAML (called M0, M1a, M2a, M3, M7 and M8; see [63, 64]) failed to show any result compatible with a positive selection regime. In particular, the critical comparisons  involving either the M1a and M2a models or the M7 and the M8 models were non-significant (not shown). On the other hand, branch models, which aim to detect positive selection acting on particular genes or gene lineages, showed more interesting results. When the simplest M0 model, which considers a single ω for the whole dataset -- being ω = dN/dS, i. e., the ratio of the nonsynonymous (dN) to the synonymous (dS) nucleotide substitution rates per codon -- was compared with a “free-ratio” model in which each branch of the tree is allowed to have a different ω value, it was found that the latter significantly improved over the first one (2 Δl = 50.53; degrees of freedom = 16; p = 0.00002; see Methods for these parameters). In the free-ratio model, ω values above 1, implying positive selection, were observed in four particular branches (labelled 1–4 in Figure 10). It is often convenient to test whether the free-ratio model, which is parameter rich, can be simplified, using models in which ω values are allowed to vary only in a few particular branches. Here, two of those simpler models were tested. The first was a “five-ratio” model in which the four potentially interesting branches detected in the free-ratio model were allowed to have their own ω values, while a fifth identical value was assigned to the rest of branches. Although this five-ratio model significantly improved over M0, the free-ratio model was still however better (2 Δl = 22.73; degrees of freedom = 12; p = 0.029). However, it was observed that only two branches (labelled 3 and 4 in Figure 4) showed values ω > 1 in the five-ratio model, which led to the idea of testing a third, “three-ratio” model, in which only three ω values were allowed, one for branch 3, another for branch 4 and a third for the rest of branches of the tree. Again, this three-ratio model significantly improved over the M0 model, but was worse than the free-ratio model, albeit with a difference that was very close to the significance level (2 Δl = 24.48; degrees of freedom = 14; p = 0.041). The ω = 999 values found in some branches in the five-ratio and three-ratio models imply a very low/zero dS value. Actually, dS values were 0.0001 in branch 4 of both the five-ratio and the three-ratio models and 0 in branch 3 of the three-ratio model. These very low dS values however do not affect the likelihood comparisons, and therefore, the conclusion that the free-ratio model is the best one stands.
If we accept the free-ratio model as the best, this means that there is evidence for positive selection acting in four cases: 1) after the duplication that gave rise to the ancestor of the UBQLNL group genes; 2) after the UBQLN1/UBQLN2 duplication; 3) in the UBQLN3 genes after the rat/human split and specifically in the lineage that gave rise to humans; and, 4) in the UBQLNL genes, also in the human lineage. However, given that the improvement of the free-ratio model over the three-ratio model is not highly significant, it cannot be disregarded at present that positive selection may have acted solely on UBQLN3 and UBQLNL, i. e. along branches 3 and 4 of Figure 10.
I finally performed branch-site models  to determine whether any codon could be detected to be under positive selection in branches 1, 2, 3 or 4 of Figure 10. Each branch was examined in an independent analysis. However, perhaps not unexpectedly, no significant results were found. The most likely reason for those negative results is that branch-site models only allow testing for positively selected codons using two different ω values, one for a given branch and a second for the rest of the tree. However, the best model observed in branch models (free-ratio model) implies a different ω value in each branch and it is thus likely that the limitation of using just two of those values for the whole tree precludes the determination of the codons under positive selection. All PAML results can be found in Additional file 3.
This study establishes for the first time in the literature the main patterns of the evolution of the ubiquilin gene family. All the data obtained are compatible with ubiquilins emerging very early in eukaryotic evolution and transmitting strictly vertically. The only potential exception found concerns the divergent ubiquilin detected in the chytridiomycote Batrachochytrium dendrobatidis, whose origin is uncertain. Very few ubiquilin genes, most likely just a single one, were present before the divergence of the main eukaryotic groups. Also, single genes were most likely present in the ancestors of animals, plants (including green algae) and fungi. Ubiquilin evolution has been in general very conservative, given that most eukaryotes have just 1 – 2 ubiquilins. The only exceptions are vertebrates and green plants, in which up to 4 – 6 ubiquilins are detected in some species. As already indicated, most plant duplications are relatively recent, often genus-specific. Just a single ancient duplication must be postulated to explain the data obtained for some monocots (poaceae) and a second one for dicots. This explains why many viridiplantae species have only 1 or 2 ubiquilin genes. Given that many plant lineages have suffered whole genome duplications, this means that ubiquilin genes are relatively “resistant” to be duplicated in plants, i. e. many of the duplicates have been lost. This is reminiscent of what is found in some families of plant ubiquitin ligases [66, 67]. At present, the reason for the expansion of the ubiquilin gene family in some plant lineages is totally unknown.
Respect to all the other groups, the amplification of the ubiquilin gene family in vertebrates, and more specifically in mammals is clearly unique. Mammals have 5 – 6 ubiquilin genes, while most animal species (i. e. with the only exception of Drosophila, all non-vertebrates, including chordates, and also actinopterygian fish) have just one. This multiplication has occurred recently. In particular, four genes are not present in sauropsids and three of them, UBQLN2, UBQLN3 and UBQLN5, are apparently absent in monotremes, meaning that they may have emerged in the last 150 millions of years. Losses of ubiquilin genes in vertebrates have occurred, but only rarely (Figure 4).
Genes of the UBQLN4 group (UBQLN1, UBQLN2 and UBQLN4) have retained a general pattern of expression which must be similar to that of the single ubiquilin gene in other animals, which is most likely expressed in all tissues (e. g. the Drosophila Ubqn results mentioned above). There is good evidence for the products of UBQLN4 group genes having roles as ubiquitin receptors either to lead to proteasome degradation of ubiquitinated proteins  or to redirect ubiquitinated proteins to the autophagy pathway [68, 69]. This second role is not apparently performed by the Saccharomyces cerevisiae single ubiquilin, DSK2 . An additional facet of the role of these proteins is related to the finding that UBQLN4 can interact with both UBQLN1 and UBQLN2, and that this interaction redirects ubiquilin-interacting proteins towards the autophagy pathway . This last result suggests that the roles in autophagy may depend on the presence of multiple different ubiquilins, and thus would have appeared only after the duplications that generated UBQLN1 and UBQLN2 from UBQLN4. This interesting functional hypothesis could be quite simply tested, by determining whether ubiquilins have roles in autophagy in actinopterygians fishes, in which only UBQLN4 is present. Other important points that deserve further study are: 1) why UBQLN4 has lower levels of expression in all tissues than the other two genes (Figures 6 and 7); 2) the cause of the increased level of expression observed for these genes in particular tissues, and especially in the nervous system (see also Figures 6 and 7). This is a question that may provide strong clues about their relationship with neurodegenerative diseases; 3) whether different UBQLN4 group proteins have different affinities for different types of ubiquitinated proteins or even for particular types of ubiquitin chains, and, 4) whether there has been positive selection on the UBQLN2 gene, as suggested by the “free-ratio” model results (see above). It is noteworthy the finding that yeast or plant ubiquilins show preferential binding to Lys-48 chains  while mammalian UBQLN1 binds both Lys-48 and Lys-63 chains with quite similar affinities [71, 72]. Whether this difference is related to the presence in mammals of multiple related proteins of the UBQLN4 group, each with its own biochemical properties, is unknown.
The testis-specific roles of the second trio, the UBQLNL group genes UBQLN3, UBQLN5 and UBLQNL, have been explored here in detail. A significant finding is that their main roles seem to occur in spermatids, postmeiotic germ cells (Figures 8 and 9). This may be used as a clue to understand the origin of the UBQLNL group genes, which cannot be deduced from the phylogenetic analyses, given that none of the UBQLN4 group genes is particularly similar to any of the UBQLNL group genes (Figure 3). I think that it is significant that UBQLN1 is consistently expressed in testis (Figures 6 and 7) while the levels of expression of UBQLN2 and UBQLN4 are much lower. Also, it has been shown that UBQLN1 has specific roles in postmeiotic germ cells, colocalizing with the manchette, a structure made of actin and microtubules that is present in elongating spermatids . These results suggest that the oldest UBQLNL group gene, UBQLNL, may have originated as a duplicate of UBQLN1.
It is likely that the production of testis-specific tandem duplicates is linked to a strong selective pressure to increase ubiquilin gene expression in that organ. In this context, the fact that Drosophila CG31528, the only duplicate found in invertebrates (i. e. with an origin totally independent from that of the UBQLNL group genes) is also testis-specific seems to be more than a coincidental finding. It is also significant that some evidence for positive selection acting after the first UBQLNL group gene originated but before the UBQLN3/UBQLNL duplication occurred has been obtained (see “free-ratio” model results, above). The finding of positive selection acting on the UBQLNL and UBQLN3 genes after the split that separated rodents from primates, specifically in the human ancestors is quite robust (Figure 10), and also deserves further study.
Two alternative hypothesis for the origin of new testis-specific ubiquilins are either that UBQLNL group genes have acquired totally novel roles in that organ (neofunctionalization) or that those genes have roles in the testis that, before their emergence, were performed by other ubiquilins, for example UBQLN1 (subfunctionalization). At present, there is no way to determine which of these options is correct, the precise roles of testis-specific ubiquilins being unknown. However, the fact that UBQLNL group proteins are structurally different from their ancestors, having lost several Sti1 domains, suggests that new functions may have arisen. In any case, this need for an increase of ubiquilins in the testis must be related to the important specific roles of the ubiquitination system in the male gonad, for which there is a growing body of evidence (see reviews: [73–75]). A recent work has pointed out that perhaps 20 % of all ubiquitin ligases, which are the proteins that provide specificity to the ubiquitination machinery, may be expressed at much higher levels or even totally specifically in the mouse testis . The reasons for this particular need for a complex, testis-specific ubiquitination machinery are unknown. A hypothesis that fits well with the data obtained in this study is that this specificity may be linked to roles that are performed in male germ cells, but not in somatic cells or female germ cells. In this context, postmeiotic patterns of expression as those found here for testis-specific ubiquilins are coincidental with those of genes critical for the substitution of histones by protamines, including transition proteins and protamines themselves , suggesting that they might be linked, directly or indirectly, to this unique need for extreme chromatin compactation. However, other spermiogenesis-specific processes (alteration of cell shape, generation of the acrosome and flagellum, etc.; see ) are also candidates for involving specific roles of the ubiquitin-proteasome in which ubiquilins might be required.
Additional evolutionary considerations may contribute to discriminate among these options. It is significant that my own analyses of the data obtained by Hou et al. indicate that those testis-specific ubiquitin ligases for which phylogenetic data is available (e. g. testis-specifically expressed members of the RBR and TRIM families that I have studied before [79, 80]) originated at very different times along eukaryotic evolution (unpublished results). Therefore, the testis-specific roles of the ubiquitination system proteins that we found now in mouse or human may have arisen at very different times, making unlikely a simple, all-encompassing explanation for their functions. This means that we should look for particular, gene-specific explanations, and then it becomes relevant exactly when each particular gene originated. Thus, the fact that group UBLQNL genes are mammalian-specific, and two of them are found only in eutherians, may provide important clues about their roles. If we assume that their roles are linked to processes that only happen in eutherians, we are left with very few options, because most processes that occur in the mammalian male germ line (e. g. the histone to protamine transition indicated above) have ancient origins, being common to all vertebrates or even to all animals. Actually, the available literature indicates that the main functional differences between the male germ cells of eutherian mammals and those of reptiles, birds or even monotremes are related to the difficulty to fecundate the heavily protected eutherian oocytes (reviewed in [81, 82]). Eutherian sperm has acquired a series of specific adaptations to achieve fecundation, involving redistribution of products through complex intracellular transport processes that lead to changes both in the physiology and the shape of the cells [83, 84]. These specific adaptations are good candidates to require proteins with novel roles. Given that the product of UBQLN1 – probably the ancestral gene from which the testis-specific ubiquilins emerged and also the only one among the UBQLN4 group genes consistently expressed in the testis -- colocalizes with the manchette , a structure linked to intracellular transport in these germ cells, an attractive hypothesis is that the emergence of the UBQLNL group proteins is related to new roles of the ubiquitination system linked to these eutherian-specific sperm adaptations. This hypothesis can be tested by generating loss-of-function mutants in the testis-specific ubiquilin genes.
The ubiquilin gene family is present in all eukaryotes and has a very conservative pattern of evolution, with many eukaryotic species having a single gene. The lineage-specific amplifications observed, among which the one detected in mammals is the most extensive, are probably linked to the acquisition of specific, potentially novel functions by the newly-emerged duplicates. In mammals, three recently arisen ubiquilins are required specifically in the testis and this is also the case for a Drosophila novel ubiquilin gene. This suggests that duplications leading to the generation of genes with testis-specific roles may occasionally provide a selective advantage in animal lineages. Potential roles in intracellular transport for the mammalian testis-specific proteins are suggested by the available data.
Ubiquilin sequences were obtained by TblastN searches against the nr, est, htgs, gss, wgs and tsa databases at the National Center for Biotechnology Information (NCBI; http://www.ncbi.nlm.nih.gov/). First, general searches and then specific searches to detect ubiquilins of particular groups (animal, plants, fungi, etc.) or critical model species which could have been missed in the original searches were performed. The query sequences for those searches were selected from all main eukaryotic groups in which ubiquitins were detected (animals, plants, fungi, amoebozoans, alveolates and flagellates). Given the high similarity among all ubiquilins, those searches soon became saturated, with no additional sequences being detectable when additional query sequences were used.
From those searches, and after eliminating duplicates and partial sequences, a total of 643 sequences were found to be complete or almost complete (i. e. only missing a few amino acids, generally at the N or C termini). A total of 619 of them had full-length UBA and UBL domains. These sequences were aligned using MAFFT v6.864b  and the alignments were manually corrected editing the sequences with GeneDoc 2.7 . Phylogenetic analyses were performed following similar procedures to those in several studies focused on ubiquitination system genes ([87, 88] and related references along the text). Three different methods were used, namely Neighbor-joining (NJ), maximum-parsimony (MP) and maximum-likelihood (ML). NJ and ML trees were obtained using MEGA 5.1  and Maximum-parsimony (MP) trees were obtained using PAUP* 4.0, beta 10 version . For NJ, Kimura´s correction was used and sites with gaps were treated with the pairwise deletion option. Parameters for MP analyses based on full-length sequences were as follows: 1) all sites included, gaps treated as unknown characters; 2) randomly generated trees used as seeds; 3) maximum number of trees saved equal to 100; and, 4) heuristic search using the tree-bisection-reconnection (TBR) algorithm, with default parameters. The same methods were used for analyses of alignments of the UBA and UBL domains, except that the faster subtree-pruning-regrafting (SPR) algorithm, also with default parameters, was used instead of the TBR algorithm. The reason for this methodological change is that the analyses based on UBA and UBL domains included a very large number of sequences, making TBR-based analyses unfeasible. Finally, for ML analyses, the BioNJ tree was taken as starting point for the iterative searches using the Jones-Taylor-Thornton (JTT) model of amino acidic substitutions. A discrete Gamma (G) distribution with five categories of sites was estimated, to account for heterogeneity in evolutionary rates. This JTT + G model was chosen because it was the best, according to the ML model comparison analyses available in MEGA 5.1. Gaps were also treated as unknown characters. Here, for the same reason indicated above, while the TBR routine with 3 levels of tree interchange was used to explore the landscape of ML trees in analyses involving full-length sequences, the faster SPR algorithm, also with 3 levels of subtree interchange, was used in the analyses involving just the UBA and UBL domains. Bootstrap tests were performed to establish the reliability of the final dendrograms obtained in the NJ, MP and ML analyses. A total of 1000 replicates were generated for NJ analyses and 100 replicates were made for the MP and ML analyses, which are more computer-intensive. MEGA 5.1 was also used to draw and edit the trees in Figures 1, 2, 3, 4 and 5.
The origin of the genes and the patterns of duplications and losses were determined by reconciling the gene trees with the species trees and, when needed, integrating additional information such as the relative location of the genes in several genomes, which ones are the genes located adjacent to those encoding ubiquilins or the ubiquilin protein structures (see Results for the details). Analyses of the genomic locations of ubiquilin genes were performed at the Ensembl genome browser web page  (http://www.ensembl.org/). Structural analyses of the ubiquilin proteins were performed using the integrated tool InterProScan  (available online at http://www.ebi.ac.uk/Tools/pfa/iprscan/).
Microarray data were obtained from the public repositories in which the datasets from the studies cited along the text were deposited. These were either the Gene Expression Omnibus (GEO) database at the NCBI (http://www.ncbi.nlm.nih.gov/geo/), the ArrayExpress database (http://www.ebi.ac.uk/arrayexpress/) or the BioGPS database (http://www.biogps.org). For some samples and genes, multiple probes or experiments were available. In those cases, I chose the ones with the highest average level of expression. The raw data of the experiments has been used in this study. Given that the quantitative values of expression of different experiments were not compared, no further normalization or other kinds of data manipulation were required.
Analyses to detect positive selection in human ubiquilins were performed using PAMLX , a graphical interface for the PAML program . Methods were very similar to those described in one of my previous papers . In brief, I took the protein sequences of the full-length UBA and UBL domains (a total of 117 amino acids) of the five ubiquilin genes present in both Homo sapiens and Rattus norvegicus (UBQLN1, UBQLN2, UBQLN3, UBQLN4 and UBQLNL) and obtained the corresponding nucleotide sequences for those regions of the genes in the two species. Trees were obtained with the ten sequences that generated the expected topology (Figure 10). Then, the codon-substitution models implemented in the CODEML program of PAML  were used to estimate the synonymous (dS) and non-synonymous (dN) rates of evolution, in order to determine: 1) whether positive selection (ω = dN/dS > 1) was detectable at some codons in the whole sets of sequences (“site models”); 2) whether positive selection was detected in particular branches of the sequences tree (“branch models”); and, 3) whether positive selection was detectable at some codons in particular branches of the trees (“branch-site models”) [63, 64]. The CODEML analyses provide the log likelihood value of a given model of codon substitution, for the sequences considered and evaluating their evolutionary relationships. All the comparisons made here among alternative models corresponded to pairs of nested models. Thus, their results can be compared by using the LRT statistic. This statistic equals 2 Δl, being Δl the difference between the log likelihoods of the two models. The LRT statistic follows a chi-square distribution with a number of degrees of freedom equal to the difference in the number of parameters used in each of the two models that are compared (see details in [63, 64] and references therein).
This study was supported by grant BFU2011-30063 (Spanish government). The funding body did not have any role in the design, analysis, or interpretation of data or in the writing of the manuscript and the decision to submit the manuscript for publication.
- Vallen EA, Ho W, Winey M, Rose MD: Genetic interactions between CDC31 and KAR1, two genes required for duplication of the microtubule organizing center in Saccharomyces cerevisiae. Genetics. 1994, 137: 407-422.PubMedPubMed CentralGoogle Scholar
- Biggins S, Ivanovska I, Rose MD: Yeast ubiquitin-like genes are involved in duplication of the microtubule organizing center. J Cell Biol. 1996, 133: 1331-1346. 10.1083/jcb.133.6.1331.PubMedView ArticleGoogle Scholar
- Funakoshi M, Geley S, Hunt T, Nishimoto T, Kobayashi H: Identification of XDRP1; a Xenopus protein related to yeast Dsk2p binds to the N-terminus of cyclin A and inhibits its degradation. EMBO J. 1999, 18: 5009-5018. 10.1093/emboj/18.18.5009.PubMedPubMed CentralView ArticleGoogle Scholar
- Wu AL, Wang J, Zheleznyak A, Brown EJ: Ubiquitin-related proteins regulate interaction of vimentin intermediate filaments with the plasma membrane. Mol Cell. 1999, 4: 619-625. 10.1016/S1097-2765(00)80212-9.PubMedView ArticleGoogle Scholar
- Kaye FJ, Modi S, Ivanovska I, Koonin EV, Thress K, Kubo A, Kornbluth S, Rose MD: A family of ubiquitin-like proteins binds the ATPase domain of Hsp70-like Stch. FEBS Lett. 2000, 467: 348-355. 10.1016/S0014-5793(00)01135-2.PubMedView ArticleGoogle Scholar
- Mah AL, Perry G, Smith MA, Monteiro MJ: Identification of ubiquilin, a novel presenilin interactor that increases presenilin protein accumulation. J Cell Biol. 2000, 151: 847-862. 10.1083/jcb.151.4.847.PubMedPubMed CentralView ArticleGoogle Scholar
- Kleijnen MF, Shih AH, Zhou P, Kumar S, Soccio RE, Kedersha NL, Gill G, Howley PM: The hPLIC proteins may provide a link between the ubiquitination machinery and the proteasome. Mol Cell. 2000, 6: 409-419. 10.1016/S1097-2765(00)00040-X.PubMedView ArticleGoogle Scholar
- Conklin D, Holderman S, Whitmore TE, Maurer M, Feldhaus AL: Molecular cloning, chromosome mapping and characterization of UBQLN3 a testis-specific gene that contains an ubiquitin-like domain. Gene. 2000, 249: 91-98. 10.1016/S0378-1119(00)00122-0.PubMedView ArticleGoogle Scholar
- Davidson JD, Riley B, Burright EN, Duvick LA, Zoghbi HY, Orr HT: Identification and characterization of an ataxin-1-interacting protein: A1Up, a ubiquitin-like nuclear protein. Hum Mol Genet. 2000, 9: 2305-2312. 10.1093/oxfordjournals.hmg.a018922.PubMedView ArticleGoogle Scholar
- Ozaki T, Hishiki T, Toyama Y, Yuasa S, Nakagawara A, Sakiyama S: Identification of a new cellular protein that can interact specifically with DAN. DNA Cell Biol. 1997, 16: 985-991. 10.1089/dna.1997.16.985.PubMedView ArticleGoogle Scholar
- Hanaoka E, Ozaki T, Ohira M, Nakamura Y, Suzuki M, Takahashi E, Moriya H, Nakagawara A, Sakiyama S: Molecular cloning and expression analysis of the human DA41 gene and its mapping to chromosome 9q21.2-q21.3. J Hum Genet. 2000, 45: 188-191. 10.1007/s100380050209.PubMedView ArticleGoogle Scholar
- Matsuda M, Koide T, Yorihuzi T, Hosokawa N, Nagata K: Molecular cloning of a novel ubiquitin-like protein, UBIN, that binds to ER targeting signal sequences. Biochem Biophys Res Commun. 2001, 280: 535-540. 10.1006/bbrc.2000.4149.PubMedView ArticleGoogle Scholar
- Hayes MG, Pluzhnikov A, Miyake K, Sun Y, Ng MC, Roe CA, Below JE, Nicolae RI, Konkashbaev A, Bell GI, Cox NJ, Hanis CL: Identification of type 2 diabetes genes in Mexican Americans through genome-wide association studies. Diabetes. 2007, 56: 3033-3044. 10.2337/db07-0482.PubMedView ArticleGoogle Scholar
- Funakoshi M, Sasaki T, Nishimoto T, Kobayashi H: Budding yeast Dsk2p is a polyubiquitin-binding protein that can interact with the proteasome. Proc Natl Acad Sci U S A. 2002, 99: 745-750. 10.1073/pnas.012585199.PubMedPubMed CentralView ArticleGoogle Scholar
- Wilkinson CR, Seeger M, Hartmann-Petersen R, Stone M, Wallace M, Semple C, Gordon C: Proteins containing the UBA domain are able to bind to multi-ubiquitin chains. Nat Cell Biol. 2001, 3: 939-943. 10.1038/ncb1001-939.PubMedView ArticleGoogle Scholar
- Elsasser S, Gali RR, Schwickart M, Larsen CN, Leggett DS, Müller B, Feng MT, Tübing F, Dittmar GA, Finley D: Proteasome subunit Rpn1 binds ubiquitin-like protein domains. Nat Cell Biol. 2002, 4: 725-730. 10.1038/ncb845.PubMedView ArticleGoogle Scholar
- Heessen S, Masucci MG, Dantuma NP: The UBA2 domain functions as an intrinsic stabilization signal that protects Rad23 from proteasomal degradation. Mol Cell. 2005, 18: 225-235. 10.1016/j.molcel.2005.03.015.PubMedView ArticleGoogle Scholar
- Regan-Klapisz E, Sorokina I, Voortman J, de Keizer P, Roovers RC, Verheesen P, Urbé S, Fallon L, Fon EA, Verkleij A, Benmerah A, van Bergen en Henegouwen PM: Ubiquilin recruits Eps15 into ubiquitin-rich cytoplasmic aggregates via a UIM-UBL interaction. J Cell Sci. 2005, 118: 4437-4450. 10.1242/jcs.02571.PubMedView ArticleGoogle Scholar
- Heir R, Ablasou C, Dumontier E, Elliott M, Fagotto-Kaufmann C, Bedford FK: The UBL domain of PLIC-1 regulates aggresome formation. EMBO Rep. 2006, 7: 1252-1258. 10.1038/sj.embor.7400823.PubMedPubMed CentralView ArticleGoogle Scholar
- Heinen C, Acs K, Hoogstraten D, Dantuma NP: C-terminal UBA domains protect ubiquitin receptors by preventing initiation of protein degradation. Nat Commun. 2011, 2: 191-PubMedPubMed CentralView ArticleGoogle Scholar
- Kang Y, Zhang N, Koepp DM, Walters KJ: Ubiquitin receptor proteins hHR23a and hPLIC2 interact. J Mol Biol. 2007, 365: 1093-1101. 10.1016/j.jmb.2006.10.056.PubMedPubMed CentralView ArticleGoogle Scholar
- Di Fiore PP, Polo S, Hofmann K: When ubiquitin meets ubiquitin receptors: a signalling connection. Nat Rev Mol Cell Biol. 2003, 4: 491-497. 10.1038/nrm1124.PubMedView ArticleGoogle Scholar
- Elsasser S, Finley D: Delivery of ubiquitinated substrates to protein-unfolding machines. Nat Cell Biol. 2005, 7: 742-749. 10.1038/ncb0805-742.PubMedView ArticleGoogle Scholar
- Su V, Lau AF: Ubiquitin-like and ubiquitin-associated domain proteins: significance in proteasomal degradation. Cell Mol Life Sci. 2009, 66: 2819-2833. 10.1007/s00018-009-0048-9.PubMedPubMed CentralView ArticleGoogle Scholar
- Lee DY, Brown EJ: Ubiquilins in the crosstalk among proteolytic pathways. Biol Chem. 2012, 393: 441-447.PubMedView ArticleGoogle Scholar
- Husnjak K, Dikic I: Ubiquitin-binding proteins: decoders of ubiquitin-mediated cellular functions. Annu Rev Biochem. 2012, 81: 291-322. 10.1146/annurev-biochem-051810-094654.PubMedView ArticleGoogle Scholar
- Haapasalo A, Viswanathan J, Bertram L, Soininen H, Tanzi RE, Hiltunen M: Emerging role of Alzheimer's disease-associated ubiquilin-1 in protein aggregation. Biochem Soc Trans. 2010, 38: 150-155. 10.1042/BST0380150.PubMedView ArticleGoogle Scholar
- El Ayadi A, Stieren ES, Barral JM, Boehning D: Ubiquilin-1 and protein quality control in Alzheimer disease. Prion. 2013, 7: 164-169. 10.4161/pri.23711.PubMedPubMed CentralView ArticleGoogle Scholar
- Massey LK, Mah AL, Ford DL, Miller J, Liang J, Doong H, Monteiro MJ: Overexpression of ubiquilin decreases ubiquitination and degradation of presenilin proteins. J Alzheimers Dis. 2004, 6: 79-92.PubMedGoogle Scholar
- Bertram L, Hiltunen M, Parkinson M, Ingelsson M, Lange C, Ramasamy K, Mullin K, Menon R, Sampson AJ, Hsiao MY, Elliott KJ, Velicelebi G, Moscarillo T, Hyman BT, Wagner SL, Becker KD, Blacker D, Tanzi RE: Family-based association between Alzheimer's disease and variants in UBQLN1. N Engl J Med. 2005, 352: 884-894. 10.1056/NEJMoa042765.PubMedView ArticleGoogle Scholar
- Bensemain F, Chapuis J, Tian J, Shi J, Thaker U, Lendon C, Iwatsubo T, Amouyel P, Mann D, Lambert JC: Association study of the Ubiquilin gene with Alzheimer's disease. Neurobiol Dis. 2006, 22: 691-693. 10.1016/j.nbd.2006.01.007.PubMedView ArticleGoogle Scholar
- Brouwers N, Sleegers K, Engelborghs S, Bogaerts V, van Duijn CM, De Deyn PP, Van Broeckhoven C, Dermaut B: The UBQLN1 polymorphism, UBQ-8i, at 9q22 is not associated with Alzheimer's disease with onset before 70 years. Neurosci Lett. 2006, 392: 72-74. 10.1016/j.neulet.2005.08.064.PubMedView ArticleGoogle Scholar
- Kamboh MI, Minster RL, Feingold E, DeKosky ST: Genetic association of ubiquilin with Alzheimer's disease and related quantitative measures. Mol Psychiatry. 2006, 11: 273-279. 10.1038/sj.mp.4001775.PubMedView ArticleGoogle Scholar
- Slifer MA, Martin ER, Bronson PG, Browning-Large C, Doraiswamy PM, Welsh-Bohmer KA, Gilbert JR, Haines JL, Pericak-Vance MA: Lack of association between UBQLN1 and Alzheimer disease. Am J Med Genet B Neuropsychiatr Genet. 2006, 141B: 208-213. 10.1002/ajmg.b.30298.PubMedView ArticleGoogle Scholar
- Smemo S, Nowotny P, Hinrichs AL, Kauwe JS, Cherny S, Erickson K, Myers AJ, Kaleem M, Marlowe L, Gibson AM, Hollingworth P, O'Donovan MC, Morris CM, Holmans P, Lovestone S, Morris JC, Thal L, Li Y, Grupe A, Hardy J, Owen MJ, Williams J, Goate A: Ubiquilin 1 polymorphisms are not associated with late-onset Alzheimer's disease. Ann Neurol. 2006, 59: 21-26. 10.1002/ana.20673.PubMedView ArticleGoogle Scholar
- Arias-Vásquez A, de Lau L, Pardo L, Liu F, Feng BJ, Bertoli-Avella A, Isaacs A, Aulchenko Y, Hofman A, Oostra B, Breteler M, van Duijn C: Relationship of the Ubiquilin 1 gene with Alzheimer's and Parkinson's disease and cognitive function. Neurosci Lett. 2007, 424: 1-5. 10.1016/j.neulet.2007.07.015.PubMedView ArticleGoogle Scholar
- Golan MP, Melquist S, Safranow K, Styczyńska M, Słowik A, Kobryś M, Zekanowski C, Barcikowska M: Analysis of UBQLN1 variants in a Polish Alzheimer's disease patient: control series. Dement Geriatr Cogn Disord. 2008, 25: 366-371. 10.1159/000121006.PubMedView ArticleGoogle Scholar
- Stieren ES, El Ayadi A, Xiao Y, Siller E, Landsverk ML, Oberhauser AF, Barral JM, Boehning D: Ubiquilin-1 is a molecular chaperone for the amyloid precursor protein. J Biol Chem. 2011, 286: 35689-35698. 10.1074/jbc.M111.243147.PubMedPubMed CentralView ArticleGoogle Scholar
- Viswanathan J, Haapasalo A, Böttcher C, Miettinen R, Kurkinen KM, Lu A, Thomas A, Maynard CJ, Romano D, Hyman BT, Berezovska O, Bertram L, Soininen H, Dantuma NP, Tanzi RE, Hiltunen M: Alzheimer's disease-associated ubiquilin-1 regulates presenilin-1 accumulation and aggresome formation. Traffic. 2011, 12: 330-348. 10.1111/j.1600-0854.2010.01149.x.PubMedPubMed CentralView ArticleGoogle Scholar
- El Ayadi A, Stieren ES, Barral JM, Boehning D: Ubiquilin-1 regulates amyloid precursor protein maturation and degradation by stimulating K63-linked polyubiquitination of lysine 688. Proc Natl Acad Sci U S A. 2012, 109: 13416-13421. 10.1073/pnas.1206786109.PubMedPubMed CentralView ArticleGoogle Scholar
- Li A, Xie Z, Dong Y, McKay KM, McKee ML, Tanzi RE: Isolation and characterization of the Drosophila ubiquilin ortholog dUbqln: in vivo interaction with early-onset Alzheimer disease genes. Hum Mol Genet. 2007, 16: 2626-2639. 10.1093/hmg/ddm219.PubMedView ArticleGoogle Scholar
- Ganguly A, Feldman RM, Guo M: Ubiquilin antagonizes presenilin and promotes neurodegeneration in Drosophila. Hum Mol Genet. 2008, 17: 293-302.PubMedView ArticleGoogle Scholar
- Deng HX, Chen W, Hong ST, Boycott KM, Gorrie GH, Siddique N, Yang Y, Fecto F, Shi Y, Zhai H, Jiang H, Hirano M, Rampersaud E, Jansen GH, Donkervoort S, Bigio EH, Brooks BR, Ajroud K, Sufit RL, Haines JL, Mugnaini E, Pericak-Vance MA, Siddique T: Mutations in UBQLN2 cause dominant X-linked juvenile and adult-onset ALS and ALS/dementia. Nature. 2011, 477: 211-215. 10.1038/nature10353.PubMedPubMed CentralView ArticleGoogle Scholar
- Daoud H, Suhail H, Szuto A, Camu W, Salachas F, Meininger V, Bouchard JP, Dupré N, Dion PA, Rouleau GA: UBQLN2 mutations are rare in French and French-Canadian amyotrophic lateral sclerosis. Neurobiol Aging. 2012, 33: 2230.e1-2230.e5.View ArticleGoogle Scholar
- Synofzik M, Maetzler W, Grehl T, Prudlo J, Vom Hagen JM, Haack T, Rebassoo P, Munz M, Schöls L, Biskup S: Screening in ALS and FTD patients reveals 3 novel UBQLN2 mutations outside the PXX domain and a pure FTD phenotype. Neurobiol Aging. 2012, 33: 2949.e13-2949.e17.View ArticleGoogle Scholar
- Williams KL, Warraich ST, Yang S, Solski JA, Fernando R, Rouleau GA, Nicholson GA, Blair IP: UBQLN2/ubiquilin 2 mutation and pathology in familial amyotrophic lateral sclerosis. Neurobiol Aging. 2012, 33: 2527.e3-2527.e10.View ArticleGoogle Scholar
- Kim SH, Shi Y, Hanson KA, Williams LM, Sakasai R, Bowler MJ, Tibbetts RS: Potentiation of amyotrophic lateral sclerosis (ALS)-associated TDP-43 aggregation by the proteasome-targeting factor, ubiquilin 1. J Biol Chem. 2009, 284: 8083-8092. 10.1074/jbc.M808064200.PubMedPubMed CentralView ArticleGoogle Scholar
- Hanson KA, Kim SH, Wassarman DA, Tibbetts RS: Ubiquilin modifies TDP-43 toxicity in a Drosophila model of amyotrophic lateral sclerosis (ALS). J Biol Chem. 2010, 285: 11068-11072. 10.1074/jbc.C109.078527.PubMedPubMed CentralView ArticleGoogle Scholar
- Brettschneider J, Van Deerlin VM, Robinson JL, Kwong L, Lee EB, Ali YO, Safren N, Monteiro MJ, Toledo JB, Elman L, McCluskey L, Irwin DJ, Grossman M, Molina-Porcel L, Lee VM, Trojanowski JQ: Pattern of ubiquilin pathology in ALS and FTLD indicates presence of C9ORF72 hexanucleotide expansion. Acta Neuropathol. 2012, 123: 825-839. 10.1007/s00401-012-0970-z.PubMedPubMed CentralView ArticleGoogle Scholar
- Riley BE, Xu Y, Zoghbi HY, Orr HT: The effects of the polyglutamine repeat protein ataxin-1 on the UbL-UBA protein A1Up. J Biol Chem. 2004, 279: 42290-42301. 10.1074/jbc.M406284200.PubMedView ArticleGoogle Scholar
- Wang H, Lim PJ, Yin C, Rieckher M, Vogel BE, Monteiro MJ: Suppression of polyglutamine-induced toxicity in cell and animal models of Huntington's disease by ubiquilin. Hum Mol Genet. 2006, 15: 1025-1041. 10.1093/hmg/ddl017.PubMedView ArticleGoogle Scholar
- Doi H, Mitsui K, Kurosawa M, Machida Y, Kuroiwa Y, Nukina N: Identification of ubiquitin-interacting proteins in purified polyglutamine aggregates. FEBS Letters. 2004, 571: 171-176. 10.1016/j.febslet.2004.06.077.PubMedView ArticleGoogle Scholar
- Mori F, Tanji K, Odagiri S, Toyoshima Y, Yoshida M, Ikeda T, Sasaki H, Kakita A, Takahashi H, Wakabayashi K: Ubiquilin immunoreactivity in cytoplasmic and nuclear inclusions in synucleinopathies, polyglutamine diseases and intranuclear inclusión body disease. Acta Neuropathol. 2012, 124: 149-151. 10.1007/s00401-012-0999-z.PubMedView ArticleGoogle Scholar
- Rutherford NJ, Lewis J, Clippinger AK, Thomas MA, Adamson J, Cruz PE, Cannon A, Xu G, Golde TE, Shaw G, Borchelt DR, Giasson BI: Unbiased screen reveals ubiquilin-a and −2 highly associated with huntingtin inclusions. Brain Res. 2013, 1524: 62-73.PubMedPubMed CentralView ArticleGoogle Scholar
- González-Pérez P, Lu Y, Chian RJ, Sapp PC, Tanzi RE, Bertram L, McKenna-Yasek D, Gao FB, Brown RH: Association of UBQLN1 mutation with Brown-Vialetto-Van Laere syndrome but not typical ALS. Neurobiol Dis. 2012, 48: 391-398. 10.1016/j.nbd.2012.06.018.PubMedPubMed CentralView ArticleGoogle Scholar
- Lipinszki Z, Pál M, Nagy O, Deák P, Hunyadi-Gulyas E, Udvardy A: Overexpression of Dsk2/dUbqln results in severe developmental defects and lethality in Drosophila melanogaster that can be rescued by overexpression of the p54/Rpn10/S5a proteasomal subunit. FEBS J. 2011, 278: 4833-4844. 10.1111/j.1742-4658.2011.08383.x.PubMedView ArticleGoogle Scholar
- Wu C, Macleod I, Su AI: BioGPS and MyGene.info: organizing online, gene-centric information. Nucleic Acids Res. 2013, 41 (Database issue): D561-D565.PubMedPubMed CentralView ArticleGoogle Scholar
- Schultz N, Hamra FK, Garbers DL: A multitude of genes expressed solely in meiotic or postmeiotic spermatogenic cells offers a myriad of contraceptive targets. Proc Natl Acad Sci U S A. 2003, 100: 12201-12206. 10.1073/pnas.1635054100.PubMedPubMed CentralView ArticleGoogle Scholar
- Shima JE, McLean DJ, McCarrey JR, Griswold MD: The murine testicular transcriptome: characterizing gene expression in the testis during the progression of spermatogenesis. Biol Reprod. 2004, 71: 319-330. 10.1095/biolreprod.103.026880.PubMedView ArticleGoogle Scholar
- Laiho A, Kotaja N, Gyenesei A, Sironen A: Transcriptome profiling of the murine testis during the first wave of spermatogenesis. PLoS One. 2013, 8: e61558-10.1371/journal.pone.0061558.PubMedPubMed CentralView ArticleGoogle Scholar
- Bao J, Zhang J, Zheng H, Xu C, Yan W: UBQLN1 interacts with SPEM1 and participates in spermiogenesis. Mol Cell Endocrinol. 2010, 327: 89-97. 10.1016/j.mce.2010.06.006.PubMedPubMed CentralView ArticleGoogle Scholar
- Chalmel F, Rolland AD, Niederhauser-Wiederkehr C, Chung SS, Demougin P, Gattiker A, Moore J, Patard JJ, Wolgemuth DJ, Jégou B, Primig M: The conserved transcriptome in human and rodent male gametogenesis. Proc Natl Acad Sci U S A. 2007, 104: 8346-8351. 10.1073/pnas.0701883104.PubMedPubMed CentralView ArticleGoogle Scholar
- Yang Z: PAML 4: Phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007, 24: 1586-1591. 10.1093/molbev/msm088.PubMedView ArticleGoogle Scholar
- Fares MA, Bezemer D, Moya A, Marín I: Selection on coding regions determined Hox7 genes evolution. Mol. Biol. Evol. 2003, 20: 2104-2112. 10.1093/molbev/msg222.PubMedView ArticleGoogle Scholar
- Xu B, Yang Z: PAMLX: A graphical user interface for PAML. Mol. Biol Evol. 2013, 30: 2723-2724. 10.1093/molbev/mst179.PubMedView ArticleGoogle Scholar
- Marín I: Diversification and specialization of plant RBR ubiquitin ligases. PLoS One. 2010, 5: e11579-10.1371/journal.pone.0011579.PubMedPubMed CentralView ArticleGoogle Scholar
- Marín I: Evolution of plant HECT ubiquitin ligases. PloS One. 2013, 8: e68536-10.1371/journal.pone.0068536.PubMedPubMed CentralView ArticleGoogle Scholar
- N'Diaye EN, Kajihara KK, Hsieh I, Morisaki H, Debnath J, Brown EJ: PLIC proteins or ubiquilins regulate autophagy-dependent cell survival during nutrient starvation. EMBO Rep. 2009, 10: 173-179. 10.1038/embor.2008.238.PubMedPubMed CentralView ArticleGoogle Scholar
- Rothenberg C, Srinivasan D, Mah L, Kaushik S, Peterhoff CM, Ugolino J, Fang S, Cuervo AM, Nixon RA, Monteiro MJ: Ubiquilin functions in autophagy and is degraded by chaperone-mediated autophagy. Hum Mol Genet. 2010, 19: 3219-3232. 10.1093/hmg/ddq231.PubMedPubMed CentralView ArticleGoogle Scholar
- Lee DY, Arnott D, Brown EJ: Ubiquilin4 is an adaptor protein that recruits Ubiquilin1 to the autophagy pathway. EMBO Rep. 2013, 14: 373-38. 10.1038/embor.2013.22.PubMedView ArticleGoogle Scholar
- Fatimababy AS, Lin YL, Usharani R, Radjacommare R, Wang HT, Tsai HL, Lee Y, Fu H: Cross-species divergence of the major recognition pathways of ubiquitylated substrates for ubiquitin/26S proteasome-mediated proteolysis. FEBS J. 2010, 277: 796-816. 10.1111/j.1742-4658.2009.07531.x.PubMedView ArticleGoogle Scholar
- Zhang D, Raasi S, Fushman D: Affinity makes the difference: nonselective interaction of the UBA domain of Ubiquilin-1 with monomeric ubiquitin and polyubiquitin chains. J Mol Biol. 2008, 377: 162-180. 10.1016/j.jmb.2007.12.029.PubMedPubMed CentralView ArticleGoogle Scholar
- Sutovsky P: Ubiquitin-dependent proteolysis in mammalian spermatogenesis, fertilization, and sperm quality control: killing three birds with one stone. Microsc Res Tech. 2003, 61: 88-102. 10.1002/jemt.10319.PubMedView ArticleGoogle Scholar
- Hermo L, Pelletier RM, Cyr DG, Smith CE: Surfing the wave, cycle, life history, and genes/proteins expressed by testicular germ cells. Part 4: intercellular bridges, mitochondria, nuclear envelope, apoptosis, ubiquitination, membrane/voltage-gated channels, methylation/acetylation, and transcription factors. Microsc Res Tech. 2010, 73: 364-408.PubMedView ArticleGoogle Scholar
- Hou CC, Yang WX: New insights to the ubiquitin-proteasome pathway (UPP) mechanism during spermatogenesis. Mol Biol Rep. 2013, 40: 3213-3230. 10.1007/s11033-012-2397-y.PubMedView ArticleGoogle Scholar
- Hou X, Zhang W, Xiao Z, Gan H, Lin X, Hou X, Zhang W, Xiao Z, Gan H, Lin X, Liao S, Han C: Mining and characterization of ubiquitin E3 ligases expressed in the mouse testis. BMC Genomics. 2012, 13: 495-10.1186/1471-2164-13-495.PubMedPubMed CentralView ArticleGoogle Scholar
- Tanaka H, Baba T: Gene expression in spermiogenesis. Cell Mol Life Sci. 2005, 62: 344-354. 10.1007/s00018-004-4394-y.PubMedView ArticleGoogle Scholar
- Hermo L, Pelletier RM, Cyr DG, Smith CE: Surfing the wave, cycle, life history, and genes/proteins expressed by testicular germ cells. Part 2: changes in spermatid organelles associated with development of spermatozoa. Microsc Res Tech. 2010, 73: 279-319.PubMedView ArticleGoogle Scholar
- Marín I: RBR ubiquitin ligases: diversification and streamlining in animal lineages. J Mol Evol. 2009, 69: 54-64. 10.1007/s00239-009-9252-3.PubMedView ArticleGoogle Scholar
- Marín I: Origin and diversification of TRIM ubiquitin ligases. PLoS One. 2012, 7: e50030-10.1371/journal.pone.0050030.PubMedPubMed CentralView ArticleGoogle Scholar
- Bedford JM: Puzzles of mammalian fertilization–and beyond. Int J Dev Biol. 2008, 52: 415-426. 10.1387/ijdb.072551jb.PubMedView ArticleGoogle Scholar
- Nixon B, Ecroyd HW, Dacheux JL, Jones RC: Monotremes provide a key to understanding the evolutionary significance of epididymal sperm maturation. J Androl. 2011, 32: 665-671. 10.2164/jandrol.110.012716.PubMedView ArticleGoogle Scholar
- Sun X, Kovacs T, Hu YJ, Yang WX: The role of actin and myosin during spermatogenesis. Mol Biol Rep. 2011, 38: 3993-4001. 10.1007/s11033-010-0517-0.PubMedView ArticleGoogle Scholar
- Sperry AO: The dynamic cytoskeleton of the developing male germ cell. Biol Cell. 2012, 104: 297-305. 10.1111/boc.201100102.PubMedView ArticleGoogle Scholar
- Katoh K, Toh H: Recent developments in the MAFFT multiple sequence alignment program. Brief Bioinform. 2008, 9: 286-298. 10.1093/bib/bbn013.PubMedView ArticleGoogle Scholar
- Nicholas KB, Nicholas HB: GeneDoc: a tool for editing and annotating multiple sequence alignments. 1997, Distributed by the authorGoogle Scholar
- Marín I: Ancient origin of animal U-box ubiquitin ligases. BMC Evol Biol. 2010, 10: 331-10.1186/1471-2148-10-331.PubMedPubMed CentralView ArticleGoogle Scholar
- Marín I: Animal HECT ubiquitin ligases: evolution and functional implications. BMC Evol Biol. 2010, 10: 56-10.1186/1471-2148-10-56.PubMedPubMed CentralView ArticleGoogle Scholar
- Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011, 28: 2731-2739. 10.1093/molbev/msr121.PubMedPubMed CentralView ArticleGoogle Scholar
- Swofford DL: PAUP*. Phylogenetic Analysis Using Parsimony (*and Other Methods). Version 4. 2003, Sunderland, Massachusetts: Sinauer AssociatesGoogle Scholar
- Flicek P, Ahmed I, Amode MR, Barrell D, Beal K, Brent S, Carvalho-Silva D, Clapham P, Coates G, Fairley S, Fitzgerald S, Gil L, García-Girón C, Gordon L, Hourlier T, Hunt S, Juettemann T, Kähäri AK, Keenan S, Komorowska M, Kulesha E, Longden I, Maurel T, McLaren WM, Muffato M, Nag R, Overduin B, Pignatelli M, Pritchard B, Pritchard E, et al: Ensembl. Nucleic Acids Res. 2013, 41 (Database issue): D48-D55.PubMedPubMed CentralView ArticleGoogle Scholar
- Zdobnov EM, Apweiler R: InterProScan - an integration platform for the signature-recognition methods in InterPro. Bioinformatics. 2001, 17: 847-848. 10.1093/bioinformatics/17.9.847.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.