Comprehensive analysis of stop codon usage in bacteria and its correlation with release factor abundance

Fluorescent reporters give new insights into antibiotics-induced nonsense and frameshift mistranslation

We developed a reporter system based on simultaneous expression of two fluorescent proteins: GFP as a reporter of the capacity of protein synthesis and mutated mScarlet-I as a reporter of translational errors. Because of the unique stop codons or frameshift mutations introduced into the mScarlet-I gene, red fluorescence was produced only after a mistranslation event. These reporters allowed us to estimate mistranslation at a single cell level using either flow cytometry or fluorescence microscopy. We found that laboratory strains of Escherichia coli are more prone to mistranslation compared to the clinical isolates. As relevant for uropathogenic E. coli, growth in human urine elevated translational frameshifting compared to standard laboratory media, whereas different standard media had a small effect on translational fidelity. Antibiotic-induced mistranslation was studied by using amikacin (aminoglycoside family) and azithromycin (macrolide family). Bactericidal amikacin induced preferably stop-codon readthrough at a moderate level. Bacteriostatic azithromycin on the other hand induced both frameshifting and stop-codon readthrough at much higher level. Single cell analysis revealed that fluorescent reporter-protein signal can be lost due to leakage from a fraction of bacteria in the presence of antibiotics, demonstrating the complexity of the antimicrobial activity.

Transient disome complex formation in native polysomes during ongoing protein synthesis captured by cryo-EM

Structural studies of translating ribosomes traditionally rely on in vitro assembly and stalling of ribosomes in defined states. To comprehensively visualize bacterial translation, we reactivated ex vivo-derived E. coli polysomes in the PURE in vitro translation system and analyzed the actively elongating polysomes by cryo-EM. We find that 31% of 70S ribosomes assemble into disome complexes that represent eight distinct functional states including decoding and termination intermediates, and a pre-nucleophilic attack state. The functional diversity of disome complexes together with RNase digest experiments suggests that paused disome complexes transiently form during ongoing elongation. Structural analysis revealed five disome interfaces between leading and queueing ribosomes that undergo rearrangements as the leading ribosome traverses through the elongation cycle. Our findings reveal at the molecular level how bL9's CTD obstructs the factor binding site of queueing ribosomes to thwart harmful collisions and illustrate how translation dynamics reshape inter-ribosomal contacts.

Analysis of programmed frameshifting during translation of prfB in Flavobacterium johnsoniae

Ribosomes of Bacteroidia fail to recognize Shine-Dalgarno (SD) sequences due to sequestration of the 3' tail of the 16S rRNA on the 30S platform. Yet in these organisms, the prfB gene typically contains the programmed +1 frameshift site with its characteristic SD sequence. Here, we investigate prfB autoregulation in Flavobacterium johnsoniae, a member of the Bacteroidia. We find that the efficiency of prfB frameshifting in F. johnsoniae is low (∼7%) relative to that in Escherichia coli (∼50%). Mutation or truncation of bS21 in F. johnsoniae increases frameshifting substantially, suggesting that anti-SD (ASD) sequestration is responsible for the reduced efficiency. The frameshift site of certain Flavobacteriales, such as Winogradskyella psychrotolerans, has no SD. In F. johnsoniae, this W. psychrotolerans sequence supports frameshifting as well as the native sequence, and mutation of bS21 causes no enhancement. These data suggest that prfB frameshifting normally occurs without SD-ASD pairing, at least under optimal laboratory growth conditions. Chromosomal mutations that remove the frameshift or ablate the SD confer subtle growth defects in the presence of paraquat or streptomycin, respectively, indicating that both the autoregulatory mechanism and the SD element contribute to F. johnsoniae cell fitness. Analysis of prfB frameshift sites across 2686 representative bacteria shows loss of the SD sequence in many clades, with no obvious relationship to genome-wide SD usage. These data reveal unexpected variation in the mechanism of frameshifting and identify another group of organisms, the Verrucomicrobiales, that globally lack SD sequences.

Host translation machinery is not a barrier to phages that interact with both CPR and non-CPR bacteria

IMPORTANCE

Here, we profiled putative phages of Saccharibacteria, which are of particular importance as Saccharibacteria influence some human oral diseases. We additionally profiled putative phages of Gracilibacteria and Absconditabacteria, two Candidate Phyla Radiation (CPR) lineages of interest given their use of an alternative genetic code. Among the phages identified in this study, some are targeted by spacers from both CPR and non-CPR bacteria and others by both bacteria that use the standard genetic code as well as bacteria that use an alternative genetic code. These findings represent new insights into possible phage replication strategies and have relevance for phage therapies that seek to manipulate microbiomes containing CPR bacteria.

Genes for highly abundant proteins in Escherichia coli avoid 5' codons that promote ribosomal initiation

In many species highly expressed genes (HEGs) over-employ the synonymous codons that match the more abundant iso-acceptor tRNAs. Bacterial transgene codon randomization experiments report, however, that enrichment with such "translationally optimal" codons has little to no effect on the resultant protein level. By contrast, consistent with the view that ribosomal initiation is rate limiting, synonymous codon usage following the 5' ATG greatly influences protein levels, at least in part by modifying RNA stability. For the design of bacterial transgenes, for simple codon based in silico inference of protein levels and for understanding selection on synonymous mutations, it would be valuable to computationally determine initiation optimality (IO) scores for codons for any given species. One attractive approach is to characterize the 5' codon enrichment of HEGs compared with the most lowly expressed genes, just as translational optimality scores of codons have been similarly defined employing the full gene body. Here we determine the viability of this approach employing a unique opportunity: for Escherichia coli there is both the most extensive protein abundance data for native genes and a unique large-scale transgene codon randomization experiment enabling objective definition of the 5' codons that cause, rather than just correlate with, high protein abundance (that we equate with initiation optimality, broadly defined). Surprisingly, the 5' ends of native genes that specify highly abundant proteins avoid such initiation optimal codons. We find that this is probably owing to conflicting selection pressures particular to native HEGs, including selection favouring low initiation rates, this potentially enabling high efficiency of ribosomal usage and low noise. While the classical HEG enrichment approach does not work, rendering simple prediction of native protein abundance from 5' codon content futile, we report evidence that initiation optimality scores derived from the transgene experiment may hold relevance for in silico transgene design for a broad spectrum of bacteria.

Evolution of termination codons of proteins and the TAG-TGA paradox

