Lon M. Chubiz

FREQ-Seq: A Rapid, Cost-Effective, Sequencing-Based Method to Determine Allele Frequencies Directly from Mixed Populations

Understanding evolutionary dynamics within microbial populations requires the ability to accurately follow allele frequencies through time. Here we present a rapid, cost-effective method (FREQ-Seq) that leverages Illumina next-generation sequencing for localized, quantitative allele frequency detection. Analogous to RNA-Seq, FREQ-Seq relies upon counts from the >105 reads generated per locus per time-point to determine allele frequencies. Loci of interest are directly amplified from a mixed population via two rounds of PCR using inexpensive, user-designed oligonucleotides and a bar-coded bridging primer system that can be regenerated in-house. The resulting bar-coded PCR products contain the adapters needed for Illumina sequencing, eliminating further library preparation. We demonstrate the utility of FREQ-Seq by determining the order and dynamics of beneficial alleles that arose as a microbial population, founded with an engineered strain of Methylobacterium, evolved to grow on methanol. Quantifying allele frequencies with minimal bias down to 1% abundance allowed effective analysis of SNPs, small in-dels and insertions of transposable elements. Our data reveal large-scale clonal interference during the early stages of adaptation and illustrate the utility of FREQ-Seq as a cost-effective tool for tracking allele frequencies in populations.

Aromatic Acid Metabolites of Escherichia coli K-12 Can Induce the marRAB Operon

MarR is a key regulator of the marRAB operon involved in antibiotic resistance and solvent stress tolerance in Escherichia coli. We show that two metabolic intermediates, 2,3-dihydroxybenzoate and anthranilate, involved in enterobactin and tryptophan biosynthesis, respectively, can activate marRAB transcription. We also found that a third intermediate involved in ubiquinone biosynthesis, 4-hydroxybenzoate, activates marRAB transcription in the absence of TolC. Of the three, however, only 2,3-dihydroxybenzoate directly binds MarR and affects its activity.

Large-Effect Beneficial Synonymous Mutations Mediate Rapid and Parallel Adaptation in a Bacterium

Contrary to previous understanding, recent evidence indicates that synonymous codon changes may sometimes face strong selection. However, it remains difficult to generalize the nature, strength, and mechanism(s) of such selection. Previously, we showed that synonymous variants of a key enzyme-coding gene (fae) of Methylobacterium extorquens AM1 decreased enzyme production and reduced fitness dramatically. We now show that during laboratory evolution, these variants rapidly regained fitness via parallel yet variant-specific, highly beneficial point mutations in the N-terminal region of fae. These mutations (including four synonymous mutations) had weak but consistently positive impacts on transcript levels, enzyme production, or enzyme activity. However, none of the proposed mechanisms (including internal ribosome pause sites or mRNA structure) predicted the fitness impact of evolved or additional, engineered point mutations. This study shows that synonymous mutations can be fixed through strong positive selection, but the mechanism for their benefit varies depending on the local sequence context.

Computational design of orthogonal ribosomes

Orthogonal ribosomes (o-ribosomes), also known as specialized ribosomes, are able to selectively translate mRNA not recognized by host ribosomes. As a result, they are powerful tools for investigating translational regulation and probing ribosome structure. To date, efforts directed towards engineering o-ribosomes have involved random mutagenesis-based approaches. As an alternative, we present here a computational method for rationally designing o-ribosomes in bacteria. Working under the assumption that base-pair interactions between the 16S rRNA and mRNA serve as the primary mode for ribosome binding and translational initiation, the algorithm enumerates all possible extended recognition sequences for 16S rRNA and then chooses those candidates that: (i) have a similar binding strength to their target mRNA as the canonical, wild-type ribosome/mRNA pair; (ii) do not bind mRNA with the wild-type, canonical Shine-Dalgarno (SD) sequence and (iii) minimally interact with host mRNA irrespective of whether a recognizable SD sequence is present. In order to test the algorithm, we experimentally characterized a number of computationally designed o-ribosomes in Escherichia coli.

Transcriptional Cross Talk within the mar-sox-rob Regulon in Escherichia coli Is Limited to the rob and marRAB Operons

