Codon replacement in the PGK1 gene of Saccharomyces cerevisiae: experimental approach to study the role of biased codon usage in gene expression

Roles of the CCR4-Not complex in translation and dynamics of co-translation events

The Ccr4-Not complex is a global regulator of mRNA metabolism in eukaryotic cells that is most well-known to repress gene expression. Delivery of the complex to mRNAs through a multitude of distinct mechanisms accelerates their decay, yet Ccr4-Not also plays an important role in co-translational processes, such as co-translational association of proteins and delivery of translating mRNAs to organelles. The recent structure of Not5 interacting with the translated ribosome has brought to light that embedded information within the codon sequence can be monitored by recruitment of the Ccr4-Not complex to elongating ribosomes. Thereby, the Ccr4-Not complex is empowered with regulatory decisions determining the fate of proteins being synthesized and their encoding mRNAs. This review will focus on the roles of the complex in translation and dynamics of co-translation events. This article is categorized under: Translation > Mechanisms Translation > Regulation.

Regulation of eukaryotic mRNA deadenylation and degradation by the Ccr4-Not complex

Accurate and precise regulation of gene expression programmes in eukaryotes involves the coordinated control of transcription, mRNA stability and translation. In recent years, significant progress has been made about the role of sequence elements in the 3' untranslated region for the regulation of mRNA degradation, and a model has emerged in which recruitment of the Ccr4-Not complex is the critical step in the regulation of mRNA decay. Recruitment of the Ccr4-Not complex to a target mRNA results in deadenylation mediated by the Caf1 and Ccr4 catalytic subunits of the complex. Following deadenylation, the 5' cap structure is removed, and the mRNA subjected to 5'-3' degradation. Here, the role of the human Ccr4-Not complex in cytoplasmic deadenylation of mRNA is reviewed, with a particular focus on mechanisms of its recruitment to mRNA by sequence motifs in the 3' untranslated region, codon usage, as well as general mechanisms involving the poly(A) tail.

Synonymous mutations in the phosphoglycerate kinase 1 gene induce an altered response to protein misfolding in Schizosaccharomyces pombe

Background

Proteostasis refers to the processes that regulate the biogenesis, folding, trafficking, and degradation of proteins. Any alteration in these processes can lead to cell malfunction. Protein synthesis, a key proteostatic process, is highly-regulated at multiple levels to ensure adequate adaptation to environmental and physiological challenges such as different stressors, proteotoxic conditions and aging, among other factors. Because alterations in protein translation can lead to protein misfolding, examining how protein translation is regulated may also help to elucidate in part how proteostasis is controlled. Codon usage bias has been implicated in the fine-tuning of translation rate, as more-frequent codons might be read faster than their less-frequent counterparts. Thus, alterations in codon usage due to synonymous mutations may alter translation kinetics and thereby affect the folding of the nascent polypeptide, without altering its primary structure. To date, it has been difficult to predict the effect of synonymous mutations on protein folding and cellular fitness due to a scarcity of relevant data. Thus, the purpose of this work was to assess the effect of synonymous mutations in discrete regions of the gene that encodes the highly-expressed enzyme 3-phosphoglycerate kinase 1 (pgk1) in the fission yeast Schizosaccharomyces pombe.

Results

By means of systematic replacement of synonymous codons along pgk1, we found slightly-altered protein folding and activity in a region-specific manner. However, alterations in protein aggregation, heat stress as well as changes in proteasome activity occurred independently of the mutated region. Concomitantly, reduced mRNA levels of the chaperones Hsp9 and Hsp16 were observed.

Conclusion

Taken together, these data suggest that codon usage bias of the gene encoding this highly-expressed protein is an important regulator of protein function and proteostasis.

LHX9, a p53-binding protein, inhibits the progression of glioma by suppressing glycolysis

PURPOSE

LHX9 methylation has been reported in many tumors, but its functions and related mechanisms in glioma are still unknown and need to be verified.

METHODS

The protein level of LHX9 in glioma tissues was examined using western blotting and immunohistochemistry, and the functions of LHX9 in glioma cell lines were investigated using MTT and colony formation assays. In addition, the interaction between LHX9 and P53 was analyzed by immunoprecipitation, and the roles of LHX9 in cancer metabolism were explored by measuring metabolites.

RESULTS

In this study, we found that the LHX9 expression level was decreased in glioma specimens, and the upregulation of LHX9 expression inhibited the growth of glioma cells in liquid medium and on soft agar. Regarding the molecular mechanism, we found that LHX9 interacted with p53, and downregulation of LHX9 promoted the expression of the glycolysis-related enzyme PGK1 and increased the lactic acid content. By interfering with the expression of LHX9, the tumorigenicity of glioma cells was promoted, an outcome blocked by further interference with PGK1 expression.

CONCLUSION

In summary, the decreased expression of LHX9 in gliomas activates the expression of the glycolysis-related enzyme PGK1, thereby promoting the development of gliomas, suggesting that the LHX9-PGK1 signaling axis can be used as a target for the treatment of glioma.

Ribosome states signal RNA quality control

Eukaryotic cells integrate multiple quality control (QC) responses during protein synthesis in the cytoplasm. These QC responses are signaled by slow or stalled elongating ribosomes. Depending on the nature of the delay, the signal may lead to translational repression, messenger RNA decay, ribosome rescue, and/or nascent protein degradation. Here, we discuss how the structure and composition of an elongating ribosome in a troubled state determine the downstream quality control pathway(s) that ensue. We highlight the intersecting pathways involved in RNA decay and the crosstalk that occurs between RNA decay and ribosome rescue.

A stress-induced tyrosine-tRNA depletion response mediates codon-based translational repression and growth suppression

Eukaryotic transfer RNAs can become selectively fragmented upon various stresses, generating tRNA-derived small RNA fragments. Such fragmentation has been reported to impact a small fraction of the tRNA pool and thus presumed to not directly impact translation. We report that oxidative stress can rapidly generate tyrosine-tRNAGUA fragments in human cells-causing significant depletion of the precursor tRNA. Tyrosine-tRNAGUA depletion impaired translation of growth and metabolic genes enriched in cognate tyrosine codons. Depletion of tyrosine tRNAGUA or its translationally regulated targets USP3 and SCD repressed proliferation-revealing a dedicated tRNA-regulated growth-suppressive pathway for oxidative stress response. Tyrosine fragments are generated in a DIS3L2 exoribonuclease-dependent manner and inhibit hnRNPA1-mediated transcript destabilization. Moreover, tyrosine fragmentation is conserved in C. elegans. Thus, tRNA fragmentation can coordinately generate trans-acting small RNAs and functionally deplete a tRNA. Our findings reveal the existence of an underlying adaptive codon-based regulatory response inherent to the genetic code.

Synonymous but Not Silent: The Codon Usage Code for Gene Expression and Protein Folding

Codon usage bias, the preference for certain synonymous codons, is found in all genomes. Although synonymous mutations were previously thought to be silent, a large body of evidence has demonstrated that codon usage can play major roles in determining gene expression levels and protein structures. Codon usage influences translation elongation speed and regulates translation efficiency and accuracy. Adaptation of codon usage to tRNA expression determines the proteome landscape. In addition, codon usage biases result in nonuniform ribosome decoding rates on mRNAs, which in turn influence the cotranslational protein folding process that is critical for protein function in diverse biological processes. Conserved genome-wide correlations have also been found between codon usage and protein structures. Furthermore, codon usage is a major determinant of mRNA levels through translation-dependent effects on mRNA decay and translation-independent effects on transcriptional and posttranscriptional processes. Here, we discuss the multifaceted roles and mechanisms of codon usage in different gene regulatory processes.

A rare codon-based translational program of cell proliferation

BACKGROUND

