Autoregulation of the tufB operon in Salmonella

Fitness benefits of a synonymous substitution in an ancient EF-Tu gene depend on the genetic background

Synonymous mutations are changes to DNA sequence, which occur within translated genes but which do not affect the protein sequence. Although often referred to as silent mutations, evidence suggests that synonymous mutations can affect gene expression, mRNA stability, and even translation efficiency. A collection of both experimental and bioinformatic data has shown that synonymous mutations can impact cell phenotype, yet less is known about the molecular mechanisms and potential of beneficial or adaptive effects of such changes within evolved populations. Here, we report a beneficial synonymous mutation acquired via experimental evolution in an essential gene variant encoding the translation elongation factor protein EF-Tu. We demonstrate that this particular synonymous mutation increases EF-Tu mRNA and protein levels as well as global polysome abundance on RNA transcripts. Although presence of the synonymous mutation is clearly causative of such changes, we also demonstrate that fitness benefits are highly contingent on other potentiating mutations present within the genetic background in which the mutation arose. Our results underscore the importance of beneficial synonymous mutations, especially those that affect levels of proteins that are key for cellular processes.IMPORTANCEThis study explores the degree to which synonymous mutations in essential genes can influence adaptation in bacteria. An experimental system whereby an Escherichia coli strain harboring an engineered translation protein elongation factor-Tu (EF-Tu) was subjected to laboratory evolution. We find that a synonymous mutation acquired on the gene encoding for EF-Tu is conditionally beneficial for bacterial fitness. Our findings provide insight into the importance of the genetic background when a synonymous substitution is favored by natural selection and how such changes have the potential to impact evolution when critical cellular processes are involved.

A beneficial synonymous substitution in EF-Tu is contingent on genetic background

Synonymous mutations are changes to DNA sequence that occur within translated genes but which do not affect the protein sequence. Although often referred to as silent mutations, evidence suggests that synonymous mutations can affect gene expression, mRNA stability, and even translation efficiency. A collection of both experimental and bioinformatic data has shown that synonymous mutations can impact cell phenotype, yet less is known about the molecular mechanisms and potential of beneficial or adaptive effects of such changes within evolved populations. Here, we report a beneficial synonymous mutation acquired via experimental evolution in an essential gene variant encoding the translation Elongation Factor protein EF-Tu. We demonstrate that this particular synonymous mutation increases EF-Tu mRNA and protein levels, as well as the polysome abundance on global transcripts. Although presence of the synonymous mutation is clearly causative of such changes, we also demonstrate that fitness benefits are highly contingent on other potentiating mutations present within the genetic background in which the mutation arose. Our results underscore the importance of beneficial synonymous mutations, especially those that affect levels of proteins that are key for cellular processes.

Regulatory interplay between small RNAs and transcription termination factor Rho

The largest and best studied group of regulatory small RNAs (sRNAs) in bacteria act by modulating translation or turnover of messenger RNAs (mRNAs) through base-pairing interactions that typically take place near the 5' end of the mRNA. This allows the sRNA to bind the complementary target sequence while the remainder of the mRNA is still being made, creating conditions whereby the action of the sRNA can extend to transcriptional steps, most notably transcription termination. Increasing evidence corroborates the existence of a functional interplay between sRNAs and termination factor Rho. Two general mechanisms have emerged. One mechanism operates in translated regions subjected to sRNA repression. By inhibiting ribosome binding co-transcriptionally, the sRNA uncouples translation from transcription, allowing Rho to bind the nascent RNA and promote termination. In the second mechanism, which functions in 5' untranslated regions, the sRNA antagonizes termination directly by interfering with Rho binding to the RNA or the subsequent translocation along the RNA. Here, we review the above literature in the context of other mechanisms that underlie the participation of Rho-dependent transcription termination in gene regulation. This article is part of a Special Issue entitled: RNA and gene control in bacteria edited by Dr. M. Guillier and F. Repoila.

Operon Concatenation Is an Ancient Feature That Restricts the Potential to Rearrange Bacterial Chromosomes

