Causes and effects of N-terminal codon bias in bacterial genes

Mechanisms governing codon usage bias and the implications for protein expression in the chloroplast of Chlamydomonas reinhardtii

Chloroplasts possess a considerably reduced genome that is decoded via an almost minimal set of tRNAs. These features make an excellent platform for gaining insights into fundamental mechanisms that govern protein expression. Here, we present a comprehensive and revised perspective of the mechanisms that drive codon selection in the chloroplast of Chlamydomonas reinhardtii and the functional consequences for protein expression. In order to extract this information, we applied several codon usage descriptors to genes with different expression levels. We show that highly expressed genes strongly favor translationally optimal codons, while genes with lower functional importance are rather affected by directional mutational bias. We demonstrate that codon optimality can be deduced from codon-anticodon pairing affinity and, for a small number of amino acids (leucine, arginine, serine, and isoleucine), tRNA concentrations. Finally, we review, analyze, and expand on the impact of codon usage on protein yield, secondary structures of mRNA, translation initiation and termination, and amino acid composition of proteins, as well as cotranslational protein folding. The comprehensive analysis of codon choice provides crucial insights into heterologous gene expression in the chloroplast of C. reinhardtii, which may also be applicable to other chloroplast-containing organisms and bacteria.

Flow-Seq Evaluation of Translation Driven by a Set of Natural Escherichia coli 5'-UTR of Variable Length

Flow-seq is a method that combines fluorescently activated cell sorting and next-generation sequencing to deduce a large amount of data about translation efficiency from a single experiment. Here, we constructed a library of fluorescent protein-based reporters preceded by a set of 648 natural 5'-untranslated regions (5'-UTRs) of Escherichia coli genes. Usually, Flow-seq libraries are constructed using uniform-length sequence elements, in contrast to natural situations, where functional elements are of heterogenous lengths. Here, we demonstrated that a 5'-UTR library of variable length could be created and analyzed with Flow-seq. In line with previous Flow-seq experiments with randomized 5'-UTRs, we observed the influence of an RNA secondary structure and Shine-Dalgarno sequences on translation efficiency; however, the variability of these parameters for natural 5'-UTRs in our library was smaller in comparison with randomized libraries. In line with this, we only observed a 30-fold difference in translation efficiency between the best and worst bins sorted with this factor. The results correlated with those obtained with ribosome profiling.

The landscape of transcriptional and translational changes over 22 years of bacterial adaptation

Organisms can adapt to an environment by taking multiple mutational paths. This redundancy at the genetic level, where many mutations have similar phenotypic and fitness effects, can make untangling the molecular mechanisms of complex adaptations difficult. Here, we use the Escherichia coli long-term evolution experiment (LTEE) as a model to address this challenge. To understand how different genomic changes could lead to parallel fitness gains, we characterize the landscape of transcriptional and translational changes across 12 replicate populations evolving in parallel for 50,000 generations. By quantifying absolute changes in mRNA abundances, we show that not only do all evolved lines have more mRNAs but that this increase in mRNA abundance scales with cell size. We also find that despite few shared mutations at the genetic level, clones from replicate populations in the LTEE are remarkably similar in their gene expression patterns at both the transcriptional and translational levels. Furthermore, we show that the majority of the expression changes are due to changes at the transcriptional level with very few translational changes. Finally, we show how mutations in transcriptional regulators lead to consistent and parallel changes in the expression levels of downstream genes. These results deepen our understanding of the molecular mechanisms underlying complex adaptations and provide insights into the repeatability of evolution.

Bioinformatic Assessment of Factors Affecting the Correlation between Protein Abundance and Elongation Efficiency in Prokaryotes

Protein abundance is crucial for the majority of genetically regulated cell functions to act properly in prokaryotic organisms. Therefore, developing bioinformatic methods for assessing the efficiency of different stages of gene expression is of great importance for predicting the actual protein abundance. One of these steps is the evaluation of translation elongation efficiency based on mRNA sequence features, such as codon usage bias and mRNA secondary structure properties. In this study, we have evaluated correlation coefficients between experimentally measured protein abundance and predicted elongation efficiency characteristics for 26 prokaryotes, including non-model organisms, belonging to diverse taxonomic groups The algorithm for assessing elongation efficiency takes into account not only codon bias, but also number and energy of secondary structures in mRNA if those demonstrate an impact on predicted elongation efficiency of the ribosomal protein genes. The results show that, for a number of organisms, secondary structures are a better predictor of protein abundance than codon usage bias. The bioinformatic analysis has revealed several factors associated with the value of the correlation coefficient. The first factor is the elongation efficiency optimization type-the organisms whose genomes are optimized for codon usage only have significantly higher correlation coefficients. The second factor is taxonomical identity-bacteria that belong to the class Bacilli tend to have higher correlation coefficients among the analyzed set. The third is growth rate, which is shown to be higher for the organisms with higher correlation coefficients between protein abundance and predicted translation elongation efficiency. The obtained results can be useful for further improvement of methods for protein abundance prediction.

Relative specificity as an important consideration in the big data era

Technological breakthroughs such as high-throughput methods, genomics, single-cell studies, and machine learning have fundamentally transformed research and ushered in the big data era of biology. Nevertheless, current data collections, analyses, and modeling frequently overlook relative specificity, a crucial property of molecular interactions in biochemical systems. Relative specificity describes how, for example, an enzyme reacts with its many substrates at different rates, and how this discriminatory action alone is sufficient to modulate the substrates and downstream events. As a corollary, it is not only important to comprehensively identify an enzyme’s substrates, but also critical to quantitatively determine how the enzyme interacts with the substrates and to evaluate how it shapes subsequent biological outcomes. Genomics and high-throughput techniques have greatly facilitated the studies of relative specificity in the 21st century, and its functional significance has been demonstrated in complex biochemical systems including transcription, translation, protein kinases, RNA-binding proteins, and animal microRNAs (miRNAs), although it remains ignored in most work. Here we analyze recent findings in big data and relative specificity studies and explain how the incorporation of relative specificity concept might enhance our mechanistic understanding of gene functions, biological phenomena, and human diseases.

Flow-Seq Method: Features and Application in Bacterial Translation Studies

The Flow-seq method is based on using reporter construct libraries, where a certain element regulating the gene expression of fluorescent reporter proteins is represented in many thousands of variants. Reporter construct libraries are introduced into cells, sorted according to their fluorescence level, and then subjected to next-generation sequencing. Therefore, it turns out to be possible to identify patterns that determine the expression efficiency, based on tens and hundreds of thousands of reporter constructs in one experiment. This method has become common in evaluating the efficiency of protein synthesis simultaneously by multiple mRNA variants. However, its potential is not confined to this area. In the presented review, a comparative analysis of the Flow-seq method and other alternative approaches used for translation efficiency evaluation of mRNA was carried out; the features of its application and the results obtained by Flow-seq were also considered.

