Essential factors determining codon usage in ubiquitin genes

Summary Ubiquitin is ubiquitous in all eukaryotes and its amino acid sequence shows extreme conservation. Ubiquitin genes comprise direct repeats of the ubiquitin coding unit with no spacers. The nucleotide sequences coding for 13 ubiquitin genes from 11 species reported so far have been compiled and analyzed. The G+C content of codon third base reveals a positive linear correlation with the genome G+C content of the corresponding species. The slope strongly suggests that the overall G+C content of codons of polyubiquitin genes clearly reflects the genome G+C content by AT/GC substitutions at the codon third position. The G+C content of ubiquitin codon third base also shows a positive linear correlation with the overall G+C content of coding regions of compiled genes, indicating the codon choices among synonymous codons reflect the average codon usage pattern of corresponding species. On the other hand, the monoubiquitin gene, which is different from the polyubiquitin gene in gene organization, gene expression, and function of the encoding protein, shows a different codon usage pattern compared with that of the polyubiquitin gene. From comparisons of the levels of synonymous substitutions among ubiquitin repeats and the homology of the amino acid sequence of the tail of monomeric ubiquitin genes, we propose that the molecular evolution of ubiquitin genes occurred as follows: Plural primitive ubiquitin sequences were dispersed on genome in ancestral eukaryotes. Some of them situated in a particular environment fused with the tail sequence to produce monomeric ubiquitin genes that were maintained across species. After divergence of species, polyubiquitin genes were formed by duplication of the other primitive ubiquitin sequences on different chromosomes. Differences in the environments in which ubiquitin genes are embedded reflect the differences in codon choice and in gene expression pattern between poly- and monomeric ubiquitin genes.     

Concerted evolution in protists: Recent homogenization of a polyubiquitin gene in Trichomonas vaginalis

Ubiquitin is a 76-amino-acid protein with a remarkably high degree of conservation between all known sequences. Ubiquitin genes are almost always multicopy in eukaryotes, and often are found as polyubiquitin genes—fused tandem repeats which are coexpressed. Seventeen ubiquitin sequences from the amitochondrial protist Trichomonas vaginalis have been examined here, including an 11-repeat fragment of a polyubiquitin gene. These sequences reveal a number of interesting features that are not seen in other eukaryotes. The predicted amino acid sequences lack several universally conserved residues, and individual units do not always encode identical peptides as is usually the case. On the nucleotide level, these repeats are in general highly variable, but one region in the polyubiquitin is extremely homogeneous, with seven repeats absolutely identical. Such extended stretches of homogeneity have never been observed in ubiquitin genes and since substitutions are common in other coding units, it is likely that these repeats are the product of a very recent homogenization or amplification. Correspondence to: P.J. Keeling     

Ubiquitin genes as a paradigm of concerted evolution of tandem repeats

Summary Ubiquitin is remarkable for its ubiquitous distribution and its extreme protein sequence conservation. Ubiquitin genes comprise direct repeats of the ubiquitin coding unit with no spacers. The nucleotide sequences of several ubiquitin repeats from each of humans, chicken,Xenopus, Drosophila, barley, and yeast have recently been determined. By analysis of these data we show that ubiquitin is evolving more slowly than any other known protein, and that this (together with its gene organization) contributes to an ideal situation for the occurrence of concerted evolution of tandem repeats. By contrast, there is little evidence of between-cluster concerted evolution. We deduce that in ubiquitin genes, concerted evolution involves both unequal crossover and gene conversion, and that the average time since two repeated units within the polyubiquitin locus most recently shared a common ancestor is approximately 38 million years (Myr) in mammals, but perhaps only 11 Myr inDrosophila. The extreme conservatism of ubiquitin evolution also allows the inference that certain synonymous serine codons differing at the first two positions were probably mutated at single steps.     

Pea polyubiquitin genes: (I) structure and genomic organization