The last common ancestor of the Gammaproteobacteria carried an important 40-kb chromosome section encoding 51 proteins of the transcriptional and translational machinery. These genes were organized into eight contiguous operons (rrnB-tufB-secE-rpoBC-str-S10-spc-alpha). Over 2 Gy of evolution, in different lineages, some of the operons became separated by multigene insertions. Surprisingly, in many Enterobacteriaceae, much of the ancient organization is conserved, indicating a strong selective force on the operons to remain colinear. Here, we show for one operon pair, tufB-secE in Salmonella, that an interruption of contiguity significantly reduces growth rate. Our data show that the tufB-secE operons are concatenated by an interoperon terminator–promoter overlap that plays a significant role regulating gene expression. Interrupting operon contiguity interferes with this regulation, reducing cellular fitness. Six operons of the ancestral chromosome section remain contiguous in Salmonella (tufB-secE-rpoBC and S10-spc-alpha) and, strikingly, each of these operon pairs is also connected by an interoperon terminator–promoter overlap. Accordingly, we propose that operon concatenation is an ancient feature that restricts the potential to rearrange bacterial chromosomes and can select for the maintenance of a colinear operon organization over billions of years.

A multivariate prediction model for Rho-dependent termination of transcription

Bacterial transcription termination proceeds via two main mechanisms triggered either by simple, well-conserved (intrinsic) nucleic acid motifs or by the motor protein Rho. Although bacterial genomes can harbor hundreds of termination signals of either type, only intrinsic terminators are reliably predicted. Computational tools to detect the more complex and diversiform Rho-dependent terminators are lacking. To tackle this issue, we devised a prediction method based on Orthogonal Projections to Latent Structures Discriminant Analysis [OPLS-DA] of a large set of in vitro termination data. Using previously uncharacterized genomic sequences for biochemical evaluation and OPLS-DA, we identified new Rho-dependent signals and quantitative sequence descriptors with significant predictive value. Most relevant descriptors specify features of transcript C>G skewness, secondary structure, and richness in regularly-spaced 5'CC/UC dinucleotides that are consistent with known principles for Rho-RNA interaction. Descriptors collectively warrant OPLS-DA predictions of Rho-dependent termination with a ∼85% success rate. Scanning of the Escherichia coli genome with the OPLS-DA model identifies significantly more termination-competent regions than anticipated from transcriptomics and predicts that regions intrinsically refractory to Rho are primarily located in open reading frames. Altogether, this work delineates features important for Rho activity and describes the first method able to predict Rho-dependent terminators in bacterial genomes.

A Fluorogenic Assay To Monitor Rho-Dependent Termination of Transcription

Transcription termination mediated by the ring-shaped, ATP-dependent Rho motor is a multipurpose regulatory mechanism specific to bacteria and constitutes an interesting target for the development of new antibiotics. Although Rho-dependent termination can punctuate gene expression or contribute to the protection of the genome at hundreds of sites within a given bacterium, its exact perimeter and site- or species-specific features remain insufficiently characterized. New advanced approaches are required to explore thoroughly the diversity of Rho-dependent terminators and the complexity of associated mechanisms. Current in vitro analyses of Rho-dependent termination rely on radiolabeling, gel electrophoresis, and phosphorimaging of transcription reaction products and are thus hazardous, inconvenient, and low-throughput. To address these limitations, we have developed the first in vitro assay using a fluorescence detection modality to study Rho-dependent transcription termination. This powerful experimental tool accurately estimates terminator strengths in a matter of minutes and is optimized for a microplate reader format allowing multiplexed characterization of putative terminator sequences and mechanisms or high-throughput screening of new drugs targeting Rho-dependent termination.

Relation of the pdxB-usg-truA-dedA Operon and the truA Gene to the Intracellular Survival of Salmonella enterica Serovar Typhimurium