In most eukaryotes and prokaryotes TGA is used at a significantly higher frequency than TAG as termination codon of protein-coding genes. Although this phenomenon has been recognized several years ago, there is no generally accepted explanation for the TAG-TGA paradox. Our analyses of human mutation data revealed that out of the eighteen sense codons that can give rise to a nonsense codon by single base substitution, the CGA codon is exceptional: it gives rise to the TGA stop codon at an order of magnitude higher rate than the other codons. Here we propose that the TAG-TGA paradox is due to methylation and hypermutabilty of CpG dinucleotides. In harmony with this explanation, we show that the coding genomes of organisms with strong CpG methylation have a significant bias for TGA whereas those from organisms that lack CpG methylation use TGA and TAG termination codons with similar probability.

Epistasis between promoter activity and coding mutations shapes gene evolvability

The evolution of protein-coding genes proceeds as mutations act on two main dimensions: regulation of transcription level and the coding sequence. The extent and impact of the connection between these two dimensions are largely unknown because they have generally been studied independently. By measuring the fitness effects of all possible mutations on a protein complex at various levels of promoter activity, we show that promoter activity at the optimal level for the wild-type protein masks the effects of both deleterious and beneficial coding mutations. Mutations that are deleterious at low activity but masked at optimal activity are slightly destabilizing for individual subunits and binding interfaces. Coding mutations that increase protein abundance are beneficial at low expression but could potentially incur a cost at high promoter activity. We thereby demonstrate that promoter activity in interaction with protein properties can dictate which coding mutations are beneficial, neutral, or deleterious.

The Influence of the Nucleotide Composition of Genes and Gene Regulatory Elements on the Efficiency of Protein Expression in Escherichia coli

Recombinant proteins expressed in Escherichia coli are widely used in biochemical research and industrial processes. At the same time, achieving higher protein expression levels and correct protein folding still remains the key problem, since optimization of nutrient media, growth conditions, and methods for induction of protein synthesis do not always lead to the desired result. Often, low protein expression is determined by the sequences of the expressed genes and their regulatory regions. The genetic code is degenerated; 18 out of 20 amino acids are encoded by more than one codon. Choosing between synonymous codons in the coding sequence can significantly affect the level of protein expression and protein folding due to the influence of the gene nucleotide composition on the probability of formation of secondary mRNA structures that affect the ribosome binding at the translation initiation phase, as well as the ribosome movement along the mRNA during elongation, which, in turn, influences the mRNA degradation and the folding of the nascent protein. The nucleotide composition of the mRNA untranslated regions, in particular the promoter and Shine-Dalgarno sequences, also affects the efficiency of mRNA transcription, translation, and degradation. In this review, we describe the genetic principles that determine the efficiency of protein production in Escherichia coli.

Bioprocessing of recombinant proteins from Escherichia coli inclusion bodies: insights from structure-function relationship for novel applications

The formation of inclusion bodies (IBs) during expression of recombinant therapeutic proteins using E. coli is a significant hurdle in producing high-quality, safe, and efficacious medicines. The improved understanding of the structure-function relationship of the IBs has resulted in the development of novel biotechnologies that have streamlined the isolation, solubilization, refolding, and purification of the active functional proteins from the bacterial IBs. Together, this overall effort promises to radically improve the scope of experimental biology of therapeutic protein production and expand new prospects in IBs usage. Notably, the IBs are increasingly used for applications in more pristine areas such as drug delivery and material sciences. In this review, we intend to provide a comprehensive picture of the bio-processing of bacterial IBs, including assessing critical gaps that still need to be addressed and potential solutions to overcome them. We expect this review to be a useful resource for those working in the area of protein refolding and therapeutic protein production.

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.

Complete Mitochondrial Genome Sequences of Korean Phytophthora infestans Isolates and Comparative Analysis of Mitochondrial Haplotypes

Potato late blight caused by Phytophthora infestans is a destructive disease in Korea. To elucidate the genomic variation of the mitochondrial (mt) genome, we assembled its complete mt genome and compared its sequence among different haplotypes. The mt genome sequences of four Korean P. infestans isolates were revealed by Illumina HiSeq. The size of the circular mt genome of the four major genotypes, KR_1_A1, KR_2_A2, SIB-1, and US-11, was 39,872, 39,836, 39,872, and 39,840 bp, respectively. All genotypes contained the same 61 genes in the same order, comprising two RNA-encoding genes, 16 ribosomal genes, 25 transfer RNA, 17 genes encoding electron transport and ATP synthesis, 11 open reading frames of unknown function, and one protein import-related gene, tatC. The coding region comprised 91% of the genome, and GC content was 22.3%. The haplotypes were further analyzed based on sequence polymorphism at two hypervariable regions (HVRi), carrying a 2 kb insertion/deletion sequence, and HVRii, carrying 36 bp variable number tandem repeats (VNTRs). All four genotypes carried the 2 kb insertion/deletion sequence in HVRi, whereas HVRii had two VNTRs in KR_1_A1 and SIB-1 but three VNTRs in US-11 and KR_2_A2. Minimal spanning network and phylogenetic analysis based on 5,814 bp of mtDNA sequences from five loci, KR_1_A1 and SIB-1 were classified as IIa-6 haplotype, and isolates KR_1_A2 and US-11 as haplotypes IIa-5 and IIb-2, respectively. mtDNA sequences of KR_1_A1 and SIB-1 shared 100% sequence identity, and both were 99.9% similar to those of KR_2_A2 and US-11.

Stop codon usage as a window into genome evolution: mutation, selection, biased gene conversion and the TAG paradox

Protein coding genes terminate with one of three stop codons (TAA, TGA, or TAG) that, like synonymous codons, are not employed equally. With TGA and TAG having identical nucleotide content, analysis of their differential usage provides an unusual window into the forces operating on what are ostensibly functionally identical residues. Across genomes and between isochores within the human genome, TGA usage increases with G + C content but, with a common G + C → A + T mutation bias, this cannot be explained by mutation bias-drift equilibrium. Increased usage of TGA in G + C-rich genomes or genomic regions is also unlikely to reflect selection for the optimal stop codon, as TAA appears to be universally optimal, probably because it has the lowest read-through rate. Despite TAA being favored by selection and mutation bias, as with codon usage bias G + C pressure is the prime determinant of between-species TGA usage trends. In species with strong G + C-biased gene conversion (gBGC), such as mammals and birds, the high usage and conservation of TGA is best explained by an A + T → G + C repair bias. How to explain TGA enrichment in other G + C-rich genomes is less clear. Enigmatically, across bacterial and archaeal species and between human isochores TAG usage is mostly unresponsive to G + C pressure. This unresponsiveness we dub the TAG paradox as currently no mutational, selective, or gBGC model provides a well-supported explanation. That TAG does increase with G + C usage across eukaryotes makes the usage elsewhere yet more enigmatic. We suggest resolution of the TAG paradox may provide insights into either an unknown but common selective preference (probably at the DNA/RNA level) or an unrecognized complexity to the action of gBGC.

Unusual mammalian usage of TGA stop codons reveals that sequence conservation need not imply purifying selection

