A role for mRNA secondary structure in the control of translation initiation

Understanding bacteriocin heterologous expression: A review

Bacteriocins are a diverse group of ribosomally synthesised antimicrobial peptides/proteins that play an important role in self-defence. They are widely used as bio-preservatives and effective substitutes for disease eradication. They can be used in conjunction with or as an alternative to antibiotics to minimize the risk of resistance development. There are remarkably few reports indicating resistance to bacteriocins. Although there are many research reports that emphasise heterologous expression of bacteriocin, there are no convincing reports on the significant role that intrinsic and extrinsic factors play in overexpression. A coordinated and cooperative expression system works in concert with multiple genetic elements encoding native proteins, immunoproteins, exporters, transporters and enzymes involved in the post-translational modification of bacteriocins. The simplest way could be to utilise the existing E. coli expression system, which is conventional, widely used for heterologous expression and has been further extended for bacteriocin expression. In this article, we will review the intrinsic and extrinsic factors, advantages, disadvantages and major problems associated with bacteriocin overexpression in E. coli. Finally, we recommend the most effective strategies as well as numerous bacteriocin expression systems from E. coli, Lactococcus, Kluveromyces lactis, Saccharomyces cerevisiae and Pichia pastoris for their suitability for successful overexpression.

Enrichment of rare codons at 5' ends of genes is a spandrel caused by evolutionary sequence turnover and does not improve translation

