Codon usage biases co-evolve with transcription termination machinery to suppress premature cleavage and polyadenylation

DNA Sequence Changes Resulting from Codon Optimization Affect Gene Expression in Pichia pastoris by Altering Chromatin Accessibility

Codon optimization is a widely employed strategy to enhance protein expression. However, it occasionally leads to unexpected transcriptional repression despite preserving amino acid sequences. This study investigates the mechanistic basis of such transcriptional attenuation by analyzing two gene candidates (0432 and Fluc) in the common expression chassis P. pastoris. Both genes experienced severe mRNA reduction following codon optimization. Evidenced by histone H3 chromatin immunoprecipitation (ChIP) and a DNase I hypersensitivity assay, gene sequences with transcriptional repression displayed elevated nucleosome occupancy and reduced chromatin accessibility. The above change was caused by an ORF sequence change independent of the promoter, since transcriptional attenuation and compromised chromatin accessibility were still observed after replacing the strong promoter PGAP with Ppor1 or Prps8b. Our findings challenge the conventional view of codon optimization as solely translation-centric, revealing its capacity to preemptively modulate transcription through chromatin accessibility. This work underscores the necessity of integrating chromatin-level considerations into synthetic gene design to avoid unintended transcriptional silencing and optimize expression outcomes.

Histone deacetylase-1 is required for epigenome stability in Neurospora crassa

Polycomb group (PcG) proteins form chromatin modifying complexes that stably repress lineage- or context-specific genes in animals, plants, and some fungi. Polycomb Repressive Complex 2 (PRC2) catalyzes trimethylation of lysine 27 on histone H3 (H3K27me3) to assemble repressive chromatin. In the model fungus Neurospora crassa, H3K27me3 deposition is controlled by the H3K36 methyltransferase ASH1 and components of constitutive heterochromatin including the H3K9me3-binding protein HETEROCHROMATIN PROTEIN 1 (HP1). Hypoacetylated histones are a defining feature of both constitutive heterochromatin and PcG-repressed chromatin, but how histone deacetylases (HDACs) contribute to normal H3K27me3 and transcriptional repression within PcG-repressed chromatin is poorly understood. We performed a genetic screen to identify HDACs required for repression of PRC2-methylated genes. In the absence of HISTONE DEACETYLASE-1 (HDA-1), PRC2-methylated genes were activated and H3K27me3 was depleted from typical PRC2-targeted regions. At constitutive heterochromatin, HDA-1 deficient cells displayed reduced H3K9me3, hyperacetylation, and aberrant enrichment of H3K27me3 and H3K36me3. CHROMODOMAIN PROTEIN-2 (CDP-2) is required to target HDA-1 to constitutive heterochromatin and was also required for normal H3K27me3 patterns. Patterns of aberrant H3K27me3 were distinct in isogenic Δhda-1 strains, suggesting that loss of HDA-1 causes stochastic or progressive epigenome dysfunction. To test this, we constructed a new Δhda-1 strain and performed a laboratory evolution experiment. Deletion of hda-1 led to progressive epigenome decay over hundreds of nuclear divisions. Together, our data indicate that HDA-1 is a critical regulator of epigenome stability in N. crassa.

SARS-CoV-2 Displays a Suboptimal Codon Usage Bias for Efficient Translation in Human Cells Diverted by Hijacking the tRNA Epitranscriptome

Codon bias analysis of SARS-CoV-2 reveals suboptimal adaptation for translation in human cells it infects. The detailed examination of the codons preferentially used by SARS-CoV-2 shows a strong preference for LysAAA, GlnCAA, GluGAA, and ArgAGA, which are infrequently used in human genes. In the absence of an adapted tRNA pool, efficient decoding of these codons requires a 5-methoxycarbonylmethyl-2-thiouridine (mcm5s2) modification at the U34 wobble position of the corresponding tRNAs (tLysUUU; tGlnUUG; tGluUUC; tArgUCU). The optimal translation of SARS-CoV-2 open reading frames (ORFs) may therefore require several adjustments to the host's translation machinery, enabling the highly biased viral genome to achieve a more favorable "Ready-to-Translate" state in human cells. Experimental approaches based on LC-MS/MS quantification of tRNA modifications and on alteration of enzymatic tRNA modification pathways provide strong evidence to support the hypothesis that SARS-CoV-2 induces U34 tRNA modifications and relies on these modifications for its lifecycle. The conclusions emphasize the need for future studies on the evolution of SARS-CoV-2 codon bias and its ability to alter the host tRNA pool through the manipulation of RNA modifications.

CbCyp51-Mediated Demethylation Inhibitor Resistance Is Modulated by Codon Bias

Cercospora leaf spot, caused by the fungus Cercospora beticola, is the most destructive foliar disease of sugarbeet worldwide. Resistance to the sterol demethylation inhibitor (DMI) fungicide tetraconazole has been previously correlated with synonymous and nonsynonymous mutations in CbCyp51. Here, we extend these analyses to the DMI fungicides prothioconazole, difenoconazole, and mefentrifluconazole in addition to tetraconazole to confirm whether the synonymous and nonsynonymous mutations at amino acid positions 144 and 170 are associated with resistance to these fungicides. Nearly half of the 593 isolates of C. beticola collected in the Red River Valley of North Dakota and Minnesota in 2021 were resistant to all four DMIs. Another 20% were resistant to tetraconazole and prothioconazole but sensitive to difenoconazole and mefentrifluconazole. A total of 13% of isolates were sensitive to all DMIs tested. We found five CbCyp51 haplotypes and associated them with phenotypes to the four DMIs. The most predominant haplotype (E170_A/L144F_C) correlated with resistance to all four DMIs with up to 97.6% accuracy. The second most common haplotype (E170_A/L144) consisted of isolates associated with resistance phenotypes to tetraconazole and prothioconazole while also exhibiting sensitive phenotypes to difenoconazole and mefentrifluconazole with up to 98.4% accuracy. Quantitative PCR did not identify differences in CbCyp51 expression between haplotypes. This study offers an understanding of the importance of codon usage in fungicide resistance and provides crop management acuity for fungicide application decision-making.

Effective Synthesis of High-Integrity mRNA Using In Vitro Transcription

mRNA vaccines are entering a period of rapid development. However, their synthesis is still plagued by challenges related to mRNA impurities and fragments (incomplete mRNA). Most impurities of mRNA products transcribed in vitro are mRNA fragments. Only full-length mRNA transcripts containing both a 5'-cap and a 3'-poly(A) structure are viable for in vivo expression. Therefore, RNA fragments are the primary product-related impurities that significantly hinder mRNA efficacy and must be effectively controlled; these species are believed to originate from either mRNA hydrolysis or premature transcriptional termination. In the manufacturing of commercial mRNA vaccines, T7 RNA polymerase-catalyzed in vitro transcription (IVT) synthesis is a well-established method for synthesizing long RNA transcripts. This study identified a pivotal domain on the T7 RNA polymerase that is associated with erroneous mRNA release. By leveraging the advantageous properties of a T7 RNA polymerase mutant and precisely optimized IVT process parameters, we successfully achieved an mRNA integrity exceeding 91%, thereby further unlocking the immense potential of mRNA therapeutics.

Alter codon bias of the P. pastoris genome to overcome a bottleneck in codon optimization strategy development and improve protein expression