The assumption that conservation of sequence implies the action of purifying selection is central to diverse methodologies to infer functional importance. GC-biased gene conversion (gBGC), a meiotic mismatch repair bias strongly favouring GC over AT, can in principle mimic the action of selection, this being thought to be especially important in mammals. As mutation is GC→AT biased, to demonstrate that gBGC does indeed cause false signals requires evidence that an AT-rich residue is selectively optimal compared to its more GC-rich allele, while showing also that the GC-rich alternative is conserved. We propose that mammalian stop codon evolution provides a robust test case. Although in most taxa TAA is the optimal stop codon, TGA is both abundant and conserved in mammalian genomes. We show that this mammalian exceptionalism is well explained by gBGC mimicking purifying selection and that TAA is the selectively optimal codon. Supportive of gBGC, we observe (i) TGA usage trends are consistent at the focal stop codon and elsewhere (in UTR sequences); (ii) that higher TGA usage and higher TAA→TGA substitution rates are predicted by a high recombination rate; and (iii) across species the difference in TAA <-> TGA substitution rates between GC-rich and GC-poor genes is largest in genomes that possess higher between-gene GC variation. TAA optimality is supported both by enrichment in highly expressed genes and trends associated with effective population size. High TGA usage and high TAA→TGA rates in mammals are thus consistent with gBGC’s predicted ability to “drive” deleterious mutations and supports the hypothesis that sequence conservation need not be indicative of purifying selection. A general trend for GC-rich trinucleotides to reside at frequencies far above their mutational equilibrium in high recombining domains supports the generality of these results.

Genetic Determinants of Tigecycline Resistance in Mycobacteroides abscessus

Mycobacteroides abscessus (formerly Mycobacterium abscessus) is a clinically important, rapid-growing non-tuberculous mycobacterium notoriously known for its multidrug-resistance phenotype. The intrinsic resistance of M. abscessus towards first- and second-generation tetracyclines is mainly due to the over-expression of a tetracycline-degrading enzyme known as MabTetX (MAB_1496c). Tigecycline, a third-generation tetracycline, is a poor substrate for the MabTetX and does not induce the expression of this enzyme. Although tigecycline-resistant strains of M. abscessus have been documented in different parts of the world, their resistance determinants remain largely elusive. Recent work on tigecycline resistance or reduced susceptibility in M. abscessus revealed the involvement of the gene MAB_3508c which encodes the transcriptional activator WhiB7, as well as mutations in the sigH-rshA genes which control heat shock and oxidative-stress responses. The deletion of whiB7 has been observed to cause a 4-fold decrease in the minimum inhibitory concentration of tigecycline. In the absence of environmental stress, the SigH sigma factor (MAB_3543c) interacts with and is inhibited by the anti-sigma factor RshA (MAB_3542c). The disruption of the SigH-RshA interaction resulting from mutations and the subsequent up-regulation of SigH have been hypothesized to lead to tigecycline resistance in M. abscessus. In this review, the evidence for different genetic determinants reported to be linked to tigecycline resistance in M. abscessus was examined and discussed.

Why Is the UAG (Amber) Stop Codon Almost Absent in Highly Expressed Bacterial Genes?

The genome hypothesis postulates that genes in a genome tend to conform to their species’ usage of the codon catalog and the GC content of the DNA. Thus, codon frequencies differ across organisms, including the three termination codons in the standard genetic code. Here, we analyze the frequencies of stop codons in a group of highly expressed genes from 196 prokaryotes under strong translational selection. The occurrence of the three translation termination codons is highly biased, with UAA (ochre) being the most prevalent in almost all bacteria. In contrast, UAG (amber) is the least frequent termination codon, e.g., only 321 occurrences (7.4%) in E. coli K-12 substr. W3110. Of the 253 highly expressed genes, only two end with an UAG codon. The strength of the selective bias against UAG in highly expressed genes varies among bacterial genomes, but it is not affected by the GC content of these genomes. In contrast, increased GC content results in a decrease in UAA abundance with a concomitant increase in UGA abundance. We propose that readthrough efficiency and context effects could explain the prevalence of UAA over UAG, particularly in highly expressed genes. Findings from this communication can be utilized for the optimization of gene expression.

Variation in Release Factor Abundance Is Not Needed to Explain Trends in Bacterial Stop Codon Usage

In bacteria stop codons are recognized by one of two class I release factors (RF1) recognizing TAG, RF2 recognizing TGA, and TAA being recognized by both. Variation across bacteria in the relative abundance of RF1 and RF2 is thus hypothesized to select for different TGA/TAG usage. This has been supported by correlations between TAG:TGA ratios and RF1:RF2 ratios across multiple bacterial species, potentially also explaining why TAG usage is approximately constant despite extensive variation in GC content. It is, however, possible that stop codon trends are determined by other forces and that RF ratios adapt to stop codon usage, rather than vice versa. Here, we determine which direction of the causal arrow is the more parsimonious. Our results support the notion that RF1/RF2 ratios become adapted to stop codon usage as the same trends, notably the anomalous TAG behavior, are seen in contexts where RF1:RF2 ratios cannot be, or are unlikely to be, causative, that is, at 3'untranslated sites never used for translation termination, in intragenomic analyses, and across archaeal species (that possess only one RF1). We conclude that specifics of RF biology are unlikely to fully explain TGA/TAG relative usage. We discuss why the causal relationships for the evolution of synonymous stop codon usage might be different from those affecting synonymous sense codon usage, noting that transitions between TGA and TAG require two-point mutations one of which is likely to be deleterious.

In-depth analysis of amino acid and nucleotide sequences of Hsp60: how conserved is this protein?

Chaperonin Hsp60, as a protein found in all organisms, is of great interest in medicine, since it is present in many tissues and can be used both as a drug and as an object of targeted therapy. Hence, Hsp60 deserves a fundamental comparative analysis to assess its evolutionary characteristics. It was found that the percent identity of Hsp60 amino acid sequences both within and between phyla was not high enough to identify Hsp60s as highly conserved proteins. However, their ATP binding sites are largely conserved. The amino acid composition of Hsp60s remained relatively constant. At the same time, the analysis of the nucleotide sequences showed that GC content in the Hsp60 genes was comparable to or greater than the genomic values, which may indicate a high resistance to mutations due to tight control of the nucleotide composition by DNA repair systems. Natural selection plays a dominant role in the evolution of Hsp60 genes. The degree of mutational pressure affecting the Hsp60 genes is quite low, and its direction does not depend on taxonomy. Interestingly, for the Hsp60 genes from Chordata, Arthropoda, and Proteobacteria the exact direction of mutational pressure could not be determined. However, upon further division into classes, it was found that the direction of the mutational pressure for Hsp60 genes from Fish differs from that for other chordates. The direction of the mutational pressure affects the synonymous codon usage bias. The number of high and low represented codons increases with increasing GC content, which can improve codon usage. Special server has been created for bioinformatics analysis of Hsp60: http://oka.protres.ru:4202/.