Previously, Tuller et al. found that the first 30-50 codons of the genes of yeast and other eukaryotes are slightly enriched for rare codons. They argued that this slowed translation, and was adaptive because it queued ribosomes to prevent collisions. Today, the translational speeds of different codons are known, and indeed rare codons are translated slowly. We re-examined this 5' slow translation 'ramp.' We confirm that 5' regions are slightly enriched for rare codons; in addition, they are depleted for downstream Start codons (which are fast), with both effects contributing to slow 5' translation. However, we also find that the 5' (and 3') ends of yeast genes are poorly conserved in evolution, suggesting that they are unstable and turnover relatively rapidly. When a new 5' end forms de novo, it is likely to include codons that would otherwise be rare. Because evolution has had a relatively short time to select against these codons, 5' ends are typically slightly enriched for rare, slow codons. Opposite to the expectation of Tuller et al., we show by direct experiment that genes with slowly translated codons at the 5' end are expressed relatively poorly, and that substituting faster synonymous codons improves expression. Direct experiment shows that slow codons do not prevent downstream ribosome collisions. Further informatic studies suggest that for natural genes, slow 5' ends are correlated with poor gene expression, opposite to the expectation of Tuller et al. Thus, we conclude that slow 5' translation is a 'spandrel'--a non-adaptive consequence of something else, in this case, the turnover of 5' ends in evolution, and it does not improve translation.

Bicistronic expression of nuclear transgenes in Chlamydomonas reinhardtii

In eukaryotic organisms, proteins are typically translated from monocistronic messenger RNAs containing a single coding sequence (CDS). However, recent long transcript sequencing identified 87 nuclear polycistronic mRNAs in Chlamydomonas reinhardtii natively carrying multiple co-expressed CDSs. In this study, we investigated the dynamics of 22 short intergenic sequences derived from these native polycistronic loci by their application in genetic constructs for synthetic transgene expression. A promising candidate sequence was identified based on the quantification of transformation efficiency and expression strength of a fluorescence reporter protein. Subsequently, the expression of independent proteins from one mRNA was verified by cDNA amplification and protein molecular mass characterization. We demonstrated engineered bicistronic expression in vivo to drive successful co-expression of several terpene synthases with the selection marker aphVIII. Bicistronic transgene design resulted in significantly increased (E)-α-bisabolene production of 7.95 mg L-1 from a single open reading frame, 18.1× fold higher than previous reports. Use of this strategy simplifies screening procedures for identification of high-level expressing transformants, does not require the application of additional fluorescence reporters, and reduces the nucleotide footprint compared to classical monocistronic expression cassettes. Although clear advantages for bicistronic transgene expression were observed, this strategy was found to be limited to the aphVIII marker, and further studies are necessary to gain insights into the underlying mechanism that uniquely permits this co-expression from the algal nuclear genome.

An Integer Linear Programming Model to Optimize Coding DNA Sequences By Joint Control of Transcript Indicators

A Coding DNA Sequence (CDS) is a fraction of DNA whose nucleotides are grouped into consecutive triplets called codons, each one encoding an amino acid. Because most amino acids can be encoded by more than one codon, the same amino acid chain can be obtained by a very large number of different CDSs. These synonymous CDSs show different features that, also depending on the organism the transcript is expressed in, could affect translational efficiency and yield. The identification of optimal CDSs with respect to given transcript indicators is in general a challenging task, but it has been observed in recent literature that integer linear programming (ILP) can be a very flexible and efficient way to achieve it. In this article, we add evidence to this observation by proposing a new ILP model that simultaneously optimizes different well-grounded indicators. With this model, we efficiently find solutions that dominate those returned by six existing codon optimization heuristics.

Synonymous edits in the Escherichia coli genome have substantial and condition-dependent effects on fitness

CRISPR-Cas-based genome editing is widely used in bacteria at scales ranging from construction of individual mutants to massively parallel libraries. This procedure relies on guide RNA-directed cleavage of the genome followed by repair with a template that introduces a desired mutation along with synonymous "immunizing" mutations to prevent re-cleavage of the genome after editing. Because the immunizing mutations do not change the protein sequence, they are often assumed to be neutral. However, synonymous mutations can change mRNA structures in ways that alter levels of the encoded proteins. We have tested the assumption that immunizing mutations are neutral by constructing a library of over 50,000 edits that consist of only synonymous mutations in Escherichia coli. Thousands of edits had substantial effects on fitness during growth of E. coli on acetate, a poor carbon source that is toxic at high concentrations. The percentage of high-impact edits varied considerably between genes and at different positions within genes. We reconstructed clones with high-impact edits and found that 69% indeed had significant effects on growth in acetate. Interestingly, fewer edits affected fitness during growth in glucose, a preferred carbon source, suggesting that changes in protein expression caused by synonymous mutations may be most important when an organism encounters challenging conditions. Finally, we showed that synonymous edits can have widespread effects; a synonymous edit at the 5' end of ptsI altered expression of hundreds of genes. Our results suggest that the synonymous immunizing edits introduced during CRISPR-Cas-based genome editing should not be assumed to be innocuous.

Molecular Mechanisms and the Significance of Synonymous Mutations

Synonymous mutations result from the degeneracy of the genetic code. Most amino acids are encoded by two or more codons, and mutations that change a codon to another synonymous codon do not change the amino acid in the gene product. Historically, such mutations have been considered silent because they were assumed to have no to very little impact. However, research in the last few decades has produced several examples where synonymous mutations play important roles. These include optimizing expression by enhancing translation initiation and accelerating or decelerating translation elongation via codon usage and mRNA secondary structures, stabilizing mRNA molecules and preventing their breakdown before translation, and faulty protein folding or increased degradation due to enhanced ubiquitination and suboptimal secretion of proteins into the appropriate cell compartments. Some consequences of synonymous mutations, such as mRNA stability, can lead to different outcomes in prokaryotes and eukaryotes. Despite these examples, the significance of synonymous mutations in evolution and in causing disease in comparison to nonsynonymous mutations that do change amino acid residues in proteins remains controversial. Whether the molecular mechanisms described by which synonymous mutations affect organisms can be generalized remains poorly understood and warrants future research in this area.

The distinct translational landscapes of gram-negative Salmonella and gram-positive Listeria

Translational control in pathogenic bacteria is fundamental to gene expression and affects virulence and other infection phenotypes. We used an enhanced ribosome profiling protocol coupled with parallel transcriptomics to capture accurately the global translatome of two evolutionarily distant pathogenic bacteria-the Gram-negative bacterium Salmonella and the Gram-positive bacterium Listeria. We find that the two bacteria use different mechanisms to translationally regulate protein synthesis. In Salmonella, in addition to the expected correlation between translational efficiency and cis-regulatory features such as Shine-Dalgarno (SD) strength and RNA secondary structure around the initiation codon, our data reveal an effect of the 2nd and 3rd codons, where the presence of tandem lysine codons (AAA-AAA) enhances translation in both Salmonella and E. coli. Strikingly, none of these features are seen in efficiently translated Listeria transcripts. Instead, approximately 20% of efficiently translated Listeria genes exhibit 70 S footprints seven nt upstream of the authentic start codon, suggesting that these genes may be subject to a novel translational initiation mechanism. Our results show that SD strength is not a direct hallmark of translational efficiency in all bacteria. Instead, Listeria has evolved additional mechanisms to control gene expression level that are distinct from those utilised by Salmonella and E. coli.

Identification of over ten thousand candidate structured RNAs in viruses and phages

Structured RNAs play crucial roles in viruses, exerting influence over both viral and host gene expression. However, the extensive diversity of structured RNAs and their ability to act in cis or trans positions pose challenges for predicting and assigning their functions. While comparative genomics approaches have successfully predicted candidate structured RNAs in microbes on a large scale, similar efforts for viruses have been lacking. In this study, we screened over 5 million DNA and RNA viral sequences, resulting in the prediction of 10,006 novel candidate structured RNAs. These predictions are widely distributed across taxonomy and ecosystem. We found transcriptional evidence for 206 of these candidate structured RNAs in the human fecal microbiome. These candidate RNAs exhibited evidence of nucleotide covariation, indicative of selective pressure maintaining the predicted secondary structures. Our analysis revealed a diverse repertoire of candidate structured RNAs, encompassing a substantial number of putative tRNAs or tRNA-like structures, Rho-independent transcription terminators, and potentially cis-regulatory structures consistently positioned upstream of genes. In summary, our findings shed light on the extensive diversity of structured RNAs in viruses, offering a valuable resource for further investigations into their functional roles and implications in viral gene expression and pave the way for a deeper understanding of the intricate interplay between viruses and their hosts at the molecular level.

Ribosome heterogeneity results in leader sequence-mediated regulation of protein synthesis in Francisella tularensis

Although ribosomes are generally examined in aggregate, ribosomes can be heterogenous in composition. Evidence is accumulating that changes in ribosome composition may result in altered function, such that ribosome heterogeneity may provide a mechanism to regulate protein synthesis. Ribosome heterogeneity in the human pathogen Francisella tularensis results from incorporation of one of three homologs of bS21, a small ribosomal subunit protein demonstrated to regulate protein synthesis in other bacteria. Loss of one homolog, bS21-2, results in genome-wide post-transcriptional changes in protein abundance. This suggests that bS21-2 can, either directly or indirectly, lead to preferential translation of particular mRNAs. Here, we examine the potential of bS21-2 to function in a leader sequence-dependent manner and to function indirectly, via Hfq. We found that the 5´ untranslated region (UTR) of some bS21-2-responsive genes, including key virulence genes, is sufficient to alter translation in cells lacking bS21-2. We further identify features of a 5´ UTR that allow responsiveness to bS21-2. These include an imperfect Shine-Dalgarno sequence and a particular six nucleotide sequence. Our results are consistent with a model in which a bS21 homolog increases the efficiency of translation initiation through interactions with specific leader sequences. With respect to bS21-2 indirectly regulating translation via the RNA-binding protein Hfq, we found that Hfq controls transcript abundance rather than protein synthesis, impacting virulence gene expression via a distinct mechanism. Together, we determined that ribosome composition in F. tularensis regulates translation in a leader sequence-dependent manner, a regulatory mechanism which may be used in other bacteria. IMPORTANCE Ribosome heterogeneity is common in bacteria, and there is mounting evidence that ribosome composition plays a regulatory role in protein synthesis. However, mechanisms of ribosome-driven gene regulation are not well understood. In the human pathogen Francisella tularensis, which encodes multiple homologs for the ribosomal protein bS21, loss of one homolog impacts protein synthesis and virulence. Here, we explore the mechanism behind bS21-mediated changes in protein synthesis, finding that they can be linked to altered translation initiation and are dependent on specific sequences in the leaders of transcripts. Our data support a model in which ribosome composition regulates gene expression through translation, a strategy that may be conserved in diverse organisms with various sources of ribosome heterogeneity.

Ribosome profiling reveals the translational landscape and allele-specific translational efficiency in rice

Translational regulation is a critical step in the process of gene expression and governs the synthesis of proteins from mRNAs. Many studies have revealed translational regulation in plants in response to various environmental stimuli. However, there have been no studies documenting the comprehensive landscape of translational regulation and allele-specific translational efficiency in multiple plant tissues, especially those of rice, a main staple crop that feeds nearly half of the world's population. Here we used RNA sequencing and ribosome profiling data to analyze the transcriptome and translatome of an elite hybrid rice, Shanyou 63 (SY63), and its parental varieties Zhenshan 97 and Minghui 63. The results revealed that gene expression patterns varied more among tissues than among varieties at the transcriptional and translational levels. We identified 3392 upstream open reading frames (uORFs), and the uORF-containing genes were enriched in transcription factors. Only 668 of 13 492 long non-coding RNAs could be translated into peptides. Finally, we discovered numerous genes with allele-specific translational efficiency in SY63 and demonstrated that some cis-regulatory elements may contribute to allelic divergence in translational efficiency. Overall, these findings may improve our understanding of translational regulation in rice and provide information for molecular breeding research.

Random and Natural Non-Coding RNA Have Similar Structural Motif Patterns but Differ in Bulge, Loop, and Bond Counts

An important question in evolutionary biology is whether (and in what ways) genotype–phenotype (GP) map biases can influence evolutionary trajectories. Untangling the relative roles of natural selection and biases (and other factors) in shaping phenotypes can be difficult. Because the RNA secondary structure (SS) can be analyzed in detail mathematically and computationally, is biologically relevant, and a wealth of bioinformatic data are available, it offers a good model system for studying the role of bias. For quite short RNA (length L≤126), it has recently been shown that natural and random RNA types are structurally very similar, suggesting that bias strongly constrains evolutionary dynamics. Here, we extend these results with emphasis on much larger RNA with lengths up to 3000 nucleotides. By examining both abstract shapes and structural motif frequencies (i.e., the number of helices, bonds, bulges, junctions, and loops), we find that large natural and random structures are also very similar, especially when contrasted to typical structures sampled from the spaces of all possible RNA structures. Our motif frequency study yields another result, where the frequencies of different motifs can be used in machine learning algorithms to classify random and natural RNA with high accuracy, especially for longer RNA (e.g., ROC AUC 0.86 for L = 1000). The most important motifs for classification are the number of bulges, loops, and bonds. This finding may be useful in using SS to detect candidates for functional RNA within ‘junk’ DNA regions.

Identification and in silico characterization of CSRP3 synonymous variants in dilated cardiomyopathy

BACKGROUND

Synonymous variations have always been ignored while studying the underlying genetic mechanisms for most of the human diseases. However, recent studies have suggested that these silent changes in the genome can alter the protein expression and folding.

METHODS AND RESULTS

CSRP3, which is a well-known candidate gene associated with dilated cardiomyopathy (DCM) and hypertrophic cardiomyopathy (HCM), was screened for 100 idiopathic DCM cases and 100 controls. Three synonymous variations were identified viz., c.96G > A, p.K32=; c.336G > A, p.A112=; c.354G > A, p.E118=. A comprehensive in silico analysis was performed using various web based widely accepted tools, Mfold, Codon Usage, HSF3.1 and RNA22. Mfold predicted structural changes in all the variants except c.96 G > A (p.K32=), however it predicted changes in the stability of mRNA due to all the synonymous variants. Codon bias was observed as evident by the Relative Synonymous Codon Usage and Log Ratio of Codon Usage Frequencies. The Human Splicing Finder also predicted remarkable changes in the regulatory elements in the variants c.336G > A and c.354 G > A. The miRNA target prediction using varied modes available in RNA22 revealed that 70.6% of the target sites of miRNAs in CSRP3 were altered due to variant c.336G > A while 29.41% sites were completely lost.

CONCLUSION

Findings of the present study suggest that synonymous variants revealed striking deviations in the structural conformation of mRNA, stability of mRNA, relative synonymous codon usage, splicing and miRNA binding sites from the wild type suggesting their possible role in the pathogenesis of DCM, either by destabilizing the mRNA structure, or codon usage bias or else altering the cis-acting regulatory elements during splicing.

IRES-mediated translation in bacteria

Despite the similarity in fundamental goals of translation initiation between different domains of life, it is one of the most phylogenetically diverse steps of the central dogma of molecular biology. In a classical view, the translation signals for prokaryotes and eukaryotes are distinct from each other. This idea was challenged by the finding that the Internal Ribosome Entry Site (IRES) belonging to Plautia stali intestine virus (PSIV) could bypass the domain-specific boundaries and effectively initiate translation in E. coli. This finding led us to investigate whether the ability of PSIV IRES to initiate translation in E. coli is specific to this IRES and also to study features that allow this viral IRES to mediate prokaryotic translation initiation. We observed that certain IRESs may also possess the ability to initiate E. coli translation. Our results also indicated that the structural integrity of the PSIV IRES in translation in prokaryotes does not appear to be as critical as it is in eukaryotes. We also demonstrated that two regions of the PSIV IRES with complementarity to 16S ribosomal RNA are important for the ability of this IRES to initiate translation in E. coli.

INRI-seq enables global cell-free analysis of translation initiation and off-target effects of antisense inhibitors

Ribosome profiling (Ribo-seq) is a powerful method for the transcriptome-wide assessment of protein synthesis rates and the study of translational control mechanisms. Yet, Ribo-seq also has limitations. These include difficulties with the analysis of translation-modulating molecules such as antibiotics, which are often toxic or challenging to deliver into living cells. Here, we have developed in vitro Ribo-seq (INRI-seq), a cell-free method to analyze the translational landscape of a fully customizable synthetic transcriptome. Using Escherichia coli as an example, we show how INRI-seq can be used to analyze the translation initiation sites of a transcriptome of interest. We also study the global impact of direct translation inhibition by antisense peptide nucleic acid (PNA) to analyze PNA off-target effects. Overall, INRI-seq presents a scalable, sensitive method to study translation initiation in a transcriptome-wide manner without the potentially confounding effects of extracting ribosomes from living cells.

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.

SARS-CoV-2 competes with host mRNAs for efficient translation by maintaining the mutations favorable for translation initiation

During SARS-CoV-2 proliferation, the translation of viral RNAs is usually the rate-limiting step. Understanding the molecular details of this step is beneficial for uncovering the origin and evolution of SARS-CoV-2 and even for controlling the pandemic. To date, it is unclear how SARS-CoV-2 competes with host mRNAs for ribosome binding and efficient translation. We retrieved the coding sequences of all human genes and SARS-CoV-2 genes. We systematically profiled the GC content and folding energy of each CDS. Considering that some fixed or polymorphic mutations exist in SARS-CoV-2 and human genomes, all algorithms and analyses were applied to both pre-mutate and post-mutate versions. In SARS-CoV-2 but not human, the 5-prime end of CDS had lower GC content and less RNA structure than the 3-prime part, which was favorable for ribosome binding and efficient translation initiation. Globally, the fixed and polymorphic mutations in SARS-CoV-2 had created an even lower GC content at the 5-prime end of CDS. In contrast, no similar patterns were observed for the fixed and polymorphic mutations in human genome. Compared with human RNAs, the SARS-CoV-2 RNAs have less RNA structure in the 5-prime end and thus are more favorable of fast translation initiation. The fixed and polymorphic mutations in SARS-CoV-2 are further amplifying this advantage. This might serve as a strategy for SARS-CoV-2 to adapt to the human host.

Stem Region of tRNA Genes Favors Transition Substitution Towards Keto Bases in Bacteria

Transversion and transition mutations have variable effects on the stability of RNA secondary structure considering that the former destabilizes the double helix geometry to a greater extent by introducing purine:purine (R:R) or pyrimidine:pyrimidine (Y:Y) base pairs. Therefore, transversion frequency is likely to be lower than that of transition in the secondary structure regions of RNA genes. Here, we performed an analysis of transition and transversion frequencies in tRNA genes defined well with secondary structure and compared with the intergenic regions in five bacterial species namely Escherichia coli, Klebsiella pneumoniae, Salmonella enterica, Staphylococcus aureus and Streptococcus pneumoniae using a large genome sequence data set. In general, the transversion frequency was observed to be lower than that of transition in both tRNA genes and intergenic regions. The transition to transversion ratio was observed to be greater in tRNA genes than that in the intergenic regions in all the five bacteria that we studied. Interestingly, the intraspecies base substitution analysis in tRNA genes revealed that non-compensatory substitutions were more frequent than compensatory substitutions in the stem region. Further, transition to transversion ratio in the loop region was observed to be significantly lesser than that among the non-compensatory substitutions in the stem region. This indicated that the transversion is more deleterious than transition in the stem regions. In addition, substitutions from amino bases (A/C) to keto bases (G/T) were also observed to be more than the reverse substitutions in the stem region. Substitution from amino bases to keto bases are likely to facilitate the stable G:U pairing unlike the reverse substitution that facilitates the unstable A:C pairing in the stem region of tRNA. This work provides additional support that the secondary structure of tRNA molecule is what drives the different substitutions in its gene sequence.

Evidence for selection on SARS-CoV-2 RNA translation revealed by the evolutionary dynamics of mutations in UTRs and CDSs

RNA translation is the rate-limiting step when cells synthesize proteins. Elevating translation efficiency (TE) is intuitively beneficial. Particularly, when viruses invade host cells, how to compete with endogenous RNAs for efficient translation is a major issue to be resolved. We collected millions of worldwide SARS-CoV-2 sequences during the past year and traced the dynamics of allele frequency of every mutation. We defined adaptive and deleterious mutations according to the rise and fall of their frequencies along time. For 5ʹUTR and synonymous mutations in SARS-CoV-2, the selection on TE is evident near start codons. Adaptive mutations generally decrease GC content while deleterious mutations increase GC content. This trend fades away with increasing distance to start codons. Mutations decreasing GC content near start codons would unravel the complex RNA structure and facilitate translation initiation, which are beneficial to SARS-CoV-2, and vice versa. During this evolutionary arms race between human and virus, SARS-CoV-2 tries to improve its cis elements to compete with host RNAs for rapid translation.

Engineering gene overlaps to sustain genetic constructs in vivo

Evolution is often an obstacle to the engineering of stable biological systems due to the selection of mutations inactivating costly gene circuits. Gene overlaps induce important constraints on sequences and their evolution. We show that these constraints can be harnessed to increase the stability of costly genes by purging loss-of-function mutations. We combine computational and synthetic biology approaches to rationally design an overlapping reading frame expressing an essential gene within an existing gene to protect. Our algorithm succeeded in creating overlapping reading frames in 80% of E. coli genes. Experimentally, scoring mutations in both genes of such overlapping construct, we found that a significant fraction of mutations impacting the gene to protect have a deleterious effect on the essential gene. Such an overlap thus protects a costly gene from removal by natural selection by associating the benefit of this removal with a larger or even lethal cost. In our synthetic constructs, the overlap converts many of the possible mutants into evolutionary dead-ends, reducing the evolutionary potential of the system and thus increasing its stability over time.

Genomes are translated by triplets of nucleotides on two different strands, allowing for six different reading frames. This permits the existence of gene overlaps, often observed in microbial genomes, where two different proteins are encoded on the same piece of DNA, but in different reading frames. Gene overlaps are classically considered an obstacle for both evolution and genetic engineering, as mutations in overlapping regions likely have pleitrotropic effects on several genes. In 2013, we identified specific evolutionary scenarios where the decrease in evolutionary potential caused by gene overlaps could instead be advantageous and selected for. In this work, we demonstrate the use of gene overlaps in another context where reducing evolutionary potential can be useful: preventing evolution from inactivating synthetic circuits. We show that gene overlaps can be engineered to increase the evolutionary stability of genes that are costly to their hosts, by entangling these costly genes with essential genes.

Direct in vivo observation of the effect of codon usage bias on gene expression in Arabidopsis hybrids

Hybrids are the perfect materials to study cis regulatory elements because the two parental alleles are subjected to identical trans environments. There has been a debate on whether synonymous codon usage could affect gene expression. In vitro experiments found that luciferase genes with enhanced codon optimality showed elevated mRNA expression. However, the underlying mechanism is still unclear, and no direct evidence is observed to support this notion. By mapping the RNA-seq data of hybrids of Arabidopsis thaliana and Arabidopsis lyrata, we quantified the allele-specific reads and estimated the relative expression of orthologous genes. We focused on orthologous genes with dN = 0 and dS > 0, which means that they only differ in synonymous codon usage. We found that orthologous genes with higher codon optimality in A. thaliana tend to have higher expression levels of the A. thaliana allele. Codon usage bias could influence gene expression. This phenomenon is not only found in in vitro experiments but also supported by in vivo observations. Therefore, synonymous mutations could have a broad impact from multiple aspects and should not be automatically ignored in genomic studies. KEY MESSAGE: In Arabidopsis hybrids, alleles with higher codon optimality tend to have higher expression levels.

Directing evolution of novel ligands by mRNA display

mRNA display is a powerful biological display platform for the directed evolution of proteins and peptides. mRNA display libraries covalently link the displayed peptide or protein (phenotype) with the encoding genetic information (genotype) through the biochemical activity of the small molecule puromycin. Selection for peptide/protein function is followed by amplification of the linked genetic material and generation of a library enriched in functional sequences. Iterative selection cycles are then performed until the desired level of function is achieved, at which time the identity of candidate peptides can be obtained by sequencing the genetic material. The purpose of this review is to discuss the development of mRNA display technology since its inception in 1997 and to comprehensively review its use in the selection of novel peptides and proteins. We begin with an overview of the biochemical mechanism of mRNA display and its variants with a particular focus on its advantages and disadvantages relative to other biological display technologies. We then discuss the importance of scaffold choice in mRNA display selections and review the results of selection experiments with biological (e.g., fibronectin) and linear peptide library architectures. We then explore recent progress in the development of "drug-like" peptides by mRNA display through the post-translational covalent macrocyclization and incorporation of non-proteogenic functionalities. We conclude with an examination of enabling technologies that increase the speed of selection experiments, enhance the information obtained in post-selection sequence analysis, and facilitate high-throughput characterization of lead compounds. We hope to provide the reader with a comprehensive view of current state and future trajectory of mRNA display and its broad utility as a peptide and protein design tool.

Natural selection on synonymous mutations in SARS-CoV-2 and the impact on estimating divergence time

To adapt to human host environment, synonymous mutations in SARS-CoV-2 are shaped by tRNA selection, energy cost and RNA structure.

Engineering Chloroplasts for High-Level Constitutive or Inducible Transgene Expression

Expression of transgenes from the plastid genome offers a number of attractions to biotechnologists, with the potential to attain very high protein accumulation levels arguably being the most attractive one. High-level transgene expression is of particular importance in resistance engineering (e.g., for expression of insecticidal proteins) and molecular farming (e.g., for expression of pharmaceutical proteins and industrial enzymes). Over the past decades, the production of many commercially valuable proteins in chloroplast-transgenic (transplastomic) plants has been attempted, including pharmaceutical proteins (e.g., subunit vaccines and protein antibiotics) and industrial enzymes. Although in some cases, spectacularly high foreign protein accumulation levels have been obtained, expression levels were disappointingly poor in other cases. In this review, I summarize our current knowledge about the factors influencing the efficiency of plastid transgene expression, and highlight possible optimization strategies to alleviate problems with poor expression levels. I also discuss available techniques for inducible expression of chloroplast transgenes.

A rare codon-based translational program of cell proliferation

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Translational initiation in E. coli occurs at the correct sites genome-wide in the absence of mRNA-rRNA base-pairing

Shine-Dalgarno (SD) motifs are thought to play an important role in translational initiation in bacteria. Paradoxically, ribosome profiling studies in E. coli show no correlation between the strength of an mRNA’s SD motif and how efficiently it is translated. Performing profiling on ribosomes with altered anti-Shine-Dalgarno sequences, we reveal a genome-wide correlation between SD strength and ribosome occupancy that was previously masked by other contributing factors. Using the antibiotic retapamulin to trap initiation complexes at start codons, we find that the mutant ribosomes select start sites correctly, arguing that start sites are hard-wired for initiation through the action of other mRNA features. We show that A-rich sequences upstream of start codons promote initiation. Taken together, our genome-wide study reveals that SD motifs are not necessary for ribosomes to determine where initiation occurs, though they do affect how efficiently initiation occurs.

Mutations that improve efficiency of a weak-link enzyme are rare compared to adaptive mutations elsewhere in the genome

New enzymes often evolve by gene amplification and divergence. Previous experimental studies have followed the evolutionary trajectory of an amplified gene, but have not considered mutations elsewhere in the genome when fitness is limited by an evolving gene. We have evolved a strain of Escherichia coli in which a secondary promiscuous activity has been recruited to serve an essential function. The gene encoding the ‘weak-link’ enzyme amplified in all eight populations, but mutations improving the newly needed activity occurred in only one. Most adaptive mutations occurred elsewhere in the genome. Some mutations increase expression of the enzyme upstream of the weak-link enzyme, pushing material through the dysfunctional metabolic pathway. Others enhance production of a co-substrate for a downstream enzyme, thereby pulling material through the pathway. Most of these latter mutations are detrimental in wild-type E. coli, and thus would require reversion or compensation once a sufficient new activity has evolved.

Synonymous mutations make dramatic contributions to fitness when growth is limited by a weak-link enzyme

Synonymous mutations do not alter the specified amino acid but may alter the structure or function of an mRNA in ways that impact fitness. There are few examples in the literature, however, in which the effects of synonymous mutations on microbial growth rates have been measured, and even fewer for which the underlying mechanism is understood. We evolved four populations of a strain of Salmonella enterica in which a promiscuous enzyme has been recruited to replace an essential enzyme. A previously identified point mutation increases the enzyme’s ability to catalyze the newly needed reaction (required for arginine biosynthesis) but decreases its ability to catalyze its native reaction (required for proline biosynthesis). The poor performance of this enzyme limits growth rate on glucose. After 260 generations, we identified two synonymous mutations in the first six codons of the gene encoding the weak-link enzyme that increase growth rate by 41 and 67%. We introduced all possible synonymous mutations into the first six codons and found substantial effects on growth rate; one doubles growth rate, and another completely abolishes growth. Computational analyses suggest that these mutations affect either the stability of a stem-loop structure that sequesters the start codon or the accessibility of the region between the Shine-Dalgarno sequence and the start codon. Thus, these mutations would be predicted to affect translational efficiency and thereby indirectly affect mRNA stability because translating ribosomes protect mRNA from degradation. Experimental data support these hypotheses. We conclude that the effects of the synonymous mutations are due to a combination of effects on mRNA stability and translation efficiency that alter levels of the weak-link enzyme. These findings suggest that synonymous mutations can have profound effects on fitness under strong selection and that their importance in evolution may be under-appreciated.

When a new enzyme is needed, microbes often recruit a pre-existing enzyme with a promiscuous activity corresponding to the newly needed activity. Such enzymes are often the “weak-link” in metabolism because they have not evolved to efficiently catalyze the new reaction. Under these circumstances, increasing the level of the weak-link enzyme can improve fitness. We evolved a strain of S. enterica in which a weak-link enzyme–E383A ProA–serves essential functions in synthesis of proline and arginine for 260 generations and then sequenced the genomes of several evolved strains. A mutation in the promoter of the operon encoding E383A ProA increased growth rate 9-fold. More surprisingly, a mutation upstream of the start codon and two synonymous mutations within the first six codons also increased growth rate by up to 68%. Introduction of all possible synonymous mutations in the first six codons showed that some doubled growth rate, while others slowed or even prevented growth. Computational and experimental data suggest that these effects were due to enhanced translational efficiency of the weak-link enzyme. These results show that synonymous mutations, once assumed to be selectively neutral, can have strong impacts on fitness when growth rate is limited by a weak-link enzyme.

Prokaryotic coding regions have little if any specific depletion of Shine-Dalgarno motifs

The Shine-Dalgarno motif occurs in front of prokaryotic start codons, and is complementary to the 3’ end of the 16S ribosomal RNA. Hybridization between the Shine-Dalgarno sequence and the anti-Shine-Dalgarno region of the16S rRNA (CCUCCU) directs the ribosome to the start AUG of the mRNA for translation. Shine-Dalgarno-like motifs (AGGAGG in E. coli) are depleted from open reading frames of most prokaryotes. This may be because hybridization of the 16S rRNA at Shine-Dalgarnos inside genes would slow translation or induce internal initiation. However, we analyzed 128 species from diverse phyla where the 16S rRNA gene(s) lack the anti-Shine-Dalgarno sequence, and so the 16S rRNA is incapable of interacting with Shine-Dalgarno-like sequences. Despite this lack of an anti-Shine-Dalgarno, half of these species still displayed depletion of Shine-Dalgarno-like sequences when analyzed by previous methods. Depletion of the same G-rich sequences was seen by these methods even in eukaryotes, which do not use the Shine-Dalgarno mechanism. We suggest previous methods are partly detecting a non-specific depletion of G-rich sequences. Alternative informatics approaches show that most prokaryotes have only slight, if any, specific depletion of Shine-Dalgarno-like sequences from open reading frames. Together with recent evidence that ribosomes do not pause at ORF-internal Shine-Dalgarno motifs, these results suggest the presence of ORF-internal Shine-Dalgarno-like motifs may be inconsequential, perhaps because internal regions of prokaryotic mRNAs may be structurally “shielded” from translation initiation.

Re-annotation of 12,495 prokaryotic 16S rRNA 3' ends and analysis of Shine-Dalgarno and anti-Shine-Dalgarno sequences

We examined 20,648 prokaryotic unique taxids with respect to the annotation of the 3' end of the 16S rRNA, which contains the anti-Shine-Dalgarno sequence. We used the sequence of highly conserved helix 45 of the 16S rRNA as a guide. By this criterion, 8,153 annotated 3' ends correctly included the anti-Shine-Dalgarno sequence, but 12,495 were foreshortened or otherwise mis-annotated, missing part or all of the anti-Shine-Dalgarno sequence, which immediately follows helix 45. We re-annotated, giving a total of 20,648 16S rRNA 3' ends. The vast majority indeed contained a consensus anti-Shine-Dalgarno sequence, embedded in a highly conserved 13 base "tail". However, 128 exceptional organisms had either a variant anti-Shine-Dalgarno, or no recognizable anti-Shine-Dalgarno, in their 16S rRNA(s). For organisms both with and without an anti-Shine-Dalgarno, we identified the Shine-Dalgarno motifs actually enriched in front of each organism's open reading frames. This showed to what extent the Shine-Dalgarno motifs correlated with anti-Shine Dalgarno motifs. In general, organisms whose rRNAs lacked a perfect anti-Shine-Dalgarno motif also lacked a recognizable Shine-Dalgarno. For organisms whose 16S rRNAs contained a perfect anti-Shine-Dalgarno motif, a variety of results were obtained. We found one genus, Alteromonas, where several taxids apparently maintain two different types of 16S rRNA genes, with different, but conserved, antiSDs. The fact that some organisms do not seem to have or use Shine-Dalgarno motifs supports the idea that prokaryotes have other robust mechanisms for recognizing start codons for translation.

Accessibility of the Shine-Dalgarno Sequence Dictates N-Terminal Codon Bias in E. coli

Despite considerable efforts, no physical mechanism has been shown to explain N-terminal codon bias in prokaryotic genomes. Using a systematic study of synonymous substitutions in two endogenous E. coli genes, we show that interactions between the coding region and the upstream Shine-Dalgarno (SD) sequence modulate the efficiency of translation initiation, affecting both intracellular mRNA and protein levels due to the inherent coupling of transcription and translation in E. coli. We further demonstrate that far-downstream mutations can also modulate mRNA levels by occluding the SD sequence through the formation of non-equilibrium secondary structures. By contrast, a non-endogenous RNA polymerase that decouples transcription and translation largely alleviates the effects of synonymous substitutions on mRNA levels. Finally, a complementary statistical analysis of the E. coli genome specifically implicates avoidance of intra-molecular base pairing with the SD sequence. Our results provide general physical insights into the coding-level features that optimize protein expression in prokaryotes.

Trans-splicing enhances translational efficiency in C. elegans

Translational efficiency is subject to extensive regulation. However, the factors influencing such regulation are poorly understood. In Caenorhabditis elegans, 62% of genes are trans-spliced to a specific spliced leader (SL1), which replaces part of the native 5' untranslated region (5' UTR). Given the pivotal role the 5' UTR plays in the regulation of translational efficiency, we hypothesized that SL1 trans-splicing functions to regulate translational efficiency. With genome-wide analysis on Ribo-seq data, polysome profiling experiments, and CRISPR-Cas9-based genetic manipulation of trans-splicing sites, we found four lines of evidence in support of this hypothesis. First, SL1 trans-spliced genes have higher translational efficiencies than non-trans-spliced genes. Second, SL1 trans-spliced genes have higher translational efficiencies than non-trans-spliced orthologous genes in other nematode species. Third, an SL1 trans-spliced isoform has higher translational efficiency than the non-trans-spliced isoform of the same gene. Fourth, deletion of trans-splicing sites of endogenous genes leads to reduced translational efficiency. Importantly, we demonstrated that SL1 trans-splicing plays a key role in enhancing translational efficiencies of essential genes. We further discovered that SL1 trans-splicing likely enhances translational efficiency by shortening the native 5' UTRs, hence reducing the presence of upstream start codons (uAUG) and weakening mRNA secondary structures. Taken together, our study elucidates the global function of trans-splicing in enhancing translational efficiency in nematodes, paving the way for further understanding the genomic mechanisms of translational regulation.