Apart from its role in translation, codon bias is also an important mechanism to regulate mRNA levels. The traditional frequency-based codon optimization strategy is rather efficient in organisms such as N. crassa, but much less in yeast P. pastoris which is a popular host for heterologous protein expression. This is because that unlike N. crassa, the preferred codons of P. pastoris are actually AU-rich and hence codon optimization for extremely low GC content comes with issues of pre-mature transcriptional termination or low RNA stability in spite of translational advantages. To overcome this bottleneck, we focused on three reporter genes in P. pastoris first and confirmed the great advantage of GC-prone codon optimization on mRNA levels. Then we altered the codon bias profile of P. pastoris by introducing additional rare tRNA gene copies. Prior to that we constructed IPTG-regulated tRNA species to enable chassis cells to switch between different codon bias status. As demonstrated again with reporter genes, protein yield of luc and 0788 was successfully increased by 4-5 folds in chassis cells. In summary, here we provide an alternative codon optimization strategy for genes with unsatisfactory performance under traditional codon frequency-based optimization.

Transcriptional and post-transcriptional regulation of genes encoding secretory proteins in Aspergillus oryzae

Aspergillus oryzae, also known as the yellow koji mold, produces various hydrolytic enzymes that are widely used in different industries. Its high capacity to produce secretory proteins makes this filamentous fungus a suitable host for heterologous protein production. Amylolytic gene promoter is widely used to express heterologous genes in A. oryzae. The expression of this promoter is strictly regulated by several transcription factors, whose activation involves various factors. Furthermore, the expression levels of amylolytic and heterologous genes are post-transcriptionally regulated by mRNA degradation mechanisms in response to aberrant transcriptional termination or endoplasmic reticulum stress. This review discusses the transcriptional and post-transcriptional regulatory mechanisms controlling the expression of genes encoding secretory proteins in A. oryzae.

An analysis of codon utilization patterns in the chloroplast genomes of three species of Coffea

BACKGROUND

The chloroplast genome of plants is known for its small size and low mutation and recombination rates, making it a valuable tool in plant phylogeny, molecular evolution, and population genetics studies. Codon usage bias, an important evolutionary feature, provides insights into species evolution, gene function, and the expression of exogenous genes. Coffee, a key crop in the global tropical agricultural economy, trade, and daily life, warrants investigation into its codon usage bias to guide future research, including the selection of efficient heterologous expression systems for coffee genetic transformation.

RESULTS

Analysis of the codon utilization patterns in the chloroplast genomes of three Coffea species revealed a high degree of similarity among them. All three species exhibited similar base compositions, with high A/T content and low G/C content and a preference for A/T-ending codons. Among the 30 high-frequency codons identified, 96.67% had A/T endings. Fourteen codons were identified as ideal. Multiple mechanisms, including natural selection, were found to influence the codon usage patterns in the three coffee species, as indicated by ENc-GC3s mapping, PR2 analysis, and neutral analysis. Nicotiana tabacum and Saccharomyces cerevisiae have potential value as the heterologous expression host for three species of coffee genes.

CONCLUSION

This study highlights the remarkable similarity in codon usage patterns among the three coffee genomes, primarily driven by natural selection. Understanding the gene expression characteristics of coffee and elucidating the laws governing its genetic evolution are facilitated by investigating the codon preferences in these species. The findings can enhance the efficacy of exogenous gene expression and serve as a basis for future studies on coffee evolution.

Translation and mRNA Stability Control

Messenger RNA (mRNA) stability and translational efficiency are two crucial aspects of the post-transcriptional process that profoundly impact protein production in a cell. While it is widely known that ribosomes produce proteins, studies during the past decade have surprisingly revealed that ribosomes also control mRNA stability in a codon-dependent manner, a process referred to as codon optimality. Therefore, codons, the three-nucleotide words read by the ribosome, have a potent effect on mRNA stability and provide cis-regulatory information that extends beyond the amino acids they encode. While the codon optimality molecular mechanism is still unclear, the translation elongation rate appears to trigger mRNA decay. Thus, transfer RNAs emerge as potential master gene regulators affecting mRNA stability. Furthermore, while few factors related to codon optimality have been identified in yeast, the orthologous genes in vertebrates do not necessary share the same functions. Here, we discuss codon optimality findings and gene regulation layers related to codon composition in different eukaryotic species.

Ribosome biogenesis and function in development and disease

Although differential transcription drives the development of multicellular organisms, the ultimate readout of a protein-coding gene is ribosome-dependent mRNA translation. Ribosomes were once thought of as uniform molecular machines, but emerging evidence indicates that the complexity and diversity of ribosome biogenesis and function should be given a fresh look in the context of development. This Review begins with a discussion of different developmental disorders that have been linked with perturbations in ribosome production and function. We then highlight recent studies that reveal how different cells and tissues exhibit variable levels of ribosome production and protein synthesis, and how changes in protein synthesis capacity can influence specific cell fate decisions. We finish by touching upon ribosome heterogeneity in stress responses and development. These discussions highlight the importance of considering both ribosome levels and functional specialization in the context of development and disease.

Evolutionary diversification reveals distinct somatic versus germline cytoskeletal functions of the Arp2 branched actin nucleator protein

Branched actin networks are critical in many cellular processes, including cell motility and division. Arp2, a protein within the 7-membered Arp2/3 complex, is responsible for generating branched actin. Given its essential roles, Arp2 evolves under stringent sequence conservation throughout eukaryotic evolution. We unexpectedly discovered recurrent evolutionary diversification of Arp2 in Drosophila, yielding independently arising paralogs Arp2D in obscura species and Arp2D2 in montium species. Both paralogs are unusually testis-enriched in expression relative to Arp2. We investigated whether their sequence divergence from canonical Arp2 led to functional specialization by replacing Arp2 in D. melanogaster with either Arp2D or Arp2D2. Despite their divergence, we surprisingly found both complement Arp2's essential function in the soma, suggesting they have preserved the ability to polymerize branched actin even in a non-native species. However, we found that Arp2D-expressing males are subfertile and display many defects throughout sperm development. We pinpointed two highly diverged structural regions in Arp2D that contribute to these defects: subdomain 2 and the C-terminus. We expected that germline function would be rescued by replacing Arp2D's long and charged C-terminus with Arp2's short C-terminus, yet surprisingly, the essential somatic function of Arp2D was lost. Therefore, while Arp2D's structural divergence is incompatible with D. melanogaster sperm development, its unique C-terminus has evolved a critical role in actin polymerization. Our findings suggest canonical Arp2's function differs between somatic versus germline contexts, and Arp2 paralogs have recurrently evolved and specialized for actin branching in the testis.

Codon usage bias and genetic diversity in chloroplast genomes of Elaeagnus species (Myrtiflorae: Elaeagnaceae)

Codon usage bias (CUB) reveals the characteristics of species and can be utilized to understand their evolutionary relationship, increase the target genes' expression in the heterologous receptor plants, and further provide theoretic assistance for correlative study on molecular biology and genetic breeding. The chief aim of this work was to analyze the CUB in chloroplast (cp.) genes in nine Elaeagnus species to provide references for subsequent studies. The codons of Elaeagnus cp. genes preferred to end with A/T bases rather than with G/C bases. Most of the cp. genes were prone to mutation, while the rps7 genes were identical in sequences. Natural selection was inferred to have a powerful impact on the CUB in Elaeagnus cp. genomes, and their CUB was extremely strong. In addition, the optimal codons were identified in the nine cp. genomes based on the relative synonymous codon usage (RSCU) values, and the optimal codon numbers were between 15 and 19. The clustering analyses based on RSCU were contrasted with the maximum likelihood (ML)-based phylogenetic tree derived from coding sequences, suggesting that the t-distributed Stochastic Neighbor Embedding clustering method was more appropriate for evolutionary relationship analysis than the complete linkage method. Moreover, the ML-based phylogenetic tree based on the conservative matK genes and the whole cp. genomes had visible differences, indicating that the sequences of specific cp. genes were profoundly affected by their surroundings. Following the clustering analysis, Arabidopsis thaliana was considered the optimal heterologous expression receptor plant for the Elaeagnus cp. genes.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-023-01289-6.

