Compensating the Fitness Costs of Synonymous Mutations

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.

Whole-genome sequencing analysis of two heat-evolved Escherichia coli strains

BACKGROUND

High temperatures cause a suite of problems for cells, including protein unfolding and aggregation; increased membrane fluidity; and changes in DNA supercoiling, RNA stability, transcription and translation. Consequently, enhanced thermotolerance can evolve through an unknown number of genetic mechanisms even in the simple model bacterium Escherichia coli. To date, each E. coli study exploring this question resulted in a different set of mutations. To understand the changes that can arise when an organism evolves to grow at higher temperatures, we sequenced and analyzed two previously described E. coli strains, BM28 and BM28 ΔlysU, that have been laboratory adapted to the highest E. coli growth temperature reported to date.

RESULTS

We found three large deletions in the BM28 and BM28 ΔlysU strains of 123, 15 and 8.5 kb in length and an expansion of IS10 elements. We found that BM28 and BM28 ΔlysU have considerably different genomes, suggesting that the BM28 culture that gave rise to BM28 and BM28 ΔlysU was a mixed population of genetically different cells. Consistent with published findings of high GroESL expression in BM28, we found that BM28 inexplicitly carries the groESL bearing plasmid pOF39 that was maintained simply by high-temperature selection pressure. We identified over 200 smaller insertions, deletions, single nucleotide polymorphisms and other mutations, including changes in master regulators such as the RNA polymerase and the transcriptional termination factor Rho. Importantly, this genome analysis demonstrates that the commonly cited findings that LysU plays a crucial role in thermotolerance and that GroESL hyper-expression is brought about by chromosomal mutations are based on a previous misinterpretation of the genotype of BM28.

CONCLUSIONS

This whole-genome sequencing study describes genetically distinct mechanisms of thermotolerance evolution from those found in other heat-evolved E. coli strains. Studying adaptive laboratory evolution to heat in simple model organisms is important in the context of climate change. It is important to better understand genetic mechanisms of enhancing thermotolerance in bacteria and other organisms, both in terms of optimizing laboratory evolution methods for various organisms and in terms of potential genetic engineering of organisms most at risk or most important to our societies and ecosystems.

Detection by metagenomic functional analysis and improvement by experimental evolution of β-lactams resistance genes present in oil contaminated soils

The spread of antibiotic resistance genes has become a global health concern identified by the World Health Organization as one of the greatest threats to health. Many of antimicrobial resistance determinants found in bacterial pathogens originate from environmental bacteria, so identifying the genes that confer resistance to antibiotics in different habitats is mandatory to better understand resistance mechanisms. Soil is one of the most diverse environments considered reservoir of antimicrobial resistance genes. The aim of this work is to study the presence of genes that provide resistance to antibiotics used in clinical settings in two oil contaminated soils by metagenomic functional analysis. Using fosmid vectors that efficiently transcribe metagenomic DNA, we have selected 12 fosmids coding for two class A β-lactamases, two subclass B1 and two subclass B3 metallo-β-lactamases, one class D β-lactamase and three efflux pumps that confer resistance to cefexime, ceftriaxone, meropenem and/or imipenem. In some of them, detection of the resistance required heterologous expression from the fosmid promoter. Although initially, these environmental genes only provide resistance to low concentrations of antibiotics, we have obtained, by experimental evolution, fosmid derivatives containing β-lactamase ORFs with a single base substitution, which substantially increase their β-lactamase activity and resistance level. None of the mutations affect β-lactamase coding sequences and are all located upstream of them. These results demonstrate the presence of enzymes that confer resistance to relevant β-lactams in these soils and their capacity to rapidly adapt to provide higher resistance levels.

Dup-In and DIRex: Techniques for Single-Step, Scar-Free Mutagenesis with Marker Recycling

This chapter describes two related recombineering-based techniques: "Duplication Insertion" (Dup-In) and "Direct- and Inverted Repeat stimulated excision" (DIRex). Dup-In is used for transferring existing mutations between strains, and DIRex for generating almost any type of mutation. Both techniques use intermediate insertions with counter-selectable cassettes, flanked by directly repeated sequences that enable exact and spontaneous excision of the cassettes. These constructs can be transferred to other strains using generalized transductions, and the final intended mutation is obtained following selection for spontaneous loss of the counter-selectable cassette, which leaves only the intended mutation behind in the final strain. The techniques have been used in several strains of Escherichia coli and Salmonella enterica, and should be readily adaptable to other organisms where λ Red recombineering or similar methods are available.