A bicistronic vector with destabilized mRNA secondary structure yields scalable higher titer expression of human neurturin in E. coli

Human neurturin (NTN) is a cystine knot growth factor with potential therapeutic use in diseases such as Parkinson's and diabetes. Scalable high titer production of native NTN is particularly challenging because of the cystine knot structure which consists of an embedded ring comprised of at least three disulfide bonds. We sought to pursue enhanced scalable production of NTN in Escherichia coli. Our initial efforts focused on codon optimization of the first two codons following AUG, but these studies resulted in only a marginal increase in NTN expression. Therefore, we pursued an alternative strategy of using a bicistronic vector for NTN expression designed to reduce mRNA secondary structure to achieve increased ribosome binding and re-initiation. The first cistron was designed to prevent sequestration of the translation initiation region in a secondary conformation. The second cistron, which contained the NTN coding sequence itself, was engineered to disrupt double bonded base pairs and destabilize the secondary structure for ribosome re-initiation. The ensemble approach of reducing NTN's mRNA secondary structure and using the bicistronic vector had an additive effect resulting in significantly increased NTN expression. Our strain selection studies were conducted in a miniaturized bioreactor. An optimized strain was selected and scaled up to a 100 L fermentor, which yielded an inclusion body titer of 2 g/L. The inclusion bodies were refolded to yield active NTN. We believe that our strategy is applicable to other candidate proteins that are difficult-to-express due to stable mRNA secondary structures. Biotechnol. Bioeng. 2017;114: 1753-1761. © 2017 Wiley Periodicals, Inc.