The Influence of the Selection at the Amino Acid Level on Synonymous Codon Usage from the Viewpoint of Alternative Genetic Codes

Synonymous codon usage can be influenced by mutations and/or selection, e.g., for speed of protein translation and correct folding. However, this codon bias can also be affected by a general selection at the amino acid level due to differences in the acceptance of the loss and generation of these codons. To assess the importance of this effect, we constructed a mutation-selection model model, in which we generated almost 90,000 stationary nucleotide distributions produced by mutational processes and applied a selection based on differences in physicochemical properties of amino acids. Under these conditions, we calculated the usage of fourfold degenerated (4FD) codons and compared it with the usage characteristic of the pure mutations. We considered both the standard genetic code (SGC) and alternative genetic codes (AGCs). The analyses showed that a majority of AGCs produced a greater 4FD codon bias than the SGC. The mutations producing more thymine or adenine than guanine and cytosine increased the differences in usage. On the other hand, the mutational pressures generating a lot of cytosine or guanine with a low content of adenine and thymine decreased this bias because the nucleotide content of most 4FD codons stayed in the compositional equilibrium with these pressures. The comparison of the theoretical results with those for real protein coding sequences showed that the influence of selection at the amino acid level on the synonymous codon usage cannot be neglected. The analyses indicate that the effect of amino acid selection cannot be disregarded and that it can interfere with other selection factors influencing codon usage, especially in AT-rich genomes, in which AGCs are usually used.

Codon Usage Provide Insights into the Adaptation of Rice Genes under Stress Condition

Plants experience different stresses, i.e., abiotic, or biotic, and to combat them, plants re-program the expression of growth-, metabolism-, and resistance-related genes. These genes differ in their synonymous codon usage frequency and show codon usage bias. Here, we investigated the correlation among codon usage bias, gene expression, and underlying mechanisms in rice under abiotic and biotic stress conditions. The results indicated that genes with higher expression (up- or downregulated) levels had high GC content (≥60%), a low effective number of codon usage (≤40), and exhibited strong biases towards the codons with C/G at the third nucleotide position, irrespective of stress received. TTC, ATC, and CTC were the most preferred codons, while TAC, CAC, AAC, GAC, and TGC were moderately preferred under any stress (abiotic or biotic) condition. Additionally, downregulated genes are under mutational pressure (R² ≥ 0.5) while upregulated genes are under natural selection pressure (R² ≤ 0.5). Based on these results, we also identified the possible target codons that can be used to design an optimized set of genes with specific codons to develop climate-resilient varieties. Conclusively, under stress, rice has a bias towards codon usage which is correlated with GC content, gene expression level, and gene length.

Nutritional compensation of the circadian clock is a conserved process influenced by gene expression regulation and mRNA stability

Compensation is a defining principle of a true circadian clock, where its approximately 24-hour period length is relatively unchanged across environmental conditions. Known compensation effectors directly regulate core clock factors to buffer the oscillator's period length from variables in the environment. Temperature Compensation mechanisms have been experimentally addressed across circadian model systems, but much less is known about the related process of Nutritional Compensation, where circadian period length is maintained across physiologically relevant nutrient levels. Using the filamentous fungus Neurospora crassa, we performed a genetic screen under glucose and amino acid starvation conditions to identify new regulators of Nutritional Compensation. Our screen uncovered 16 novel mutants, and together with 4 mutants characterized in prior work, a model emerges where Nutritional Compensation of the fungal clock is achieved at the levels of transcription, chromatin regulation, and mRNA stability. However, eukaryotic circadian Nutritional Compensation is completely unstudied outside of Neurospora. To test for conservation in cultured human cells, we selected top hits from our fungal genetic screen, performed siRNA knockdown experiments of the mammalian orthologs, and characterized the cell lines with respect to compensation. We find that the wild-type mammalian clock is also compensated across a large range of external glucose concentrations, as observed in Neurospora, and that knocking down the mammalian orthologs of the Neurospora compensation-associated genes CPSF6 or SETD2 in human cells also results in nutrient-dependent period length changes. We conclude that, like Temperature Compensation, Nutritional Compensation is a conserved circadian process in fungal and mammalian clocks and that it may share common molecular determinants.

Composition and Codon Usage Pattern Results in Divergence of the Zinc Binuclear Cluster (Zn(II)2Cys6) Sequences among Ascomycetes Plant Pathogenic Fungi

Zinc binuclear cluster proteins (ZBC; Zn(II)2Cys6) are unique to the fungi kingdom and associated with a series of functions, viz., the utilization of macromolecules, stress tolerance, and most importantly, host-pathogen interactions by imparting virulence to the pathogen. Codon usage bias (CUB) is the phenomenon of using synonymous codons in a non-uniform fashion during the translation event, which has arisen because of interactions among evolutionary forces. The Zn(II)2Cys6 coding sequences from nine Ascomycetes plant pathogenic species and model system yeast were analysed for compositional and codon usage bias patterns. The clustering analysis diverged the Ascomycetes fungi into two clusters. The nucleotide compositional and relative synonymous codon usage (RSCU) analysis indicated GC biasness toward Ascomycetes fungi compared with the model system S. cerevisiae, which tends to be AT-rich. Further, plant pathogenic Ascomycetes fungi belonging to cluster-2 showed a higher number of GC-rich high-frequency codons than cluster-1 and was exclusively AT-rich in S. cerevisiae. The current investigation also showed the mutual effect of the two evolutionary forces, viz. natural selection and compositional constraints, on the CUB of Zn(II)2Cys6 genes. The perseverance of GC-rich codons of Zn(II)2Cys6 in Ascomycetes could facilitate the invasion process. The findings of the current investigation show the role of CUB and nucleotide composition in the evolutionary divergence of Ascomycetes plant pathogens and paves the way to target specific codons and sequences to modulate host-pathogen interactions through genome editing and functional genomics tools.

Mitogenome-wise codon usage pattern from comparative analysis of the first mitogenome of Blepharipa sp. (Muga uzifly) with other Oestroid flies