Intergenic ORFs as elementary structural modules of de novo gene birth and protein evolution

The noncoding genome plays an important role in de novo gene birth and in the emergence of genetic novelty. Nevertheless, how noncoding sequences' properties could promote the birth of novel genes and shape the evolution and the structural diversity of proteins remains unclear. Therefore, by combining different bioinformatic approaches, we characterized the fold potential diversity of the amino acid sequences encoded by all intergenic open reading frames (ORFs) of S. cerevisiae with the aim of (1) exploring whether the structural states' diversity of proteomes is already present in noncoding sequences, and (2) estimating the potential of the noncoding genome to produce novel protein bricks that could either give rise to novel genes or be integrated into pre-existing proteins, thus participating in protein structure diversity and evolution. We showed that amino acid sequences encoded by most yeast intergenic ORFs contain the elementary building blocks of protein structures. Moreover, they encompass the large structural state diversity of canonical proteins, with the majority predicted as foldable. Then, we investigated the early stages of de novo gene birth by reconstructing the ancestral sequences of 70 yeast de novo genes and characterized the sequence and structural properties of intergenic ORFs with a strong translation signal. This enabled us to highlight sequence and structural factors determining de novo gene emergence. Finally, we showed a strong correlation between the fold potential of de novo proteins and one of their ancestral amino acid sequences, reflecting the relationship between the noncoding genome and the protein structure universe.

ReScribe: An Unrestrained Tool Combining Multiplex Recombineering and Minimal-PAM ScCas9 for Genome Recoding Pseudomonas putida

Genome recoding enables incorporating new functions into the DNA of microorganisms. By reassigning codons to noncanonical amino acids, the generation of new-to-nature proteins offers countless opportunities for bioproduction and biocontainment in industrial chassis. A key bottleneck in genome recoding efforts, however, is the low efficiency of recombineering, which hinders large-scale applications at acceptable speed and cost. To relieve this bottleneck, we developed ReScribe, a highly optimized recombineering tool enhanced by CRISPR-Cas9-mediated counterselection built upon the minimal PAM 5'-NNG-3' of the Streptococcus canis Cas9 (ScCas9). As a proof of concept, we used ReScribe to generate a minimally recoded strain of the industrial chassis Pseudomonas putida by replacing TAG stop codons (functioning as PAMs) of essential metabolic genes with the synonymous TAA. We showed that ReScribe enables nearly 100% engineering efficiency of multiple loci in P. putida, opening promising avenues for genome editing and applications thereof in this bacterium and beyond.

Sense codon reassignment enables viral resistance and encoded polymer synthesis

It is widely hypothesized that removing cellular transfer RNAs (tRNAs)-making their cognate codons unreadable-might create a genetic firewall to viral infection and enable sense codon reassignment. However, it has been impossible to test these hypotheses. In this work, following synonymous codon compression and laboratory evolution in Escherichia coli, we deleted the tRNAs and release factor 1, which normally decode two sense codons and a stop codon; the resulting cells could not read the canonical genetic code and were completely resistant to a cocktail of viruses. We reassigned these codons to enable the efficient synthesis of proteins containing three distinct noncanonical amino acids. Notably, we demonstrate the facile reprogramming of our cells for the encoded translation of diverse noncanonical heteropolymers and macrocycles.

Effective Population Size Predicts Local Rates but Not Local Mitigation of Read-through Errors

In correctly predicting that selection efficiency is positively correlated with the effective population size (Ne), the nearly neutral theory provides a coherent understanding of between-species variation in numerous genomic parameters, including heritable error (germline mutation) rates. Does the same theory also explain variation in phenotypic error rates and in abundance of error mitigation mechanisms? Translational read-through provides a model to investigate both issues as it is common, mostly nonadaptive, and has good proxy for rate (TAA being the least leaky stop codon) and potential error mitigation via "fail-safe" 3' additional stop codons (ASCs). Prior theory of translational read-through has suggested that when population sizes are high, weak selection for local mitigation can be effective thus predicting a positive correlation between ASC enrichment and Ne. Contra to prediction, we find that ASC enrichment is not correlated with Ne. ASC enrichment, although highly phylogenetically patchy, is, however, more common both in unicellular species and in genes expressed in unicellular modes in multicellular species. By contrast, Ne does positively correlate with TAA enrichment. These results imply that local phenotypic error rates, not local mitigation rates, are consistent with a drift barrier/nearly neutral model.

Engineering molecular translation systems

Molecular translation systems provide a genetically encoded framework for protein synthesis, which is essential for all life. Engineering these systems to incorporate non-canonical amino acids (ncAAs) into peptides and proteins has opened many exciting opportunities in chemical and synthetic biology. Here, we review recent advances that are transforming our ability to engineer molecular translation systems. In cell-based systems, new processes to synthesize recoded genomes, tether ribosomal subunits, and engineer orthogonality with high-throughput workflows have emerged. In cell-free systems, adoption of flexizyme technology and cell-free ribosome synthesis and evolution platforms are expanding the limits of chemistry at the ribosome's RNA-based active site. Looking forward, innovations will deepen understanding of molecular translation and provide a path to polymers with previously unimaginable structures and functions.

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

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

Spurious regulatory connections dictate the expression-fitness landscape of translation factors

During steady-state cell growth, individual enzymatic fluxes can be directly inferred from growth rate by mass conservation, but the inverse problem remains unsolved. Perturbing the flux and expression of a single enzyme could have pleiotropic effects that may or may not dominate the impact on cell fitness. Here, we quantitatively dissect the molecular and global responses to varied expression of translation termination factors (peptide release factors, RFs) in the bacterium Bacillus subtilis. While endogenous RF expression maximizes proliferation, deviations in expression lead to unexpected distal regulatory responses that dictate fitness reduction. Molecularly, RF depletion causes expression imbalance at specific operons, which activates master regulators and detrimentally overrides the transcriptome. Through these spurious connections, RF abundances are thus entrenched by focal points within the regulatory network, in one case located at a single stop codon. Such regulatory entrenchment suggests that predictive bottom-up models of expression-fitness landscapes will require near-exhaustive characterization of parts.

Evolutionary epidemiology of Streptococcus iniae: Linking mutation rate dynamics with adaptation to novel immunological landscapes