Salmonella is the genus of Gram-negative, facultative intracellular pathogens that have the ability to infect large numbers of animal or human hosts. The S. enterica usg gene is associated with intracellular survival based on ortholog screening and identification. In this study, the λ-Red recombination system was used to construct gene deletion strains and to investigate whether the identified operon was related to intracellular survival. The pdxB-usg-truA-dedA operon enhanced the intracellular survival of S. enterica by resisting the oxidative environment and the usg and truA gene expression was induced by H₂O₂. Moreover, the genes in this operon (except for dedA) contributed to virulence in mice. These findings indicate that the pdxB-usg-truA-dedA operon functions in resistance to oxidative environments during intracellular survival and is required for in vivo S. enterica virulence. This study provides insight toward a better understand of the characteristics of intracellular pathogens and explores the gene modules involved in their intracellular survival.

Proteolytic Surface-Shaving and Serotype-Dependent Expression of SPI-1 Invasion Proteins in Salmonella enterica Subspecies enterica

We performed proteolytic surface-shaving with trypsin on three strains/sevovars of Salmonella enterica enterica (SEE): Newport, Kentucky, and Thompson. Surfaced-exposed proteins of live bacterial cells were digested for 15 min. A separate 20 h re-digestion was also performed on the supernatant of each shaving experiment to more completely digest protein fragments into detectable peptides for proteomic analysis by nano-liquid chromatography-electrospray ionization-Orbitrap mass spectrometry. Control samples (i.e., no trypsin during surface-shaving step) were also performed in parallel. We detected peptides of flagella proteins: FliC (filament), FliD (cap), and FlgL (hook-filament junction) as well as peptides of FlgM (anti-σ²⁸ factor), i.e., the negative regulator of flagella synthesis. For SEE Newport and Thompson, we detected Salmonella pathogenicity island 1 (SPI-1) secreted effector/invasion proteins: SipA, SipB, SipC, and SipD, whereas no Sip proteins were detected in control samples. No Sip proteins were detected for SEE Kentucky (or its control) although sip genes were confirmed to be present. Our results may suggest a biological response (<15 min) to proteolysis of live cells for these SEE strains and, in the case of Newport and Thompson, a possible invasion response.

Transcriptional Regulation Buffers Gene Dosage Effects on a Highly Expressed Operon in Salmonella

A feature of bacterial chromosomes is that highly expressed essential genes are usually located close to the origin of replication. Because bacteria have overlapping cycles of replication, genes located close to the origin will often be present in multiple copies, and this is thought to be of selective benefit where high levels of expression support high growth rate. However, the magnitude of this selective effect and whether other forces could be at play are poorly understood. To study this, we translocated a highly expressed essential operon, tufB, to different locations and measured growth fitness. We found that transcriptional regulation buffered the effects of translocation and that even under conditions where growth rate was reduced, genetic changes that increased the expression of tufB were easily and rapidly selected. We conclude, at least for tufB, that forces other than gene dosage may be significant in selecting for chromosomal location.

Highly expressed genes are commonly located close to the origin of replication of bacterial chromosomes (OriC). This location skew is thought to reflect selective advantages associated with gene dosage effects during the replication cycle. The expression of constitutively expressed genes can vary up to fivefold based on chromosomal location, but it is not clear what level of variation would occur in naturally regulated operons. We tested the magnitude of the chromosome location effect using EF-Tu (tufA, tufB), an abundant protein whose cellular level correlates with, and limits, the maximum growth rate. We translocated the Salmonella tufB operon to four locations across the chromosome. The distance from OriC had only a small effect on growth rate, consistent with this operon having the natural ability to upregulate expression and compensate for reduced gene dosage. In contrast, when the total EF-Tu concentration was limiting for the growth rate (tufA deleted), we observed a strong gene dosage effect when tufB was located further from OriC. However, only a short period of experimental evolution was required before the bacteria adapted to this EF-Tu starvation situation by acquiring genetic changes that increased expression levels from the translocated tufB gene, restoring growth rates. Our findings demonstrate that, at least for the tufB operon, gene dosage is probably not the dominant force selecting for a chromosomal location close to OriC. We suggest that the colocation of highly expressed genes close to OriC might instead be selected because it enhances their coregulation during various growth states, with gene dosage being a secondary benefit.