Bacteria possess multiple mechanisms to survive exposure to various chemical stresses and antimicrobial compounds. In the enteric bacterium Escherichia coli, three homologous transcription factors-MarA, SoxS, and Rob-play a central role in coordinating this response. Three separate systems are known to regulate the expression and activities of MarA, SoxS, and Rob. However, a number of studies have shown that the three do not function in isolation but rather are coregulated through transcriptional cross talk. In this work, we systematically investigated the extent of transcriptional cross talk in the mar-sox-rob regulon. While the three transcription factors were found to have the potential to regulate each others expression when ectopically expressed, the only significant interactions observed under physiological conditions were between mar and rob systems. MarA, SoxS, and Rob all activate the marRAB promoter, more so when they are induced by their respective inducers: salicylate, paraquat, and decanoate. None of the three proteins affects the soxS promoter, though unexpectedly, it was mildly repressed by decanoate by an unknown mechanism. SoxS is the only one of the three proteins to repress the rob promoter. Surprisingly, salicylate somewhat activates transcription of rob, while decanoate represses it a bit. Rob, in turn, activates not only its downstream promoters in response to salicylate but also the marRAB promoter. These results demonstrate that the mar and rob systems function together in response to salicylate.

Role of the mar-sox-rob Regulon in Regulating Outer Membrane Porin Expression

Multiple factors control the expression of the outer membrane porins OmpF and OmpC in Escherichia coli. In this work, we investigated the role of the mar-sox-rob regulon in regulating outer membrane porin expression in response to salicylate. We provide both genetic and physiological evidence that MarA and Rob can independently activate micF transcription in response to salicylate, leading to reduced OmpF expression. MarA was also found to repress OmpF expression through a MicF-independent pathway. In the case of OmpC, we found that its transcription was moderately increased in response to salicylate. However, this increase was independent of MarA and Rob. Finally, we found that the reduction in OmpF expression in a tolC mutant is due primarily to Rob. Collectively, this work further clarifies the coordinated role of MarA and Rob in regulating the expression of the outer membrane porins.

A novel pair of inducible expression vectors for use in Methylobacterium extorquens

BackgroundDue to the ever increasing use of diverse microbial taxa in basic research and industrial settings, there is a growing need for genetic tools to alter the physiology of these organisms. In particular, there is a dearth of inducible expression systems available for bacteria outside commonly used γ-proteobacteria, such as Escherichia coli or Pseudomonas species. To this end, we have sought to develop a pair of inducible expression vectors for use in the α-proteobacterium Methylobacterium extorquens, a model methylotroph.FindingsWe found that the P R promoter from rhizobial phage 16-3 was active in M. extorquens and engineered the promoter to be inducible by either p-isopropyl benzoate (cumate) or anhydrotetracycline. These hybrid promoters, P R/cmtO and P R/tetO, were found to have high levels of expression in M. extorquens with a regulatory range of 10-fold and 30-fold, respectively. Compared to an existing cumate-inducible (10-fold range), high-level expression system for M. extorquens, P R/cmtO and P R/tetO have 33% of the maximal activity but were able to repress gene expression 3 and 8-fold greater, respectively. Both promoters were observed to exhibit homogeneous, titratable activation dynamics rather than on-off, switch-like behavior. The utility of these promoters was further demonstrated by complementing loss of function of ftfL - essential for growth on methanol - where we show P R/tetO is capable of not only fully complementing function but also producing a conditional null phenotype. These promoters have been incorporated into a broad-host-range backbone allowing for potential use in a variety of bacterial hosts.ConclusionsWe have developed two novel expression systems for use in M. extorquens. The expression range of these vectors should allow for increased ability to explore cellular physiology in M. extorquens. Further, the P R/tetO promoter is capable of producing conditional null phenotypes, previously unattainable in M. extorquens. As both expression systems rely on the use of membrane permeable inducers, we suspect these expression vectors will be useful for ectopic gene expression in numerous proteobacteria.

Parallel and Divergent Evolutionary Solutions for the Optimization of an Engineered Central Metabolism in Methylobacterium extorquens AM1