Four polyubiquitin genes, PUB1, PUB2, PUB3 and PUB4, were isolated from a pea (Pisum sativum L. cv Alaska) genomic library and completely sequenced. They represent all of the four polyubiquitin genes of the ubiquitin gene family in pea. The coding regions of the genes contain five or six coding units arranged as tandem repeats. The different coding repeats of the four genes share homologies between 75 and 97%, encoding the same protein of 76 amino acids identical to those from other higher plants. The open reading frames of PUB1, PUB2 and PUB4 terminate in the additional amino acid, phenylalanine (F), and PUB3 terminates in isoleucine (I). The polyubiquitin genes all contain intron sequences ranging from 584 to 1114 bp immediately 5′ to the ATG initiation codon of the first coding sequence. Of the four genes, three are associated with long AT-rich (85.4–89.4% A+T) sequences ranging from about 331 to 478 bp at their 5′ or 3′ ends. The PUB4 gene was found to be linked to a moderate to highly repetitive DNA at its 5′ flanking sequence. The greater sequence homology between different genes than among individual repeating units of a gene suggests that the polyubiquitin genes may have arisen by gene duplication of a single gene sequence.     

Structure and expression of the polyubiquitin gene in sea urchin embryos.

A cloned Lytechinus pictus cDNA has been identified, which includes seven direct repeats of a 228 bp sequence encoding ubiquitin and about 450 bp of 3' noncoding sequence. The deduced amino acid sequence is identical to that of ubiquitins of other animals (though repeats 3 and 5 each have single amino acid substitutions at different positions). Southern blot analysis revealed that the sea urchin genome contains a single copy of the polyubiquitin gene, and the number of 228 bp repeat units appears to vary from seven to ten among different alleles; no other ubiquitin coding sequences were detected. The size distribution of polyubiquitin mRNA is polymorphic among different individuals, probably corresponding to the differences in copy number of the repetitive coding sequence. The abundance of cytoplasmic polyubiquitin RNA is constant throughout embryogenesis and is similar in ectoderm, endoderm, and mesoderm cells. The constant prevalence of polyubiquitin mRNA apparently results from a balance between ontogenetic changes in its rate of synthesis and its stability in the presence of actinomycin D. Accumulation of polyubiquitin RNA was not heat shock-inducible during embryogenesis.     

Interspecific Comparison in the Frequency of Concerted Evolution at the Polyubiquitin Gene Locus

The polyubiquitin gene, encoding tandemly repeated multiple ubiquitins, constitutes a uniquitin gene subfamily. It has been demonstrated that polyubiquitin genes are subject to concerted evolution; namely, the individual ubiquitin coding units contained within a polyubiquitin gene are more similar to one another than they are to the ubiquitin coding units in the orthologous gene from other species. However there has been no comprehensive study on the concerted evolution of polyubiquitin genes in a wide range of species, because the relationships (orthologous or paralogous) among multiple polyubiquitin genes from different species have not been extensively analyzed yet. In this report, we present the results of analyzing the nucleotide sequence of polyubiquitin genes of mammals, available in the DDBJ/EMBL/GenBank nucleotide sequence databases, in which we found that there are two groups of polyubiquitin genes in an orthologous relationship. Based on this result, we analyzed the concerted evolution of the polyubiquitin gene in various species and compared the frequency of concerted evolutionary events interspecifically by taking into consideration that the rate of synonymous substitution at the polyubiquitin gene locus may vary depending on species. We found that the concerted evolutionary events in polyubiquitin genes have been more frequent in rats and Chinese hamsters than those in humans, cows, and sheep. The guinea pig polyubiquitin gene was an intermediate example. The frequency of concerted evolution in the mouse gene was unexpectedly low compared to that of other rodent genes. Received: 18 January 2000 / Accepted: 26 April 2000     

Cloning of a polyubiquitin gene fromNicotiana tabacum and comparison to other polyubiquitin genes