The cost of evolved constitutive lac gene expression is usually, but not always, maintained during evolution of generalist populations

Beneficial mutations can become costly following an environmental change. Compensatory mutations can relieve these costs, while not affecting the selected function, so that the benefits are retained if the environment shifts back to be similar to the one in which the beneficial mutation was originally selected. Compensatory mutations have been extensively studied in the context of antibiotic resistance, responses to specific genetic perturbations, and in the determination of interacting gene network components. Few studies have focused on the role of compensatory mutations during more general adaptation, especially as the result of selection in fluctuating environments where adaptations to different environment components may often involve trade-offs. We examine whether costs of a mutation in lacI, which deregulated the expression of the lac operon in evolving populations of Escherichia coli bacteria, were compensated. This mutation occurred in multiple replicate populations selected in environments that fluctuated between growth on lactose, where the mutation was beneficial, and on glucose, where it was deleterious. We found that compensation for the cost of the lacI mutation was rare, but, when it did occur, it did not negatively affect the selected benefit. Compensation was not more likely to occur in a particular evolution environment. Compensation has the potential to remove pleiotropic costs of adaptation, but its rarity indicates that the circumstances to bring about the phenomenon may be peculiar to each individual or impeded by other selected mutations.

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

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

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.

Reconstructing the Evolutionary History of a Highly Conserved Operon Cluster in Gammaproteobacteria and Bacilli

The evolution of gene order rearrangements within bacterial chromosomes is a fast process. Closely related species can have almost no conservation in long-range gene order. A prominent exception to this rule is a >40 kb long cluster of five core operons (secE-rpoBC-str-S10-spc-alpha) and three variable adjacent operons (cysS, tufB, and ecf) that together contain 57 genes of the transcriptional and translational machinery. Previous studies have indicated that at least part of this operon cluster might have been present in the last common ancestor of bacteria and archaea. Using 204 whole genome sequences, ∼2 Gy of evolution of the operon cluster were reconstructed back to the last common ancestors of the Gammaproteobacteria and of the Bacilli. A total of 163 independent evolutionary events were identified in which the operon cluster was altered. Further examination showed that the process of disconnecting two operons generally follows the same pattern. Initially, a small number of genes is inserted between the operons breaking the concatenation followed by a second event that fully disconnects the operons. While there is a general trend for loss of gene synteny over time, there are examples of increased alteration rates at specific branch points or within specific bacterial orders. This indicates the recurrence of relaxed selection on the gene order within bacterial chromosomes. The analysis of the alternation events indicates that segmental genome duplications and/or transposon-directed recombination play a crucial role in rearrangements of the operon cluster.

Genome-Based Analyses of Fitness Effects and Compensatory Changes Associated with Acquisition of bla CMY-, bla CTX-M-, and bla OXA-48/VIM-1-Containing Plasmids in Escherichia coli

(1) Background: Resistance plasmids are under selective conditions beneficial for the bacterial host, but in the absence of selective pressure, this carriage may cause fitness costs. Compensation of this fitness burden is important to obtain competitive ability under antibiotic-free conditions. In this study, we investigated fitness effects after a conjugative transfer of plasmids containing various beta-lactamase genes transferred into Escherichia coli. (2) Methods: Fourteen beta-lactamase-encoding plasmids were transferred from clinical donor strains to E. coli J53. Growth rates were compared for all transconjugants and the recipient. Selected transconjugants were challenged in long-term growth experiments. Growth rates were assessed at different time points during growth for 500 generations. Whole-genome sequencing (WGS) of initial and evolved transconjugants was determined. Results: Most plasmid acquisitions resulted in growth differences, ranging from -4.5% to 7.2%. Transfer of a single bla _CMY-16-carrying plasmid resulted in a growth burden and a growth benefit in independent mating. Long-term growth led to a compensation of fitness burdens and benefits. Analyzing WGS revealed genomic changes caused by Single Nucleotide Polymorphisms (SNPs) and insertion sequences over time. Conclusions: Fitness effects associated with plasmid acquisitions were variable. Potential compensatory mutations identified in transconjugants' genomes after 500 generations give interesting insights into aspects of plasmid-host adaptations.