Ovine RAP1GAP and rBAT gene polymorphisms and their association with tail fat deposition in Hu sheep

Excessive fat deposition in the tail of sheep will affect its feed efficiency, which will increase the feeding cost. The purpose of this study was to identify the single nucleotide polymorphisms (SNPs) of RAP1GAP and rBAT genes by PCR amplification and Sanger sequencing, the SNPs were genotyped by KASP genotyping assays to evaluate their association with tail fat deposition traits. The results showed that two intronic mutations of g.13561 G > A and g.1460 T > C were found in RAP1GAP and rBAT, respectively. There were three genotypes of GG, AG, AA and CC, CT and TT at these two loci, respectively. Association analysis showed that g.13561 G > A of RAP1GAP was associated with tail width, tail fat weight and relative tail fat weight (P < 0.05). The g.1460 T > C of rBAT was associated with tail width and tail fat weight (P < 0.05). Different combinations of genotypes also differed significantly with tail fat deposition traits. In the tail fat tissue, the expression levels of RAP1GAP gene was significantly higher in small-tailed sheep than in big-tailed sheep, and the expression levels of rBAT gene was significantly higher in big-tailed sheep than in small-tailed sheep. In the liver, the expression levels of RAP1GAP and rBAT gene was significantly higher at 6 months than at 0 and 3 months. In conclusion, RAP1GAP and rBAT polymorphisms can be used as a candidate molecular marker to reduce tail fat deposition in sheep.

Deciphering the Rules of Ribosome Binding Site Differentiation in Context Dependence

The ribosome binding site (RBS) is a crucial element regulating translation. However, the activity of RBS is poorly predictable, because it is strongly affected by the local possible secondary structure, that is, context dependence. By the Flowseq technique, over 20 000 RBS variants were sorted and sequenced, and the translation of multiple genes under the same RBS was quantitatively characterized to evaluate the context dependence of each RBS variant in E. coli. Two regions, (-7 to -2) and (-17 to -12), of RBS were predicted with a higher possibility to pair with each other to slow down the translation initiation. Associations between phenotypes and the intrinsic factors suspected to affect translation efficiency and context dependence of the RBS, including nucleotide bias at each position, free energy, and conservation, were disentangled. The results showed that translation efficiency was influenced more significantly by conservation of the SD region (-16 to -8), while an AC-rich spacer region (-7 to -1) was associated with low context dependence. We confirmed these characteristics using a series of synthesized RBSs. The average correlation between multiple reporters was significantly higher for RBSs with an AC-rich spacer (0.714) compared with a GU-rich spacer (0.286). Overall, we proposed general design criteria to improve programmability and minimize context dependence of RBS. The characteristics unraveled here can be adapted to other bacteria for fine-tuning target-gene expression.

Analysis of evolutionary and genetic patterns in structural genes of primate lentiviruses

Background

Primate lentiviruses (HIV1, HIV2, and Simian immunodeficiency virus [SIV]) cause immune deficiency, encephalitis, and infectious anemia in mammals such as cattle, cat, goat, sheep, horse, and puma.

Objective

This study was designed and conducted with the main purpose of confirming the overall codon usage pattern of primate lentiviruses and exploring the evolutionary and genetic characteristics commonly or specifically expressed in HIV1, HIV2, and SIV.

Methods

The gag, pol, and env gene sequences of HIV1, HIV2, and SIV were analyzed to determine their evolutionary relationships, nucleotide compositions, codon usage patterns, neutrality, selection pressure (influence of mutational pressure and natural selection), and viral adaptation to human codon usage.

Results

A strong ‘A’ bias was confirmed in all three structural genes, consistent with previous findings regarding HIV. Notably, the ENC-GC3s plot and neutral evolution analysis showed that all primate lentiviruses were more affected by selection pressure than by mutation caused by the GC composition of the gene, consistent with prior reports regarding HIV1. The overall codon usage bias of pol was highest among the structural genes, while the codon usage bias of env was lowest. The virus groups showing high codon bias in all three genes were HIV1 and SIVcolobus. The codon adaptation index (CAI) and similarity D(A, B) values indicated that although there was a high degree of similarity to human codon usage in all three structural genes of HIV, this similarity was not caused by translation pressure. In addition, compared with HIV1, the codon usage of HIV2 is more similar to the human codon usage, but the overall codon usage bias is lower.

Conclusion

The origin viruses of HIV (SIVcpz_gor and SIVsmm) exhibit greater similarity to human codon usage in the gag gene, confirming their robust adaptability to human codon usage. Therefore, HIV1 and HIV2 may have evolved to avoid human codon usage by selection pressure in the gag gene after interspecies transmission from SIV hosts to humans. By overcoming safety and stability issues, information from codon usage analysis will be useful for attenuated HIV1 vaccine development. A recoded HIV1 variant can be used as a vaccine vector or in immunotherapy to induce specific innate immune responses. Further research regarding HIV1 dinucleotide usage and codon pair usage will facilitate new approaches to the treatment of AIDS.

Enhanced extracellular raw starch-degrading α-amylase production in Bacillus subtilis by promoter engineering and translation initiation efficiency optimization

Background

A raw starch-degrading α-amylase from Pontibacillus sp. ZY (AmyZ1), previously screened by our laboratory, showed a promising application potential for starch-processing industries. However, the AmyZ1 secretory production still under investigation, which seriously restricts its application in the starch-processing industry. On the other hand, Bacillus subtilis is widely used to achieve the extracellular expression of target proteins.

Results

AmyZ1 secretory production was achieved in B. subtilis and was enhanced by promoter engineering and translation initiation efficiency optimization. First, based on the different phase-dependent promoters, the dual-promoter P_spoVG–P_spoVG142 was constructed by combining dual-promoter engineering and promoter modification. The corresponding strain BZd34 showed an extracellular AmyZ1 activity of 1437.6 U/mL during shake flask cultivation, which was 3.11-fold higher than that of the original strain BZ1 (P_groE). Then, based on translation initiation efficiency optimization, the best strain BZd343 containing optimized 5'-proximal coding sequence (opt3) produced the highest extracellular α-amylase activity of 1691.1 U/mL, which was 3.65-fold higher than that of the strain BZ1. Finally, cultivation of BZd343 in 3-L fermenter exhibited an extracellular AmyZ1 activity of 14,012 U/mL at 48 h, with productivity of 291.9 U/mL·h.