Using a tobacco cDNA clone as a probe, a genomic clone named TUQG-4, coding for a tobacco polyubiquitin protein with the five head-to-tail repeats of ubiquitin monomer was isolated. The five ubiquitin units were completely conserved except for the extra phenylalanine at the carboxy terminus of the last ubiquitin monomer. The putative open reading frame identified from the nucleotide sequence showed two possible intron sequences in the coding region for the first ubiquitin monomer. When the amino acid sequence deduced from the nucleotide sequence of TUQG-4 was compared to the amino acid sequences coded by other polyubiquitin genes of tobacco, there were three or four amino acid differences in the sequence. When the nucleotide sequences coding for the ubiquitin monomers were compared for various species origins, the degree of identity was at the highest between the ubiquitin monomers in one polyubiquitin and did not reflect the distance of the phylogenetic relationship.     

The yeast ubiquitin genes: a family of natural gene fusions.

Ubiquitin is a 76-residue protein highly conserved among eukaryotes. Conjugation of ubiquitin to intracellular proteins mediates their selective degradation in vivo. We describe a family of four ubiquitin-coding loci in the yeast Saccharomyces cerevisiae. UB11, UB12 and UB13 encode hybrid proteins in which ubiquitin is fused to unrelated ('tail') amino acid sequences. The ubiquitin coding elements of UB11 and UB12 are interrupted at identical positions by non-homologous introns. UB11 and UB12 encode identical 52-residue tails, whereas UB13 encodes a different 76-residue tail. The tail amino acid sequences are highly conserved between yeast and mammals. Each tail contains a putative metal-binding, nucleic acid-binding domain of the form Cys-X2-4-Cys-X2-15-Cys-X2-4-Cys, suggesting that these proteins may function by binding to DNA. The fourth gene, UB14, encodes a polyubiquitin precursor protein containing five ubiquitin repeats in a head-to-tail, spacerless arrangement. All four ubiquitin genes are expressed in exponentially growing cells, while in stationary-phase cells the expression of UB11 and UB12 is repressed. The UB14 gene, which is strongly inducible by starvation, high temperatures and other stresses, contains in its upstream region strong homologies to the consensus 'heat shock box' nucleotide sequence. Elsewhere we show that the essential function of the UB14 gene is to provide ubiquitin to cells under stress.     

Recent stable insertion of mitochondrial DNA into an Arabidopsis polyubiquitin gene by nonhomologous recombination.

Sequence analysis of a newly identified polyubiquitin gene (UBQ13) from the Columbia ecotype of Arabidopsis thaliana revealed that the gene contained a 3.9-kb insertion in the coding region. All subclones of the 3.9-kb insert hybridized to isolated mitochondrial DNA. The insert was found to consist of at least two, possibly three, distinct DNA segments from the mitochondrial genome. A 590-bp region of the insert is nearly identical to the Arabidopsis mitochondrial nad1 gene. UBQ13 restriction fragments in total cellular DNA from ecotypes Ler, No-0, Be-0, WS, and RLD were identified and, with the exception of Be-0, their sizes were equivalent to that predicted from the corresponding ecotype Columbia UBQ13 restriction fragment without the mitochondrial insert. Isolation by polymerase chain reaction and sequence determination of UBQ13 sequences from the other ecotypes showed that all lacked the mitochondrial insert. All ecotypes examined, except Columbia, contain intact open reading frames in the region of the insert, including four ubiquitin codons which Columbia lacks. This indicates that the mitochondrial DNA in UBQ13 in ecotype Columbia is the result of an integration event that occurred after speciation of Arabidopsis rather than a deletion event that occurred in all ecotypes except Columbia. This stable movement of mitochondrial DNA to the nucleus is so recent that there are few nucleotide changes subsequent to the transfer event. This allows for precise analysis of the sequences involved and elucidation of the possible mechanism. The presence of intron sequences in the transferred nucleic acid indicates that DNA was the transfer intermediate. The lack of sequence identity between the integrating sequence and the target site, represented by the other Arabidopsis ecotypes, suggests that integration occurred via nonhomologus recombination. This nuclear/organellar gene transfer event is strikingly similar to the experimentally accessible process of nuclear integration of introduced heterologous DNA.     

Variations and evolution of polyubiquitin genes from ciliates