Operon mRNAs are organized into ORF-centric structures that predict translation efficiency

Bacterial mRNAs are organized into operons consisting of discrete open reading frames (ORFs) in a single polycistronic mRNA. Individual ORFs on the mRNA are differentially translated, with rates varying as much as 100-fold. The signals controlling differential translation are poorly understood. Our genome-wide mRNA secondary structure analysis indicated that operonic mRNAs are comprised of ORF-wide units of secondary structure that vary across ORF boundaries such that adjacent ORFs on the same mRNA molecule are structurally distinct. ORF translation rate is strongly correlated with its mRNA structure in vivo, and correlation persists, albeit in a reduced form, with its structure when translation is inhibited and with that of in vitro refolded mRNA. These data suggest that intrinsic ORF mRNA structure encodes a rough blueprint for translation efficiency. This structure is then amplified by translation, in a self-reinforcing loop, to provide the structure that ultimately specifies the translation of each ORF.

DOI:http://dx.doi.org/10.7554/eLife.22037.001

Proteins make up much of the biological machinery inside cells and perform the essential tasks needed to keep each cell alive. Cells contain thousands of different proteins and the instructions needed to build each protein are encoded in genes. However, these instructions cannot be used directly to manufacture the proteins. Instead, a messenger molecule called mRNA is needed to carry the information stored within genes to the parts of the cell where proteins are made.

