Reprogrammed genetic decoding in cellular gene expression

Reprogrammed genetic decoding signals in mRNAs productively overwrite the normal decoding rules of translation. These "recoding" signals are associated with sites of programmed ribosomal frameshifting, hopping, termination codon suppression, and the incorporation of the unusual amino acids selenocysteine and pyrrolysine. This review summarizes current knowledge of the structure and function of recoding signals in cellular genes, the biological importance of recoding in gene regulation, and ways to identify new recoded genes.     

Recoding in archaea

Standard decoding of the genetic information into polypeptides is performed by one of the most sophisticated cell machineries, the translating ribosome, which, by following the genetic code, ensures the correspondence between the mature mRNA and the protein sequence. However, the expression of a minority of genes requires programmed deviations from the standard decoding rules, globally named recoding. This includes ribosome programmed -/+1 frameshifting, ribosome hopping, and stop codon readthrough. Recoding in Archaea was unequivocally demonstrated only for the translation of the UGA stop codon into the amino acid selenocysteine. However, a new recoding event leading to the 22nd amino acid pyrrolysine and the preliminary reports on a gene regulated by programmed -1 frameshifting have been recently described in Archaea. Therefore, it appears that the study of this phenomenon in Archaea is still at its dawn and that most of the genes whose expression is regulated by recoding are still uncharacterized.     

The gene of an archaeal alpha-L-fucosidase is expressed by translational frameshifting

The standard rules of genetic translational decoding are altered in specific genes by different events that are globally termed recoding. In Archaea recoding has been unequivocally determined so far only for termination codon readthrough events. We study here the mechanism of expression of a gene encoding for a alpha-l-fucosidase from the archaeon Sulfolobus solfataricus (fucA1), which is split in two open reading frames separated by a -1 frameshifting. The expression in Escherichia coli of the wild-type split gene led to the production by frameshifting of full-length polypeptides with an efficiency of 5%. Mutations in the regulatory site where the shift takes place demonstrate that the expression in vivo occurs in a programmed way. Further, we identify a full-length product of fucA1 in S.solfataricus extracts, which translate this gene in vitro by following programmed -1 frameshifting. This is the first experimental demonstration that this kind of recoding is present in Archaea.     

Efficient incorporation of multiple selenocysteines involves an inefficient decoding step serving as a potential translational checkpoint and ribosome bottleneck

Selenocysteine is incorporated into proteins via "recoding" of UGA from a stop codon to a sense codon, a process that requires specific secondary structures in the 3' untranslated region, termed selenocysteine incorporation sequence (SECIS) elements, and the protein factors that they recruit. Whereas most selenoprotein mRNAs contain a single UGA codon and a single SECIS element, selenoprotein P genes encode multiple UGAs and two SECIS elements. We have identified evolutionary adaptations in selenoprotein P genes that contribute to the efficiency of incorporating multiple selenocysteine residues in this protein. The first is a conserved, inefficiently decoded UGA codon in the N-terminal region, which appears to serve both as a checkpoint for the presence of factors required for selenocysteine incorporation and as a "bottleneck," slowing down the progress of elongating ribosomes. The second adaptation involves the presence of introns downstream of this inefficiently decoded UGA which confer the potential for nonsense-mediated decay when factors required for selenocysteine incorporation are limiting. Third, the two SECIS elements in selenoprotein P mRNA function with differing efficiencies, affecting both the rate and the efficiency of decoding different UGAs. The implications for how these factors contribute to the decoding of multiple selenocysteine residues are discussed.     

In silico identification of novel selenoproteins in the Drosophila melanogaster genome

In selenoproteins, incorporation of the amino acid selenocysteine is specified by the UGA codon, usually a stop signal. The alternative decoding of UGA is conferred by an mRNA structure, the SECIS element, located in the 3′-untranslated region of the selenoprotein mRNA. Because of the non-standard use of the UGA codon, current computational gene prediction methods are unable to identify selenoproteins in the sequence of the eukaryotic genomes. Here we describe a method to predict selenoproteins in genomic sequences, which relies on the prediction of SECIS elements in coordination with the prediction of genes in which the strong codon bias characteristic of protein coding regions extends beyond a TGA codon interrupting the open reading frame. We applied the method to the Drosophila melanogaster genome, and predicted four potential selenoprotein genes. One of them belongs to a known family of selenoproteins, and we have tested experimentally two other predictions with positive results. Finally, we have characterized the expression pattern of these two novel selenoprotein genes.     