Uziflies (Family: Tachinidae) are dipteran endoparasites of sericigenous insects which cause major economic loss in the silk industry globally. Here, we are presenting the first full mitogenome of Blepharipa sp. (Acc: KY644698, 15,080 bp, A + T = 78.41%), a dipteran parasitoid of Muga silkworm (Antheraea assamensis) found in the Indian states of Assam and Meghalaya. This study has confirmed that Blepharipa sp. mitogenome gene content and arrangement is similar to other Tachinidae and Sarcophagidae flies of Oestroidea superfamily, typical of ancestral Diptera. Although, Calliphoridae and Oestridae flies have undergone tRNA translocation and insertion, forming unique intergenic spacers (IGS) and overlapping regions (OL) and a few of them (IGS, OL) have been conserved across Oestroidea flies. The Tachinidae mitogenomes exhibit more AT content and AT biased codons in their protein-coding genes (PCGs) than the Oestroidea counterpart. About 92.07% of all (3722) codons in PCGs of this new species have A/T in their 3rd codon position. The high proportion of AT and repeats in the control region (CR) affects sequence coverage, resulting in a short CR (Blepharipa sp.: 168 bp) and a smaller tachinid mitogenome. Our research unveils those genes with a high AT content had a reduced effective number of codons, leading to high codon usage bias. The neutrality test shows that natural selection has a stronger influence on codon usage bias than directed mutational pressure. This study also reveals that longer PCGs (e.g., nad5, cox1) have a higher codon usage bias than shorter PCGs (e.g., atp8, nad4l). The divergence rates increase nonlinearly as AT content at the 3rd codon position increases and higher rate of synonymous divergence than nonsynonymous divergence causes strong purifying selection. The phylogenetic analysis explains that Blepharipa sp. is well suited in the family of insectivorous tachinid maggots. It's possible that biased codon usage in the Tachinidae family reduces the effective number of codons, and purifying selection retains the core functions in their mitogenome, which could help with efficient metabolism in their endo-parasitic life style and survival strategy.

Strong negative correlation between codon usage bias and protein structural disorder impedes protein expression after codon optimization

As a common phenomenon existing in almost all genomes, codon usage bias has been studied for a long time. Codon optimization is a frequently used strategy to accelerate protein synthesis rate. Besides regulating protein translation speed, codon usage bias has also be reported to affect co-translation folding and transcription. P. pastoris is a well-developed expression system, whose efficiency is tightly correlated with commercial value. However, few studies focus on the role of codon usage bias in affecting protein expression in P. pastoris. Besides, many genes in P. pastoris genome show significant negative correlation between codon usage bias and protein structural disorder tendency. It's not known whether this feature is important for their expression. In order to answer these questions, we picked 4P. pastoris gene candidates with strong negative correlation between codon usage bias and protein structural disorder. We then performed full-length codon optimization which completely eliminated the correlation. Protein and RNA assays were then used to compare protein and mRNA levels before and after codon optimization. As a result, codon optimization failed to elevate their protein expression levels, and even resulted in a decrease. As represented by the trypsin sensitivity assays, codon optimization also altered the protein structure of 0616 and 0788. Besides protein, codon optimization also affected mRNA levels. Shown by in vitro and in vivo RNA degradation assays, the mRNA stability of 0616, 0788 and 0135 were also altered by codon optimization. For each gene, the detailed effect may be related with its specific sequence and protein structure. Our results suggest that codon usage bias is an important factor to regulate gene expression level, as well as mRNA and protein stabilities in P. pastoris. "Extreme" codon optimization in genes with strong negative correlation between codon usage bias and protein structural disorder tendency may not be favored. Compromised strategies should be tried if expression is not successful. Besides, codon optimization may affect protein structural conformation more severely in structural disordered proteins.

Codon usage and protein length-dependent feedback from translation elongation regulates translation initiation and elongation speed

Essential cellular functions require efficient production of many large proteins but synthesis of large proteins encounters many obstacles in cells. Translational control is mostly known to be regulated at the initiation step. Whether translation elongation process can feedback to regulate initiation efficiency is unclear. Codon usage bias, a universal feature of all genomes, plays an important role in determining gene expression levels. Here, we discovered that there is a conserved but codon usage-dependent genome-wide negative correlation between protein abundance and CDS length. The codon usage effects on protein expression and ribosome flux on mRNAs are influenced by CDS length; optimal codon usage preferentially promotes production of large proteins. Translation of mRNAs with long CDS and non-optimal codon usage preferentially induces phosphorylation of initiation factor eIF2α, which inhibits translation initiation efficiency. Deletion of the eIF2α kinase CPC-3 (GCN2 homolog) in Neurospora preferentially up-regulates large proteins encoded by non-optimal codons. Surprisingly, CPC-3 also inhibits translation elongation rate in a codon usage and CDS length-dependent manner, resulting in slow elongation rates for long CDS mRNAs. Together, these results revealed a codon usage and CDS length-dependent feedback mechanism from translation elongation to regulate both translation initiation and elongation kinetics.

Genome-Wide Association and Selective Sweep Studies Reveal the Complex Genetic Architecture of DMI Fungicide Resistance in Cercospora beticola

The rapid and widespread evolution of fungicide resistance remains a challenge for crop disease management. The demethylation inhibitor (DMI) class of fungicides is a widely used chemistry for managing disease, but there has been a gradual decline in efficacy in many crop pathosystems. Reliance on DMI fungicides has increased resistance in populations of the plant pathogenic fungus Cercospora beticola worldwide. To better understand the genetic and evolutionary basis for DMI resistance in C. beticola, a genome-wide association study (GWAS) and selective sweep analysis were conducted for the first time in this species. We performed whole-genome resequencing of 190 C. beticola isolates infecting sugar beet (Beta vulgaris ssp. vulgaris). All isolates were phenotyped for sensitivity to the DMI tetraconazole. Intragenic markers on chromosomes 1, 4, and 9 were significantly associated with DMI fungicide resistance, including a polyketide synthase gene and the gene encoding the DMI target CbCYP51. Haplotype analysis of CbCYP51 identified a synonymous mutation (E170) and nonsynonymous mutations (L144F, I387M, and Y464S) associated with DMI resistance. Genome-wide scans of selection showed that several of the GWAS mutations for fungicide resistance resided in regions that have recently undergone a selective sweep. Using radial plate growth on selected media as a fitness proxy, we did not find a trade-off associated with DMI fungicide resistance. Taken together, we show that population genomic data from a crop pathogen can allow the identification of mutations conferring fungicide resistance and inform about their origins in the pathogen population.

Genome-wide role of codon usage on transcription and identification of potential regulators

Codon usage bias is a fundamental feature of all genomes and plays an important role in determining gene expression levels. The codon usage was thought to influence gene expression mainly due to its impact on translation. Recently, however, codon usage was shown to affect transcription of fungal and mammalian genes, indicating the existence of a gene regulatory phenomenon with unknown mechanism. In Neurospora, codon usage biases strongly correlate with mRNA levels genome-wide, and here we show that the correlation between codon usage and RNA levels is maintained in the nucleus. In addition, codon optimality is tightly correlated with both total and nuclear RNA levels, suggesting that codon usage broadly influences mRNA levels through transcription in a translation-independent manner. A large-scale RNA sequencing-based genetic screen in Neurospora identified 18 candidate factors that when deleted decreased the genome-wide correlation between codon usage and RNA levels and reduced the codon usage effect on gene expression. Most of these factors, such as the H3K36 methyltransferase, are chromatin regulators or transcription factors. Together, our results suggest that the transcriptional effect of codon usage is mediated by multiple transcriptional regulatory mechanisms.

Analysis of codon usage of Horseshoe Bat Hepatitis B virus and its host

In the present analysis, codon usage strategies and base distribution of Horseshoe bat hepatitis B virus (HBHBV) were analyzed and compared with its host Rhinolophus sinicus, as no work was yet reported. The magnitude of synonymous codon usage bias (CUB) in the virus and its host was low with higher proportion of the base C. Notably, 21 more frequently used codons, 19 less frequently used codons and 3 underrepresented codons (TCG, ACG and GCG) were found to be similar in both virus and its host coding sequences. Neutrality plot analysis reported greater role of natural selection in HBHBV (67.84%) and R. sinicus (76.90%) over mutation pressure. Base skewness and protein properties also influenced the CUB of genes. Further, codon usage analysis depicted, HBHBV and R. sinicus had many similarities in codon usage patterns that might reflect viral adaptation to its host.