In bacteria, one mRNA molecule can include information from several genes. This group of genes is called an operon and produces a set of proteins that perform a shared task. Although these proteins work together, some of them are needed in greater numbers than others. Because they are all made using information from the same mRNA, some instructions on the mRNA must be read more times than others. It is unclear how bacterial cells control how many proteins are produced from each part of one mRNA but it is thought to relate to the three-dimensional shape of the molecule itself.

Burkhardt, Rouskin, Zhang et al. have now examined the production of proteins from mRNAs in the commonly studied bacterium, Escherichia coli. The results showed that each set of instructions on the mRNA formed a three-dimensional structure that corresponds to the amount of protein produced from that portion of the mRNA. When this three-dimensional structure is more stable or rigid, the corresponding instructions tended to produce fewer proteins than if the structure was relatively simple and unstable.

Further investigation showed that these three-dimensional mRNA structures could form spontaneously outside of cells, suggesting that molecules other than the mRNA itself have a relatively small role in controlling the number of proteins produced. This also suggests that the entire structure of each mRNA is important and is likely to be essential for cell survival. The next step is to understand why bacteria organise their genes in this way and how the different mRNA structures control how proteins are produced. Moreover, because many bacteria are used like biological factories to produce a variety of commercially useful molecules, these new insights have the potential to enhance a number of manufacturing processes.