Codon bias at the 3'-side of the initiation codon is correlated with translation initiation efficiency in Escherichia coli

The codon that follows the AUG initiation triplet (+2 codon) affects gene expression in Escherichia coli. We have extended this analysis using two model genes lacking any apparent Shine-Dalgarno sequence. Depending on the identity of the +2 codon a difference in gene expression up to 20-fold could be obtained. The effects did not correlate with the levels of intracellular pools of cognate tRNA for the +2 codon, with putative secondary mRNA structures, or with mRNA stability. However, most +2 iso-codons that were decoded by the same species of tRNA gave pairwise similar effects, suggesting that the effect on gene expression was associated with the decoding tRNA. High adenine content of the +2 codon was associated with high gene expression. Of the fourteen +2 codons that mediated the highest efficiency, all except two had an adenine as the first base of the codon. Analysis of the 3540 E. coli genes from the TransTerm database revealed that codons associated with high gene expression in the two expression systems are over-represented at the +2 position in natural genes. Codons that are associated with low gene expression are under-represented. The data suggest that evolution has favored codons at the +2 position that give high translation initiation.     

Highly expressed and alien genes of the Synechocystis genome

Comparisons of codon frequencies of genes to several gene classes are used to characterize highly expressed and alien genes on the SYNECHOCYSTIS: PCC6803 genome. The primary gene classes include the ensemble of all genes (average gene), ribosomal protein (RP) genes, translation processing factors (TF) and genes encoding chaperone/degradation proteins (CH). A gene is predicted highly expressed (PHX) if its codon usage is close to that of the RP/TF/CH standards but strongly deviant from the average gene. Putative alien (PA) genes are those for which codon usage is significantly different from all four classes of gene standards. In SYNECHOCYSTIS:, 380 genes were identified as PHX. The genes with the highest predicted expression levels include many that encode proteins vital for photosynthesis. Nearly all of the genes of the RP/TF/CH gene classes are PHX. The principal glycolysis enzymes, which may also function in CO(2) fixation, are PHX, while none of the genes encoding TCA cycle enzymes are PHX. The PA genes are mostly of unknown function or encode transposases. Several PA genes encode polypeptides that function in lipopolysaccharide biosynthesis. Both PHX and PA genes often form significant clusters (operons). The proteins encoded by PHX and PA genes are described with respect to functional classifications, their organization in the genome and their stoichiometry in multi-subunit complexes.     

The CUG codon is decoded in vivo as serine and not leucine in Candida albicans.

Previous studies have shown that the yeast Candida albicans encodes a unique seryl-tRNA(CAG) that should decode the leucine codon CUG as serine. However, in vitro translation of several different CUG-containing mRNAs in the presence of this unusual seryl-tRNA(CAG) result in an apparent increase in the molecular weight of the encoded polypeptides as judged by SDS-PAGE even though the molecular weight of serine is lower than that of leucine. A possible explanation for this altered electrophoretic mobility is that the CUG codon is decoded as modified serine in vitro. To elucidate the nature of CUG decoding in vivo, a reporter system based on the C. albicans gene (RBP1) encoding rapamycin-binding protein (RBP), coupled to the promoter of the C. albicans TEF3 gene, was utilized. Sequencing and mass-spectrometry analysis of the recombinant RBP expressed in C. albicans demonstrated that the CUG codon was decoded exclusively as serine while the related CUU codon was translated as leucine. A database search revealed that 32 out of the 65 C. albicans gene sequences available have CUG codons in their open reading frames. The CUG-containing genes do not belong to any particular gene family. Thus the amino acid specified by the CUG codon has been reassigned within the mRNAs of C. albicans. We argue here that this unique genetic code change in cellular mRNAs cannot be explained by the 'Codon Reassignment Theory'.     

Nucleotide sequence and transcript analysis of three photosystem II genes from the cyanobacterium Synechococcus sp. PCC7942