Experimental Determination and Prediction of the Fitness Effects of Random Point Mutations in the Biosynthetic Enzyme HisA

The distribution of fitness effects of mutations is a factor of fundamental importance in evolutionary biology. We determined the distribution of fitness effects of 510 mutants that each carried between 1 and 10 mutations (synonymous and nonsynonymous) in the hisA gene, encoding an essential enzyme in the l-histidine biosynthesis pathway of Salmonella enterica. For the full set of mutants, the distribution was bimodal with many apparently neutral mutations and many lethal mutations. For a subset of 81 single, nonsynonymous mutants most mutations appeared neutral at high expression levels, whereas at low expression levels only a few mutations were neutral. Furthermore, we examined how the magnitude of the observed fitness effects was correlated to several measures of biophysical properties and phylogenetic conservation.We conclude that for HisA: (i) The effect of mutations can be masked by high expression levels, such that mutations that are deleterious to the function of the protein can still be neutral with regard to organism fitness if the protein is expressed at a sufficiently high level; (ii) the shape of the fitness distribution is dependent on the extent to which the protein is rate-limiting for growth; (iii) negative epistatic interactions, on an average, amplified the combined effect of nonsynonymous mutations; and (iv) no single sequence-based predictor could confidently predict the fitness effects of mutations in HisA, but a combination of multiple predictors could predict the effect with a SD of 0.04 resulting in 80% of the mutations predicted within 12% of their observed selection coefficients.

Functional Constraints on Replacing an Essential Gene with Its Ancient and Modern Homologs

Genes encoding proteins that carry out essential informational tasks in the cell, in particular where multiple interaction partners are involved, are less likely to be transferable to a foreign organism. Here, we investigated the constraints on transfer of a gene encoding a highly conserved informational protein, translation elongation factor Tu (EF-Tu), by systematically replacing the endogenous tufA gene in the Escherichia coli genome with its extant and ancestral homologs. The extant homologs represented tuf variants from both near and distant homologous organisms. The ancestral homologs represented phylogenetically resurrected tuf sequences dating from 0.7 to 3.6 billion years ago (bya). Our results demonstrate that all of the foreign tuf genes are transferable to the E. coli genome, provided that an additional copy of the EF-Tu gene, tufB, remains present in the E. coli genome. However, when the tufB gene was removed, only the variants obtained from the gammaproteobacterial family (extant and ancestral) supported growth which demonstrates the limited functional interchangeability of E. coli tuf with its homologs. Relative bacterial fitness correlated with the evolutionary distance of the extant tuf homologs inserted into the E. coli genome. This reduced fitness was associated with reduced levels of EF-Tu and reduced rates of protein synthesis. Increasing the expression of tuf partially ameliorated these fitness costs. In summary, our analysis suggests that the functional conservation of protein activity, the amount of protein expressed, and its network connectivity act to constrain the successful transfer of this essential gene into foreign bacteria.IMPORTANCE Horizontal gene transfer (HGT) is a fundamental driving force in bacterial evolution. However, whether essential genes can be acquired by HGT and whether they can be acquired from distant organisms are very poorly understood. By systematically replacing tuf with ancestral homologs and homologs from distantly related organisms, we investigated the constraints on HGT of a highly conserved gene with multiple interaction partners. The ancestral homologs represented phylogenetically resurrected tuf sequences dating from 0.7 to 3.6 bya. Only variants obtained from the gammaproteobacterial family (extant and ancestral) supported growth, demonstrating the limited functional interchangeability of E. coli tuf with its homologs. Our analysis suggests that the functional conservation of protein activity, the amount of protein expressed, and its network connectivity act to constrain the successful transfer of this essential gene into foreign bacteria.

Elfamycins: inhibitors of elongation factor-Tu