DOI:http://dx.doi.org/10.7554/eLife.22037.002

A Codon Deletion at the Beginning of Green Fluorescent Protein Genes Enhances Protein Expression

Recombinant protein expression is one of the key issues in protein engineering and biotechnology. Among the different models for assessing protein production and structure-function studies, green fluorescent protein (GFP) is one of the preferred models because of its importance as a reporter in cellular and molecular studies. In this research we analyze the effect of codon deletions near the amino terminus of different GFP proteins on fluorescence. Our study includes Gly4 deletions in the enhanced GFP (EGFP), the red-shifted GFP and the red-shifted EGFP. The Gly4 deletion mutants and their corresponding wild-type counterparts were transcribed under the control of the T7 or Trc promoters and their expression patterns were analyzed. Different fluorescent outcomes were observed depending on the type of fluorescent gene versions. In silico analysis of the RNA secondary structures near the ribosome binding site revealed a direct relationship between their minimum free energy and GFP production. Integrative analysis of these results, including SDS-PAGE analysis, led us to conclude that the fluorescence improvement of cells expressing different versions of GFPs with Gly4 deleted is due to an enhancement of the accessibility of the ribosome binding site by reducing the stability of the RNA secondary structures at their mRNA leader regions.

Co-expression of the protease furin in Nicotiana benthamiana leads to efficient processing of latent transforming growth factor-β1 into a biologically active protein