The genome of the cyanobacterium Synechococcus sp. PCC7942 contains two genes encoding the D2 polypeptide of photosystem II (PSII), which are designated here as psbDI and psbDII. The psbDI gene, like the psbD gene of plant chloroplasts, is cotranscribed with and overlaps the open reading frame of the psbC gene, encoding the PSII protein CP43. The psbDII gene is not linked to psbC, and appears to be transcribed as a monocistronic message. The two psbD genes encode identical polypeptides of 352 amino acids, which are 86% conserved with the D2 polypeptide of spinach. In plants, the translational start codon of the psbC gene has been reported to be an ATG codon 50 bp upstream from the end of the psbD gene. This triplet is not present in the psbDI sequence of Synechococcus sp., but is replaced by ACG, a codon which is very unlikely to initiate translation. Translation of the psbC gene may begin at a GTG codon which overlap the psbDI open reading frame by 14 bp and is preceded by a block of homology to the 3' end of the 16S ribosomal RNA, a potential ribosome-binding site. There are only two bp differences between the sequences of the two psbD genes; one of these results in substitution in psbDII of GCG for the presumed GTG start codon in psbDI.     

Compositional dynamics and codon usage pattern of BRCA1 gene across nine mammalian species

The BRCA1 gene is located on the human chromosome 17q21.31 and plays important role in biological processes. The aminoacyl-tRNA synthetases (AARS) are a family of heterogenous enzymes responsible protein synthesis and whose secondary functions include a role in autoimmune myositis. Our findings reveal that the compositional constraint and the preference of more A/T –ending codons determine the codon usage patterns in BRCA1 gene while more G/C-ending codons influence the codon usage pattern of AARS gene among mammals. The codon usage bias in BRCA1 and AARS genes is low. The codon CGC encoding arginine amino acid and the codon TTA encoding leucine were uniformly distributed in BRCA1 and AARS genes, respectively in mammals including human. Natural selection might have played a major role while mutation pressure might have played a minor role in shaping the codon usage pattern of BRCA1 and AARS genes.     

RECODE: a database of frameshifting, bypassing and codon redefinition utilized for gene expression

The RECODE database is a compilation of 'programmed' translational recoding events taken from the scientific literature and personal communications. The database deals with programmed ribosomal frameshifting, codon redefinition and translational bypass occurring in a variety of organisms. The entries for each event include the sequences of the corresponding genes, their encoded proteins for both the normal and alternate decoding, the types of the recoding events involved, trans-factors and cis-elements that influence recoding. The database is freely available at http://recode.genetics. utah.edu/.     

Factors that determine cell-specific gene expression in pancreatic endocrine tumor cells

Phenotypically distinct islet tumor cell lines may recapitulate certain of the developmental pathways of normal islet cell differentiation by expressing a combinatorial set of positively and negatively acting DNA-binding proteins to allow for the programmed expression of genes encoding polypeptide hormones. The structure of one of these DNA-binding proteins, a cyclic AMP-responsive protein (CREB) that binds specific DNA regulatory elements in the somatostatin gene, has been deduced from the sequence of a cloned cDNA. The CREB protein contains a DNA-binding domain separate from a cAMP-dependent protein kinase A activation domain. Further characterizations of the genes encoding the DNA-binding proteins should help to elucidate the cellular processes involved in islet cell differentiation and the genesis of tumors.     

Effect of codon adaptation on codon-level and gene-level translation efficiency in vivo

Background

There is a significant difference between synonymous codon usage in many organisms, and it is known that codons used more frequently generally showed efficient decoding rate. At the gene level, however, there are conflicting reports on the existence of a correlation between codon adaptation and translation efficiency, even in the same organism.

Results

To resolve this issue, we cultured Escherichia coli under conditions designed to maintain constant levels of mRNA and protein and subjected the cells to ribosome profiling (RP) and mRNA-seq analyses. We showed that the RP results correlated more closely with protein levels generated under similar culture conditions than with the mRNA abundance from the mRNA-seq. Our result indicated that RP/mRNA ratio could be used as a measure of translation efficiency at gene level. On the other hand, the RP data showed that codon-specific ribosome density at the decoding site negatively correlated with codon usage, consistent with the hypothesis that preferred codons display lower ribosome densities due to their faster decoding rate. However, highly codon-adapted genes showed higher ribosome densities at the gene level, indicating that the efficiency of translation initiation, rather than higher elongation efficiency of preferred codons, exerted a greater effect on ribosome density and thus translation efficiency.

Conclusions

These findings indicate that evolutionary pressure on highly expressed genes influenced both codon bias and translation initiation efficiency and therefore explains contradictory findings that codon usage bias correlates with translation efficiency of native genes, but not with the artificially created gene pool, which was not subjected to evolution pressure.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1115) contains supplementary material, which is available to authorized users.