Depletion of Shine-Dalgarno Sequences Within Bacterial Coding Regions Is Expression Dependent

CryptKeeper: a negative design tool for reducing unintentional gene expression in bacteria

Foundational techniques in molecular biology-such as cloning genes, tagging biomolecules for purification or identification, and overexpressing recombinant proteins-rely on introducing non-native or synthetic DNA sequences into organisms. These sequences may be recognized by the transcription and translation machinery in their new context in unintended ways. The cryptic gene expression that sometimes results has been shown to produce genetic instability and mask experimental signals. Computational tools have been developed to predict individual types of gene expression elements, but it can be difficult for researchers to contextualize their collective output. Here, we introduce CryptKeeper, a software pipeline that visualizes predictions of bacterial gene expression signals and estimates the translational burden possible from a DNA sequence. We investigate several published examples where cryptic gene expression in E. coli interfered with experiments. CryptKeeper accurately postdicts unwanted gene expression from both eukaryotic virus infectious clones and individual proteins that led to genetic instability. It also identifies off-target gene expression elements that resulted in truncations that confounded protein purification. Incorporating negative design using CryptKeeper into reverse genetics and synthetic biology workflows can help to mitigate cloning challenges and avoid unexplained failures and complications that arise from unintentional gene expression.

Initiator AUGs Are Discriminated from Elongator AUGs Predominantly through mRNA Accessibility in C. crescentus

The initiation of translation in bacteria is thought to occur upon base pairing between the Shine-Dalgarno (SD) site in the mRNA and the anti-SD site in the rRNA. However, in many bacterial species, such as Caulobacter crescentus, a minority of mRNAs have SD sites. To examine the functional importance of SD sites in C. crescentus, we analyzed the transcriptome and found that more SD sites exist in the coding sequence than in the preceding start codons. To examine the function of SD sites in initiation, we designed a series of mutants with altered ribosome accessibility and SD content in translation initiation regions (TIRs) and in elongator AUG regions (EARs). A lack of mRNA structure content is required for initiation in TIRs, and, when introduced into EARs, can stimulate initiation, thereby suggesting that low mRNA structure content is a major feature that is required for initiation. SD sites appear to stimulate initiation in TIRs, which generally lack structure content, but SD sites only stimulate initiation in EARs if RNA secondary structures are destabilized. Taken together, these results suggest that the difference in secondary structure between TIRs and EARs directs ribosomes to start codons where SD base pairing can tune the efficiency of initiation, but SDs in EARs do not stimulate initiation, as they are blocked by stable secondary structures. This highlights the importance of studying translation initiation mechanisms in diverse bacterial species. IMPORTANCE Start codon selection is an essential process that is thought to occur via the base pairing of the rRNA to the SD site in the mRNA. This model is based on studies in E. coli, yet whole-genome sequencing revealed that SD sites are absent at start codons in many species. By examining the transcriptome of C. crescentus, we found more SD-AUG pairs in the CDS of mRNAs than preceding start codons, yet these internal sites do not initiate. Instead, start codon regions have lower mRNA secondary structure content than do internal SD-AUG regions. Therefore, we find that start codon selection is not controlled by the presence of SD sites, which tune initiation efficiency, but by lower mRNA structure content surrounding the start codon.

Prokaryotic rRNA-mRNA interactions are involved in all translation steps and shape bacterial transcripts

The well-established Shine-Dalgarno model suggests that translation initiation in bacteria is regulated via base-pairing between ribosomal RNA (rRNA) and mRNA. We used novel computational analyses and modelling of 823 bacterial genomes coupled with experiments to demonstrate that rRNA-mRNA interactions are diverse and regulate all translation steps from pre-initiation to termination. Previous research has reported the significant influence of rRNA-mRNA interactions, mainly in the initiation phase of translation. The results reported in this paper suggest that, in addition to the rRNA-mRNA interactions near the start codon that trigger initiation in bacteria, rRNA-mRNA interactions affect all sub-stages of the translation process (pre-initiation, initiation, elongation, termination). As these interactions dictate translation efficiency, they serve as an evolutionary driving force for shaping transcripts in bacteria while considering trade-offs between the effects of different interactions across different transcript regions on translation efficacy and efficiency. We observed selection for strong interactions in regions where such interactions are likely to enhance initiation, regulate early elongation, and ensure translation termination fidelity. We discovered selection against strong interactions and for intermediate interactions in coding regions and presented evidence that these patterns maximize elongation efficiency while also enhancing initiation. These finding are relevant to all biomedical disciplines due to the centrality of the translation process and the effect of rRNA-mRNA interactions on transcript evolution.