ΦX174 Attenuation by Whole-Genome Codon Deoptimization

Natural selection acting on synonymous mutations in protein-coding genes influences genome composition and evolution. In viruses, introducing synonymous mutations in genes encoding structural proteins can drastically reduce viral growth, providing a means to generate potent, live-attenuated vaccine candidates. However, an improved understanding of what compositional features are under selection and how combinations of synonymous mutations affect viral growth is needed to predictably attenuate viruses and make them resistant to reversion. We systematically recoded all nonoverlapping genes of the bacteriophage ΦX174 with codons rarely used in its Escherichia coli host. The fitness of recombinant viruses decreases as additional deoptimizing mutations are made to the genome, although not always linearly, and not consistently across genes. Combining deoptimizing mutations may reduce viral fitness more or less than expected from the effect size of the constituent mutations and we point out difficulties in untangling correlated compositional features. We test our model by optimizing the same genes and find that the relationship between codon usage and fitness does not hold for optimization, suggesting that wild-type ΦX174 is at a fitness optimum. This work highlights the need to better understand how selection acts on patterns of synonymous codon usage across the genome and provides a convenient system to investigate the genetic determinants of virulence.

Towards understanding the evolutionary dynamics of mtDNA

Historically, mtDNA was considered a selectively neutral marker that was useful for estimating the population genetic history of the maternal lineage. Over time there has been an increasing appreciation of mtDNA and mitochondria in maintaining cellular and organismal health. Beyond energy production, mtDNA and mitochondria have critical cellular roles in signalling. Here we briefly review the structure of mtDNA and the role of the mitochondrion in energy production. We then discuss the predictions that can be obtained from quaternary structure modelling and focus on mitochondrial complex I. Complex I is the primary entry point for electrons into the electron transport system is the largest respiratory complex of the chain and produces about 40% of the proton flux used to synthesize ATP. A focus of the review is Drosophila's utility as a model organism to study the selective advantage of specific mutations. However, we note that the incorporation of insights from a multitude of systems is necessary to fully understand the range of roles that mtDNA has in organismal fitness. We speculate that dietary changes can illicit stress responses that influence the selective advantage of specific mtDNA mutations and cause spatial and temporal fluctuations in the frequencies of mutations. We conclude that developing our understanding of the roles mtDNA has in determining organismal fitness will enable increased evolutionary insight and propose we can no longer assume it is evolving as a strictly neutral marker without testing this hypothesis.

Population Bottlenecks Strongly Influence the Evolutionary Trajectory to Fluoroquinolone Resistance in Escherichia coli

Experimental evolution is a powerful tool to study genetic trajectories to antibiotic resistance under selection. A confounding factor is that outcomes may be heavily influenced by the choice of experimental parameters. For practical purposes (minimizing culture volumes), most experimental evolution studies with bacteria use transmission bottleneck sizes of 5 × 106 cfu. We currently have a poor understanding of how the choice of transmission bottleneck size affects the accumulation of deleterious versus high-fitness mutations when resistance requires multiple mutations, and how this relates outcome to clinical resistance. We addressed this using experimental evolution of resistance to ciprofloxacin in Escherichia coli. Populations were passaged with three different transmission bottlenecks, including single cell (to maximize genetic drift) and bottlenecks spanning the reciprocal of the frequency of drug target mutations (108 and 1010). The 1010 bottlenecks selected overwhelmingly mutations in drug target genes, and the resulting genotypes corresponded closely to those found in resistant clinical isolates. In contrast, both the 108 and single-cell bottlenecks selected mutations in three different gene classes: 1) drug targets, 2) efflux pump repressors, and 3) transcription-translation genes, including many mutations with low fitness. Accordingly, bottlenecks smaller than the average nucleotide substitution rate significantly altered the experimental outcome away from genotypes observed in resistant clinical isolates. These data could be applied in designing experimental evolution studies to increase their predictive power and to explore the interplay between different environmental conditions, where transmission bottlenecks might vary, and resulting evolutionary trajectories.

Synonymous genome recoding: a tool to explore microbial biology and new therapeutic strategies