The speed of translation elongation is primarily determined by the abundance of tRNAs. Thus, the codon usage influences the rate with which individual mRNAs are translated. As the nature of tRNA pools and modifications can vary across biological conditions, codon elongation rates may also vary, leading to fluctuations in the protein production from individual mRNAs. Although it has been observed that functionally related mRNAs exhibit similar codon usage, presumably to provide an effective way to coordinate expression of multiple proteins, experimental evidence for codon-mediated translation efficiency modulation of functionally related mRNAs in specific conditions is scarce and the associated mechanisms are still debated.

RESULTS

Here, we reveal that mRNAs whose expression increases during cell proliferation are enriched in rare codons, poorly adapted to tRNA pools. Ribosome occupancy profiling and proteomics measurements show that upon increased cell proliferation, transcripts enriched in rare codons undergo a higher translation boost than transcripts with common codons. Re-coding of a fluorescent reporter with rare codons increased protein output by ~ 30% relative to a reporter re-coded with common codons. Although the translation capacity of proliferating cells was higher compared to resting cells, we did not find evidence for the regulation of individual tRNAs. Among the models that were proposed so far to account for codon-mediated translational regulation upon changing conditions, the one that seems most consistent with our data involves a global upregulation of ready-to-translate tRNAs, which we show can lead to a higher increase in the elongation velocity at rare codons compared to common codons.

CONCLUSIONS

We propose that the alleviation of translation bottlenecks in rapidly dividing cells enables preferential upregulation of pro-proliferation proteins, encoded by mRNAs that are enriched in rare codons.

Quantification and discovery of sequence determinants of protein-per-mRNA amount in 29 human tissues

Despite their importance in determining protein abundance, a comprehensive catalogue of sequence features controlling protein-to-mRNA (PTR) ratios and a quantification of their effects are still lacking. Here, we quantified PTR ratios for 11,575 proteins across 29 human tissues using matched transcriptomes and proteomes. We estimated by regression the contribution of known sequence determinants of protein synthesis and degradation in addition to 45 mRNA and 3 protein sequence motifs that we found by association testing. While PTR ratios span more than 2 orders of magnitude, our integrative model predicts PTR ratios at a median precision of 3.2-fold. A reporter assay provided functional support for two novel UTR motifs, and an immobilized mRNA affinity competition-binding assay identified motif-specific bound proteins for one motif. Moreover, our integrative model led to a new metric of codon optimality that captures the effects of codon frequency on protein synthesis and degradation. Altogether, this study shows that a large fraction of PTR ratio variation in human tissues can be predicted from sequence, and it identifies many new candidate post-transcriptional regulatory elements.

Non-invasive measurement of mRNA decay reveals translation initiation as the major determinant of mRNA stability

The cytoplasmic abundance of mRNAs is strictly controlled through a balance of production and degradation. Whereas the control of mRNA synthesis through transcription has been well characterized, less is known about the regulation of mRNA turnover, and a consensus model explaining the wide variations in mRNA decay rates remains elusive. Here, we combine non-invasive transcriptome-wide mRNA production and stability measurements with selective and acute perturbations to demonstrate that mRNA degradation is tightly coupled to the regulation of translation, and that a competition between translation initiation and mRNA decay -but not codon optimality or elongation- is the major determinant of mRNA stability in yeast. Our refined measurements also reveal a remarkably dynamic transcriptome with an average mRNA half-life of only 4.8 min - much shorter than previously thought. Furthermore, global mRNA destabilization by inhibition of translation initiation induces a dose-dependent formation of processing bodies in which mRNAs can decay over time.

Codon optimality, bias and usage in translation and mRNA decay

The advent of ribosome profiling and other tools to probe mRNA translation has revealed that codon bias - the uneven use of synonymous codons in the transcriptome - serves as a secondary genetic code: a code that guides the efficiency of protein production, the fidelity of translation and the metabolism of mRNAs. Recent advancements in our understanding of mRNA decay have revealed a tight coupling between ribosome dynamics and the stability of mRNA transcripts; this coupling integrates codon bias into the concept of codon optimality, or the effects that specific codons and tRNA concentrations have on the efficiency and fidelity of the translation machinery. In this Review, we first discuss the evidence for codon-dependent effects on translation, beginning with the basic mechanisms through which translation perturbation can affect translation efficiency, protein folding and transcript stability. We then discuss how codon effects are leveraged by the cell to tailor the proteome to maintain homeostasis, execute specific gene expression programmes of growth or differentiation and optimize the efficiency of protein production.

The link between adjacent codon pairs and mRNA stability

Background

Evidence in diverse organisms suggests that codon optimality is a major determinant of mRNA translation and degradation. Codon optimality is thought to act by modulating the efficiency of ribosome elongation. In Saccharomyces cerevisiae, a recent study has identified 17 adjacent codon pairs that mediate strong inhibition of translation elongation. However, relationships between the inhibitory codon pairs and other aspects of gene expression are unknown.

Results

To gain insights into how the inhibitory codon pairs may affect aspects of gene expression, we utilized existing datasets to conduct genome-scale analyses in S. cerevisiae. Our analysis revealed the following points. First, the inhibitory codon pairs are significantly associated with faster mRNA decay. The association is not solely due to the content of nucleotides, individual codons, or dipeptides encoded by the inhibitory codon pairs. Second, the inhibitory codon pairs cannot fully explain the previously known relationship of codon optimality with mRNA stability, suggesting that optimality of individual codons and properties of adjacent codon pairs both contribute to gene regulation. Finally, although the inhibitory codon pairs are associated with slower mRNA synthesis and protein instability, the associations can be attributed to usage bias in individual codons.

Conclusions

This study suggests an association of inhibitory codon pairs with mRNA stability and thus another layer of complexity in the codon-mediated gene regulation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3749-8) contains supplementary material, which is available to authorized users.

Synonymous Codons: Choose Wisely for Expression

The genetic code, which defines the amino acid sequence of a protein, also contains information that influences the rate and efficiency of translation. Neither the mechanisms nor functions of codon-mediated regulation were well understood. The prevailing model was that the slow translation of codons decoded by rare tRNAs reduces efficiency. Recent genome-wide analyses have clarified several issues. Specific codons and codon combinations modulate ribosome speed and facilitate protein folding. However, tRNA availability is not the sole determinant of rate; rather, interactions between adjacent codons and wobble base pairing are key. One mechanism linking translation efficiency and codon use is that slower decoding is coupled to reduced mRNA stability. Changes in tRNA supply mediate biological regulationfor instance,, changes in tRNA amounts facilitate cancer metastasis.

Analysis of the association between codon optimality and mRNA stability in Schizosaccharomyces pombe

Background

Recent experiments have shown that codon optimality is a major determinant of mRNA stability in Saccharomyces cerevisiae and that this phenomenon may be conserved in Escherichia coli and some metazoans, although work in Neurospora crassa is not consistent with this model.

Results

We examined the association between codon optimality and mRNA stability in the fission yeast Schizosaccharomyces pombe. Our analysis revealed the following points. First, we observe a genome-wide association between codon optimality and mRNA stability also in S. pombe, suggesting evolutionary conservation of the phenomenon. Second, in both S. pombe and S. cerevisiae, mRNA synthesis rates are also correlated at the genome-wide analysis with codon optimality, suggesting that the long-appreciated association between codon optimality and mRNA abundance is due to regulation of both mRNA synthesis and degradation. However, when we examined correlation of codon optimality and either mRNA half-lives or synthesis rates controlling for mRNA abundance, codon optimality was still positively correlated with mRNA half-lives in S. cerevisiae, but the association was no longer significant for mRNA half-lives in S. pombe or for synthesis rates in either organism. This illustrates how only the pairwise analysis of multiple correlating variables may limit these types of analyses. Finally, in S. pombe, codon optimality is associated with known DNA/RNA sequence motifs that are associated with mRNA production/stability, suggesting these two features have been under similar selective pressures for optimal gene expression.