Compensating Complete Loss of Signal Recognition Particle During Co-translational Protein Targeting by the Translation Speed and Accuracy

Signal recognition particle (SRP) is critical for delivering co-translational proteins to the bacterial inner membrane. Previously, we identified SRP suppressors in Escherichia coli that inhibit translation initiation and elongation, which provided insights into the mechanism of bypassing the requirement of SRP. Suppressor mutations tended to be located in regions that govern protein translation under evolutionary pressure. To test this hypothesis, we re-executed the suppressor screening of SRP. Here, we isolated a novel SRP suppressor mutation located in the Shine–Dalgarno sequence of the S10 operon, which partially offset the targeting defects of SRP-dependent proteins. We found that the suppressor mutation decreased the protein translation rate, which extended the time window of protein targeting. This increased the possibility of the correct localization of inner membrane proteins. Furthermore, the fidelity of translation was decreased in suppressor cells, suggesting that the quality control of translation was inactivated to provide an advantage in tolerating toxicity caused by the loss of SRP. Our results demonstrated that the inefficient protein targeting due to SRP deletion can be rescued through modulating translational speed and accuracy.

A novel framework for evaluating the performance of codon usage bias metrics

The unequal utilization of synonymous codons affects numerous cellular processes including translation rates, protein folding and mRNA degradation. In order to understand the biological impact of variable codon usage bias (CUB) between genes and genomes, it is crucial to be able to accurately measure CUB for a given sequence. A large number of metrics have been developed for this purpose, but there is currently no way of systematically testing the accuracy of individual metrics or knowing whether metrics provide consistent results. This lack of standardization can result in false-positive and false-negative findings if underpowered or inaccurate metrics are applied as tools for discovery. Here, we show that the choice of CUB metric impacts both the significance and measured effect sizes in numerous empirical datasets, raising questions about the generality of findings in published research. To bring about standardization, we developed a novel method to create synthetic protein-coding DNA sequences according to different models of codon usage. We use these benchmark sequences to identify the most accurate and robust metrics with regard to sequence length, GC content and amino acid heterogeneity. Finally, we show how our benchmark can aid the development of new metrics by providing feedback on its performance compared to the state of the art.

Within-Gene Shine-Dalgarno Sequences Are Not Selected for Function

The Shine-Dalgarno (SD) sequence motif facilitates translation initiation and is frequently found upstream of bacterial start codons. However, thousands of instances of this motif occur throughout the middle of protein coding genes in a typical bacterial genome. Here, we use comparative evolutionary analysis to test whether SD sequences located within genes are functionally constrained. We measure the conservation of SD sequences across Enterobacteriales, and find that they are significantly less conserved than expected. Further, the strongest SD sequences are the least conserved whereas we find evidence of conservation for the weakest possible SD sequences given amino acid constraints. Our findings indicate that most SD sequences within genes are likely to be deleterious and removed via selection. To illustrate the origin of these deleterious costs, we show that ATG start codons are significantly depleted downstream of SD sequences within genes, highlighting the constraint that these sequences impose on the surrounding nucleotides to minimize the potential for erroneous translation initiation.

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.

Adenine Enrichment at the Fourth CDS Residue in Bacterial Genes Is Consistent with Error Proofing for +1 Frameshifts

Beyond selection for optimal protein functioning, coding sequences (CDSs) are under selection at the RNA and DNA levels. Here, we identify a possible signature of “dual-coding,” namely extensive adenine (A) enrichment at bacterial CDS fourth sites. In 99.07% of studied bacterial genomes, fourth site A use is greater than expected given genomic A-starting codon use. Arguing for nucleotide level selection, A-starting serine and arginine second codons are heavily utilized when compared with their non-A starting synonyms. Several models have the ability to explain some of this trend. In part, A-enrichment likely reduces 5′ mRNA stability, promoting translation initiation. However T/U, which may also reduce stability, is avoided. Further, +1 frameshifts on the initiating ATG encode a stop codon (TGA) provided A is the fourth residue, acting either as a frameshift “catch and destroy” or a frameshift stop and adjust mechanism and hence implicated in translation initiation. Consistent with both, genomes lacking TGA stop codons exhibit weaker fourth site A-enrichment. Sequences lacking a Shine–Dalgarno sequence and those without upstream leader genes, that may be more error prone during initiation, have greater utilization of A, again suggesting a role in initiation. The frameshift correction model is consistent with the notion that many genomic features are error-mitigation factors and provides the first evidence for site-specific out of frame stop codon selection. We conjecture that the NTG universal start codon may have evolved as a consequence of TGA being a stop codon and the ability of NTGA to rapidly terminate or adjust a ribosome.