Pathogens continuously adapt to changing host environments where variation in their virulence and antigenicity is critical to their long-term evolutionary success. The emergence of novel variants is accelerated in microbial mutator strains (mutators) deficient in DNA repair genes, most often from mismatch repair and oxidized-guanine repair systems (MMR and OG respectively). Bacterial MMR/OG mutants are abundant in clinical samples and show increased adaptive potential in experimental infection models, yet the role of mutators in the epidemiology and evolution of infectious disease is not well understood. Here we investigated the role of mutation rate dynamics in the evolution of a broad host range pathogen, Streptococcus iniae, using a set of 80 strains isolated globally over 40 years. We have resolved phylogenetic relationships using non-recombinant core genome variants, measured in vivo mutation rates by fluctuation analysis, identified variation in major MMR/OG genes and their regulatory regions, and phenotyped the major traits determining virulence in streptococci. We found that both mutation rate and MMR/OG genotype are remarkably conserved within phylogenetic clades but significantly differ between major phylogenetic lineages. Further, variation in MMR/OG loci correlates with occurrence of atypical virulence-associated phenotypes, infection in atypical hosts (mammals), and atypical (osseous) tissue of a vaccinated primary host. These findings suggest that mutators are likely to facilitate adaptations preceding major diversification events and may promote emergence of variation permitting colonization of a novel host tissue, novel host taxa (host jumps), and immune-escape in the vaccinated host.

Collateral Toxicity Limits the Evolution of Bacterial Release Factor 2 toward Total Omnipotence

When new genes evolve through modification of existing genes, there are often tradeoffs between the new and original functions, making gene duplication and amplification necessary to buffer deleterious effects on the original function. We have used experimental evolution of a bacterial strain lacking peptide release factor 1 (RF1) in order to study how peptide release factor 2 (RF2) evolves to compensate the loss of RF1. As expected, amplification of the RF2-encoding gene prfB to high copy number was a rapid initial response, followed by the appearance of mutations in RF2 and other components of the translation machinery. Characterization of the evolved RF2 variants by their effects on bacterial growth rate, reporter gene expression, and in vitro translation termination reveals a complex picture of reduced discrimination between the cognate and near-cognate stop codons and highlights a functional tradeoff that we term "collateral toxicity." We suggest that this type of tradeoff may be a more serious obstacle in new gene evolution than the more commonly discussed evolutionary tradeoffs between "old" and "new" functions of a gene, as it cannot be overcome by gene copy number changes. Further, we suggest a model for how RF2 autoregulation responds to alterations in the demand not only for RF2 activity but also for RF1 activity.

Gene autoregulation by 3' UTR-derived bacterial small RNAs

Negative feedback regulation, that is the ability of a gene to repress its own synthesis, is the most abundant regulatory motif known to biology. Frequently reported for transcriptional regulators, negative feedback control relies on binding of a transcription factor to its own promoter. Here, we report a novel mechanism for gene autoregulation in bacteria relying on small regulatory RNA (sRNA) and the major endoribonuclease, RNase E. TIER-seq analysis (transiently-inactivating-an-endoribonuclease-followed-by-RNA-seq) revealed ~25,000 RNase E-dependent cleavage sites in Vibrio cholerae, several of which resulted in the accumulation of stable sRNAs. Focusing on two examples, OppZ and CarZ, we discovered that these sRNAs are processed from the 3' untranslated region (3' UTR) of the oppABCDF and carAB operons, respectively, and base-pair with their own transcripts to inhibit translation. For OppZ, this process also triggers Rho-dependent transcription termination. Our data show that sRNAs from 3' UTRs serve as autoregulatory elements allowing negative feedback control at the post-transcriptional level.

In eubacteria, unlike eukaryotes, there is no evidence for selection favouring fail-safe 3' additional stop codons

Errors throughout gene expression are likely deleterious, hence genomes are under selection to ameliorate their consequences. Additional stop codons (ASCs) are in-frame nonsense ‘codons’ downstream of the primary stop which may be read by translational machinery should the primary stop have been accidentally read through. Prior evidence in several eukaryotes suggests that ASCs are selected to prevent potentially-deleterious consequences of read-through. We extend this evidence showing that enrichment of ASCs is common but not universal for single cell eukaryotes. By contrast, there is limited evidence as to whether the same is true in other taxa. Here, we provide the first systematic test of the hypothesis that ASCs act as a fail-safe mechanism in eubacteria, a group with high read-through rates. Contra to the predictions of the hypothesis we find: there is paucity, not enrichment, of ASCs downstream; substitutions that degrade stops are more frequent in-frame than out-of-frame in 3’ sequence; highly expressed genes are no more likely to have ASCs than lowly expressed genes; usage of the leakiest primary stop (TGA) in highly expressed genes does not predict ASC enrichment even compared to usage of non-leaky stops (TAA) in lowly expressed genes, beyond downstream codon +1. Any effect at the codon immediately proximal to the primary stop can be accounted for by a preference for a T/U residue immediately following the stop, although if anything, TT- and TC- starting codons are preferred. We conclude that there is no compelling evidence for ASC selection in eubacteria. This presents an unusual case in which the same error could be solved by the same mechanism in eukaryotes and prokaryotes but is not. We discuss two possible explanations: that, owing to the absence of nonsense mediated decay, bacteria may solve read-through via gene truncation and in eukaryotes certain prion states cause raised read-through rates.

In all organisms, gene expression is error-prone. One such error, translational read-through, occurs where the primary stop codon of an expressed gene is missed by the translational machinery. Failure to terminate is likely to be costly, hence genomes are under selection to prevent this from happening. One proposed error-proofing strategy involves in-frame proximal additional stop codons (ASCs) which may act as a ‘fail-safe’ mechanism by providing another opportunity for translation to terminate. There is evidence for ASC enrichment in several eukaryotes. We extend this evidence showing it to be common but not universal in single celled eukaryotes. However, the situation in bacteria is poorly understood, despite bacteria having high read-through rates. Here, we test the fail-safe hypothesis within a broad range of bacteria. To our surprise, we find that not only are ASCs not enriched, but they may even be selected against. This provides evidence for an unusual circumstance where eukaryotes and prokaryotes could solve the same problem the same way but don’t. What are we to make of this? We suggest that if read-through is the problem, ASCs are not necessarily the expected solution. Owing to the absence of nonsense-mediated decay, a process that makes gene truncation in eukaryotes less viable, we propose bacteria may rescue a leaky stop by mutation that creates a new stop upstream. Alternatively, raised read-through rates in some particular conditions in eukaryotes might explain the difference.

Understanding the Genetic Code

The universal triple-nucleotide genetic code is often viewed as a given, randomly selected through evolution. However, as summarized in this article, many observations and deductions within structural and thermodynamic frameworks help to explain the forces that must have shaped the code during the early evolution of life on Earth.

Forces acting on codon bias in malaria parasites