Synthetic genome recoding is a new means of generating designed organisms with altered phenotypes. Synonymous mutations introduced into the protein coding region tolerate modifications in DNA or mRNA without modifying the encoded proteins. Synonymous genome-wide recoding has allowed the synthetic generation of different small-genome viruses with modified phenotypes and biological properties. Recently, a decreased cost of chemically synthesizing DNA and improved methods for assembling DNA fragments (e.g. lambda red recombination and CRISPR-based editing) have enabled the construction of an Escherichia coli variant with a 4-Mb synthetic synonymously recoded genome with a reduced number of sense codons (n = 59) encoding the 20 canonical amino acids. Synonymous genome recoding is increasing our knowledge of microbial interactions with innate immune responses, identifying functional genome structures, and strategically ameliorating cis-inhibitory signaling sequences related to splicing, replication (in eukaryotes), and complex microbe functions, unraveling the relevance of codon usage for the temporal regulation of gene expression and the microbe mutant spectrum and adaptability. New biotechnological and therapeutic applications of this methodology can easily be envisaged. In this review, we discuss how synonymous genome recoding may impact our knowledge of microbial biology and the development of new and better therapeutic methodologies.

Synonymous Mutations in rpsT Lead to Ribosomal Assembly Defects That Can Be Compensated by Mutations in fis and rpoA

We previously described how four deleterious synonymous mutations in the Salmonella enterica rpsT gene (encoding ribosomal protein S20) result in low S20 levels that can be compensated by mutations that restore [S20]. Here, we have further studied the cause for the deleterious effects of S20 deficiency and found that the S20 mutants were also deficient in four other 30S proteins (S1, S2, S12, and S21), which is likely due to an assembly defect of the S20 deficient 30S subunits. We examined the compensatory effect by six additional mutations affecting the global regulator Fis and the C-terminal domain of the α subunit of RNA polymerase (encoded by rpoA). The fis and rpoA mutations restored the S20 levels, concomitantly restoring the assembly defect and the levels of S1, S2, S12, and S21. These results illustrate the complexity of compensatory evolution and how the negative effects of deleterious mutations can be suppressed by a multitude of mechanisms. Additionally, we found that the mutations in fis and rpoA caused reduced expression of other ribosomal components. Notably, some of the fis mutations and the rpoA mutation corrected the fitness of the rpsT mutants to wild-type levels, although expression of other ribosomal components was reduced compared to wild-type. This finding raises new questions regarding the relation between translation capacity and growth rate.

Repeated Phenotypic Evolution by Different Genetic Routes in Pseudomonas fluorescens SBW25

Repeated evolution of functionally similar phenotypes is observed throughout the tree of life. The extent to which the underlying genetics are conserved remains an area of considerable interest. Previously, we reported the evolution of colony switching in two independent lineages of Pseudomonas fluorescens SBW25. The phenotypic and genotypic bases of colony switching in the first lineage (Line 1) have been described elsewhere. Here, we deconstruct the evolution of colony switching in the second lineage (Line 6). We show that, as for Line 1, Line 6 colony switching results from an increase in the expression of a colanic acid-like polymer (CAP). At the genetic level, nine mutations occur in Line 6. Only one of these—a nonsynonymous point mutation in the housekeeping sigma factor rpoD—is required for colony switching. In contrast, the genetic basis of colony switching in Line 1 is a mutation in the metabolic gene carB. A molecular model has recently been proposed whereby the carB mutation increases capsulation by redressing the intracellular balance of positive (ribosomes) and negative (RsmAE/CsrA) regulators of a positive feedback loop in capsule expression. We show that Line 6 colony switching is consistent with this model; the rpoD mutation generates an increase in ribosomal gene expression, and ultimately an increase in CAP expression.

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.

Codon usage influences fitness through RNA toxicity

Many organisms are subject to selective pressure that gives rise to unequal usage of synonymous codons, known as codon bias. To experimentally dissect the mechanisms of selection on synonymous sites, we expressed several hundred synonymous variants of the GFP gene in Escherichia coli, and used quantitative growth and viability assays to estimate bacterial fitness. Unexpectedly, we found many synonymous variants whose expression was toxic to E. coli Unlike previously studied effects of synonymous mutations, the effect that we discovered is independent of translation, but it depends on the production of toxic mRNA molecules. We identified RNA sequence determinants of toxicity and evolved suppressor strains that can tolerate the expression of toxic GFP variants. Genome sequencing of these suppressor strains revealed a cluster of promoter mutations that prevented toxicity by reducing mRNA levels. We conclude that translation-independent RNA toxicity is a previously unrecognized obstacle in bacterial gene expression.