Conclusions

Consistent with the emerging body of studies, this study suggests that the association between codon optimality and mRNA stability may be a broadly conserved phenomenon. It also suggests that the association can be explained at least in part by independent adaptations of codon optimality and other transcript features for elevated expression during evolution.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3237-6) contains supplementary material, which is available to authorized users.

The Yin and Yang of codon usage

The genetic code is degenerate. With the exception of two amino acids (Met and Trp), all other amino acid residues are each encoded by multiple, so-called synonymous codons. Synonymous codons were initially presumed to have entirely equivalent functions, however, the finding that synonymous codons are not present at equal frequencies in genes/genomes suggested that codon choice might have functional implications beyond amino acid coding. The pattern of non-uniform codon use (known as codon usage bias) varies between organisms and represents a unique feature of an organism. Organism-specific codon choice is related to organism-specific differences in populations of cognate tRNAs. This implies that, in a given organism, frequently used codons will be translated more rapidly than infrequently used ones and vice versa A theory of codon-tRNA co-evolution (necessary to balance accurate and efficient protein production) was put forward to explain the existence of codon usage bias. This model suggests that selection favours preferred (frequent) over un-preferred (rare) codons in order to sustain efficient protein production in cells and that a given un-preferred codon will have the same effect on an organism's fitness regardless of its position within an mRNA's open reading frame. However, many recent studies refute this prediction. Un-preferred codons have been found to have important functional roles and their effects appeared to be position-dependent. Synonymous codon usage affects the efficiency/stringency of mRNA decoding, mRNA biogenesis/stability, and protein secretion and folding. This review summarizes recent developments in the field that have identified novel functions of synonymous codons and their usage.

Codon optimality is a major determinant of mRNA stability

mRNA degradation represents a critical regulated step in gene expression. Although the major pathways in turnover have been identified, accounting for disparate half-lives has been elusive. We show that codon optimality is one feature that contributes greatly to mRNA stability. Genome-wide RNA decay analysis revealed that stable mRNAs are enriched in codons designated optimal, whereas unstable mRNAs contain predominately non-optimal codons. Substitution of optimal codons with synonymous, non-optimal codons results in dramatic mRNA destabilization, whereas the converse substitution significantly increases stability. Further, we demonstrate that codon optimality impacts ribosome translocation, connecting the processes of translation elongation and decay through codon optimality. Finally, we show that optimal codon content accounts for the similar stabilities observed in mRNAs encoding proteins with coordinated physiological function. This work demonstrates that codon optimization exists as a mechanism to finely tune levels of mRNAs and, ultimately, proteins.

A critical analysis of codon optimization in human therapeutics

Codon optimization describes gene engineering approaches that use synonymous codon changes to increase protein production. Applications for codon optimization include recombinant protein drugs and nucleic acid therapies, including gene therapy, mRNA therapy, and DNA/RNA vaccines. However, recent reports indicate that codon optimization can affect protein conformation and function, increase immunogenicity, and reduce efficacy. We critically review this subject, identifying additional potential hazards including some unique to nucleic acid therapies. This analysis highlights the evolved complexity of codon usage and challenges the scientific bases for codon optimization. Consequently, codon optimization may not provide the optimal strategy for increasing protein production and may decrease the safety and efficacy of biotech therapeutics. We suggest that the use of this approach is reconsidered, particularly for in vivo applications.

Selection on codon bias in yeast: a transcriptional hypothesis

Codons that code for the same amino acid are often used with unequal frequencies. This phenomenon is termed codon bias. Here, we report a computational analysis of codon bias in yeast using experimental and theoretical genome-wide data. We show that the most used codons in highly expressed genes can be predicted by mRNA structural data and that the codon choice at each synonymous site within an mRNA is not random with respect to the local secondary structure. Because we also found that the folding stability of intron sequences is strongly correlated with codon bias and mRNA level, our results suggest that codon bias is linked to mRNA folding structure through a mechanism that, at least partially, operates before pre-mRNA splicing. Consistent with this, we report evidence supporting the adaptation of the tRNA pool to the codon profile of the most expressed genes rather than vice versa. We show that the correlation of codon usage with the gene expression level also includes the stop codons that are normally not decoded by aminoacyl-tRNAs. The results reported here are consistent with a role for transcriptional forces in driving codon usage bias via a mechanism that improves gene expression by optimizing mRNA folding structures.

Influence of codon bias on heterologous production of human papillomavirus type 16 major structural protein L1 in yeast

Heterologous gene expression is dependent on multistep processes involving regulation at the level of transcription, mRNA turnover, protein translation, and posttranslational modifications. Codon bias has a significant influence on protein yields. However, sometimes it is not clear which parameter causes observed differences in heterologous gene expression as codon adaptation typically optimizes many sequence properties at once. In the current study, we evaluated the influence of codon bias on heterologous production of human papillomavirus type 16 (HPV-16) major structural protein L1 in yeast by expressing five variants of codon-modified open reading frames (OFRs) encoding HPV-16 L1 protein. Our results showed that despite the high toleration of various codons used throughout the length of the sequence of heterologously expressed genes in transformed yeast, there was a significant positive correlation between the gene's expression level and the degree of its codon bias towards the favorable codon usage. The HPV-16 L1 protein expression in yeast can be optimized by adjusting codon composition towards the most preferred codon adaptation, and this effect most probably is dependent on the improved translational elongation.

Mapping codon usage of the translation initiation region in porcine reproductive and respiratory syndrome virus genome

BACKGROUND

Porcine reproductive and respitatory syndrome virus (PRRSV) is a recently emerged pathogen and severely affects swine populations worldwide. The replication of PRRSV is tightly controlled by viral gene expression and the codon usage of translation initiation region within each gene could potentially regulate the translation rate. Therefore, a better understanding of the codon usage pattern of the initiation translation region would shed light on the regulation of PRRSV gene expression.

RESULTS

In this study, the codon usage in the translation initiation region and in the whole coding sequence was compared in PRRSV ORF1a and ORFs2-7. To investigate the potential role of codon usage in affecting the translation initiation rate, we established a codon usage model for PRRSV translation initiation region. We observed that some non-preferential codons are preferentially used in the translation initiation region in particular ORFs. Although some positions vary with codons, they intend to use codons with negative CUB. Furthermore, our model of codon usage showed that the conserved pattern of CUB is not directly consensus with the conserved sequence, but shaped under the translation selection.

CONCLUSIONS

The non-variation pattern with negative CUB in the PRRSV translation initiation region scanned by ribosomes is considered the rate-limiting step in the translation process.

Synonymous genes explore different evolutionary landscapes

The evolutionary potential of a gene is constrained not only by the amino acid sequence of its product, but by its DNA sequence as well. The topology of the genetic code is such that half of the amino acids exhibit synonymous codons that can reach different subsets of amino acids from each other through single mutation. Thus, synonymous DNA sequences should access different regions of the protein sequence space through a limited number of mutations, and this may deeply influence the evolution of natural proteins. Here, we demonstrate that this feature can be of value for manipulating protein evolvability. We designed an algorithm that, starting from an input gene, constructs a synonymous sequence that systematically includes the codons with the most different evolutionary perspectives; i.e., codons that maximize accessibility to amino acids previously unreachable from the template by point mutation. A synonymous version of a bacterial antibiotic resistance gene was computed and synthesized. When concurrently submitted to identical directed evolution protocols, both the wild type and the recoded sequence led to the isolation of specific, advantageous phenotypic variants. Simulations based on a mutation isolated only from the synthetic gene libraries were conducted to assess the impact of sub-functional selective constraints, such as codon usage, on natural adaptation. Our data demonstrate that rational design of synonymous synthetic genes stands as an affordable improvement to any directed evolution protocol. We show that using two synonymous DNA sequences improves the overall yield of the procedure by increasing the diversity of mutants generated. These results provide conclusive evidence that synonymous coding sequences do experience different areas of the corresponding protein adaptive landscape, and that a sequence's codon usage effectively constrains the evolution of the encoded protein.