IMITATION SWITCH is required for normal chromatin structure and gene repression in PRC2 target domains

Polycomb Group (PcG) proteins are part of an epigenetic cell memory system that plays essential roles in multicellular development, stem cell biology, X chromosome inactivation, and cancer. In animals, plants, and many fungi, Polycomb Repressive Complex 2 (PRC2) catalyzes trimethylation of histone H3 lysine 27 (H3K27me3) to assemble transcriptionally repressed facultative heterochromatin. PRC2 is structurally and functionally conserved in the model fungus Neurospora crassa, and recent work in this organism has generated insights into PRC2 control and function. To identify components of the facultative heterochromatin pathway, we performed a targeted screen of Neurospora deletion strains lacking individual ATP-dependent chromatin remodeling enzymes. We found the Neurospora homolog of IMITATION SWITCH (ISW) is critical for normal transcriptional repression, nucleosome organization, and establishment of typical histone methylation patterns in facultative heterochromatin domains. We also found that stable interaction between PRC2 and chromatin depends on ISW. A functional ISW ATPase domain is required for gene repression and normal H3K27 methylation. ISW homologs interact with accessory proteins to form multiple complexes with distinct functions. Using proteomics and molecular approaches, we identified three distinct Neurospora ISW-containing complexes. A triple mutant lacking three ISW accessory factors and disrupting multiple ISW complexes led to widespread up-regulation of PRC2 target genes and altered H3K27 methylation patterns, similar to an ISW-deficient strain. Taken together, our data show that ISW is a key component of the facultative heterochromatin pathway in Neurospora, and that distinct ISW complexes perform an apparently overlapping role to regulate chromatin structure and gene repression at PRC2 target domains.

Learning the Regulatory Code of Gene Expression

Data-driven machine learning is the method of choice for predicting molecular phenotypes from nucleotide sequence, modeling gene expression events including protein-DNA binding, chromatin states as well as mRNA and protein levels. Deep neural networks automatically learn informative sequence representations and interpreting them enables us to improve our understanding of the regulatory code governing gene expression. Here, we review the latest developments that apply shallow or deep learning to quantify molecular phenotypes and decode the cis-regulatory grammar from prokaryotic and eukaryotic sequencing data. Our approach is to build from the ground up, first focusing on the initiating protein-DNA interactions, then specific coding and non-coding regions, and finally on advances that combine multiple parts of the gene and mRNA regulatory structures, achieving unprecedented performance. We thus provide a quantitative view of gene expression regulation from nucleotide sequence, concluding with an information-centric overview of the central dogma of molecular biology.

Transcription, mRNA export and immune evasion shape the codon usage of viruses

The nucleotide composition, dinucleotide composition, and codon usage of many viruses differs from their hosts. These differences arise because viruses are subject to unique mutation and selection pressures that do not apply to host genomes; however, the molecular mechanisms that underlie these evolutionary forces are unclear. Here, we analysed the patterns of codon usage in 1,520 vertebrate-infecting viruses, focusing on parameters known to be under selection and associated with gene regulation. We find that GC content, dinucleotide content, and splicing and m6A modification-related sequence motifs are associated with the type of genetic material (DNA or RNA), strandedness, and replication compartment of viruses. In an experimental follow-up, we find that the effects of GC content on gene expression depend on whether the genetic material is delivered to the cell as DNA or mRNA, whether it is transcribed by endogenous or exogenous RNA polymerase, and whether transcription takes place in the nucleus or cytoplasm. Our results suggest that viral codon usage cannot be explained by a simple adaptation to the codon usage of the host - instead, it reflects the combination of multiple selective and mutational pressures, including the need for efficient transcription, export, and immune evasion.

Read between the lines: Diversity of non-translational selection pressures on local codon usage

Protein coding genes can contain specific motifs within their nucleotide sequence that function as a signal for various biological pathways. The presence of such sequence motifs within a gene can have beneficial or detrimental effects on the phenotype and fitness of an organism, and this can lead to the enrichment or avoidance of this sequence motif. The degeneracy of the genetic code allows for the existence of alternative synonymous sequences that exclude or include these motifs, while keeping the encoded amino acid sequence intact. This implies that locally, there can be a selective pressure for preferentially using a codon over its synonymous alternative in order to avoid or enrich a specific sequence motif. This selective pressure could -in addition to mutation, drift and selection for translation efficiency and accuracy- contribute to shape the codon usage bias. In this review, we discuss patterns of avoidance of (or enrichment for) the various biological signals contained in specific nucleotide sequence motifs: transcription and translation initiation and termination signals, mRNA maturation signals, and antiviral immune system targets. Experimental data on the phenotypic or fitness effects of synonymous mutations in these sequence motifs confirm that they can be targets of local selection pressures on codon usage. We also formulate the hypothesis that transposable elements could have a similar impact on codon usage through their preferred integration sequences. Overall, selection on codon usage appears to be a combination of a global selection pressure imposed by the translation machinery, and a patchwork of local selection pressures related to biological signals contained in specific sequence motifs.

Protein products of nonstop mRNA disrupt nucleolar homeostasis

Stalled mRNA translation results in the production of incompletely synthesized proteins that are targeted for degradation by ribosome-associated quality control (RQC). Here we investigated the fate of defective proteins translated from stall-inducing, nonstop mRNA that escape ubiquitylation by the RQC protein LTN1. We found that nonstop protein products accumulated in nucleoli and this localization was driven by polylysine tracts produced by translation of the poly(A) tails of nonstop mRNA. Nucleolar sequestration increased the solubility of invading proteins but disrupted nucleoli, altering their dynamics, morphology, and resistance to stress in cell culture and intact flies. Our work elucidates how stalled translation may affect distal cellular processes and may inform studies on the pathology of diseases caused by failures in RQC and characterized by nucleolar stress.

Codon usage bias regulates gene expression and protein conformation in yeast expression system P. pastoris

Background

Protein synthesis is one of the extremely important anabolic pathways in the yeast expression system Pichia pastoris. Codon optimization is a commonly adopted strategy for improved protein expression, although unexpected failures did appear sometimes waiting for further exploration. Recently codon bias has been studied to regulate protein folding and activity in many other organisms.

Results

Here the codon bias profile of P. pastoris genome was examined first and a direct correlation between codon translation efficiency and usage frequency was identified. By manipulating the codon choices of both endogenous and heterologous signal peptides, secretion abilities of N-terminal signal peptides were shown to be tolerant towards codon changes. Then two gene candidates with different levels of structural disorder were studied, and full-length codon optimization was found to affect their expression profiles differentially. Finally, more evidences were provided to support possible protein conformation change brought by codon optimization in structurally disordered proteins.

Conclusion

Our results suggest that codon bias regulates gene expression by modulating several factors including transcription and translation efficiency, protein folding and activity. Because of sequences difference, the extent of affection may be gene specific. For some genes, special codon optimization strategy should be adopted to ensure appropriate expression and conformation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12934-021-01580-9.

Effects of codon usage on gene expression are promoter context dependent