The fitness landscape of the codon space across environments

Fitness landscapes map the relationship between genotypes and fitness. However, most fitness landscape studies ignore the genetic architecture imposed by the codon table and thereby neglect the potential role of synonymous mutations. To quantify the fitness effects of synonymous mutations and their potential impact on adaptation on a fitness landscape, we use a new software based on Bayesian Monte Carlo Markov Chain methods and re-estimate selection coefficients of all possible codon mutations across 9 amino acid positions in Saccharomyces cerevisiae Hsp90 across 6 environments. We quantify the distribution of fitness effects of synonymous mutations and show that it is dominated by many mutations of small or no effect and few mutations of larger effect. We then compare the shape of the codon fitness landscape across amino acid positions and environments, and quantify how the consideration of synonymous fitness effects changes the evolutionary dynamics on these fitness landscapes. Together these results highlight a possible role of synonymous mutations in adaptation and indicate the potential mis-inference when they are neglected in fitness landscape studies.

Genetic Adaptation to Growth Under Laboratory Conditions in Escherichia coli and Salmonella enterica

Experimental evolution under controlled laboratory conditions is becoming increasingly important to address various evolutionary questions, including, for example, the dynamics and mechanisms of genetic adaptation to different growth and stress conditions. In such experiments, mutations typically appear that increase the fitness under the conditions tested (medium adaptation), but that are not necessarily of interest for the specific research question. Here, we have identified mutations that appeared during serial passage of E. coli and S. enterica in four different and commonly used laboratory media and measured the relative competitive fitness and maximum growth rate of 111 genetically re-constituted strains, carrying different single and multiple mutations. Little overlap was found between the mutations that were selected in the two species and the different media, implying that adaptation occurs via different genetic pathways. Furthermore, we show that commonly occurring adaptive mutations can generate undesired genetic variation in a population and reduce the accuracy of competition experiments. However, by introducing media adaptation mutations with large effects into the parental strain that was used for the evolution experiment, the variation (standard deviation) was decreased 10-fold, and it was possible to measure fitness differences between two competitors as small as |s| < 0.001.

The Codon Usage of Lowly Expressed Genes Is Subject to Natural Selection

Codon usage bias affects the genomes of organisms from all kingdoms of life and results from both background substitution biases and natural selection. Natural selection on codon usage to increase translation accuracy and efficiency has long been known to affect gene sequences. Such selection is stronger on highly, compared with lowly expressed genes, resulting in higher levels of codon bias within genes with higher expression levels. Additionally, selection on translation accuracy affects more strongly codons encoding conserved amino acids, since these will more often affect protein folding and/or function. By applying tests of selection on the gene sequences of the bacterium Escherichia coli, we demonstrate that both highly and lowly expressed genes display signals of selection on codon usage. Such signals are found for both conserved and less conserved amino acid positions, even within the 10% of E. coli genes expressed at the lowest levels. We further demonstrate experimentally that single synonymous codon replacements within a lowly expressed, essential gene can carry substantial effects on bacterial fitness. Combined, our results demonstrate that even within genes expressed at relatively low levels there is substantial selection on codon usage and that single synonymous codon replacements within such genes can have a marked effect on bacterial fitness.

Direct and Inverted Repeat stimulated excision (DIRex): Simple, single-step, and scar-free mutagenesis of bacterial genes

The need for generating precisely designed mutations is common in genetics, biochemistry, and molecular biology. Here, I describe a new λ Red recombineering method (Direct and Inverted Repeat stimulated excision; DIRex) for fast and easy generation of single point mutations, small insertions or replacements as well as deletions of any size, in bacterial genes. The method does not leave any resistance marker or scar sequence and requires only one transformation to generate a semi-stable intermediate insertion mutant. Spontaneous excision of the intermediate efficiently and accurately generates the final mutant. In addition, the intermediate is transferable between strains by generalized transductions, enabling transfer of the mutation into multiple strains without repeating the recombineering step. Existing methods that can be used to accomplish similar results are either (i) more complicated to design, (ii) more limited in what mutation types can be made, or (iii) require expression of extrinsic factors in addition to λ Red. I demonstrate the utility of the method by generating several deletions, small insertions/replacements, and single nucleotide exchanges in Escherichia coli and Salmonella enterica. Furthermore, the design parameters that influence the excision frequency and the success rate of generating desired point mutations have been examined to determine design guidelines for optimal efficiency.

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.