Hepatitis A virus mutant spectra under the selective pressure of monoclonal antibodies: codon usage constraints limit capsid variability

Severe structural constraints in the hepatitis A virus (HAV) capsid have been suggested as the reason for the lack of emergence of new serotypes in spite of the occurrence of complex distributions of mutants or quasispecies. Analysis of the HAV mutant spectra under immune pressure by the monoclonal antibodies (MAbs) K34C8 (immunodominant site) and H7C27 (glycophorin binding site) has revealed different evolutionary dynamics. Populations composed of complex ensembles of mutants with very low fitness or single dominant mutants with high fitness permit the acquisition of resistance to each of the MAbs, respectively. Deletion mutants were detected as components of the mutant spectra: up to 61 residues, with an average of 19, and up to 83 residues, with an average of 45, in VP3 and VP1 proteins, respectively. A clear negative selection of those replacements affecting the residues encoded by rare codons of the capsid surface has been detected through the present quasispecies analysis, confirming a certain beneficial role of such clusters. Since these clusters are located near or at the epitope regions, the need to maintain such clusters might prevent the emergence of new serotypes.

Modulation of poliovirus replicative fitness in HeLa cells by deoptimization of synonymous codon usage in the capsid region

We replaced degenerate codons for nine amino acids within the capsid region of the Sabin type 2 oral poliovirus vaccine strain with corresponding nonpreferred synonymous codons. Codon replacements were introduced into four contiguous intervals spanning 97% of the capsid region. In the capsid region of the most highly modified virus construct, the effective number of codons used (N(C)) fell from 56.2 to 29.8, the number of CG dinucleotides rose from 97 to 302, and the G+C content increased from 48.4% to 56.4%. Replicative fitness in HeLa cells, measured by plaque areas and virus yields in single-step growth experiments, decreased in proportion to the number of replacement codons. Plaque areas decreased over an approximately 10-fold range, and virus yields decreased over an approximately 65-fold range. Perhaps unexpectedly, the synthesis and processing of viral proteins appeared to be largely unaltered by the restriction in codon usage. In contrast, total yields of viral RNA in infected cells were reduced approximately 3-fold and specific infectivities of purified virions (measured by particle/PFU ratios) decreased approximately 18-fold in the most highly modified virus. The replicative fitness of both codon replacement viruses and unmodified viruses increased with the passage number in HeLa cells. After 25 serial passages (approximately 50 replication cycles), most codon replacements were retained, and the relative fitness of the modified viruses remained well below that of the unmodified virus. The increased replicative fitness of high-passage modified virus was associated with the elimination of several CG dinucleotides. Potential applications for the systematic modulation of poliovirus replicative fitness by deoptimization of codon usage are discussed.

tRNASer(CGA) differentially regulates expression of wild-type and codon-modified papillomavirus L1 genes

Exogenous transfer RNAs (tRNAs) favor translation of bovine papillomavirus 1 wild-type (wt) L1 mRNA in in vitro translation systems (Zhou et al. 1999, J. Virol., 73, 4972-4982). We, therefore, investigated whether papillomavirus (PV) wt L1 protein expression could be enhanced in eukaryotic cells following exogenous tRNA supplementation. Both Chinese hamster ovary (CHO) and Cos1 cells, transfected with PV1 wt L1 genes, effectively transcribed the genes but did not translate them. However, L1 protein translation was demonstrated following co-transfection with the L1 gene and a gene expressing tRNA(Ser)(CGA). Cell lines, stably transfected with a bovine papillomavirus 1 (BPV1) wt L1 expression construct, produced L1 protein after the transfection of the tRNA(Ser)(CGA) gene, but not following the transfection with basal vectors, suggesting that tRNA(Ser)(CGA) gene enhanced wt L1 translation as a result of endogenous tRNA alterations and phosphorylation of translation initiation factors elF4E and elF2alpha in the tRNA(Ser)(CGA) transfected L1 cell lines. The tRNA(Ser)(CGA) gene expression significantly reduced translation of L1 proteins expressed from codon-modified (HB) PV L1 genes utilizing mammalian preferred codons, but had variable effects on translation of green fluorescent proteins (GFPs) expressed from six serine GFP variants. The changes of tRNA pools appear to match the codon composition of PV wt and HB L1 genes and serine GFP variants to regulate translation of their mRNAs. These findings demonstrate for the first time in eukaryotic cells that translation of the target genes can be differentially influenced by the provision of a single tRNA expression construct.

Mammalian mutation pressure, synonymous codon choice, and mRNA degradation

The usage of synonymous codons (SCs) in mammalian genes is highly correlated with local base composition and is therefore thought to be determined by mutation pressure. The usage is nonetheless structured. For instance, mammals share with Saccharomyces and Drosophila most preferences for the C-ending over the G-ending codon (or vice versa) within each fourfold-degenerate SC family and the fact that their SCs are placed along coding regions in ways that minimize the number of T|A and C|G dinucleotides ("|" being the codon boundary). TA and CG underrepresentations are observed everywhere in the mammalian genome affecting the SC usage, the amino acid composition of proteins, and the primary structure of introns and noncoding DNA. While the rarity of CG is ascribed to the high mutability of this dinucleotide, the rarity of TA in coding regions is considered adaptive because UA dinucleotides are cleaved by endoribonucleases. Here we present in vivo experimental evidence indicating that the number of T|A and/or C|G dinucleotides of a human gene can affect strongly the expression level and degradation of its mRNA. Our results are consistent with indirect evidence produced by other workers and with the detailed work that has been devoted to characterize UA cleavage in vitro and in vivo. We conclude that SC choice can influence strongly mRNA function and gene expression through effects not directly related to the codon-anticodon interaction. These effects should constrain heavily the nucleotide motif composition of the most abundant mRNAs in the transcriptome, in particular, their SC usage, a usage that must be reflected by cellular tRNA concentrations and thus defines for all other genes which SCs are translated fastest and most accurately. Furthermore, the need to avoid such effects genome-wide appears serious enough to have favored the evolution of biases in context-dependent mutation that reduce the occurrence of intrinsically unfavorable motifs, and/or, when possible, to have induced the molecular machinery mediating such effects to rely opportunistically on already existing motif rarities and abundances. This may explain why nucleotide motif preferences are very similar in transcribed and nontranscribed mammalian DNA even though the preferences appear to be adaptive only in transcribed DNA.

Intercodon dinucleotides affect codon choice in plant genes

In this work, 710 CDSs corresponding to over 290 000 codons equally distributed between Brassica napus, Arabidopsis thaliana, Lycopersicon esculentum, Nicotiana tabacum, Pisum sativum, Glycine max, Oryza sativa, Triticum aestivum, Hordeum vulgare and Zea mays were considered. For each amino acid, synonymous codon choice was determined in the presence of A, G, C or T as the initial nucleotide of the subsequent triplet; data were statistically analysed under the hypothesis of an independent assortment of codons. In 33.4% of cases, a frequency significantly (P: = 0.01) different from that expected was recorded. This was mainly due to a pervasive intercodon TpA and CpG deficiency. As a general rule, intercodon TpAs and CpGs were preferably replaced by CpAs and TpGs, respectively. In several instances, codon frequencies were also modified to avoid homotetramer and homotrimer formation, to reduce intercodon ApCs downstream (1,2) GG or AG dinucleotides, as well as to increase GpA or ApG intercodons under certain contexts. Since TpA, CpG and homotetra(tri)mer deficiency directly or indirectly accounted for 77% of significant variation in the codon frequency, it can be concluded that codon usage mirrors precise needs at the DNA structure level. Plant species exhibited a phylogenetically-related adaptation to structural constraints. Codon usage flexibility was reflected in strikingly different arrays of optimum codons for probe design.