Elfamycins are a relatively understudied group of antibiotics that target the essential process of translation through impairment of EF-Tu function. For the most part, the utility of these compounds has been as laboratory tools for the study of EF-Tu and the ribosome, as their poor pharmacokinetic profile and solubility has prevented implementation as therapeutic agents. However, due to the slowing of the antibiotic pipeline and the rapid emergence of resistance to approved antibiotics, this group is being reconsidered. Some researchers are using screens for novel naturally produced variants, while others are making directed, systematic chemical improvements on publically disclosed compounds. As an example of the latter approach, a GE2270 A derivative, LFF571, has completed phase 2 clinical trials, thus demonstrating the potential for elfamycins to become more prominent antibiotics in the future.

Experimental Evolution of Escherichia coli Harboring an Ancient Translation Protein

The ability to design synthetic genes and engineer biological systems at the genome scale opens new means by which to characterize phenotypic states and the responses of biological systems to perturbations. One emerging method involves inserting artificial genes into bacterial genomes and examining how the genome and its new genes adapt to each other. Here we report the development and implementation of a modified approach to this method, in which phylogenetically inferred genes are inserted into a microbial genome, and laboratory evolution is then used to examine the adaptive potential of the resulting hybrid genome. Specifically, we engineered an approximately 700-million-year-old inferred ancestral variant of tufB, an essential gene encoding elongation factor Tu, and inserted it in a modern Escherichia coli genome in place of the native tufB gene. While the ancient homolog was not lethal to the cell, it did cause a twofold decrease in organismal fitness, mainly due to reduced protein dosage. We subsequently evolved replicate hybrid bacterial populations for 2000 generations in the laboratory and examined the adaptive response via fitness assays, whole genome sequencing, proteomics, and biochemical assays. Hybrid lineages exhibit a general adaptive strategy in which the fitness cost of the ancient gene was ameliorated in part by upregulation of protein production. Our results suggest that an ancient-modern recombinant method may pave the way for the synthesis of organisms that exhibit ancient phenotypes, and that laboratory evolution of these organisms may prove useful in elucidating insights into historical adaptive processes.

Learning from the Leaders: Gene Regulation by the Transcription Termination Factor Rho

The RNA helicase Rho triggers 20-30% of transcription termination events in bacteria. While Rho is associated with most transcription elongation complexes, it only promotes termination of a subset. Recent studies of individual Rho-dependent terminators located within the 5' leader regions of bacterial mRNAs have identified novel mechanisms that govern Rho target specificity and have revealed unanticipated physiological functions for Rho. In particular, the multistep nature of Rho-dependent termination enables regulatory input from determinants beyond the sequence of the Rho loading site, and allows a given Rho-dependent terminator to respond to multiple signals. Further, the unique position of Rho as a sensor of cellular translation has been exploited to regulate the transcription of genes required for protein synthesis, including those specifying Mg(2+) transporters.

Compensating the Fitness Costs of Synonymous Mutations

Synonymous mutations do not change the sequence of the polypeptide but they may still influence fitness. We investigated in Salmonella enterica how four synonymous mutations in the rpsT gene (encoding ribosomal protein S20) reduce fitness (i.e., growth rate) and the mechanisms by which this cost can be genetically compensated. The reduced growth rates of the synonymous mutants were correlated with reduced levels of the rpsT transcript and S20 protein. In an adaptive evolution experiment, these fitness impairments could be compensated by mutations that either caused up-regulation of S20 through increased gene dosage (due to duplications), increased transcription of the rpsT gene (due to an rpoD mutation or mutations in rpsT), or increased translation from the rpsT transcript (due to rpsT mutations). We suggest that the reduced levels of S20 in the synonymous mutants result in production of a defective subpopulation of 30S subunits lacking S20 that reduce protein synthesis and bacterial growth and that the compensatory mutations restore S20 levels and the number of functional ribosomes. Our results demonstrate how specific synonymous mutations can cause substantial fitness reductions and that many different types of intra- and extragenic compensatory mutations can efficiently restore fitness. Furthermore, this study highlights that also synonymous sites can be under strong selection, which may have implications for the use of dN/dS ratios as signature for selection.