Malaria parasite genomes have a range of codon biases, with Plasmodium falciparum one of the most AT-biased genomes known. We examined the make up of synonymous coding sites and stop codons in the core genomes of representative malaria parasites, showing first that local DNA context influences codon bias similarly across P. falciparum, P. vivax and P. berghei, with suppression of CpG dinucleotides and enhancement of CpC dinucleotides, both within and aross codons. Intense asexual phase gene expression in P. falciparum and P. berghei is associated with increased A3:G3 bias but reduced T3:C3 bias at 2-fold sites, consistent with adaptation of codons to tRNA pools and avoidance of wobble tRNA interactions that potentially slow down translation. In highly expressed genes, the A3:G3 ratio can exceed 30-fold while the T3:C3 ratio can be less than 1, according to the encoded amino acid and subsequent base. Lysine codons (AAA/G) show distinctive behaviour with substantially reduced A3:G3 bias in highly expressed genes, perhaps because of selection against frameshifting when the AAA codon is followed by another adenine. Intense expression is also associated with a strong bias towards TAA stop codons (found in 94% and 89% of highly expressed P. falciparum and P. berghei genes respectively) and a proportional rise in the TAAA stop ‘tetranucleotide’. The presence of these expression-linked effects in the relatively AT-rich malaria parasite species adds weight to the suggestion that AT-richness in the Plasmodium genus might be a fitness adaptation. Potential explanations for the relative lack of codon bias in P. vivax include the distinct features of its lifecycle and its effective population size over evolutionary time.

Refining the Ambush Hypothesis: Evidence That GC- and AT-Rich Bacteria Employ Different Frameshift Defence Strategies

Stop codons are frequently selected for beyond their regular termination function for error control. The “ambush hypothesis” proposes out-of-frame stop codons (OSCs) terminating frameshifted translations are selected for. Although early indirect evidence was partially supportive, recent evidence suggests OSC frequencies are not exceptional when considering underlying nucleotide content. However, prior null tests fail to control amino acid/codon usages or possible local mutational biases. We therefore return to the issue using bacterial genomes, considering several tests defining and testing against a null. We employ simulation approaches preserving amino acid order but shuffling synonymous codons or preserving codons while shuffling amino acid order. Additionally, we compare codon usage in amino acid pairs, where one codon can but the next, otherwise identical codon, cannot encode an OSC. OSC frequencies exceed expectations typically in AT-rich genomes, the +1 frame and for TGA/TAA but not TAG. With this complex evidence, simply rejecting or accepting the ambush hypothesis is not warranted. We propose a refined post hoc model, whereby AT-rich genomes have more accidental frameshifts, handled by RF2–RF3 complexes (associated with TGA/TAA) and are mostly +1 (or −2) slips. Supporting this, excesses positively correlate with in silico predicted frameshift probabilities. Thus, we propose a more viable framework, whereby genomes broadly adopt one of the two strategies to combat frameshifts: preventing frameshifting (GC-rich) or permitting frameshifts but minimizing impacts when most are caught early (AT-rich). Our refined framework holds promise yet some features, such as the bias of out-of-frame sense codons, remain unexplained.

Bioinformatics and Translation Elongation

Codon usage depends on mutation bias, tRNA-mediated selection, and the need for high efficiency and accuracy in translation. One codon in a synonymous codon family is often strongly over-used, especially in highly expressed genes, which often leads to a high dN/dS ratio because dS is very small. Many different codon usage indices have been proposed to measure codon usage and codon adaptation. Sense codon could be misread by release factors and stop codons misread by tRNAs, which also contribute to codon usage in rare cases. This chapter outlines the conceptual framework on codon evolution, illustrates codon-specific and gene-specific codon usage indices, and presents their applications. A new index for codon adaptation that accounts for background mutation bias (Index of Translation Elongation) is presented and contrasted with codon adaptation index (CAI) which does not consider background mutation bias. They are used to re-analyze data from a recent paper claiming that translation elongation efficiency matters little in protein production. The reanalysis disproves the claim.

Purifying and positive selection in the evolution of stop codons

Modes of evolution of stop codons in protein-coding genes, especially the conservation of UAA, have been debated for many years. We reconstructed the evolution of stop codons in 40 groups of closely related prokaryotic and eukaryotic genomes. The results indicate that the UAA codons are maintained by purifying selection in all domains of life. In contrast, positive selection appears to drive switches from UAG to other stop codons in prokaryotes but not in eukaryotes. Changes in stop codons are significantly associated with increased substitution frequency immediately downstream of the stop. These positions are otherwise more strongly conserved in evolution compared to sites farther downstream, suggesting that such substitutions are compensatory. Although GC content has a major impact on stop codon frequencies, its contribution to the decreased frequency of UAA differs between bacteria and archaea, presumably, due to differences in their translation termination mechanisms.

Operative Binding of Class I Release Factors and YaeJ Stabilizes the Ribosome in the Nonrotated State

During translation, the small subunit of the ribosome rotates with respect to the large subunit primarily between two states as mRNA is being translated into a protein. At the termination of bacterial translation, class I release factors (RFs) bind to a stop codon in the A-site and catalyze the release of the peptide chain from the ribosome. Periodically, mRNA is truncated prematurely, and the translating ribosome stalls at the end of the mRNA forming a nonstop complex requiring one of several ribosome rescue factors to intervene. One factor, YaeJ, is structurally homologous with the catalytic region of RFs but differs by binding to the ribosome directly through its C-terminal tail. Structures of the ribosome show that the ribosome adopts the nonrotated state conformation when these factors are bound. However, these studies do not elucidate the influence of binding to cognate or noncognate codons on the dynamics of intersubunit rotation. Here, we investigate the effects of wild-type and mutant forms of RF1, RF2, and YaeJ binding on ribosome intersubunit rotation using single-molecule Förster resonance energy transfer. We show that both RF1 binding and RF2 binding are sufficient to shift the population of posthydrolysis ribosome complexes from primarily the rotated to the nonrotated state only when a cognate stop codon is present in the A-site. Similarly, YaeJ binding stabilizes nonstop ribosomal complexes in the nonrotated state. Along with previous studies, these results are consistent with the idea that directed conformational changes and binding of subsequent factors to the ribosome are requisite for efficient termination and ribosome recycling.

The production of SPusp45-MSP-1₁₉ gene construct and its recombinant protein in <em>Lactococcus lactis</em> to be used as a malaria vaccine