Stp1p, Stp2p and Abf1p are involved in regulation of expression of the amino acid transporter gene BAP3 of Saccharomyces cerevisiae

Expression of the BAP3 gene of Saccharomyces cerevisiae, encoding a branched chain amino acid permease, is induced in response to the availability of several naturally occurring amino acids in the medium. This induction is mediated via an upstream activating sequence (called UAS(aa)) in the BAP3 promoter, and dependent on Stp1p, a nuclear protein with zinc finger domains, suggesting that Stp1p is a transcription factor involved in BAP3 expression. In this paper, we show that Stp2p, a protein with considerable similarity to Stp1p, is also involved in the induction of BAP3 expression. To gain more insight into the roles of STP1 and STP2, we have overexpressed both Stp1p and Stp2p in yeast cells. Gel shift assays with the UAS(aa)as a probe show that the UAS(aa)can form two major complexes. One complex is dependent on Stp2p overexpression and the other is formed independently of STP1 or STP2, suggesting that the UAS(aa)is also bound by another factor. Here we show that the other factor is Abf1p, which binds specifically to the UAS(aa)of BAP3.

Expression of giant silkmoth cecropin B genes in tobacco

Cecropin B is a small antibacterial peptide from the giant silkmoth Hyalophora cecropia. To reveal the potential of this peptide for engineering bacterial disease resistance into crops, several cecropin B gene constructs were made either for expression in the cytosol or for secretion. All constructs were cloned in a plant expression vector and introduced in tobacco via Agrobacterium tumefaciens. A cDNA-derived cecropin B gene construct lacking the amino-terminal signal peptide was poorly expressed in transgenic plants at the mRNA level, whereas plants harbouring a full-length cDNA-derived construct containing the insect signal peptide, showed increased cecropin B-mRNA levels. Highest expression was found in plants harbouring a construct with a plant-gene-derived signal peptide. In none of the transgenic plants could the cecropin B peptide be detected. This is most likely caused by breakdown of the peptide by plant endogenous proteases, since a chemically synthesized cecropin B peptide was degraded within seconds in various plant cell extracts. This degradation could be prevented by the addition of specific protease inhibitors and by boiling the extract prior to adding the peptide. In addition, anionic detergents, in contrast to cationic, zwitter-ionic or non-ionic detergents, could prevent this degradation. Nevertheless, transgenic tobacco plants were evaluated for resistance to Pseudomonas solanacearum, the causal agent of bacterial wilt of many crops, and P. syringae pv. tabaci, the causal agent of bacterial wildfire, which are highly susceptible to cecropin B in vitro. No resistance was found. These experiments indicate that introduction and expression of cecropin B genes in tobacco does not result in detectable cecropin B protein levels and resistance to bacterial infections, most likely due to degradation of the protein by endogenous proteases.

Initiation of translation can occur only in a restricted region of the CYC1 mRNA of Saccharomyces cerevisiae

The steady-state levels and half-lives of CYC1 mRNAs were estimated in a series of mutant strains of Saccharomyces cerevisiae containing (i) TAA nonsense codons, (ii) ATG initiator codons, or (iii) the sequence ATA ATG ACT TAA (denoted ATG-TAA) at various positions along the CYC1 gene, which encodes iso-1-cytochrome c. These mutational alterations were made in backgrounds lacking all internal in-frame and out-of-frame ATG triplets or containing only one ATG initiator codon at the normal position. The results revealed a "sensitive" region encompassing approximately the first half of the CYC1 mRNA, in which nonsense codons caused Upf1-dependent degradation. This result and the stability of CYC1 mRNAs lacking all ATG triplets, as well as other results, suggested that degradation occurs unless elements associated with this sensitive region are covered with 80S ribosomes, 40S ribosomal subunits, or ribonucleoprotein particle proteins. While elongation by 80S ribosomes could be prematurely terminated by TAA codons, the scanning of 40S ribosomal units could not be terminated solely by TAA codons but could be disrupted by the ATG-TAA sequence, which caused the formation and subsequent prompt release of 80S ribosomes. The ATG-TAA sequence caused degradation of the CYC1 mRNA only when it was in the region spanning nucleotide positions -27 to +37 but not in the remaining 3' distal region, suggesting that translation could initiate only in this restricted initiation region. CYC1 mRNA distribution on polyribosomes confirmed that only ATG codons within the initiation region were translated at high efficiency. This initiation region was not entirely dependent on the distance from the 5' cap site and was not obviously dependent on the short-range secondary structure but may simply reflect an open structural requirement for initiation of translation of the CYC1 mRNA.

Expression enhancement of the Tn5 neomycin-resistance gene by removal of upstream ATG sequences and its use for probing heterologous upstream activating sequences in yeast

We have constructed a series of promoter or upstream activating sequence (UAS)-probe plasmids carrying the Tn5-derived neomycin resistance gene whose seven additional ATG codons in the 5'-untranslated region were completely or partially removed. When the deleted version of the neo sequence retaining only one additional ATG (NeoD) was expressed under the control of a TDH3 promoter whose UAS was deleted, the transformed cells were unable to grow at a low concentration of the antibiotic G418. In contrast with this, yeast cells expressing the NeoC sequence and having no additional ATG exhibited a high level of G418-resistance. Moreover, the UAS-probe system using NeoD has been successfully applied for the identification of several E. coli DNA sequences that clearly function as UASs in yeast cells. Two of these prokaryotic sequences with UAS activity were identified as a part of the coding region of the tgt and the hydG gene, respectively.

Identification, sequence analysis, and expression of a Corynebacterium glutamicum gene cluster encoding the three glycolytic enzymes glyceraldehyde-3-phosphate dehydrogenase, 3-phosphoglycerate kinase, and triosephosphate isomerase

To investigate a possible chromosomal clustering of glycolytic enzyme genes in Corynebacterium glutamicum, a 6.4-kb DNA fragment located 5' adjacent to the structural phosphoenolpyruvate carboxylase (PEPCx) gene ppc was isolated. Sequence analysis of the ppc-proximal part of this fragment identified a cluster of three glycolytic genes, namely, the glyceraldehyde-3-phosphate dehydrogenase (GAPDH) gene gap, the 3-phosphoglycerate kinase (PGK) gene pgk, and the triosephosphate isomerase (TPI) gene tpi. The four genes are organized in the order gap-pgk-tpi-ppc and are separated by 215 bp (gap and pgk), 78 bp (pgk and tpi), and 185 bp (tpi and ppc). The predicted gene product of gap consists of 336 amino acids (M(r) of 36,204), that of pgk consists of 403 amino acids (M(r) of 42,654), and that of tpi consists of 259 amino acids (M(r) of 27,198). The amino acid sequences of the three enzymes show up to 62% (GAPDH), 48% (PGK), and 44% (TPI) identity in comparison with respective enzymes from other organisms. The gap, pgk, tpi, and ppc genes were cloned into the C. glutamicum-Escherichia coli shuttle vector pEK0 and introduced into C. glutamicum. Relative to the wild type, the recombinant strains showed up to 20-fold-higher specific activities of the respective enzymes. On the basis of codon usage analysis of gap, pgk, tpi, and previously sequenced genes from C. glutamicum, a codon preference profile for this organism which differs significantly from those of E. coli and Bacillus subtilis is presented.

Synthetic cryIIIA gene from Bacillus thuringiensis improved for high expression in plants