Bioengineering holds great promise to provide fast and efficient biocatalysts for methanol-based biotechnology, but necessitates proven methods to optimize physiology in engineered strains. Here, we highlight experimental evolution as an effective means for optimizing an engineered Methylobacterium extorquens AM1. Replacement of the native formaldehyde oxidation pathway with a functional analog substantially decreased growth in an engineered Methylobacterium, but growth rapidly recovered after six hundred generations of evolution on methanol. We used whole-genome sequencing to identify the basis of adaptation in eight replicate evolved strains, and examined genomic changes in light of other growth and physiological data. We observed great variety in the numbers and types of mutations that occurred, including instances of parallel mutations at targets that may have been “rationalized” by the bioengineer, plus other “illogical” mutations that demonstrate the ability of evolution to expose unforeseen optimization solutions. Notably, we investigated mutations to RNA polymerase, which provided a massive growth benefit but are linked to highly aberrant transcriptional profiles. Overall, we highlight the power of experimental evolution to present genetic and physiological solutions for strain optimization, particularly in systems where the challenges of engineering are too many or too difficult to overcome via traditional engineering methods.

Parallel Mutations Result in a Wide Range of Cooperation and Community Consequences in a Two-Species Bacterial Consortium

Multi-species microbial communities play a critical role in human health, industry, and waste remediation. Recently, the evolution of synthetic consortia in the laboratory has enabled adaptation to be addressed in the context of interacting species. Using an engineered bacterial consortium, we repeatedly evolved cooperative genotypes and examined both the predictability of evolution and the phenotypes that determine community dynamics. Eight Salmonella enterica serovar Typhimurium strains evolved methionine excretion sufficient to support growth of an Escherichia coli methionine auxotroph, from whom they required excreted growth substrates. Non-synonymous mutations in metA, encoding homoserine trans-succinylase (HTS), were detected in each evolved S. enterica methionine cooperator and were shown to be necessary for cooperative consortia growth. Molecular modeling was used to predict that most of the non-synonymous mutations slightly increase the binding affinity for HTS homodimer formation. Despite this genetic parallelism and trend of increasing protein binding stability, these metA alleles gave rise to a wide range of phenotypic diversity in terms of individual versus group benefit. The cooperators with the highest methionine excretion permitted nearly two-fold faster consortia growth and supported the highest fraction of E. coli, yet also had the slowest individual growth rates compared to less cooperative strains. Thus, although the genetic basis of adaptation was quite similar across independent origins of cooperative phenotypes, quantitative measurements of metabolite production were required to predict either the individual-level growth consequences or how these propagate to community-level behavior.

Identification of the potentiating mutations and synergistic epistasis that enabled the evolution of inter-species cooperation

Microbes often engage in cooperation through releasing biosynthetic compounds required by other species to grow. Given that production of costly biosynthetic metabolites is generally subjected to multiple layers of negative feedback, single mutations may frequently be insufficient to generate cooperative phenotypes. Synergistic epistatic interactions between multiple coordinated changes may thus often underlie the evolution of cooperation through overproduction of metabolites. To test the importance of synergistic mutations in cooperation we used an engineered bacterial consortium of an Escherichia coli methionine auxotroph and Salmonella enterica. S. enterica relies on carbon by-products from E. coli if lactose is the only carbon source. Directly selecting wild-type S. enterica in an environment that favored cooperation through secretion of methionine only once led to a methionine producer, and this producer both took a long time to emerge and was not very effective at cooperating. On the other hand, when an initial selection for resistance of S. enterica to a toxic methionine analog, ethionine, was used, subsequent selection for cooperation with E. coli was rapid, and the resulting double mutants were much more effective at cooperation. We found that potentiating mutations in metJ increase expression of metA, which encodes the first step of methionine biosynthesis. This increase in expression is required for the previously identified actualizing mutations in metA to generate cooperation. This work highlights that where biosynthesis of metabolites involves multiple layers of regulation, significant secretion of those metabolites may require multiple mutations, thereby constraining the evolution of cooperation.