Codon usage bias is a universal feature of all genomes. Although codon usage has been shown to regulate mRNA and protein levels by influencing mRNA decay and transcription in eukaryotes, little or no genome-wide correlations between codon usage and mRNA levels are detected in mammalian cells, raising doubt on the significance of codon usage effect on gene expression. Here we show that gene-specific regulation reduces the genome-wide codon usage and mRNA correlations: Constitutively expressed genes exhibit much higher genome-wide correlations than differentially expressed genes from fungi to human cells. Using Drosophila S2 cells as a model system, we showed that the effect of codon usage on mRNA expression level is promoter-dependent. Regions downstream of the core promoters of differentially expressed genes can repress the codon usage effects on mRNA expression. An element in the Hsp70 promoter was identified to be necessary and sufficient for this inhibitory effect. The promoter-dependent codon usage effects on mRNA levels are regulated at the transcriptional level through modulation of histone modifications, nucleosome densities and premature termination. Together, our results demonstrate that promoters play a major role in determining whether codon usage influences gene expression and further establish the transcription-dependent codon usage effects on gene expression.

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.

Deep learning suggests that gene expression is encoded in all parts of a co-evolving interacting gene regulatory structure

Understanding the genetic regulatory code governing gene expression is an important challenge in molecular biology. However, how individual coding and non-coding regions of the gene regulatory structure interact and contribute to mRNA expression levels remains unclear. Here we apply deep learning on over 20,000 mRNA datasets to examine the genetic regulatory code controlling mRNA abundance in 7 model organisms ranging from bacteria to Human. In all organisms, we can predict mRNA abundance directly from DNA sequence, with up to 82% of the variation of transcript levels encoded in the gene regulatory structure. By searching for DNA regulatory motifs across the gene regulatory structure, we discover that motif interactions could explain the whole dynamic range of mRNA levels. Co-evolution across coding and non-coding regions suggests that it is not single motifs or regions, but the entire gene regulatory structure and specific combination of regulatory elements that define gene expression levels.

Nonoptimal Codon Usage Is Critical for Protein Structure and Function of the Master General Amino Acid Control Regulator CPC-1

Under amino acid starvation conditions, eukaryotic organisms activate a general amino acid control response. In Neurospora crassa, Cross Pathway Control Protein 1 (CPC-1), the ortholog of the Saccharomyces cerevisiae bZIP transcription factor GCN4, functions as the master regulator of the general amino acid control response. Codon usage biases are a universal feature of eukaryotic genomes and are critical for regulation of gene expression. Although codon usage has also been implicated in the regulation of protein structure and function, genetic evidence supporting this conclusion is very limited. Here, we show that Neurospora cpc-1 has a nonoptimal NNU-rich codon usage profile that contrasts with the strong NNC codon preference in the genome. Although substitution of the cpc-1 NNU codons with synonymous NNC codons elevated CPC-1 expression in Neurospora, it altered the CPC-1 degradation rate and abolished its amino acid starvation-induced protein stabilization. The codon-manipulated CPC-1 protein also exhibited different sensitivity to limited protease digestion. Furthermore, CPC-1 functions in rescuing the cell growth of the cpc-1 deletion mutant and activation of the expression of its target genes were impaired by the synonymous codon changes. Together, these results reveal the critical role of codon usage in regulation of CPC-1 expression and function and establish a genetic example of the importance of codon usage in protein folding.IMPORTANCE The general amino acid control response is critical for adaptation of organisms to amino acid starvation conditions. The preference to use certain synonymous codons is a universal feature of all genomes. Synonymous codon changes were previously thought to be silent mutations. In this study, we showed that the Neurospora cpc-1 gene has an unusual codon usage profile compared to other genes in the genome. We found that codon optimization of the cpc-1 gene without changing its amino acid sequence resulted in elevated CPC-1 expression, an altered protein degradation rate, and impaired protein functions due to changes in protein structure. Together, these results reveal the critical role of synonymous codon usage in regulation of CPC-1 expression and function and establish a genetic example of the importance of codon usage in protein structure.

A code within the genetic code: codon usage regulates co-translational protein folding

The genetic code is degenerate, and most amino acids are encoded by two to six synonymous codons. Codon usage bias, the preference for certain synonymous codons, is a universal feature of all genomes examined. Synonymous codon mutations were previously thought to be silent; however, a growing body evidence now shows that codon usage regulates protein structure and gene expression through effects on co-translational protein folding, translation efficiency and accuracy, mRNA stability, and transcription. Codon usage regulates the speed of translation elongation, resulting in non-uniform ribosome decoding rates on mRNAs during translation that is adapted to co-translational protein folding process. Biochemical and genetic evidence demonstrate that codon usage plays an important role in regulating protein folding and function in both prokaryotic and eukaryotic organisms. Certain protein structural types are more sensitive than others to the effects of codon usage on protein folding, and predicted intrinsically disordered domains are more prone to misfolding caused by codon usage changes than other domain types. Bioinformatic analyses revealed that gene codon usage correlates with different protein structures in diverse organisms, indicating the existence of a codon usage code for co-translational protein folding. This review focuses on recent literature on the role and mechanism of codon usage in regulating translation kinetics and co-translational protein folding.

Adaptation of codon usage to tRNA I34 modification controls translation kinetics and proteome landscape

Codon usage bias is a universal feature of all genomes and plays an important role in regulating protein expression levels. Modification of adenosine to inosine at the tRNA anticodon wobble position (I34) by adenosine deaminases (ADATs) is observed in all eukaryotes and has been proposed to explain the correlation between codon usage and tRNA pool. However, how the tRNA pool is affected by I34 modification to influence codon usage-dependent gene expression is unclear. Using Neurospora crassa as a model system, by combining molecular, biochemical and bioinformatics analyses, we show that silencing of adat2 expression severely impaired the I34 modification levels for the ADAT-related tRNAs, resulting in major ADAT-related tRNA profile changes and reprogramming of translation elongation kinetics on ADAT-related codons. adat2 silencing also caused genome-wide codon usage-biased ribosome pausing on mRNAs and proteome landscape changes, leading to selective translational repression or induction of different mRNAs. The induced expression of CPC-1, the Neurospora ortholog of yeast GCN4p, mediates the transcriptional response after adat2 silencing and amino acid starvation. Together, our results demonstrate that the tRNA I34 modification by ADAT plays a major role in driving codon usage-biased translation to shape proteome landscape.

Modification of transfer RNA (tRNA) can have profound impacts on gene expression by shaping cellular tRNA pool. How codon usage bias and tRNA profiles synergistically regulate gene expression is unclear. By combining molecular, biochemical and bioinformatics analyses, we showed that the correlation between genome codon usage and tRNA I34 (inosine 34) modification modulates translation elongation kinetics and proteome landscape. Inhibition of tRNA I34 modification causes codon usage-dependent ribosome pausing on mRNAs during translation and changes cellular protein contents in a codon usage biased manner. Together, our results demonstrate that the tRNA I34 modification plays a major role in driving codon usage-dependent translation to determine proteome landscape in a eukaryotic organism.

Codon Usage and Splicing Jointly Influence mRNA Localization

In the human genome, most genes undergo splicing, and patterns of codon usage are splicing dependent: guanine and cytosine (GC) content is the highest within single-exon genes and within first exons of multi-exon genes. However, the effects of codon usage on gene expression are typically characterized in unspliced model genes. Here, we measured the effects of splicing on expression in a panel of synonymous reporter genes that varied in nucleotide composition. We found that high GC content increased protein yield, mRNA yield, cytoplasmic mRNA localization, and translation of unspliced reporters. Splicing did not affect the expression of GC-rich variants. However, splicing promoted the expression of AT-rich variants by increasing their steady-state protein and mRNA levels, in part through promoting cytoplasmic localization of mRNA. We propose that splicing promotes the nuclear export of AU-rich mRNAs and that codon- and splicing-dependent effects on expression are under evolutionary pressure in the human genome.