A 1974 bp synthetic gene was constructed from chemically synthesized oligonucleotides in order to improve transgenic protein expression of the cryIIIA gene from Bacillus thuringiensis var. tenebrionis in transgenic tobacco. The crystal toxin genes (cry) from B. thuringiensis are difficult to express in plants even when under the control of efficient plant regulatory sequences. We identified and eliminated five classes of sequence found throughout the cryIIIA gene that mimic eukaryotic processing signals and which may be responsible for the low levels of transcription and translation. Furthermore, the GC content of the gene was raised from 36% to 49% and the codon usage was changed to be more plant-like. When the synthetic gene was placed behind the cauliflower mosaic virus 35S promoter and the alfalfa mosaic virus translational enhancer, up to 0.6% of the total protein in transgenic tobacco plants was cryIIIA as measured from immunoblot analysis. Bioassay data using potato beetle larvae confirmed this estimate.

The rate-limiting step in yeast PGK1 mRNA degradation is an endonucleolytic cleavage in the 3'-terminal part of the coding region

Insertion of an 18-nucleotide-long poly(G) tract into the 3'-terminal untranslated region of yeast phosphoglycerate kinase (PGK1) mRNA increases its chemical half-life by about a factor of 2 (P. Vreken, R. Van der Veen, V. C. H. F. de Regt, A. L. de Maat, R. J. Planta, and H. A. Raué, Biochimie 73:729-737, 1991). In this report, we show that this insertion also causes the accumulation of a degradation intermediate extending from the poly(G) sequence down to the transcription termination site. Reverse transcription and S1 nuclease mapping experiments demonstrated that this intermediate is the product of shorter-lived primary fragments resulting from endonucleolytic cleavage immediately downstream from the U residue of either of two 5'-GGUG-3' sequences present between positions 1100 and 1200 close to the 3' terminus (position 1251) of the coding sequence. Similar endonucleolytic cleavages appear to initiate degradation of wild-type PGK1 mRNA. Insertion of a poly(G) tract just upstream from the AUG start codon resulted in the accumulation of a 5'-terminal degradation intermediate extending from the insertion to the 1100-1200 region. RNase H degradation in the presence of oligo(dT) demonstrated that the wild-type and mutant PGK1 mRNAs are deadenylated prior to endonucleolytic cleavage and that the half-life of the poly(A) tail is three- to sixfold lower than that of the remainder of the mRNA. Thus, the endonucleolytic cleavage constitutes the rate-limiting step in degradation of both wild-type and mutant PGK1 transcripts, and the resulting fragments are degraded by a 5'----3' exonuclease, which appears to be severely retarded by a poly(G) sequence.

The rate-limiting step in yeast PGK1 mRNA degradation is an endonucleolytic cleavage in the 3'-terminal part of the coding region

Insertion of an 18-nucleotide-long poly(G) tract into the 3'-terminal untranslated region of yeast phosphoglycerate kinase (PGK1) mRNA increases its chemical half-life by about a factor of 2 (P. Vreken, R. Van der Veen, V. C. H. F. de Regt, A. L. de Maat, R. J. Planta, and H. A. Raué, Biochimie 73:729-737, 1991). In this report, we show that this insertion also causes the accumulation of a degradation intermediate extending from the poly(G) sequence down to the transcription termination site. Reverse transcription and S1 nuclease mapping experiments demonstrated that this intermediate is the product of shorter-lived primary fragments resulting from endonucleolytic cleavage immediately downstream from the U residue of either of two 5'-GGUG-3' sequences present between positions 1100 and 1200 close to the 3' terminus (position 1251) of the coding sequence. Similar endonucleolytic cleavages appear to initiate degradation of wild-type PGK1 mRNA. Insertion of a poly(G) tract just upstream from the AUG start codon resulted in the accumulation of a 5'-terminal degradation intermediate extending from the insertion to the 1100-1200 region. RNase H degradation in the presence of oligo(dT) demonstrated that the wild-type and mutant PGK1 mRNAs are deadenylated prior to endonucleolytic cleavage and that the half-life of the poly(A) tail is three- to sixfold lower than that of the remainder of the mRNA. Thus, the endonucleolytic cleavage constitutes the rate-limiting step in degradation of both wild-type and mutant PGK1 transcripts, and the resulting fragments are degraded by a 5'----3' exonuclease, which appears to be severely retarded by a poly(G) sequence.

A cell-free extract from yeast cells for studying mRNA turnover

We have isolated a cell-free extract from yeast cells that reproduces the differences observed in vivo in the rate of turnover of individual yeast mRNAs. Detailed analysis of the degradation of yeast phosphoglycerate kinase (PGK) mRNA in this system demonstrated that both natural and synthetically prepared PGK transcripts are degraded by the same pathway previously established by us in vivo, consisting of endonucleolytic cleavage at a number of 5'-GGUG-3' sequence motifs within a short target region located close to the 3'-end of the coding sequence followed by 5'-3' exonucleolytic removal of the resulting fragments. The extract, therefore, is suitable for studying the mechanistic details of mRNA turnover in yeast. As a first application of this system we have performed a limited mutational analysis of two of the GGUG motifs within the endonucleolytic target region of the PGK transcript. The results show that sequence changes in either motif abolish cleavage at the mutated site only, indicating the involvement of the residues in question in selection of the cleavage positions.

Analysis of CpG suppression in methylated and nonmethylated species

The development of nearest-neighbor analysis led to the finding that the frequency of the dinucleotide CpG is markedly depressed in vertebrates. One explanation of this suppression is that methylation of CpG found in vertebrates represents a mutational hot spot through deamination of methylcytidine to thymidine. We have examined the role of methylated CpG as a factor in CpG suppression by comparing CpG distributions in coding regions of 121 genes from six species, three with methylated DNA and three with nonmethylated DNA. Overall base composition shows that all species exhibit CpG suppression, with the methylated forms showing significantly greater suppression than nonmethylated forms. When the data are analyzed by CpG position, the mean values of the methylated forms exhibit greater suppression than nonmethylated forms at positions I-II and II-III, but there is considerable overlap of suppression scores for individual species. At position III-I, CpG suppression is marked in all methylated species, and it is reversed in all nonmethylated species. Our analysis supports the hypothesis that CpG patterns at positions II-III and III-I in methylated forms are affected by mutation acting through deamination of methylcytidine to thymidine. We speculate that the excess of CpGs at position III-I in nonmethylated forms may be related to a requirement for minimal thermal stability of the DNA.

Directional mutation pressure, selective constraints, and genetic equilibria

Rates of substitution mutations in two directions, v [from an A-T or T-A nucleotide pair (AT-pair) to a G-C or C-G nucleotide pair (GC-pair)] and u [from a GC-pair to an AT-pair], are usually not the same. The net effect, v/(u + v), has previously been defined as directional mutation pressure (mu D), which explains the wide interspecific variation and narrow intragenomic heterogeneity of DNA G + C content in bacteria. In this article, first, a theory of the evolution of DNA G + C content is presented that is based on the equilibrium among three components: directional mutation pressure, DNA G + C content, and selective constraints. According to this theory, consideration of both u and v as well as selective constraints is essential to explain the molecular evolution of the DNA base composition and sequence. Second, the theory of directional mutation pressure is applied to the analysis of the wide intragenomic heterogeneity of DNA G + C content in multicellular eukaryotes. The theory explains the extensive intragenomic heterogeneity of G + C content of higher eukaryotes primarily as the result of the intragenomic differences of directional mutation pressure and selective constraints rather than the result of positive selections for functional advantages of the DNA G + C content itself.

Cloning and sequencing of the sacA gene: characterization of a sucrase from Zymomonas mobilis