Background: Merozoite surface protein 1 (MSP-1) is a major protein used by the Plasmodium during red blood cells invasion in malaria. MSP-119, one of MSP-1 is highly conserved, and  it is a potential malaria vaccine candidate because the monoclonal antibodies are capable blocking erythrocyte invasion in vitro. The aim of this study was to produce MSP-119 gene construct and the recombinant protein in Lactococcus lactis.Methods: Usp45-MSP-119, derived from codon optimization and the synthetic gene, was inserted into the pMAT cloning vector. A vector expressing MSP-119 included usp45 has been constructed by the manipulation of recombinant DNA using restriction enzymes. The MSP-119 protein was expressed to 45% ammonium sulfate precipitation and purified using Sephadex-G50 gel filtration chromatography. The expressed protein was characterized by SDS-PAGE and dot blot.Results: usp45-MSP-119 gene was amplified using specific primers and inserted into the multiple cloning sites in the expression vector pNZ8148 with size 3,538 bp as a recombinant vector. The protein of  MSP-119 was successfully expressed in L. lactis with molecular weight of 10.45 kDa. The dot blot was tested in 3 different comparisons between the host cells, non-induced cells, and induced cells with 10 ng/ml nisin. The results showed that 10 ng/ml nisin gave a positive reaction as detected by dot blot assay.Conclusion: This study confirmed that the usp45-MSP-119 gene was successfully inserted into the multiple cloning sites of the pNZ8148 expression vector and the MSP-119 protein expressed in the NICE system of the L. lactis host cell.

R213I mutation in release factor 2 (RF2) is one step forward for engineering an omnipotent release factor in bacteria Escherichia coli

The current understanding of the specificity of the bacterial class I release factors (RFs) in decoding stop codons has evolved beyond a simple tripeptide anticodon model. A recent molecular dynamics study for deciphering the principles for specific stop codon recognition by RFs identified Arg-213 as a crucial residue on Escherichia coli RF2 for discriminating guanine in the third position (G3). Interestingly, Arg-213 is highly conserved in RF2 and substituted by Ile-196 in the corresponding position in RF1. Another similar pair is Leu-126 in RF1 and Asp-143 in RF2, which are also conserved within their respective groups. With the hypothesis that replacement of Arg-213 and Asp-143 with the corresponding RF1 residues will reduce G3 discrimination by RF2, we swapped these residues between E. coli RF1 and RF2 by site-directed mutagenesis and characterized their preference for different codons using a competitive peptide release assay. Among these, the R213I mutant of RF2 showed 5-fold improved reading of the RF1-specific UAG codon relative to UAA, the universal stop codon, compared with the wild type (WT). In-depth fast kinetic studies revealed that the gain in UAG reading by RF2 R213I is associated with a reduced efficiency of termination on the cognate UAA codon. Our work highlights the notion that stop codon recognition involves complex interactions with multiple residues beyond the PXT/SPF motifs. We propose that the R213I mutation in RF2 brings us one step forward toward engineering an omnipotent RF in bacteria, capable of reading all three stop codons.

Consensus architecture of promoters and transcription units in Escherichia coli: design principles for synthetic biology

Genetic information in genomes is ordered, arranged in such a way that it constitutes a code, the so-called cis regulatory code. The regulatory machinery of the cell, termed trans-factors, decodes and expresses this information. In this way, genomes maintain a potential repertoire of genetic programs, parts of which are executed depending on the presence of active regulators in each condition. These genetic programs, executed by the regulatory machinery, have functional units in the genome delimited by punctuation-like marks. In genetic terms, these informational phrases correspond to transcription units, which are the minimal genetic information expressed consistently from initiation to termination marks. Between the start and final punctuation marks, additional marks are present that are read by the transcriptional and translational machineries. In this work, we look at all the experimentally described and predicted genetic elements in the bacterium Escherichia coli K-12 MG1655 and define a comprehensive architectural organization of transcription units to reveal the natural genome-design and to guide the construction of synthetic genetic programs.

A Screen for Protein-Protein Interactions in Live Mycobacteria Reveals a Functional Link between the Virulence-Associated Lipid Transporter LprG and the Mycolyltransferase Antigen 85A

Outer membrane lipids in pathogenic mycobacteria are important for virulence and survival. Although the biosynthesis of these lipids has been extensively studied, mechanisms responsible for their assembly in the outer membrane are not understood. In the study of Gram-negative outer membrane assembly, protein-protein interactions define transport mechanisms, but analogous interactions have not been explored in mycobacteria. Here we identified interactions with the lipid transport protein LprG. Using site-specific photo-cross-linking in live mycobacteria, we mapped three major interaction interfaces within LprG. We went on to identify proteins that cross-link at the entrance to the lipid binding pocket, an area likely relevant to LprG transport function. We verified LprG site-specific interactions with two hits, the conserved lipoproteins LppK and LppI. We further showed that LprG interacts physically and functionally with the mycolyltransferase Ag85A, as loss of either protein leads to similar defects in cell growth and mycolylation. Overall, our results support a model in which protein interactions coordinate multiple pathways in outer membrane biogenesis and connect lipid biosynthesis to transport.

Expanding the Genetic Code of a Photoautotrophic Organism

The photoautotrophic freshwater cyanobacterium Synechococcus elongatus is widely used as a chassis for biotechnological applications as well as a photosynthetic bacterial model. In this study, a method for expanding the genetic code of this cyanobacterium has been established, thereby allowing the incorporation of unnatural amino acids into proteins. This was achieved through UAG stop codon suppression, using an archaeal pyrrolysyl orthogonal translation system. We demonstrate incorporation of unnatural amino acids into green fluorescent protein with 20 ± 3.5% suppression efficiency. The introduced components were shown to be orthogonal to the host translational machinery. In addition, we observed that no significant growth impairment resulted from the integration of the system. To interpret the observations, we modeled and investigated the competition over the UAG codon between release factor 1 and pyl-tRNA_CUA. On the basis of the model results, and the fact that 39.6% of the stop codons in the S. elongatus genome are UAG stop codons, the suppression efficiency in S. elongatus is unexpectedly high. The reason for this unexpected suppression efficiency has yet to be determined.

Application of Noncanonical Amino Acids for Protein Labeling in a Genomically Recoded Escherichia coli

Small synthetic fluorophores are in many ways superior to fluorescent proteins as labels for imaging. A major challenge is to use them for a protein-specific labeling in living cells. Here, we report on our use of noncanonical amino acids that are genetically encoded via the pyrrolysyl-tRNA/pyrrolysyl-RNA synthetase pair at artificially introduced TAG codons in a recoded E. coli strain. The strain is lacking endogenous TAG codons and the TAG-specific release factor RF1. The amino acids contain bioorthogonal groups that can be clicked to externally supplied dyes, thus enabling protein-specific labeling in live cells. We find that the noncanonical amino acid incorporation into the target protein is robust for diverse amino acids and that the usefulness of the recoded E. coli strain mainly derives from the absence of release factor RF1. However, the membrane permeable dyes display high nonspecific binding in intracellular environment and the electroporation of hydrophilic nonmembrane permeable dyes severely impairs growth of the recoded strain. In contrast, proteins exposed on the outer membrane of E. coli can be labeled with hydrophilic dyes with a high specificity as demonstrated by labeling of the osmoporin OmpC. Here, labeling can be made sufficiently specific to enable single molecule studies as exemplified by OmpC single particle tracking.