Duplication-Insertion Recombineering: a fast and scar-free method for efficient transfer of multiple mutations in bacteria

We have developed a new λ Red recombineering methodology for generating transient selection markers that can be used to transfer mutations between bacterial strains of both Escherichia coli and Salmonella enterica. The method is fast, simple and allows for the construction of strains with several mutations without any unwanted sequence changes (scar-free). The method uses λ Red recombineering to generate a marker-held tandem duplication, termed Duplication-Insertion (Dup-In). The Dup-Ins can easily be transferred between strains by generalized transduction and are subsequently rapidly lost by homologous recombination between the two copies of the duplicated sequence, leaving no scar sequence or antibiotic resistance cassette behind. We demonstrate the utility of the method by generating several Dup-Ins in E. coli and S. enterica to transfer genetically linked mutations in both essential and non-essential genes. We have successfully used this methodology to re-construct mutants found after various types of selections, and to introduce foreign genes into the two species. Furthermore, recombineering with two overlapping fragments was as efficient as recombineering with the corresponding single large fragment, allowing more complicated constructions without the need for overlap extension PCR.

Variation in Mutational Robustness between Different Proteins and the Predictability of Fitness Effects

Random mutations in genes from disparate protein classes may have different distributions of fitness effects (DFEs) depending on different structural, functional, and evolutionary constraints. We measured the fitness effects of 156 single mutations in the genes encoding AraC (transcription factor), AraD (enzyme), and AraE (transporter) used for bacterial growth on l-arabinose. Despite their different molecular functions these genes all had bimodal DFEs with most mutations either being neutral or strongly deleterious, providing a general expectation for the DFE. This contrasts with the unimodal DFEs previously obtained for ribosomal protein genes where most mutations were slightly deleterious. Based on theoretical considerations, we suggest that the 33-fold higher average mutational robustness of ribosomal proteins is due to stronger selection for reduced costs of translational and transcriptional errors. Whereas the large majority of synonymous mutations were deleterious for ribosomal proteins genes, no fitness effects could be detected for the AraCDE genes. Four mutations in AraC and AraE increased fitness, suggesting that slightly advantageous mutations make up a significant fraction of the DFE, but that they often escape detection due to the limited sensitivity of commonly used fitness assays. We show that the fitness effects of amino acid substitutions can be predicted based on evolutionary conservation, but those weakly deleterious mutations are less reliably detected. This suggests that large-effect mutations and the fraction of highly deleterious mutations can be computationally predicted, but that experiments are required to characterize the DFE close to neutrality, where many mutations ultimately fixed in a population will occur.

Lateral Gene Transfer Dynamics in the Ancient Bacterial Genus Streptomyces

Lateral gene transfer (LGT) profoundly shapes the evolution of bacterial lineages. LGT across disparate phylogenetic groups and genome content diversity between related organisms suggest a model of bacterial evolution that views LGT as rampant and promiscuous. It has even driven the argument that species concepts and tree-based phylogenetics cannot be applied to bacteria. Here, we show that acquisition and retention of genes through LGT are surprisingly rare in the ubiquitous and biomedically important bacterial genus Streptomyces Using a molecular clock, we estimate that the Streptomyces bacteria are ~380 million years old, indicating that this bacterial genus is as ancient as land vertebrates. Calibrating LGT rate to this geologic time span, we find that on average only 10 genes per million years were acquired and subsequently maintained. Over that same time span, Streptomyces accumulated thousands of point mutations. By explicitly incorporating evolutionary timescale into our analyses, we provide a dramatically different view on the dynamics of LGT and its impact on bacterial evolution.IMPORTANCE Tree-based phylogenetics and the use of species as units of diversity lie at the foundation of modern biology. In bacteria, these pillars of evolutionary theory have been called into question due to the observation of thousands of lateral gene transfer (LGT) events within and between lineages. Here, we show that acquisition and retention of genes through LGT are exceedingly rare in the bacterial genus Streptomyces, with merely one gene acquired in Streptomyces lineages every 100,000 years. These findings stand in contrast to the current assumption of rampant genetic exchange, which has become the dominant hypothesis used to explain bacterial diversity. Our results support a more nuanced understanding of genetic exchange, with LGT impacting evolution over short timescales but playing a significant role over long timescales. Deeper understanding of LGT provides new insight into the evolutionary history of life on Earth, as the vast majority of this history is microbial.