The Zymomonas mobilis gene (sacA) encoding a protein with sucrase activity has been cloned in Escherichia coli and its nucleotide sequence has been determined. Potential ribosome-binding site and promoter sequences were identified in the region upstream of the gene which were homologous to E. coli and Z. mobilis consensus sequences. Extracts from E. coli cells, containing the sacA gene, displayed a sucrose-hydrolyzing activity. However, no transfructosylation activity (exchange reaction or levan formation) could be detected. This sucrase activity was different from that observed with the purified extracellular protein B46 from Z. mobilis. These two proteins showed different electrophoretic mobilities and molecular masses and shared no immunological similarity. Thus, the product of sacA (a polypeptide of 58.4-kDa molecular mass) is a new sucrase from Z. mobilis. The amino acid sequence, deduced from the nucleotide sequence of sacA, showed strong homologies with the sucrases from Bacillus subtilis, Salmonella typhimurium, and Vibrio alginolyticus.

G418-resistance as a dominant marker and reporter for gene expression in Saccharomyces cerevisiae

Coding sequence cartridges for aminoglycoside phosphotransferase (APT) were isolated from bacterial transposon Tn903. When incorporated into a heterologous gene construction utilising the PGK1 promoter and terminator, the heterologous APT gene provided a G418-resistance determinant that functioned efficiently as a dominant marker for yeast in both multiple- and single-copy. Transformant colonies on selective medium appeared rapidly, within 36-48 h, and growth rate of the transformed cells was normal. A simple and highly sensitive radiolabelling assay for APT enzyme activity was developed for use with crude cell protein extracts. Enzyme activity units were equated to the amount of APT protein present in the cells, and the APT protein was shown to be stable in yeast. Heterologous APT expression was 130-fold reduced compared with homologous PGK1. This resulted from an estimated two-fold decrease in mRNA level and a 65-fold decrease in translation efficiency. The latter was unaffected by AUG sequence context change, but corresponded with a high frequency of minor codons in the APT-coding sequence. APT can be used as a semi-quantitative reporter of gene expression, whose useful features are in vivo detection via the G418-resistance phenotype and powerful cell-free assay.

Translation to near the distal end of the penultimate exon is required for normal levels of spliced triosephosphate isomerase mRNA

The translation of human triosephosphate isomerase (TPI) mRNA normally terminates at codon 249 within exon 7, the final exon. Frameshift and nonsense mutations within the TPI gene that cause translation to terminate prematurely at or upstream of codon 189, within exon 6, result in a decreased level of TPI mRNA (I.O. Daar and L.E. Maquat, Mol. Cell. Biol. 8:802-813, 1988). For all mutations in this group, the decrease is to the same extent, i.e., to approximately 20% of the normal level. We show here that a second group of nonsense mutations that cause translation to terminate prematurely at or downstream of codon 208, in exon 6, did not affect TPI mRNA abundance. Deletion analysis demonstrated that the abundance of translationally active TPI mRNA is a function of both the distance and the polarity of the nonsense codon relative to the final intron in TPI pre-mRNA. Our results indicate that if translating ribosomes are unable to progress to at least a certain position within the penultimate exon relative to the final intron, then the level of the corresponding mRNA will be abnormally low. Studies inhibiting RNA synthesis with dactinomycin demonstrated that a block in translation does not affect the half-life of mature TPI mRNA. The simplest interpretation of our data is that the translation of TPI mRNA in the cytoplasm facilitates the splicing of TPI pre-mRNA or the transport of TPI mRNA across the nuclear envelope or both.

Translation to near the distal end of the penultimate exon is required for normal levels of spliced triosephosphate isomerase mRNA

The translation of human triosephosphate isomerase (TPI) mRNA normally terminates at codon 249 within exon 7, the final exon. Frameshift and nonsense mutations within the TPI gene that cause translation to terminate prematurely at or upstream of codon 189, within exon 6, result in a decreased level of TPI mRNA (I.O. Daar and L.E. Maquat, Mol. Cell. Biol. 8:802-813, 1988). For all mutations in this group, the decrease is to the same extent, i.e., to approximately 20% of the normal level. We show here that a second group of nonsense mutations that cause translation to terminate prematurely at or downstream of codon 208, in exon 6, did not affect TPI mRNA abundance. Deletion analysis demonstrated that the abundance of translationally active TPI mRNA is a function of both the distance and the polarity of the nonsense codon relative to the final intron in TPI pre-mRNA. Our results indicate that if translating ribosomes are unable to progress to at least a certain position within the penultimate exon relative to the final intron, then the level of the corresponding mRNA will be abnormally low. Studies inhibiting RNA synthesis with dactinomycin demonstrated that a block in translation does not affect the half-life of mature TPI mRNA. The simplest interpretation of our data is that the translation of TPI mRNA in the cytoplasm facilitates the splicing of TPI pre-mRNA or the transport of TPI mRNA across the nuclear envelope or both.

Codon preferences in free-living microorganisms

A popular interpretation of the major codon preference is that it reflects the operation of a regulatory device that controls the expression of individual proteins. In this popular model, rapidly translated codons are thought to promote the accumulation of the highly expressed proteins and slowly translated codons are thought to retard the expression of poorly expressed proteins. However, this widely accepted model is not supported by kinetic theory or by experimental results. A less fashionable model in which the major codon preference has nothing to do with the expression level of the individual proteins is forwarded. In this model, the major codon preference is viewed as a global strategy to support the efficient function of the translation system and thereby to maximize the growth rates of cells under favorable conditions.

Identification and comparison of stable and unstable mRNAs in Saccharomyces cerevisiae

We developed a procedure to measure mRNA decay rates in the yeast Saccharomyces cerevisiae and applied it to the determination of half-lives for 20 mRNAs encoded by well-characterized genes. The procedure utilizes Northern (RNA) or dot blotting to quantitate the levels of individual mRNAs after thermal inactivation of RNA polymerase II in an rpb1-1 temperature-sensitive mutant. We compared the results of this procedure with results obtained by two other procedures (approach to steady-state labeling and inhibition of transcription with Thiolutin) and also evaluated whether heat shock alter mRNA decay rates. We found that there are no significant differences in the mRNA decay rates measured in heat-shocked and non-heat-shocked cells and that, for most mRNAs, different procedures yield comparable relative decay rates. Of the 20 mRNAs studied, 11, including those encoded by HIS3, STE2, STE3, and MAT alpha 1, were unstable (t1/2 less than 7 min) and 4, including those encoded by ACT1 and PGK1, were stable (t1/2 greater than 25 min). We have begun to assess the basis and significance of such differences in the decay rates of these two classes of mRNA. Our results indicate that (i) stable and unstable mRNAs do not differ significantly in their poly(A) metabolism; (ii) deadenylation does not destabilize stable mRNAs; (iii) there is no correlation between mRNA decay rate and mRNA size; (iv) the degradation of both stable and unstable mRNAs depends on concomitant translational elongation; and (v) the percentage of rare codons present in most unstable mRNAs is significantly higher than in stable mRNAs.

Identification and comparison of stable and unstable mRNAs in Saccharomyces cerevisiae

We developed a procedure to measure mRNA decay rates in the yeast Saccharomyces cerevisiae and applied it to the determination of half-lives for 20 mRNAs encoded by well-characterized genes. The procedure utilizes Northern (RNA) or dot blotting to quantitate the levels of individual mRNAs after thermal inactivation of RNA polymerase II in an rpb1-1 temperature-sensitive mutant. We compared the results of this procedure with results obtained by two other procedures (approach to steady-state labeling and inhibition of transcription with Thiolutin) and also evaluated whether heat shock alter mRNA decay rates. We found that there are no significant differences in the mRNA decay rates measured in heat-shocked and non-heat-shocked cells and that, for most mRNAs, different procedures yield comparable relative decay rates. Of the 20 mRNAs studied, 11, including those encoded by HIS3, STE2, STE3, and MAT alpha 1, were unstable (t1/2 less than 7 min) and 4, including those encoded by ACT1 and PGK1, were stable (t1/2 greater than 25 min). We have begun to assess the basis and significance of such differences in the decay rates of these two classes of mRNA. Our results indicate that (i) stable and unstable mRNAs do not differ significantly in their poly(A) metabolism; (ii) deadenylation does not destabilize stable mRNAs; (iii) there is no correlation between mRNA decay rate and mRNA size; (iv) the degradation of both stable and unstable mRNAs depends on concomitant translational elongation; and (v) the percentage of rare codons present in most unstable mRNAs is significantly higher than in stable mRNAs.