Diversity of Translation Initiation Mechanisms across Bacterial Species Is Driven by Environmental Conditions and Growth Demands

The Shine-Dalgarno (SD) sequence motif is frequently found upstream of protein coding genes and is thought to be the dominant mechanism of translation initiation used by bacteria. Experimental studies have shown that the SD sequence facilitates start codon recognition and enhances translation initiation by directly interacting with the highly conserved anti-SD sequence on the 30S ribosomal subunit. However, the proportion of SD-led genes within a genome varies across species and the factors governing this variation in translation initiation mechanisms remain largely unknown. Here, we conduct a phylogenetically informed analysis and find that species capable of rapid growth contain a higher proportion of SD-led genes throughout their genomes. We show that SD sequence utilization covaries with a suite of genomic features that are important for efficient translation initiation and elongation. In addition to these endogenous genomic factors, we further show that exogenous environmental factors may influence the evolution of translation initiation mechanisms by finding that thermophilic species contain significantly more SD-led genes than mesophiles. Our results demonstrate that variation in translation initiation mechanisms across bacterial species is predictable and is a consequence of differential life-history strategies related to maximum growth rate and environmental-specific constraints.

Orthogonal Ribosome Biofirewall

Biocontainment systems are crucial for preventing genetically modified organisms from escaping into natural ecosystems. Here, we describe the orthogonal ribosome biofirewall, which consists of an activation circuit and a degradation circuit. The activation circuit is a genetic AND gate based on activation of the encrypted pathway by the orthogonal ribosome in response to specific environmental signals. The degradation circuit is a genetic NOT gate with an output of I-SceI homing endonuclease, which conditionally degrades the orthogonal ribosome genes. We demonstrate that the activation circuit can be flexibly incorporated into genetic circuits and metabolic pathways for encryption. The plasmid-based encryption of the deoxychromoviridans pathway and the genome-based encryption of lacZ are tightly regulated and can decrease the expression to 7.3% and 7.8%, respectively. We validated the ability of the degradation circuit to decrease the expression levels of the target plasmids and the orthogonal rRNA (O-rRNA) plasmids to 0.8% in lab medium and 0.76% in nonsterile soil medium, respectively. Our orthogonal ribosome biofirewall is a versatile platform that can be useful in biosafety research and in the biotechnology industry.

Comparative genomic analysis of translation initiation mechanisms for genes lacking the Shine-Dalgarno sequence in prokaryotes

In prokaryotes, translation initiation is believed to occur through an interaction between the 3΄ tail of a 16S rRNA and a corresponding Shine–Dalgarno (SD) sequence in the 5΄ untranslated region (UTR) of an mRNA. However, some genes lack SD sequences (non-SD genes), and the fraction of non-SD genes in a genome varies depending on the prokaryotic species. To elucidate non-SD translation initiation mechanisms in prokaryotes from an evolutionary perspective, we statistically examined the nucleotide frequencies around the initiation codons in non-SD genes from 260 prokaryotes (235 bacteria and 25 archaea). We identified distinct nucleotide frequency biases upstream of the initiation codon in bacteria and archaea, likely because of the presence of leaderless mRNAs lacking a 5΄ UTR. Moreover, we observed overall similarities in the nucleotide patterns between upstream and downstream regions of the initiation codon in all examined phyla. Symmetric nucleotide frequency biases might facilitate translation initiation by preventing the formation of secondary structures around the initiation codon. These features are more prominent in species’ genomes that harbor large fractions of non-SD sequences, suggesting that a reduced stability around the initiation codon is important for efficient translation initiation in prokaryotes.

Synonymous codon bias as a basis for novel antibiotic design: from nucleotide wobble constraint to ribosomal garrotte

AIM

The observation of low guanine frequencies at wobble position codons over Rossmann-fold GXGXXG motifs, led to the discovery of sequence-wide synonymous codon bias, from this we propose novel ribosomal inhibitors. Methodology & results: The wobble bases of multiple sequence alignments of diverse Rossmann-fold enzymes were counted, A, C, T and G wobble frequencies of consecutive codons were displayed as wobble plots. Synonymous codon constraints were found throughout the length of proteins, particularly at Rossmann folds. Prokaryote mRNAs have disallowed nucleotide sequences with wobble bases, evolutionary tailored to safeguard translation. We propose one such editing mechanism ensures that glycine motifs avoid becoming Shine-Dalgarno motifs.

CONCLUSION

We propose using peptide nucleic acids and cyclic peptidyl scaffolds as a means of presenting disallowed chemical interactions to the anti-Shine-Dalgarno.