Codon usage bias affects α-amylase mRNA level by altering RNA stability and cytosine methylation patterns in Escherichia coli

Codon usage bias exists in almost every organism and is reported to regulate protein translation efficiency and folding. Besides translation, the preliminary role of codon usage bias on gene transcription has also been revealed in some eukaryotes such as Neurospora crassa. In this study, we took as an example the α-amylase-coding gene (amyA) and examined the role of codon usage bias in regulating gene expression in the typical prokaryote Escherichia coli. We confirmed the higher translation efficiency on codon-optimized amyA RNAs and found that the RNA level itself was also affected by codon optimization. The decreased RNA level was caused at least in part by altered mRNA stability at the post-transcriptional level. Codon optimization also altered the number of cytosine methylation sites. Examination on dcm knockouts suggested that cytosine methylation may be a minor mechanism adopted by codon bias to regulate gene RNA levels. More studies are required to verify the global effect of codon usage and to reveal its detailed mechanism on transcription.

Hidden patterns of codon usage bias across kingdoms

The genetic code is necessarily degenerate with 64 possible nucleotide triplets being translated into 20 amino acids. Eighteen out of the 20 amino acids are encoded by multiple synonymous codons. While synonymous codons are clearly equivalent in terms of the information they carry, it is now well established that they are used in a biased fashion. There is currently no consensus as to the origin of this bias. Drawing on ideas from stochastic thermodynamics we derive from first principles a mathematical model describing the statistics of codon usage bias. We show that the model accurately describes the distribution of codon usage bias of genomes in the fungal and bacterial kingdoms. Based on it, we derive a new computational measure of codon usage bias-the distance D capturing two aspects of codon usage bias: (i) differences in the genome-wide frequency of codons and (ii) apparent non-random distributions of codons across mRNAs. By means of large scale computational analysis of over 900 species across two kingdoms of life, we demonstrate that our measure provides novel biological insights. Specifically, we show that while codon usage bias is clearly based on heritable traits and closely related species show similar degrees of bias, there is considerable variation in the magnitude of D within taxonomic classes suggesting that the contribution of sequence-level selection to codon bias varies substantially within relatively confined taxonomic groups. Interestingly, commonly used model organisms are near the median for values of D for their taxonomic class, suggesting that they may not be good representative models for species with more extreme D, which comprise organisms of medical and agricultural interest. We also demonstrate that amino acid specific patterns of codon usage are themselves quite variable between branches of the tree of life, and that some of this variability correlates with organismal tRNA content.

eRF1 mediates codon usage effects on mRNA translation efficiency through premature termination at rare codons

Codon usage bias is a universal feature of eukaryotic and prokaryotic genomes and plays an important role in regulating gene expression levels. A major role of codon usage is thought to regulate protein expression levels by affecting mRNA translation efficiency, but the underlying mechanism is unclear. By analyzing ribosome profiling results, here we showed that codon usage regulates translation elongation rate and that rare codons are decoded more slowly than common codons in all codon families in Neurospora. Rare codons resulted in ribosome stalling in manners both dependent and independent of protein sequence context and caused premature translation termination. This mechanism was shown to be conserved in Drosophila cells. In both Neurospora and Drosophila cells, codon usage plays an important role in regulating mRNA translation efficiency. We found that the rare codon-dependent premature termination is mediated by the translation termination factor eRF1, which recognizes ribosomes stalled on rare sense codons. Silencing of eRF1 expression resulted in codon usage-dependent changes in protein expression. Together, these results establish a mechanism for how codon usage regulates mRNA translation efficiency.

Alternative Splicing in Heat Shock Protein Transcripts as a Mechanism of Cell Adaptation in Trichophyton rubrum

Heat shock proteins (HSPs) are involved in critical processes like host tissue invasion, resistance, and pathogenicity in dermatophytes. RNA-Seq analysis of Trichophyton rubrum exposed to undecanoic acid (UDA) revealed intron retention events in HSP transcripts. Because HSPs are modulated in response to various stimuli and as alternative splicing (AS) can result in a broad diversity in the proteome of eukaryotic cells, our objective was to confirm the aforementioned retention events, investigating their consequences and extent. Furthermore, we aimed to determine: (1) the expression profile of HSP genes in an infection-like scenario and (2) the importance of Hsp90 for the keratinolytic potential of T. rubrum. RT and qPCR analyses comparing the exposure to UDA and terbinafine (TRB) confirmed the presence of two mRNA isoforms of the hsp7-like gene, with distinct expression patterns in response to UDA and TRB. The HSP expression profile revealed two upregulated, three downregulated, and four unmodulated transcripts; Hsp90 inhibition by 17-AAG resulted in a significant decrease in keratinolytic potential at 37 °C. Altogether, these results broaden the current knowledge on the importance of HSP-mediated pathways for cell adaptation and other aspects of dermatophyte biology, indicating that HSP network proteins can be potential targets for antifungal therapy.

Codon bias confers stability to human mRNAs

Codon bias has been implicated as one of the major factors contributing to mRNA stability in several model organisms. However, the molecular mechanisms of codon bias on mRNA stability remain unclear in humans. Here, we show that human cells possess a mechanism to modulate RNA stability through a unique codon bias. Bioinformatics analysis showed that codons could be clustered into two distinct groups-codons with G or C at the third base position (GC3) and codons with either A or T at the third base position (AT3): the former stabilizing while the latter destabilizing mRNA. Quantification of codon bias showed that increased GC3-content entails proportionately higher GC-content. Through bioinformatics, ribosome profiling, and in vitro analysis, we show that decoupling the effects of codon bias reveals two modes of mRNA regulation, one GC3- and one GC-content dependent. Employing an immunoprecipitation-based strategy, we identify ILF2 and ILF3 as RNA-binding proteins that differentially regulate global mRNA abundances based on codon bias. Our results demonstrate that codon bias is a two-pronged system that governs mRNA abundance.

Rapidly evolving protointrons in Saccharomyces genomes revealed by a hungry spliceosome

Introns are a prevalent feature of eukaryotic genomes, yet their origins and contributions to genome function and evolution remain mysterious. In budding yeast, repression of the highly transcribed intron-containing ribosomal protein genes (RPGs) globally increases splicing of non-RPG transcripts through reduced competition for the spliceosome. We show that under these “hungry spliceosome” conditions, splicing occurs at more than 150 previously unannotated locations we call protointrons that do not overlap known introns. Protointrons use a less constrained set of splice sites and branchpoints than standard introns, including in one case AT-AC in place of GT-AG. Protointrons are not conserved in all closely related species, suggesting that most are not under positive selection and are fated to disappear. Some are found in non-coding RNAs (e. g. CUTs and SUTs), where they may contribute to the creation of new genes. Others are found across boundaries between noncoding and coding sequences, or within coding sequences, where they offer pathways to the creation of new protein variants, or new regulatory controls for existing genes. We define protointrons as (1) nonconserved intron-like sequences that are (2) infrequently spliced, and importantly (3) are not currently understood to contribute to gene expression or regulation in the way that standard introns function. A very few protointrons in S. cerevisiae challenge this classification by their increased splicing frequency and potential function, consistent with the proposed evolutionary process of “intronization”, whereby new standard introns are created. This snapshot of intron evolution highlights the important role of the spliceosome in the expansion of transcribed genomic sequence space, providing a pathway for the rare events that may lead to the birth of new eukaryotic genes and the refinement of existing gene function.