Polyubiquitin genes from seven ciliate species were amplified, cloned and sequenced. It is estimated that Strombidium sulcatum, Euplotes vannus, E. rariseta and Anteholosticha manca have a polyubiquitin gene of 3 repeats, and A. parawarreni, Paramecium caudatum and Pseudokeronopsis flava 4 repeats. The newly obtained ubiquitins mostly differ from that of humans by 1–5 residues in amino acid sequences. A neighbor-joining tree constructed based on monomeric ubiquitin genes supports the monophyly of an assemblage comprising the litostomateans and some oligohymenophoreans, but not the class Spirotrichea. The monomers from the same species are generally placed together and highly supported for the class Litostomatea, the genera Paramecium and Ichthyophthirius, but not for other species. The non-synonymous/synonymous rate ratio (dN/dS) at the protein level are less than 1, and the synonymous nucleotide differences per synonymous site (p_S) from intraspecific comparisons are fairly high (0.02–0.72). These results indicate that ciliates have not only the conserved, but also some quite divergent, polyubiquitin genes and confirm that the polyubiquitin genes in ciliates evolve according to the birth-and-death mode of evolution under strong purifying selection.     

Molecular cloning, chromosome mapping and characterization of UBQLN3 a testis-specific gene that contains an ubiquitin-like domain

The sequence of the ubiquitin protein is highly conserved between species and has facilitated the cloning of numerous ubiquitin-like proteins. In the present study, we report the cloning of the cDNA for human ubiquilin 3 (UBQLN3). The deduced amino acid sequence of UBQLN3 contains a UBQ domain (ubiquitin-like) in the amino terminus as well as two highly conserved domains found in several recently cloned ubiquitin-like proteins. One of these domains, termed the NP domain, is a highly conserved 93 amino acid region present in UBQLN3 and several ubiquitin-like proteins. The last conserved domain is the UBA domain (ubiquitin-associated) found in a variety of proteins of the ubiquination pathway. The human UBQLN3 gene was mapped to the 11p15 region of chromosome 11. Northern blot analysis of multiple human and mouse tissues demonstrated UBQLN3 mRNA expression specifically in testis.     

Expression Enhancement of a Rice Polyubiquitin Gene Promoter

An 808 bp promoter from a rice polyubiquitin gene, rubi3, has been isolated. The rubi3 gene contained an open reading frame of 1140 bp encoding a pentameric polyubiquitin arranged as five tandem, head-to-tail repeats of 76 aa. The 1140 bp 5′ UTR intron of the gene enhanced its promoter activity in transient expression assays by 20-fold. Translational fusion of the GUS reporter gene to the coding sequence of the ubiquitin monomer enhanced GUS enzyme activity in transient expression assays by 4.3-fold over the construct containing the original rubi3 promoter (including the 5′ UTR intron) construct. The enhancing effect residing in the ubiquitin monomer coding sequence has been narrowed down to the first 9 nt coding for the first three amino acid residues of the ubiquitin protein. Mutagenesis at the third nucleotide of this 9 nt sequence still maintains the enhancing effect, but leads to translation of the native GUS protein rather than a fusion protein. The resultant 5′ regulatory sequence, consisting of the rubi3 promoter, 5′ UTR exon and intron, and the mutated first 9 nt coding sequence, has an activity nearly 90-fold greater than the rubi3 promoter only (without the 5′ UTR intron), and 2.2-fold greater than the maize Ubi1 gene promoter (including its 5′ UTR intron). The newly created expression vector is expected to enhance transgene expression in monocot plants. Considering the high conservation of the polyubiquitin gene structure in higher plants, the observed enhancement in gene expression may apply to 5′ regulatory sequences of other plant polyubiquitin genes.     

The intron of Arabidopsis thaliana polyubiquitin genes is conserved in location and is a quantitative determinant of chimeric gene expression