Complete mitogenome mapping of potato late blight pathogen, Phytophthora infestans A2 mating type

Complete mitochondrial genome of Phytophthora infestans, A₂ mating type (MT) with a size of ≅37,767 bp was sequenced. A total of 53 protein-coding genes are predicted on both strands, including 25 tRNA, 2 rRNA, and 18 respiratory proteins. Gene order of A₂MT was consistent with that established in A₁, despite high level of polymorphism in both coding and non-coding regions. The mtDNA of A₂MT was found to have 99.5% and 99.4% homology with Ia and Ib, whereas 94.7% and 94.3% with IIa and IIb, respectively. Study of repeats revealed a dinucleotide (AT)₉ specific to A₁ and homology of cox1 gene sequence revealed the relationship among 50 Phytophthora species.

The Role of +4U as an Extended Translation Termination Signal in Bacteria

Termination efficiency of stop codons depends on the first 3' flanking (+4) base in bacteria and eukaryotes. In both Escherichia coli and Saccharomyces cerevisiae, termination read-through is reduced in the presence of +4U; however, the molecular mechanism underlying +4U function is poorly understood. Here, we perform comparative genomics analysis on 25 bacterial species (covering Actinobacteria, Bacteriodetes, Cyanobacteria, Deinococcus-Thermus, Firmicutes, Proteobacteria, and Spirochaetae) with bioinformatics approaches to examine the influence of +4U in bacterial translation termination by contrasting highly- and lowly-expressed genes (HEGs and LEGs, respectively). We estimated gene expression using the recently formulated Index of Translation Elongation, ITE, and identified stop codon near-cognate transfer RNAs (tRNAs) from well-annotated genomes. We show that +4U was consistently overrepresented in UAA-ending HEGs relative to LEGs. The result is consistent with the interpretation that +4U enhances termination mainly for UAA. Usage of +4U decreases in GC-rich species where most stop codons are UGA and UAG, with few UAA-ending genes, which is expected if UAA usage in HEGs drives up +4U usage. In HEGs, +4U usage increases significantly with abundance of UAA nc_tRNAs (near-cognate tRNAs that decode codons differing from UAA by a single nucleotide), particularly those with a mismatch at the first stop codon site. UAA is always the preferred stop codon in HEGs, and our results suggest that UAAU is the most efficient translation termination signal in bacteria.

A Nonsense Mutation in Mycobacterium marinum That Is Suppressible by a Novel Mechanism

Mycobacterial pathogens use the ESAT-6 system 1 (Esx-1) exporter to promote virulence. Previously, we used gene disruption and complementation to conclude that the MMAR_0039 gene in Mycobacterium marinum is required to promote Esx-1 export. Here we applied molecular genetics, proteomics, and whole-genome sequencing to demonstrate that the MMAR_0039 gene is not required for Esx-1 secretion or virulence. These findings suggest that we initially observed an indirect mechanism of genetic complementation. We identified a spontaneous nonsense mutation in a known Esx-1-associated gene which causes a loss of Esx-1 activity. We show that the Esx-1 function was restored by nonsense suppression. Moreover, we identified a polar mutation in the ppsC gene which reduced cellular impermeability but did not impact cytotoxicity in macrophages. Our studies reveal insight into Esx-1 export, nonsense suppression, and cell envelope lipid biogenesis.

Coevolution between Stop Codon Usage and Release Factors in Bacterial Species

Three stop codons in bacteria represent different translation termination signals, and their usage is expected to depend on their differences in translation termination efficiency, mutation bias, and relative abundance of release factors (RF1 decoding UAA and UAG, and RF2 decoding UAA and UGA). In 14 bacterial species (covering Proteobacteria, Firmicutes, Cyanobacteria, Actinobacteria and Spirochetes) with cellular RF1 and RF2 quantified, UAA is consistently over-represented in highly expressed genes (HEGs) relative to lowly expressed genes (LEGs), whereas UGA usage is the opposite even in species where RF2 is far more abundant than RF1. UGA usage relative to UAG increases significantly with P_RF2 [=RF2/(RF1 + RF2)] as expected from adaptation between stop codons and their decoders. P_RF2 is > 0.5 over a wide range of AT content (measured by P_AT3 as the proportion of AT at third codon sites), but decreases rapidly toward zero at the high range of P_AT3. This explains why bacterial lineages with high P_AT3 often have UGA reassigned because of low RF2. There is no indication that UAG is a minor stop codon in bacteria as claimed in a recent publication. The claim is invalid because of the failure to apply the two key criteria in identifying a minor codon: (1) it is least preferred by HEGs (or most preferred by LEGs) and (2) it corresponds to the least abundant decoder. Our results suggest a more plausible explanation for why UAA usage increases, and UGA usage decreases, with P_AT3, but UAG usage remains low over the entire P_AT3 range.

Augmented genetic decoding: global, local and temporal alterations of decoding processes and codon meaning

The non-universality of the genetic code is now widely appreciated. Codes differ between organisms, and certain genes are known to alter the decoding rules in a site-specific manner. Recently discovered examples of decoding plasticity are particularly spectacular. These examples include organisms and organelles with disruptions of triplet continuity during the translation of many genes, viruses that alter the entire genetic code of their hosts and organisms that adjust their genetic code in response to changing environments. In this Review, we outline various modes of alternative genetic decoding and expand existing terminology to accommodate recently discovered manifestations of this seemingly sophisticated phenomenon.

Protein Mis-Termination Initiates Genetic Diseases, Cancers, and Restricts Bacterial Genome Expansion

Protein termination is an important cellular process. Protein termination relies on the stop-codons in the mRNA interacting properly with the releasing factors on the ribosome. One third of inherited diseases, including cancers, are associated with the mutation of the stop-codons. Many pathogens and viruses are able to manipulate their stop-codons to express their virulence. The influence of stop-codons is not limited to the primary reading frame of the genes. Stop-codons in the second and third reading frames are referred as premature stop signals (PSC). Stop-codons and PSCs together are collectively referred as stop-signals. The ratios of the stop-signals (referred as translation stop-signals ratio or TSSR) of genetically related bacteria, despite their great differences in gene contents, are much alike. This nearly identical Genomic-TSSR value of genetically related bacteria may suggest that bacterial genome expansion is limited by their unique stop-signals bias. We review the protein termination process and the different types of stop-codon mutation in plants, animals, microbes, and viruses, with special emphasis on the role of PSCs in directing bacterial evolution in their natural environments. Knowing the limit of genomic boundary could facilitate the formulation of new strategies in controlling the spread of diseases and combat antibiotic-resistant bacteria.