Transforming growth factor beta (TGF-β) is a signalling molecule that plays a key role in developmental and immunological processes in mammals. Three TGF-β isoforms exist in humans, and each isoform has unique therapeutic potential. Plants offer a platform for the production of recombinant proteins, which is cheap and easy to scale up and has a low risk of contamination with human pathogens. TGF-β3 has been produced in plants before using a chloroplast expression system. However, this strategy requires chemical refolding to obtain a biologically active protein. In this study, we investigated the possibility to transiently express active human TGF-β1 in Nicotiana benthamiana plants. We successfully expressed mature TGF-β1 in the absence of the latency-associated peptide (LAP) using different strategies, but the obtained proteins were inactive. Upon expression of LAP-TGF-β1, we were able to show that processing of the latent complex by a furin-like protease does not occur in planta. The use of a chitinase signal peptide enhanced the expression and secretion of LAP-TGF-β1, and co-expression of human furin enabled the proteolytic processing of latent TGF-β1. Engineering the plant post-translational machinery by co-expressing human furin also enhanced the accumulation of biologically active TGF-β1. This engineering step is quite remarkable, as furin requires multiple processing steps and correct localization within the secretory pathway to become active. Our data demonstrate that plants can be a suitable platform for the production of complex proteins that rely on specific proteolytic processing.

A comparative genomics study on the effect of individual amino acids on ribosome stalling

BACKGROUND

During protein synthesis, the nascent peptide chain emerges from the ribosome through the ribosomal exit tunnel. Biochemical interactions between the nascent peptide and the tunnel may stall the ribosome movement and thus affect the expression level of the protein being synthesized. Earlier studies focused on one model organism (S. cerevisiae), have suggested that certain amino acid sequences may be responsible for ribosome stalling; however, the stalling effect at the individual amino acid level across many organisms has not yet been quantified.

RESULTS

By analyzing multiple ribosome profiling datasets from different organisms (including prokaryotes and eukaryotes), we report for the first time the organism-specific amino acids that significantly lead to ribosome stalling. We show that the identity of the stalling amino acids vary across the tree of life. In agreement with previous studies, we observed a remarkable stalling signal of proline and arginine in S. cerevisiae. In addition, our analysis supports the conjecture that the stalling effect of positively charged amino acids is not universal and that in certain conditions, negative charge may also induce ribosome stalling. Finally, we show that the beginning part of the tunnel tends to undergo more interactions with the translated amino acids than other positions along the tunnel.

CONCLUSIONS

The reported results support the conjecture that the ribosomal exit tunnel interacts with various amino acids and that the nature of these interactions varies among different organisms. Our findings should contribute towards better understanding of transcript and proteomic evolution and translation elongation regulation.

Enhanced Protein Production in Escherichia coli by Optimization of Cloning Scars at the Vector-Coding Sequence Junction

Protein production in Escherichia coli is a fundamental activity for a large fraction of academic, pharmaceutical, and industrial research laboratories. Maximum production is usually sought, as this reduces costs and facilitates downstream purification steps. Frustratingly, many coding sequences are poorly expressed even when they are codon-optimized and expressed from vectors with powerful genetic elements. In this study, we show that poor expression can be caused by certain nucleotide sequences (e.g., cloning scars) at the junction between the vector and the coding sequence. Since these sequences lie between the Shine-Dalgarno sequence and the start codon, they are an integral part of the translation initiation region. To identify the most optimal sequences, we devised a simple and inexpensive PCR-based step that generates sequence variants at the vector-coding sequence junction. These sequence variants modulated expression by up to 1000-fold. FACS-seq analyses indicated that low GC content and relaxed mRNA stability (ΔG) in this region were important, but not the only, determinants for high expression.

Overexpression of Arabidopsis NADPH-dependent thioredoxin reductase C (AtNTRC) confers freezing and cold shock tolerance to plants

Overexpression of AtNTRC (AtNTRC(OE)) in Arabidopsis thaliana led to a freezing and cold stress tolerance, whereas a knockout mutant (atntrc) showed a stress-sensitive phenotype. Biochemical analyses showed that the recombinant AtNTRC proteins exhibited a cryoprotective activity for malate dehydrogenase and lactic dehydrogenase. Furthermore, conclusive evidence of its interaction with nucleic acids in vitro is provided here on the basis of gel shift and electron microscopy analysis. Recombinant AtNTRC efficiently protected RNA and DNA from RNase A and metal catalyzed oxidation damage, respectively. The C-terminal thioredoxin domain is required for the nucleic acid-protein complex formation. From these results, it can be hypothesized that AtNTRC, which is known to be an electron donor of peroxiredoxin, contributes the stability of macromolecules under cold stress.

Downregulating viral gene expression: codon usage bias manipulation for the generation of novel influenza A virus vaccines

Vaccination represents the best option to protect humans against influenza virus. However, improving the effectiveness of current vaccines could better stifle the health burden caused by viral infection. Protein synthesis from individual genes can be downregulated by synthetically deoptimizing a gene's codon usage. With more rapid and affordable nucleotide synthesis, generating viruses that contain genes with deoptimized codons is now feasible. Attenuated, vaccine-candidate viruses can thus be engineered with hitherto uncharacterized properties. With eight gene segments, influenza A viruses with variably recoded genomes can produce a spectrum of attenuation that is contingent on the gene segment targeted and the number of codon changes. This review summarizes different targets and approaches to deoptimize influenza A virus codons for novel vaccine generation.

Secondary structural entropy in RNA switch (Riboswitch) identification

BACKGROUND