The protein coding information in eukaryotic genes is broken by intervening sequences called introns that are removed from RNA during transcription by a large protein-RNA complex called the spliceosome. Where introns come from and how the spliceosome contributes to genome evolution are open questions. In this study, we find more than 150 new places in the yeast genome that are recognized by the spliceosome and spliced out as introns. Since they appear to have arisen very recently in evolution by sequence drift and do not appear to contribute to gene expression or its regulation, we call these protointrons. Protointrons are found in both protein-coding and non-coding RNAs and are not efficiently removed by the splicing machinery. Although most protointrons are not conserved and will likely disappear as evolution proceeds, a few are spliced more efficiently, and are located where they might begin to play functional roles in gene expression, as predicted by the proposed process of intronization. The challenge now is to understand how spontaneously appearing splicing events like protointrons might contribute to the creation of new genes, new genetic controls, and new protein isoforms as genomes evolve.

Rapidly evolving protointrons in Saccharomyces genomes revealed by a hungry spliceosome

Introns are a prevalent feature of eukaryotic genomes, yet their origins and contributions to genome function and evolution remain mysterious. In budding yeast, repression of the highly transcribed intron-containing ribosomal protein genes (RPGs) globally increases splicing of non-RPG transcripts through reduced competition for the spliceosome. We show that under these "hungry spliceosome" conditions, splicing occurs at more than 150 previously unannotated locations we call protointrons that do not overlap known introns. Protointrons use a less constrained set of splice sites and branchpoints than standard introns, including in one case AT-AC in place of GT-AG. Protointrons are not conserved in all closely related species, suggesting that most are not under positive selection and are fated to disappear. Some are found in non-coding RNAs (e. g. CUTs and SUTs), where they may contribute to the creation of new genes. Others are found across boundaries between noncoding and coding sequences, or within coding sequences, where they offer pathways to the creation of new protein variants, or new regulatory controls for existing genes. We define protointrons as (1) nonconserved intron-like sequences that are (2) infrequently spliced, and importantly (3) are not currently understood to contribute to gene expression or regulation in the way that standard introns function. A very few protointrons in S. cerevisiae challenge this classification by their increased splicing frequency and potential function, consistent with the proposed evolutionary process of "intronization", whereby new standard introns are created. This snapshot of intron evolution highlights the important role of the spliceosome in the expansion of transcribed genomic sequence space, providing a pathway for the rare events that may lead to the birth of new eukaryotic genes and the refinement of existing gene function.

Transcriptional Control by Premature Termination: A Forgotten Mechanism

The concept of early termination as an important means of transcriptional control has long been established. Even so, its role in metazoan gene expression is underappreciated. Recent technological advances provide novel insights into premature transcription termination (PTT). This process is frequent, widespread, and can occur close to the transcription start site (TSS), or within the gene body. Stable prematurely terminated transcripts contribute to the transcriptome as instances of alternative polyadenylation (APA). Independently of transcript stability and function, premature termination opposes the formation of full-length transcripts, thereby negatively regulating gene expression, especially of transcriptional regulators. Premature termination can be beneficial or harmful, depending on its context. As a result, multiple factors have evolved to control this process.

Variation and selection on codon usage bias across an entire subphylum

Variation in synonymous codon usage is abundant across multiple levels of organization: between codons of an amino acid, between genes in a genome, and between genomes of different species. It is now well understood that variation in synonymous codon usage is influenced by mutational bias coupled with both natural selection for translational efficiency and genetic drift, but how these processes shape patterns of codon usage bias across entire lineages remains unexplored. To address this question, we used a rich genomic data set of 327 species that covers nearly one third of the known biodiversity of the budding yeast subphylum Saccharomycotina. We found that, while genome-wide relative synonymous codon usage (RSCU) for all codons was highly correlated with the GC content of the third codon position (GC3), the usage of codons for the amino acids proline, arginine, and glycine was inconsistent with the neutral expectation where mutational bias coupled with genetic drift drive codon usage. Examination between genes' effective numbers of codons and their GC3 contents in individual genomes revealed that nearly a quarter of genes (381,174/1,683,203; 23%), as well as most genomes (308/327; 94%), significantly deviate from the neutral expectation. Finally, by evaluating the imprint of translational selection on codon usage, measured as the degree to which genes' adaptiveness to the tRNA pool were correlated with selective pressure, we show that translational selection is widespread in budding yeast genomes (264/327; 81%). These results suggest that the contribution of translational selection and drift to patterns of synonymous codon usage across budding yeasts varies across codons, genes, and genomes; whereas drift is the primary driver of global codon usage across the subphylum, the codon bias of large numbers of genes in the majority of genomes is influenced by translational selection.

Codon usage regulates human KRAS expression at both transcriptional and translational levels

KRAS and HRAS are highly homologous oncogenic Ras GTPase family members that are mutated in a wide spectrum of human cancers. Despite having high amino acid identity, KRAS and HRAS have very different codon usage biases: the HRAS gene contains many common codons, and KRAS is enriched for rare codons. Rare codons in KRAS suppress its protein expression, which has been shown to affect both normal and cancer biology in mammals. Here, using HRAS or KRAS expression in different human cell lines and in vitro transcription and translation assays, we show that KRAS rare codons inhibit both translation efficiency and transcription and that the contribution of these two processes varies among different cell lines. We observed that codon usage regulates mRNA translation efficiency such that WT KRAS mRNA is poorly translated. On the other hand, common codons increased transcriptional rates by promoting activating histone modifications and recruitment of transcriptional coactivators. Finally, we found that codon usage also influences KRAS protein conformation, likely because of its effect on co-translational protein folding. Together, our results reveal that codon usage has multidimensional effects on protein expression, ranging from effects on transcription to protein folding in human cells.

Synonymous nucleotide modification of the KCNH2 gene affects both mRNA characteristics and translation of the encoded hERG ion channel

Synonymous nucleotide variation is increasingly recognized as a factor than can affect protein expression, but the underlying mechanisms are incompletely understood. Here, we investigated whether synonymous changes could affect expression of the potassium voltage-gated channel subfamily H member 2 (KCNH2) gene, encoding the human ether-a-go-go-related gene (hERG) ion channel, which is linked to hereditary cardiac arrhythmia. We examined a previously described synthetic version (hERG-codon modified (CM)) with synonymous substitutions designed to reduce GC content, rare codons, and mRNA secondary structure relative to the native construct (hERG-NT). hERG-CM exhibited lower protein expression than hERG-NT in HEK293T cells. We found that the steady-state abundance of hERG-NT mRNA was greater than hERG-CM because of an enhanced transcription rate and increased mRNA stability for hERG-NT. Translation of hERG-CM was independently reduced, contributing to the overall greater synthesis of hERG-NT channel protein. This was partially offset, however, by a higher aggregation of a newly synthesized hERG-NT channel, resulting in nonfunctional protein. Regional mRNA analyses of chimeras of hERG-NT and hERG-CM revealed that synonymous changes in the 5' segments of the coding region had the greatest influence on hERG synthesis at both the mRNA and protein levels. Taken together, these results indicate that synonymous nucleotide variations within the coding region, particularly in the 5' region of the hERG mRNA, can affect both transcription and translation. These findings support the notion that greater attention should be given to the effects of synonymous genetic variation when analyzing hERG DNA sequences in the study of hereditary cardiac disease.