Structural and functional innovations in the real-time evolution of new (βα)8 barrel enzymes

New genes can arise by duplication and divergence, but there is a fundamental gap in our understanding of the relationship between these genes, the evolving proteins they encode, and the fitness of the organism. Here we used crystallography, NMR dynamics, kinetics, and mass spectrometry to explain the molecular innovations that arose during a previous real-time evolution experiment. In that experiment, the (βα)₈ barrel enzyme HisA was under selection for two functions (HisA and TrpF), resulting in duplication and divergence of the hisA gene to encode TrpF specialists, HisA specialists, and bifunctional generalists. We found that selection affects enzyme structure and dynamics, and thus substrate preference, simultaneously and sequentially. Bifunctionality is associated with two distinct sets of loop conformations, each essential for one function. We observed two mechanisms for functional specialization: structural stabilization of each loop conformation and substrate-specific adaptation of the active site. Intracellular enzyme performance, calculated as the product of catalytic efficiency and relative expression level, was not linearly related to fitness. Instead, we observed thresholds for each activity above which further improvements in catalytic efficiency had little if any effect on growth rate. Overall, we have shown how beneficial substitutions selected during real-time evolution can lead to manifold changes in enzyme function and bacterial fitness. This work emphasizes the speed at which adaptive evolution can yield enzymes with sufficiently high activities such that they no longer limit the growth of their host organism, and confirms the (βα)₈ barrel as an inherently evolvable protein scaffold.

RNA Structural Determinants of Optimal Codons Revealed by MAGE-Seq

Synonymous codon choices at the beginning of genes optimize 5' RNA structures for enhanced translation initiation, but less is known about mechanisms that drive codon optimization downstream within the gene. To understand what determines codon choices across a gene, we generated 12,726 in situ codon mutants in the Escherichia coli essential gene infA and measured their fitness by combining multiplex automated genome engineering mutagenesis with amplicon deep sequencing (MAGE-seq). Correlating predicted 5' RNA structure with fitness revealed that codons even far from the start of the gene are deleterious if they disrupt the native 5' RNA conformation. These long-range structural interactions generate context-dependent rules that constrain codon choices beyond intrinsic codon preferences. Genome-wide RNA folding predictions confirm that natural codon choices far from the start codon are optimized in part to prevent disruption of native structures near the 5' UTR. Our results shed light on natural codon distributions and should improve engineering of gene expression for synthetic biology applications.

Unstable tandem gene amplification generates heteroresistance (variation in resistance within a population) to colistin in Salmonella enterica

Heteroresistance, a phenomenon where subpopulations of a bacterial isolate exhibit different susceptibilities to an antibiotic, is a growing clinical problem where the underlying genetic mechanisms in most cases remain unknown. We isolated colistin resistant mutants in Escherichia coli and Salmonella enterica serovar Typhimurium at different concentrations of colistin. Genetic analysis showed that genetically stable pmrAB point mutations were responsible for colistin resistance during selection at high drug concentrations for both species and at low concentrations for E. coli. In contrast, for S. Typhimurium mutants selected at low colistin concentrations, amplification of different large chromosomal regions conferred a heteroresistant phenotype. All amplifications included the pmrD gene, which encodes a positive regulator that up-regulates proteins that modify lipid A, and as a result increase colistin resistance. Inactivation and over-expression of the pmrD gene prevented and conferred resistance, respectively, demonstrating that the PmrD protein is required and sufficient to confer resistance. The heteroresistance phenotype is explained by the variable gene dosage of pmrD in a population, where sub-populations with different copy number of the pmrD gene show different levels of colistin resistance. We propose that variability in gene copy number of resistance genes can explain the heteroresistance observed in clinically isolated pathogenic bacteria.