We have isolated and determined DNA sequence for the 5-flanking regions of three Arabidopsis thaliana polyubiquitin genes, UBQ3, UBQ10, and UBQ11. Comparison to cDNA sequences revealed the presence of an intron in the 5-untranslated region at the same position immediately upstream of the initiator methionine codon in each of the three genes. An intron at this position is also present in two sunflower and two maize polyubiquitin genes. An intron is also found in the 5-untranslated regions of several animal polyubiquitin genes, although the exact intron position is not conserved among them, and none are in the same position as those in the higher plant polyubiquitin genes. Chimeric genes containing the 5-flanking regions of UBQ3, UBQ10, and UBQ11 in front of the coding regions for the reporter enzyme Escherichia coli -glucuronidase (GUS) were constructed. When introduced transiently into Arabidopsis leaves via microprojectile bombardment, all resulted in readily detectable levels of GUS activity that were quantitatively similar. The introns of UBQ3 and UBQ10 in the corresponding promoter fragments were removed by replacement with flanking cDNA sequences and chimeric genes constructed. These constructs resulted in 2.5- to 3-fold lower levels of marker enzyme activity after transient introduction into Arabidopsis leaves. The UBQ10 promoter without the 5 intron placed upstream of firefly luciferase (LUX) resulted in an average of 3-fold lower LUX activity than from an equivalent construct with the UBQ10 intron. A UBQ3 promoter cassette was constructed for the constitutive expression of open reading frames in dicot plants and it produced readily detectable levels of GUS activity in transient assays.     

Structure and expression of ubiquitin genes of Drosophila melanogaster.

We isolated and characterized two related ubiquitin genes from Drosophila melanogaster, polyubiquitin and UB3-D. The polyubiquitin gene contained 18 repeats of the 228-base-pair monomeric ubiquitin-encoding unit arranged in tandem. This gene was localized to a minor heat shock puff site, 63F, and it encoded a constitutively expressed 4.4-kilobase polyubiquitin-encoding mRNA, whose level was induced threefold by heat shock. To investigate the pattern of expression of the polyubiquitin gene in developing animals, a polyubiquitin-lacZ fusion gene was introduced into the Drosophila genome by germ line transformation. The fusion gene was expressed at high levels in a tissue-general manner at all life stages assayed. The ubiquitin-encoding gene, UB3-D, consisted of one ubiquitin-encoding unit directly fused, in frame, to a nonhomologous tail sequence. The amino acid sequence of the tail portion of the protein had 65% positional identity with that of yeast UBI3 protein, including a region that contained a potential nucleic acid-binding motif. The Drosophila UB3-D gene hybridized to a 0.9-kilobase mRNA that was constitutively expressed, and in contrast to the polyubiquitin gene, it was not inducible by heat shock.     

UbiA, the major polyubiquitin locus in Caenorhabditis elegans, has unusual structural features and is constitutively expressed.

Ubiquitin is a multifunctional 76-amino-acid protein which plays critical roles in many aspects of cellular metabolism. In Caenorhabditis elegans, the major source of ubiquitin RNA is the polyubiquitin locus, UbiA. UbiA is transcribed as a polycistronic mRNA which contains 11 tandem repeats of ubiquitin sequence and possesses a 2-amino-acid carboxy-terminal extension on the final repeat. The UbiA locus possesses several unusual features not seen in the ubiquitin genes of other organisms studied to date. Mature UbiA mRNA acquires a 22-nucleotide leader sequence via a trans-splicing reaction involving a 100-nucleotide splice leader RNA derived from a different chromosome. UbiA is also unique among known polyubiquitin genes in containing four cis-spliced introns within its coding sequence. Thus, UbiA is one of a small class of genes found in higher eucaryotes whose heterogeneous nuclear RNA undergoes both cis and trans splicing. The putative promoter region of UbiA contains a number of potential regulatory elements: (i) a cytosine-rich block, (ii) two sequences resembling the heat shock regulatory element, and (iii) a palindromic sequence with homology to the DNA-binding site of the mammalian steroid hormone receptor. The expression of the UbiA gene has been studied under various heat shock conditions and has been monitored during larval moulting and throughout the major stages of development. These studies indicate that the expression of the UbiA gene is not inducible by acute or chronic heat shock and does not appear to be under nutritional or developmental regulation in C. elegans.