RNA regulatory elements play a significant role in gene regulation. Riboswitches, a widespread group of regulatory RNAs, are vital components of many bacterial genomes. These regulatory elements generally function by forming a ligand-induced alternative fold that controls access to ribosome binding sites or other regulatory sites in RNA. Riboswitch-mediated mechanisms are ubiquitous across bacterial genomes. A typical class of riboswitch has its own unique structural and biological complexity, making de novo riboswitch identification a formidable task. Traditionally, riboswitches have been identified through comparative genomics based on sequence and structural homology. The limitations of structural-homology-based approaches, coupled with the assumption that there is a great diversity of undiscovered riboswitches, suggests the need for alternative methods for riboswitch identification, possibly based on features intrinsic to their structure. As of yet, no such reliable method has been proposed.

RESULTS

We used structural entropy of riboswitch sequences as a measure of their secondary structural dynamics. Entropy values of a diverse set of riboswitches were compared to that of their mutants, their dinucleotide shuffles, and their reverse complement sequences under different stochastic context-free grammar folding models. Significance of our results was evaluated by comparison to other approaches, such as the base-pairing entropy and energy landscapes dynamics. Classifiers based on structural entropy optimized via sequence and structural features were devised as riboswitch identifiers and tested on Bacillus subtilis, Escherichia coli, and Synechococcus elongatus as an exploration of structural entropy based approaches. The unusually long untranslated region of the cotH in Bacillus subtilis, as well as upstream regions of certain genes, such as the sucC genes were associated with significant structural entropy values in genome-wide examinations.

CONCLUSIONS

Various tests show that there is in fact a relationship between higher structural entropy and the potential for the RNA sequence to have alternative structures, within the limitations of our methodology. This relationship, though modest, is consistent across various tests. Understanding the behavior of structural entropy as a fairly new feature for RNA conformational dynamics, however, may require extensive exploratory investigation both across RNA sequences and folding models.

Adaptive evolution of an artificial RNA genome to a reduced ribosome environment

The reconstitution of an artificial system that has the same evolutionary ability as a living thing is a major challenge in the in vitro synthetic biology. In this study, we tested the adaptive evolutionary ability of an artificial RNA genome replication system, termed the translation-coupled RNA replication (TcRR) system. In a previous work, we performed a study of the long-term evolution of the genome with an excess amount of ribosome. In this study, we continued the evolution experiment in a reduced-ribosome environment and observed that the mutant genome compensated for the reduced ribosome concentration. This result demonstrated the ability of the TcRR system to adapt and may be a step toward generating living things with evolutionary ability.

The pCri System: a vector collection for recombinant protein expression and purification

A major bottleneck in structural, biochemical and biophysical studies of proteins is the need for large amounts of pure homogenous material, which is generally obtained by recombinant overexpression. Here we introduce a vector collection, the pCri System, for cytoplasmic and periplasmic/extracellular expression of heterologous proteins that allows the simultaneous assessment of prokaryotic and eukaryotic host cells (Escherichia coli, Bacillus subtilis, and Pichia pastoris). By using a single polymerase chain reaction product, genes of interest can be directionally cloned in all vectors within four different rare restriction sites at the 5'end and multiple cloning sites at the 3'end. In this way, a number of different fusion tags but also signal peptides can be incorporated at the N- and C-terminus of proteins, facilitating their expression, solubility and subsequent detection and purification. Fusion tags can be efficiently removed by treatment with site-specific peptidases, such as tobacco etch virus proteinase, thrombin, or sentrin specific peptidase 1, which leave only a few extra residues at the N-terminus of the protein. The combination of different expression systems in concert with the cloning approach in vectors that can fuse various tags makes the pCri System a valuable tool for high throughput studies.

Influence of the operon structure on poly(3-hydroxypropionate) synthesis in Shimwellia blattae

Glycerol has become a cheap and abundant carbon source due to biodiesel production at a large scale, and it is available for several biotechnological applications. We recently established poly(3-hydroxypropionate) [poly(3HP)] synthesis in a recombinant Shimwellia blattae strain (Heinrich et al. Appl Environ Microbiol 79:3582-3589, 2013). The major drawbacks of the current strains are (i) low poly(3HP) yields, (ii) low plasmid stability and (iii) insufficient conversion rates. In this study, we demonstrated the influence of alterations of the operon structure, consisting of 1,3-propanediol dehydrogenase (dhaT) and aldehyde dehydrogenase (aldD) of Pseudomonas putida KT2442, propionate:coenzyme A (propionate-CoA) transferase (pct) of Clostridium propionicum X2 and polyhydroxyalkanoate (PHA) synthase (phaC1) of Ralstonia eutropha H16. It was shown that S. blattae ATCC33430/pBBR1MCS-2::dhaT::pct::aldD::phaC1 synthesized up to 14.5 % (wtPHA/wtCDW) in a 2-L fed-batch fermentation process. Furthermore, we overcame the problem of plasmid losses during the fermentation period by engineering a carbon source-dependent plasmid addiction system in a triose phosphate isomerase knockout mutant. An assumed poly(3-hydroxyalkanoic acid) degrading activity of the lipase/esterase YbfF could not be confirmed.

Transcript level and sequence determinants of protein abundance and noise in Escherichia coli

The range over which a protein is expressed, and its cell-to-cell variability, is often thought to be linked to the demand for its activity. Steady-state protein level is determined by multiple mechanisms controlling transcription and translation, many of which are limited by DNA- and RNA-encoded signals that affect initiation, elongation and termination of polymerases and ribosomes. We performed a comprehensive analysis of >100 sequence features to derive a predictive model composed of a minimal non-redundant set of factors explaining 66% of the total variation of protein abundance observed in >800 genes in Escherichia coli. The model suggests that protein abundance is primarily determined by the transcript level (53%) and by effectors of translation elongation (12%), whereas only a small fraction of the variation is explained by translational initiation (1%). Our analyses uncover a new sequence determinant, not previously described, affecting translation initiation and suggest that elongation rate is affected by both codon biases and specific amino acid composition. We also show that transcription and translation efficiency may have an effect on expression noise, which is more similar than previously assumed.

D-Tailor: automated analysis and design of DNA sequences

Motivation: Current advances in DNA synthesis, cloning and sequencing technologies afford high-throughput implementation of artificial sequences into living cells. However, flexible computational tools for multi-objective sequence design are lacking, limiting the potential of these technologies.

Results: We developed DNA-Tailor (D-Tailor), a fully extendable software framework, for property-based design of synthetic DNA sequences. D-Tailor permits the seamless integration of multiple sequence analysis tools into a generic Monte Carlo simulation that evolves sequences toward any combination of rationally defined properties. As proof of principle, we show that D-Tailor is capable of designing sequence libraries comprising all possible combinations among three different sequence properties influencing translation efficiency in Escherichia coli. The capacity to design artificial sequences that systematically sample any given parameter space should support the implementation of more rigorous experimental designs.

Availability: Source code is available for download at https://sourceforge.net/projects/dtailor/

Contact:aparkin@lbl.gov or cambray.guillaume@gmail.com

Supplementary information:Supplementary data are available at Bioinformatics online (D-Tailor Tutorial).