Conclusions

This is the first report of recombinant expression of AmyZ1 in B. subtilis and the expression level of AmyZ1 represents the highest raw starch-degrading α-amylase level in B. subtilis to date. The high-level expression of AmyZ1 in this work provides a foundation for its industrial production. The strategies used in this study also provide a strategic reference for improving the secretory expression of other enzymes in B. subtilis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12934-022-01855-9.

The Ramp Atlas: facilitating tissue and cell-specific ramp sequence analyses through an intuitive web interface

Ramp sequences occur when the average translational efficiency of codons near the 5′ end of highly expressed genes is significantly lower than the rest of the gene sequence, which counterintuitively increases translational efficiency by decreasing downstream ribosomal collisions. Here, we show that the relative codon adaptiveness within different tissues changes the existence of a ramp sequence without altering the underlying genetic code. We present the first comprehensive analysis of tissue and cell type-specific ramp sequences and report 3108 genes with ramp sequences that change between tissues and cell types, which corresponds with increased gene expression within those tissues and cells. The Ramp Atlas (https://ramps.byu.edu/) allows researchers to query precomputed ramp sequences in 18 388 genes across 62 tissues and 66 cell types and calculate tissue-specific ramp sequences from user-uploaded FASTA files through an intuitive web interface. We used The Ramp Atlas to identify seven SARS-CoV-2 genes and seven human SARS-CoV-2 entry factor genes with tissue-specific ramp sequences that may help explain viral proliferation within those tissues. We anticipate that The Ramp Atlas will facilitate personalized and creative tissue-specific ramp sequence analyses for both human and viral genes that will increase our ability to utilize this often-overlooked regulatory region.

Random peptides rich in small and disorder-promoting amino acids are less likely to be harmful

Proteins are the workhorses of the cell, yet they carry great potential for harm via misfolding and aggregation. Despite the dangers, proteins are sometimes born de novo from non-coding DNA. Proteins are more likely to be born from non-coding regions that produce peptides that do little to no harm when translated than from regions that produce harmful peptides. To investigate which newborn proteins are most likely to "first, do no harm", we estimate fitnesses from an experiment that competed Escherichia coli lineages that each expressed a unique random peptide. A variety of peptide metrics significantly predict lineage fitness, but this predictive power stems from simple amino acid frequencies rather than the ordering of amino acids. Amino acids that are smaller and that promote intrinsic structural disorder have more benign fitness effects. We validate that the amino acids that indicate benign effects in random peptides expressed in E. coli also do so in an independent dataset of random N-terminal tags in which it is possible to control for expression level. The same amino acids are also enriched in young animal proteins.

Model-Guided Metabolic Rewiring for Gamma-Aminobutyric Acid and Butyrolactam Biosynthesis in Corynebacterium glutamicum ATCC13032

Gamma-aminobutyric acid (GABA) can be used as a bioactive component in the pharmaceutical industry and a precursor for the synthesis of butyrolactam, which functions as a monomer for the synthesis of polyamide 4 (nylon 4) with improved thermal stability and high biodegradability. The bio-based fermentation production of chemicals using microbes as a cell factory provides an alternative to replace petrochemical-based processes. Here, we performed model-guided metabolic engineering of Corynebacterium glutamicum for GABA and butyrolactam fermentation. A GABA biosynthetic pathway was constructed using a bi-cistronic expression cassette containing mutant glutamate decarboxylase. An in silico simulation showed that the increase in the flux from acetyl-CoA to α-ketoglutarate and the decrease in the flux from α-ketoglutarate to succinate drove more flux toward GABA biosynthesis. The TCA cycle was reconstructed by increasing the expression of acn and icd genes and deleting the sucCD gene. Blocking GABA catabolism and rewiring the transport system of GABA further improved GABA production. An acetyl-CoA-dependent pathway for in vivo butyrolactam biosynthesis was constructed by overexpressing act-encoding ß-alanine CoA transferase. In fed-batch fermentation, the engineered strains produced 23.07 g/L of GABA with a yield of 0.52 mol/mol from glucose and 4.58 g/L of butyrolactam. The metabolic engineering strategies can be used for genetic modification of industrial strains to produce target chemicals from α-ketoglutarate as a precursor, and the engineered strains will be useful to synthesize the bio-based monomer of polyamide 4 from renewable resources.

Translation initiation site of mRNA is selected through dynamic interaction with the ribosome

SignificanceRibosomes translate the genetic codes of messenger RNA (mRNA) to make proteins. Translation must begin at the correct initiation site; otherwise, abnormal proteins will be produced. Here, we show that a short ribosome-specific sequence in the upstream followed by an unstructured downstream sequence is a favorable initiation site. Those mRNAs lacking either of these two characteristics do not associate tightly with the ribosome. Initiator transfer RNA (tRNA) and initiation factors facilitate the binding. However, when the downstream site forms structures, initiation factor 3 triggers the dissociation of the accommodated initiator tRNA and the subsequent disassembly of the ribosome-mRNA complex. Thus, initiation factors help the ribosome distinguish unfavorable structured sequences that may not act as the mRNA translation initiation site.

Oxygen level regulates N-terminal translation elongation of selected proteins through deoxyhypusine hydroxylation

Hypusine is a post-translational modification on eukaryotic translation initiation factor 5A (eIF5A). The last step of hypusine biosynthesis, deoxyhypusine hydroxylation, is an oxygen-dependent reaction. Here we show that deletion of the deoxyhypusine hydroxylase Lia1 compromises yeast respiration through translation downregulation of selected proteins in the respiration pathway. The translation suppression, because of the lack of deoxyhypusine hydroxylation, mainly affects translation of the N termini of the proteins, independent of the presence of proline residues but likely dependent on the interaction between the N-terminal nascent peptide and the ribosomal peptide exit tunnel. Proteomics and biochemical studies reveal that Lia1 deletion decreases N-terminal translation of proteins involved in mitochondrial respiration, oxidative stress response, and protein folding. Our work uncovers functions of the hypusine modification by considering the substrate requirement of the post-translational modification, highlights the unique challenges of translating the N termini of proteins, and reveals an oxygen-sensing mechanism in eukaryotic cells.

Optimized Heterologous Expression and Efficient Purification of a New TRAIL-Based Antitumor Fusion Protein SRH-DR5-B with Dual VEGFR2 and DR5 Receptor Specificity

In the last two decades, bifunctional proteins have been created by genetic and protein engineering methods to increase therapeutic effects in various diseases, including cancer. Unlike conventional small molecule or monotargeted drugs, bifunctional proteins have increased biological activity while maintaining low systemic toxicity. The recombinant anti-cancer cytokine TRAIL has shown a limited therapeutic effect in clinical trials. To enhance the efficacy of TRAIL, we designed the HRH–DR5-B fusion protein based on the DR5-selective mutant variant of TRAIL fused to the anti-angiogenic synthetic peptide HRHTKQRHTALH. Initially low expression of HRH–DR5-B was enhanced by the substitution of E. coli-optimized codons with AT-rich codons in the DNA sequence encoding the first 7 amino acid residues of the HRH peptide. However, the HRH–DR5-B degraded during purification to form two adjacent protein bands on the SDS-PAGE gel. The replacement of His by Ser at position P2 immediately after the initiator Met dramatically minimized degradation, allowing more than 20 mg of protein to be obtained from 200 mL of cell culture. The resulting SRH–DR5-B fusion bound the VEGFR2 and DR5 receptors with high affinity and showed increased cytotoxic activity in 3D multicellular tumor spheroids. SRH–DR5-B can be considered as a promising candidate for therapeutic applications.

DeepTESR: A Deep Learning Framework to Predict the Degree of Translational Elongation Short Ramp for Gene Expression Control

Controlling translational elongation is essential for efficient protein synthesis. Ribosome profiling has revealed that the speed of ribosome movement is correlated with translational efficiency in the translational elongation ramp. In this work, we present a new deep learning model, called DeepTESR, to predict the degree of translational elongation short ramp (TESR) from mRNA sequence. The proposed deep learning model exhibited superior performance in predicting the TESR scores for 226 981 TESR sequences, resulting in the mean absolute error (MAE) of 0.285 and a coefficient of determination R² of 0.627, superior to the conventional machine learning models (e.g., MAE of 0.335 and R² of 0.571 for LightGBM). We experimentally validated that heterologous fluorescence expression of proteins with randomly selected TESR was moderately correlated with the predictions. Furthermore, a genome-wide analysis of TESR prediction in the 4305 coding sequences of Escherichia coli showed conserved TESRs over the clusters of orthologous groups. In this sense, DeepTESR can be used to predict the degree of TESR for gene expression control and to decipher the mechanism of translational control with ribosome profiling. DeepTESR is available at https://github.com/fmblab/DeepTESR.

Modeling the effect of rRNA-mRNA interactions and mRNA folding on mRNA translation in chloroplasts

The process of translation initiation in prokaryotes is mediated by the hybridization of the 16S rRNA of the small ribosomal subunit with the mRNA in a short region called the ribosomal binding site. However, translation initiation in chloroplasts, which have evolved from an ancestral bacterium, is not well understood. Some studies suggest that in many cases it differs from translation initiation in bacteria and involves various novel interactions of the mRNA structures with intracellular factors; however currently, there is no generic quantitative model related to these aspects in chloroplasts. We developed a novel computational pipeline and models that can be used for understanding and modeling translation regulation in chloroplasts. We demonstrate that local folding and co-folding energy of the rRNA and the mRNA correlates with codon usage estimators of expression levels (r = -0.63) and infer predictive models that connect these energies and codon usage to protein levels (with correlation up to 0.71). In addition, we demonstrate that the ends of the transcripts in chloroplasts are populated with various structural elements that may be functional. Furthermore, we report a database of 166 novel structures in the chloroplast transcripts that are predicted to be functional. We believe that the models reported here improve existing understandings of genomic evolution and the biophysics of translation in chloroplasts; as such, they can aid gene expression engineering in chloroplasts for various biotechnological objectives.

Active secretion of a thermostable transglutaminase variant in Escherichia coli

Background

Streptomyces mobaraenesis transglutaminase (smTG) is widely used to generate protein crosslinking or attachment of small molecules. However, the low thermostability is a main obstacle for smTG application. In addition, it is still hard to achieve the secretory expression of active smTG in E. coli, which benefits the enzyme evolution. In this study, a combined strategy was conducted to improve the thermostability and secretory expression of active smTG in E. coli.

Results

First, the thermostable S. mobaraenesis transglutaminase variant S2P-S23V-Y24N-S199A-K294L (TGm1) was intracellularly expressed in pro-enzyme form in E. coli. Fusing the pro-region of Streptomyces hygroscopicus transglutaminase (proH) and TrxA achieved a 9.78 U/mL of intracellular smTG activity, 1.37-fold higher than the TGm1 fused with its native pro-region. After in vitro activation by dispase, the TGm1 with proH yielded FRAPD-TGm1, exhibiting 0.95 ℃ and 94.25% increases in melting temperature and half-life at 60 ℃ compared to FRAP-TGm1 derived from the expression using its native pro-region, respectively. Second, the TGm1 with proH was co-expressed with transglutaminase activating protease and chaperones (DnaK, DnaJ, and GrpE) in E. coli, achieving 9.51 U/mL of intracellular FRAPD-TGm1 without in vitro activation. Third, the pelB signal peptide was used to mediate the secretory expression of active TGm in E. coli, yielding 0.54 U/mL of the extracellular FRAPD-TGm1. A script was developed to shuffle the codon of pelB and calculate the corresponding mRNA folding energy. A 1.8-fold increase in the extracellular expression of FRAPD-TGm1 was achieved by the Top-9 pelB sequence derived from the coding sequences with the lowest mRNA folding energy. Last, deleting the gene of Braun’s lipoprotein further increased the extracellular yield of FRAPD-TGm1 by 31.2%, reached 1.99 U/mL.

Conclusions

The stabilized FRAPD-smTG here could benefit the enzyme application in food and non-food sectors, while the E. coli system that enables secretory expression of active smTG will facilitate the directed evolution for further improved catalytic properties. The combined strategy (N-terminal modification, co-expression with chaperones, mRNA folding energy optimization of signal peptide, and lipoprotein deletion) may also improve the secretory expression of other functional proteins in E. coli.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12934-022-01801-9.

Cross-kingdom expression of synthetic genetic elements promotes discovery of metabolites in the human microbiome

Small molecules encoded by biosynthetic pathways mediate cross-species interactions and harbor untapped potential, which has provided valuable compounds for medicine and biotechnology. Since studying biosynthetic gene clusters in their native context is often difficult, alternative efforts rely on heterologous expression, which is limited by host-specific metabolic capacity and regulation. Here, we describe a computational-experimental technology to redesign genes and their regulatory regions with hybrid elements for cross-species expression in Gram-negative and -positive bacteria and eukaryotes, decoupling biosynthetic capacity from host-range constraints to activate silenced pathways. These synthetic genetic elements enabled the discovery of a class of microbiome-derived nucleotide metabolites-tyrocitabines-from Lactobacillus iners. Tyrocitabines feature a remarkable orthoester-phosphate, inhibit translational activity, and invoke unexpected biosynthetic machinery, including a class of "Amadori synthases" and "abortive" tRNA synthetases. Our approach establishes a general strategy for the redesign, expression, mobilization, and characterization of genetic elements in diverse organisms and communities.

A genetic toolkit and gene switches to limit Mycoplasma growth for biosafety applications

Mycoplasmas have exceptionally streamlined genomes and are strongly adapted to their many hosts, which provide them with essential nutrients. Owing to their relative genomic simplicity, Mycoplasmas have been used to develop chassis for biotechnological applications. However, the dearth of robust and precise toolkits for genomic manipulation and tight regulation has hindered any substantial advance. Herein we describe the construction of a robust genetic toolkit for M. pneumoniae, and its successful deployment to engineer synthetic gene switches that control and limit Mycoplasma growth, for biosafety containment applications. We found these synthetic gene circuits to be stable and robust in the long-term, in the context of a minimal cell. With this work, we lay a foundation to develop viable and robust biosafety systems to exploit a synthetic Mycoplasma chassis for live attenuated vectors for therapeutic applications.

QRNAstruct: a method for extracting secondary structural features of RNA via regression with biological activity

Recent technological advances have enabled the generation of large amounts of data consisting of RNA sequences and their functional activity. Here, we propose a method for extracting secondary structure features that affect the functional activity of RNA from sequence–activity data. Given pairs of RNA sequences and their corresponding bioactivity values, our method calculates position-specific structural features of the input RNA sequences, considering every possible secondary structure of each RNA. A Ridge regression model is trained using the structural features as feature vectors and the bioactivity values as response variables. Optimized model parameters indicate how secondary structure features affect bioactivity. We used our method to extract intramolecular structural features of bacterial translation initiation sites and self-cleaving ribozymes, and the intermolecular features between rRNAs and Shine–Dalgarno sequences and between U1 RNAs and splicing sites. We not only identified known structural features but also revealed more detailed insights into structure–activity relationships than previously reported. Importantly, the datasets we analyzed here were obtained from different experimental systems and differed in size, sequence length and similarity, and number of RNA molecules involved, demonstrating that our method is applicable to various types of data consisting of RNA sequences and bioactivity values.

Calculating and comparing codon usage values in rare disease genes highlights codon clustering with disease-and tissue- specific hierarchy

We designed a novel strategy to define codon usage bias (CUB) in 6 specific small cohorts of human genes. We calculated codon usage (CU) values in 29 non-disease-causing (NDC) and 31 disease-causing (DC) human genes which are highly expressed in 3 distinct tissues, kidney, muscle, and skin. We applied our strategy to the same selected genes annotated in 15 mammalian species. We obtained CUB hierarchical clusters for each gene cohort which showed tissue-specific and disease-specific CUB fingerprints. We showed that DC genes (especially those expressed in muscle) display a low CUB, well recognizable in codon hierarchical clustering. We defined the extremely biased codons as "zero codons" and found that their number is significantly higher in all DC genes, all tissues, and that this trend is conserved across mammals. Based on this calculation in different gene cohorts, we identified 5 codons which are more differentially used across genes and mammals, underlining that some genes have favorite synonymous codons in use. Since of the muscle genes clear clusters, and, among these, dystrophin gene surprisingly does not show any "zero codon" we adopted a novel approach to study CUB, we called "mapping-on-codons". We positioned 2828 dystrophin missense and nonsense pathogenic variations on their respective codon, highlighting that its frequency and occurrence is not dependent on the CU values. We conclude our strategy consents to identify a hierarchical clustering of CU values in a gene cohort-specific fingerprints, with recognizable trend across mammals. In DC muscle genes also a disease-related fingerprint can be observed, allowing discrimination between DC and NDC genes. We propose that using our strategy which studies CU in specific gene cohorts, as rare disease genes, and tissue specific genes, may provide novel information about the CUB role in human and medical genetics, with implications on synonymous variations interpretation and codon optimization algorithms.

Codon-Restrained Method for Both Eliminating and Creating Intragenic Bacterial Promoters

Future applications of synthetic biology will require refactored genetic sequences devoid of internal regulatory elements within coding sequences. These regulatory elements include cryptic and intragenic promoters, which may constitute up to a third of the predicted Escherichia coli promoters. The promoter activity is dependent on the structural interaction of core bases with a σ factor. Rational engineering can be used to alter key promoter element nucleotides interacting with σ factors and eliminate downstream transcriptional activity. In this paper, we present codon-restrained promoter silencing (CORPSE), a system for removing intragenic promoters. CORPSE exploits the DNA-σ factor structural relationship to disrupt σ⁷⁰ promoters embedded within gene coding sequences with a minimum of synonymous codon changes. Additionally, we present an inverted CORPSE system, iCORPSE, which can create highly active promoters within a gene sequence while not perturbing the function of the modified gene.

Modeling the ribosomal small subunit dynamic in Saccharomyces cerevisiae based on TCP-seq data

Translation Complex Profile Sequencing (TCP-seq), a protocol that was developed and implemented on Saccharomyces cerevisiae, provides the footprints of the small subunit (SSU) of the ribosome (with additional factors) across the entire transcriptome of the analyzed organism. In this study, based on the TCP-seq data, we developed for the first-time a predictive model of the SSU density and analyzed the effect of transcript features on the dynamics of the SSU scan in the 5′UTR. Among others, our model is based on novel tools for detecting complex statistical relations tailored to TCP-seq. We quantitatively estimated the effect of several important features, including the context of the upstream AUG, the upstream ORF length and the mRNA folding strength. Specifically, we suggest that around 50% of the variance related to the read counts (RC) distribution near a start codon can be attributed to the AUG context score. We provide the first large scale direct quantitative evidence that shows that indeed AUG context affects the small sub-unit movement. In addition, we suggest that strong folding may cause the detachment of the SSU from the mRNA. We also identified a number of novel sequence motifs that can affect the SSU scan; some of these motifs affect transcription factors and RNA binding proteins. The results presented in this study provide a better understanding of the biophysical aspects related to the SSU scan along the 5′UTR and of translation initiation in S. cerevisiae, a fundamental step toward a comprehensive modeling of initiation.

Directed evolution of the rRNA methylating enzyme Cfr reveals molecular basis of antibiotic resistance

Alteration of antibiotic binding sites through modification of ribosomal RNA (rRNA) is a common form of resistance to ribosome-targeting antibiotics. The rRNA-modifying enzyme Cfr methylates an adenosine nucleotide within the peptidyl transferase center, resulting in the C-8 methylation of A2503 (m8A2503). Acquisition of cfr results in resistance to eight classes of ribosome-targeting antibiotics. Despite the prevalence of this resistance mechanism, it is poorly understood whether and how bacteria modulate Cfr methylation to adapt to antibiotic pressure. Moreover, direct evidence for how m8A2503 alters antibiotic binding sites within the ribosome is lacking. In this study, we performed directed evolution of Cfr under antibiotic selection to generate Cfr variants that confer increased resistance by enhancing methylation of A2503 in cells. Increased rRNA methylation is achieved by improved expression and stability of Cfr through transcriptional and post-transcriptional mechanisms, which may be exploited by pathogens under antibiotic stress as suggested by natural isolates. Using a variant that achieves near-stoichiometric methylation of rRNA, we determined a 2.2 Å cryo-electron microscopy structure of the Cfr-modified ribosome. Our structure reveals the molecular basis for broad resistance to antibiotics and will inform the design of new antibiotics that overcome resistance mediated by Cfr.

Strategies to Enhance Periplasmic Recombinant Protein Production Yields in Escherichia coli

Main reasons to produce recombinant proteins in the periplasm of E. coli rather than in its cytoplasm are to -i- enable disulfide bond formation, -ii- facilitate protein isolation, -iii- control the nature of the N-terminus of the mature protein, and -iv- minimize exposure to cytoplasmic proteases. However, hampered protein targeting, translocation and folding as well as protein instability can all negatively affect periplasmic protein production yields. Strategies to enhance periplasmic protein production yields have focused on harmonizing secretory recombinant protein production rates with the capacity of the secretory apparatus by transcriptional and translational tuning, signal peptide selection and engineering, increasing the targeting, translocation and periplasmic folding capacity of the production host, preventing proteolysis, and, finally, the natural and engineered adaptation of the production host to periplasmic protein production. Here, we discuss these strategies using notable examples as a thread.

Nascent polypeptide within the exit tunnel stabilizes the ribosome to counteract risky translation

Continuous translation elongation, irrespective of amino acid sequences, is a prerequisite for living organisms to produce their proteomes. However, nascent polypeptide products bear an inherent risk of elongation abortion. For example, negatively charged sequences with occasional intermittent prolines, termed intrinsic ribosome destabilization (IRD) sequences, weaken the translating ribosomal complex, causing certain nascent chain sequences to prematurely terminate translation. Here, we show that most potential IRD sequences in the middle of open reading frames remain cryptic and do not interrupt translation, due to two features of the nascent polypeptide. Firstly, the nascent polypeptide itself spans the exit tunnel, and secondly, its bulky amino acid residues occupy the tunnel entrance region, thereby serving as a bridge and protecting the large and small ribosomal subunits from dissociation. Thus, nascent polypeptide products have an inbuilt ability to ensure elongation continuity.

Genetic fusion of P450 BM3 and formate dehydrogenase towards self-sufficient biocatalysts with enhanced activity

Fusion of multiple enzymes to multifunctional constructs has been recognized as a viable strategy to improve enzymatic properties at various levels such as stability, activity and handling. In this study, the genes coding for cytochrome P450 BM3 from B. megaterium and formate dehydrogenase from Pseudomonas sp. were fused to enable both substrate oxidation catalyzed by P450 BM3 and continuous cofactor regeneration by formate dehydrogenase within one construct. The order of the genes in the fusion as well as the linkers that bridge the enzymes were varied. The resulting constructs were compared to individual enzymes regarding substrate conversion, stability and kinetic parameters to examine whether fusion led to any substantial improvements of enzymatic properties. Most noticeably, an activity increase of up to threefold was observed for the fusion constructs with various substrates which were partly attributed to the increased diflavin reductase activity of the P450 BM3. We suggest that P450 BM3 undergoes conformational changes upon fusion which resulted in altered properties, however, no NADPH channeling was detected for the fusion constructs.

Importance of the 5' regulatory region to bacterial synthetic biology applications

The field of synthetic biology is evolving at a fast pace. It is advancing beyond single-gene alterations in single hosts to the logical design of complex circuits and the development of integrated synthetic genomes. Recent breakthroughs in deep learning, which is increasingly used in de novo assembly of DNA components with predictable effects, are also aiding the discipline. Despite advances in computing, the field is still reliant on the availability of pre-characterized DNA parts, whether natural or synthetic, to regulate gene expression in bacteria and make valuable compounds. In this review, we discuss the different bacterial synthetic biology methodologies employed in the creation of 5' regulatory regions - promoters, untranslated regions and 5'-end of coding sequences. We summarize methodologies and discuss their significance for each of the functional DNA components, and highlight the key advances made in bacterial engineering by concentrating on their flaws and strengths. We end the review by outlining the issues that the discipline may face in the near future.

A single synonymous nucleotide change impacts the male-killing phenotype of prophage WO gene wmk

Wolbachia are the most widespread bacterial endosymbionts in animals. Within arthropods, these maternally transmitted bacteria can selfishly hijack host reproductive processes to increase the relative fitness of their transmitting females. One such form of reproductive parasitism called male killing, or the selective killing of infected males, is recapitulated to degrees by transgenic expression of the prophage WO-mediated killing (wmk) gene. Here, we characterize the genotype-phenotype landscape of wmk-induced male killing in D. melanogaster using transgenic expression. While phylogenetically distant wmk homologs induce no sex-ratio bias, closely-related homologs exhibit complex phenotypes spanning no death, male death, or death of all hosts. We demonstrate that alternative start codons, synonymous codons, and notably a single synonymous nucleotide in wmk can ablate killing. These findings reveal previously unrecognized features of transgenic wmk-induced killing and establish new hypotheses for the impacts of post-transcriptional processes in male killing variation. We conclude that synonymous sequence changes are not necessarily silent in nested endosymbiotic interactions with life-or-death consequences.

A nascent polypeptide sequence modulates DnaA translation elongation in response to nutrient availability

The ability to regulate DNA replication initiation in response to changing nutrient conditions is an important feature of most cell types. In bacteria, DNA replication is triggered by the initiator protein DnaA, which has long been suggested to respond to nutritional changes; nevertheless, the underlying mechanisms remain poorly understood. Here, we report a novel mechanism that adjusts DnaA synthesis in response to nutrient availability in Caulobacter crescentus. By performing a detailed biochemical and genetic analysis of the dnaA mRNA, we identified a sequence downstream of the dnaA start codon that inhibits DnaA translation elongation upon carbon exhaustion. Our data show that the corresponding peptide sequence, but not the mRNA secondary structure or the codon choice, is critical for this response, suggesting that specific amino acids in the growing DnaA nascent chain tune translational efficiency. Our study provides new insights into DnaA regulation and highlights the importance of translation elongation as a regulatory target. We propose that translation regulation by nascent chain sequences, like the one described, might constitute a general strategy for modulating the synthesis rate of specific proteins under changing conditions.

Variability in Codon Usage in Coronaviruses Is Mainly Driven by Mutational Bias and Selective Constraints on CpG Dinucleotide

The Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the third human-emerged virus of the 21st century from the Coronaviridae family, causing the ongoing coronavirus disease 2019 (COVID-19) pandemic. Due to the high zoonotic potential of coronaviruses, it is critical to unravel their evolutionary history of host species breadth, host-switch potential, adaptation and emergence, to identify viruses posing a pandemic risk in humans. We present here a comprehensive analysis of the composition and codon usage bias of the 82 Orthocoronavirinae members, infecting 47 different avian and mammalian hosts. Our results clearly establish that synonymous codon usage varies widely among viruses, is only weakly dependent on their primary host, and is dominated by mutational bias towards AU-enrichment and by CpG avoidance. Indeed, variation in GC3 explains around 34%, while variation in CpG frequency explains around 14% of total variation in codon usage bias. Further insight on the mutational equilibrium within Orthocoronavirinae revealed that most coronavirus genomes are close to their neutral equilibrium, the exception being the three recently infecting human coronaviruses, which lie further away from the mutational equilibrium than their endemic human coronavirus counterparts. Finally, our results suggest that, while replicating in humans, SARS-CoV-2 is slowly becoming AU-richer, likely until attaining a new mutational equilibrium.

Strategies and Applications of Antigen-Binding Fragment (Fab) Production in Escherichia coli

With the advancement of genetic engineering, monoclonal antibodies (mAbs) have made far-reaching progress in the treatment of various human diseases. However, due to the high cost of production, the increasing demands for antibody-based therapies have not been fully met. Currently, mAb-derived alternatives, such as antigen-binding fragments (Fab), single-chain variable fragments, bispecifics, nanobodies, and conjugated mAbs have emerged as promising new therapeutic modalities. They can be readily prepared in bacterial systems with well-established fermentation technology and ease of manipulation, leading to the reduction of overall cost. This review aims to shed light on the strategies to improve the expression, purification, and yield of Fab fragments in Escherichia coli expression systems, as well as current advances in the applications of Fab fragments.

A 68-codon genetic code to incorporate four distinct non-canonical amino acids enabled by automated orthogonal mRNA design

Orthogonal (O) ribosome-mediated translation of O-mRNAs enables the incorporation of up to three distinct non-canonical amino acids (ncAAs) into proteins in Escherichia coli (E. coli). However, the general and efficient incorporation of multiple distinct ncAAs by O-ribosomes requires scalable strategies for both creating efficiently and specifically translated O-mRNAs, and the compact expression of multiple O-aminoacyl-tRNA synthetase (O-aaRS)/O-tRNA pairs. We automate the discovery of O-mRNAs that lead to up to 40 times more protein, and are up to 50-fold more orthogonal, than previous O-mRNAs; protein yields from our O-mRNAs match or exceed those from wild-type mRNAs. These advances enable a 33-fold increase in yield for incorporating three distinct ncAAs. We automate the creation of operons for O-tRNA genes, and develop operons for O-aaRS genes. Combining our advances creates a 68-codon, 24-amino-acid genetic code to efficiently incorporate four distinct ncAAs into a single protein in response to four distinct quadruplet codons.

Implementing Novel Designs in pET Expression Plasmids that Increase Protein Production

pET expression plasmids are widely used in the biotechnology, biopharmaceutical, and basic research sectors for the production of recombinant proteins. Typically, they are used off-the-shelf because they support high production titers; however, we have identified two design flaws in many pET plasmids that limit their production capacity. We used modern methods of DNA assembly and directed evolution to identify improved designs for these modules and demonstrated that these designs support higher protein production yields. Herein, we present two PCR protocols for implementing the designs and increasing protein production from existing pET expression plasmids. Graphic abstract: A simple workflow for implementing novel designs in pET expression plasmids.

Quantitative Control for Stoichiometric Protein Synthesis

Bacterial protein synthesis rates have evolved to maintain preferred stoichiometries at striking precision, from the components of protein complexes to constituents of entire pathways. Setting relative protein production rates to be well within a factor of two requires concerted tuning of transcription, RNA turnover, and translation, allowing many potential regulatory strategies to achieve the preferred output. The last decade has seen a greatly expanded capacity for precise interrogation of each step of the central dogma genome-wide. Here, we summarize how these technologies have shaped the current understanding of diverse bacterial regulatory architectures underpinning stoichiometric protein synthesis. We focus on the emerging expanded view of bacterial operons, which encode diverse primary and secondary mRNA structures for tuning protein stoichiometry. Emphasis is placed on how quantitative tuning is achieved. We discuss the challenges and open questions in the application of quantitative, genome-wide methodologies to the problem of precise protein production. Expected final online publication date for the Annual Review of Microbiology, Volume 75 is October 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Effects of synonymous mutations beyond codon bias: The evidence for adaptive synonymous substitutions from microbial evolution experiments

Synonymous mutations are often assumed to be neutral with respect to fitness because they do not alter the encoded amino acid and so cannot be 'seen' by natural selection. Yet a growing body of evidence suggests that synonymous mutations can have fitness effects that drive adaptive evolution through their impacts on gene expression and protein folding. Here, we review what microbial experiments have taught us about the contribution of synonymous mutations to adaptation. A survey of site-directed mutagenesis experiments reveals the distributions of fitness effects for nonsynonymous and synonymous mutations are more similar, especially for beneficial mutations, than expected if all synonymous mutations were neutral, suggesting they should drive adaptive evolution more often than is typically observed. A review of experimental evolution studies where synonymous mutations have contributed to adaptation shows they can impact fitness through a range of mechanisms including the creation of illicit RNA polymerase binding sites impacting transcription and changes to mRNA folding stability that modulate translation. We suggest that clonal interference in evolving microbial populations may be the reason synonymous mutations play a smaller role in adaptive evolution than expected based on their observed fitness effects. We finish by discussing the impacts of falsely assuming synonymous mutations are neutral and discuss directions for future work exploring the role of synonymous mutations in adaptive evolution.

The impact of COMT, BDNF and 5-HTT brain-genes on the development of anorexia nervosa: a systematic review

PURPOSE

The genetic aspect of anorexia nervosa (AN) involving specific genes of the central-nervous-system has not yet been clearly understood. The aim of this systematic review is to assess the impact of three candidate genes of the brain: catechol-O-methyltransferase, brain-derived neurotrophic factor (BDNF) and serotonin transporter protein, on the susceptibility to AN and identify whether a clear connection persists between each of the gene-polymorphisms and AN.

METHODS

A total of 21 articles were selected for this review conforming to the PRISMA guidelines. Detailed keyword combinations were implemented within specific databases such as MEDLINE, SCIENCEDIRECT and PUBMED.

RESULTS

The catechol-O-methyltransferase gene-polymorphism did not show any change in phenotypic variation between AN and control subjects; but the familial association was rather strong with an over-transmission of the H allele. The latter also correlated with several dimensions of the Temperament and Character Inventory (TCI) scale. A notable relation was indicated between BDNF gene-polymorphism and anorexia-restrictive in terms of phenotypic distribution; the Met66-allele also depicted high association with anorexic behavioral traits. The 5-HTTLPR gene-polymorphism was found to be significantly associated with AN susceptibility with an over-transmission of the S-allele from parents to offspring.

CONCLUSION

The systematic review distinctively emphasized the genetic contribution of the brain-genes on the development of AN. Despite significant study findings, no clear and standardized genetic route was determined to be the cause of AN development. Future research is needed on these specific genes to closely monitor the genetic polymorphisms and their mechanism on AN susceptibility.

LEVEL OF EVIDENCE

I, systematic review.

Comparative analysis of the relationship between translation efficiency and sequence features of endogenous proteins in multiple organisms

The translation efficiency of protein genes is known to be affected by sequence features. Previous studies have found that various sequence features based on codon usage and mRNA secondary structure contribute to translation efficiency. However, most studies have focused on a specific organism, usually a model organism such as Escherichia coli or Saccharomyces cerevisiae. Here, we investigate whether the relationship between translation efficiency and sequence features is conserved among multiple organisms using publicly available ribosome profiling data and RNA-Seq data. We analyze nine organisms from various taxa: Staphylococcus aureus, five species of Streptomyces, two strains of E. coli, and S. cerevisiae. We reveal that the relationship between translation efficiency and sequence features differs across organisms, partly reflecting their taxonomy. The codon adaptation index shows high correlation in all analyzed organisms. Our study provides an insight into the diversity and commonality of sequence determinants of protein expression in these organisms.

Production of proteins and commodity chemicals using engineered Bacillus subtilis platform strain

Currently, increasing demand of biochemicals produced from renewable resources has motivated researchers to seek microbial production strategies instead of traditional chemical methods. As a microbial platform, Bacillus subtilis possesses many advantages including the generally recognized safe status, clear metabolic networks, short growth cycle, mature genetic editing methods and efficient protein secretion systems. Engineered B. subtilis strains are being increasingly used in laboratory research and in industry for the production of valuable proteins and other chemicals. In this review, we first describe the recent advances of bioinformatics strategies during the research and applications of B. subtilis. Secondly, the applications of B. subtilis in enzymes and recombinant proteins production are summarized. Further, the recent progress in employing metabolic engineering and synthetic biology strategies in B. subtilis platform strain to produce commodity chemicals is systematically introduced and compared. Finally, the major limitations for the further development of B. subtilis platform strain and possible future directions for its research are also discussed.

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.

Sigma factor dependent translational activation in Bacillus subtilis

Sigma factors are an important class of bacterial transcription factors that lend specificity to RNA polymerases by binding to distinct promoter elements for genes in their regulons. Here we show that activation of the general stress sigma factor, σB, in Bacillus subtilis paradoxically leads to dramatic induction of translation for a subset of its regulon genes. These genes are translationally repressed when transcribed by the housekeeping sigma factor, σA, owing to extended RNA secondary structures as determined in vivo using DMS-MaPseq. Transcription from σB-dependent promoters ablates the secondary structures and activates translation, leading to dual induction. Translation efficiencies between σB- and σA-dependent RNA isoforms can vary by up to 100-fold, which in multiple cases exceeds the magnitude of transcriptional induction. These results highlight the role of long-range RNA folding in modulating translation and demonstrate that a transcription factor can regulate protein synthesis beyond its effects on transcript levels.

A role for circular code properties in translation

Circular codes represent a form of coding allowing detection/correction of frame-shift errors. Building on recent theoretical advances on circular codes, we provide evidence that protein coding sequences exhibit in-frame circular code marks, that are absent in introns and are intimately linked to the keto-amino transformation of codon bases. These properties strongly correlate with translation speed, codon influence and protein synthesis levels. Strikingly, circular code marks are absent at the beginning of coding sequences, but stably occur 40 codons after the initiator codon, hinting at the translation elongation process. Finally, we use the lens of circular codes to show that codon influence on translation correlates with the strong-weak dichotomy of the first two bases of the codon. The results can lead to defining new universal tools for sequence indicators and sequence optimization for bioinformatics and biotechnological applications, and can shed light on the molecular mechanisms behind the decoding process.

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.

Inferring Adaptive Codon Preference to Understand Sources of Selection Shaping Codon Usage Bias

Alternative synonymous codons are often used at unequal frequencies. Classically, studies of such codon usage bias (CUB) attempted to separate the impact of neutral from selective forces by assuming that deviations from a predicted neutral equilibrium capture selection. However, GC-biased gene conversion (gBGC) can also cause deviation from a neutral null. Alternatively, selection has been inferred from CUB in highly expressed genes, but the accuracy of this approach has not been extensively tested, and gBGC can interfere with such extrapolations (e.g., if expression and gene conversion rates covary). It is therefore critical to examine deviations from a mutational null in a species with no gBGC. To achieve this goal, we implement such an analysis in the highly AT rich genome of Dictyostelium discoideum, where we find no evidence of gBGC. We infer neutral CUB under mutational equilibrium to quantify "adaptive codon preference," a nontautologous genome wide quantitative measure of the relative selection strength driving CUB. We observe signatures of purifying selection consistent with selection favoring adaptive codon preference. Preferred codons are not GC rich, underscoring the independence from gBGC. Expression-associated "preference" largely matches adaptive codon preference but does not wholly capture the influence of selection shaping patterns across all genes, suggesting selective constraints associated specifically with high expression. We observe patterns consistent with effects on mRNA translation and stability shaping adaptive codon preference. Thus, our approach to quantifying adaptive codon preference provides a framework for inferring the sources of selection that shape CUB across different contexts within the genome.