Cloning of the Zymomonas mobilis structural gene encoding alcohol dehydrogenase I (adhA): sequence comparison and expression in Escherichia coli

Zymomonas mobilis ferments sugars to produce ethanol with two biochemically distinct isoenzymes of alcohol dehydrogenase. The adhA gene encoding alcohol dehydrogenase I has now been sequenced and compared with the adhB gene, which encodes the second isoenzyme. The deduced amino acid sequences for these gene products exhibited no apparent homology. Alcohol dehydrogenase I contained 337 amino acids, with a subunit molecular weight of 36,096. Based on comparisons of primary amino acid sequences, this enzyme belongs to the family of zinc alcohol dehydrogenases which have been described primarily in eucaryotes. Nearly all of the 22 strictly conserved amino acids in this group were also conserved in Z. mobilis alcohol dehydrogenase I. Alcohol dehydrogenase I is an abundant protein, although adhA lacked many of the features previously reported in four other highly expressed genes from Z. mobilis. Codon usage in adhA is not highly biased and includes many codons which were unused by pdc, adhB, gap, and pgk. The ribosomal binding region of adhA lacked the canonical Shine-Dalgarno sequence found in the other highly expressed genes from Z. mobilis. Although these features may facilitate the expression of high enzyme levels, they do not appear to be essential for the expression of Z. mobilis adhA.

Translation and a 42-nucleotide segment within the coding region of the mRNA encoded by the MAT alpha 1 gene are involved in promoting rapid mRNA decay in yeast

In yeast, the mRNA encoded by the MAT alpha 1 gene is unstable (t1/2 = 5 min) and the mRNAs encoded by the ACT1 gene (t1/2 = 30 min) and the PGK1 gene (t1/2 = 45 min) are stable. To understand the RNA structural features that dictate mRNA decay rates in yeast, we have constructed PGK1/MAT alpha 1 and ACT1/MAT alpha 1 gene fusions and analyzed the decay rates of the resultant chimeric transcripts. Fusion of a MAT alpha 1 segment containing 73% of the coding region and the 3' untranslated region to either of the stable genes is sufficient to cause rapid decay of the chimeric mRNAs (t1/2 = 6-7.5 min). Sequences required for this rapid decay are not found in the MAT alpha 1 3' untranslated region but are located within a 42-nucleotide segment of the coding region that has a high content (8 out of 14) of rare codons. Introduction of a translational stop codon upstream of this region stabilizes the hybrid mRNAs, indicating that the rapid decay promoted by these sequences is dependent on ribosomal translocation.

Suppression of the negative effect of minor arginine codons on gene expression; preferential usage of minor codons within the first 25 codons of the Escherichia coli genes

AGA and AGG codons for arginine are the least used codons in Escherichia coli, which are encoded by a rare tRNA, the product of the dnaY gene. We examined the positions of arginine residues encoded by AGA/AGG codons in 678 E. coli proteins. It was found that AGA/AGG codons appear much more frequently within the first 25 codons. This tendency becomes more significant in those proteins containing only one AGA or AGG codon. Other minor codons such as CUA, UCA, AGU, ACA, GGA, CCC and AUA are also found to be preferentially used within the first 25 codons. The effects of the AGG codon on gene expression were examined by inserting one to five AGG codons after the 10th codon from the initiation codon of the lacZ gene. The production of beta-galactosidase decreased as more AGG codons were inserted. With five AGG codons, the production of beta-galactosidase (Gal-AGG5) completely ceased after a mid-log phase of cell growth. After 22 hr induction of the lacZ gene, the overall production of Gal-AGG5 was 11% of the control production (no insertion of arginine codons). When five CGU codons, the major arginine codon were inserted instead of AGG, the production of beta-galactosidase (Gal-CGU5) continued even after stationary phase and the overall production was 66% of the control. The negative effect of the AGG codons on the Gal-AGG5 production was found to be dependent upon the distance between the site of the AGG codons and the initiation codon. As the distance was increased by inserting extra sequences between the two codons, the production of Gal-AGG5 increased almost linearly up to 8 fold. From these results, we propose that the position of the minor codons in an mRNA plays an important role in the regulation of gene expression possibly by modulating the stability of the initiation complex for protein synthesis.

Regulation of the nuclear genes encoding the cytoplasmic and mitochondrial leucyl-tRNA synthetases of Neurospora crassa

We show that the nuclear genes for the cytoplasmic and mitochondrial leucyl-tRNA synthetase (LeuRS) of Neurospora crassa are distinct in their encoded proteins, codon usage, mRNA levels, and regulation. The 4.2-kilobase-pair region representing the structural gene for cytoplasmic LeuRS and flanking regions has been sequenced. The positions of the 5' and 3' ends of mRNA and of a single 62-base-pair intron have been mapped. The methionine-initiated open reading frame encoded a protein of 1,123 amino acids and displayed a strong codon bias. Although cytoplasmic LeuRS shares with mitochondrial LeuRS some general features common to most aminoacyl-tRNA synthetases, there is little amino acid sequence similarity between them, mRNA levels for cytoplasmic LeuRS were much higher than those for mitochondrial LeuRS. This observation and the strong codon bias in the cytoplasmic LeuRS gene may contribute to a greater abundance of cytoplasmic LeuRS than mitochondrial LeuRS. The genes for cytoplasmic and mitochondrial LeuRS are regulated independently. The cytoplasmic LeuRS gene is regulated by the cross-pathway control system in N. crassa, which is analogous to general amino acid control in Saccharomyces cerevisiae. The cytoplasmic LeuRS mRNA levels are induced by amino acid starvation resulting from the addition of aminotriazole. Part of this increase is due to utilization of new transcription start sites. In contrast, the mitochondrial LeuRS gene is not induced by amino acid limitation. However, the mitochondrial LeuRS mRNA levels did increase dramatically upon inhibition of mitochondrial protein synthesis by chloramphenicol or ethidium bromide or in the temperature-sensitive strain leu-5 carrying a mutation in the mitochondrial LeuRS structural gene.

Relationship between protein synthesis and concentrations of charged and uncharged tRNATrp in Escherichia coli

We have continuously monitored Trp-tRNA(Trp) concentrations in vivo and, in the same cultures, measured rates of protein synthesis in isogenic stringent and relaxed strains. We have also manipulated cellular charged and uncharged [tRNA(Trp)] by two means: (i) the strain used contains a Trp-tRNA synthetase mutation that increases the Km for Trp; thus, varying exogenous Trp varies cellular Trp-tRNA(Trp); and (ii) we have introduced into the mutant strain a plasmid containing the tRNA(Trp) gene behind an inducible promoter; thus, total [tRNA(Trp)] also can be varied depending on length of induction. The use of these conditions, combined with a previously characterized assay system, has allowed us to demonstrate that (i) the rate of incorporation of Trp into protein is proportional to the fraction of tRNA(Trp) that is charged; for any given total [tRNA(Trp)], this rate is also proportional to the [Trp-tRNA(Trp)]; (ii) uncharged tRNA(Trp) inhibits incorporation of Trp into protein; and (iii) rates of incorporation into protein of at least two other amino acids, Lys and Cys, are also sensitive to [Trp-tRNA(Trp)] and are inhibited by uncharged tRNA(Trp). Our results are consistent with models of translational control that postulate modulating polypeptide chain elongation efficiency by varying concentrations of specific tRNAs.