Evolution of bidirectional costly mutualism from byproduct consumption.

Mutualisms are essential for life, yet it is unclear how they arise. A two-stage process has been proposed for the evolution of mutualisms that involve exchanges of two costly resources. First, costly provisioning by one species may be selected for if that species gains a benefit from costless byproducts generated by a second species, and cooperators get disproportionate access to byproducts. Selection could then drive the second species to provide costly resources in return. Previously, a synthetic consortium evolved the first stage of this scenario: Salmonella enterica evolved costly production of methionine in exchange for costless carbon byproducts generated by an auxotrophic Escherichia coli Growth on agar plates localized the benefits of cooperation around methionine-secreting S. enterica Here, we report that further evolution of these partners on plates led to hypercooperative E. coli that secrete the sugar galactose. Sugar secretion arose repeatedly across replicate communities and is costly to E. coli producers, but enhances the growth of S. enterica The tradeoff between individual costs and group benefits led to maintenance of both cooperative and efficient E. coli genotypes in this spatially structured environment. This study provides an experimental example of de novo, bidirectional costly mutualism evolving from byproduct consumption. The results validate the plausibility of costly cooperation emerging from initially costless exchange, a scenario widely used to explain the origin of the mutualistic species interactions that are central to life on Earth.

Multidrug Resistance Regulators MarA, SoxS, Rob, and RamA Repress Flagellar Gene Expression and Motility in Salmonella enterica Serovar Typhimurium.

Production of flagella is costly and subject to global multilayered regulation, which is reflected in the hierarchical control of flagellar production in many bacterial species. For Salmonella enterica serovar Typhimurium and its relatives, global regulation of flagellar production primarily occurs through the control of flhDC transcription and mRNA translation. In this study, the roles of the homologous multidrug resistance regulators MarA, SoxS, Rob, and RamA (constituting the mar-sox-rob regulon in S Typhimurium) in regulating flagellar gene expression were explored. Each of these regulators was found to inhibit flagellar gene expression, production of flagella, and motility. To different degrees, repression via these transcription factors occurred through direct interactions with the flhDC promoter, particularly for MarA and Rob. Additionally, SoxS repressed flagellar gene expression via a posttranscriptional pathway, reducing flhDC translation. The roles of these transcription factors in reducing motility in the presence of salicylic acid were also elucidated, adding a genetic regulatory element to the response of S Typhimurium to this well-characterized chemorepellent. Integration of flagellar gene expression into the mar-sox-rob regulon in S Typhimurium contrasts with findings for closely related species such as Escherichia coli, providing an example of plasticity in the mar-sox-rob regulon throughout the Enterobacteriaceae family.IMPORTANCE The mar-sox-rob regulon is a large and highly conserved stress response network in the Enterobacteriaceae family. Although it is well characterized in E. coli, the extent of this regulon in related species is unclear. Here, the control of costly flagellar gene expression is connected to the mar-sox-rob regulon of S Typhimurium, contrasting with the E. coli regulon model. These findings demonstrate the flexibility of the mar-sox-rob regulon to accommodate novel regulatory targets, and they provide evidence for its broader regulatory role within this family of diverse bacteria.

Growth Trade-Offs Accompany the Emergence of Glycolytic Metabolism in Shewanella oneidensis MR-1.

Bacteria increase their metabolic capacity via the acquisition of genetic material or by the mutation of genes already present in the genome. Here, we explore the mechanisms and trade-offs involved when Shewanella oneidensis, a bacterium that typically consumes small organic and amino acids, rapidly evolves to expand its metabolic capacity to catabolize glucose after a short period of adaptation to a glucose-rich environment. Using whole-genome sequencing and genetic approaches, we discovered that deletions in a region including the transcriptional repressor (nagR) that regulates the expression of genes associated with catabolism of N-acetylglucosamine are the common basis for evolved glucose metabolism across populations. The loss of nagR results in the constitutive expression of genes for an N-acetylglucosamine permease (nagP) and kinase (nagK). We demonstrate that promiscuous activities of both NagP and NagK toward glucose allow for the transport and phosphorylation of glucose to glucose-6-phosphate, the initial events of glycolysis otherwise thought to be absent in S. oneidensis13C-based metabolic flux analysis uncovered that subsequent utilization was mediated by the Entner-Doudoroff pathway. This is an example whereby gene loss and preexisting enzymatic promiscuity, and not gain-of-function mutations, were the drivers of increased metabolic capacity. However, we observed a significant decrease in the growth rate on lactate after adaptation to glucose catabolism, suggesting that trade-offs may explain why glycolytic function may not be readily observed in S. oneidensis in natural environments despite it being readily accessible through just a single mutational event.IMPORTANCE Gains in metabolic capacity are frequently associated with the acquisition of novel genetic material via natural or engineered horizontal gene transfer events. Here, we explored how a bacterium that typically consumes small organic acids and amino acids expands its metabolic capacity to include glucose via a loss of genetic material, a process frequently associated with a deterioration of metabolic function. Our findings highlight how the natural promiscuity of transporters and enzymes can be a key driver in expanding metabolic diversity and that many bacteria may possess a latent metabolic capacity accessible through one or a few mutations that remove regulatory functions. Our discovery of trade-offs between growth on lactate and on glucose suggests why this easily gained trait is not observed in nature.

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.

The Mar, Sox, and Rob Systems.

Environments inhabited by Enterobacteriaceae are diverse and often stressful. This is particularly true for Escherichia coli and Salmonella during host association in the gastrointestinal systems of animals. There, E. coli and Salmonella must survive exposure to various antimicrobial compounds produced or ingested by their host. A myriad of changes to cellular physiology and metabolism are required to achieve this feat. A central regulatory network responsible for sensing and responding to intracellular chemical stressors like antibiotics are the Mar, Sox, and Rob systems found throughout the Enterobacteriaceae. Each of these distinct regulatory networks controls expression of an overlapping set of downstream genes whose collective effects result in increased resistance to a wide array of antimicrobial compounds. This collection of genes is known as the mar-sox-rob regulon. This review will provide an overview of the mar-sox-rob regulon and molecular architecture of the Mar, Sox, and Rob systems.

Composition and function of the Galapagos penguin gut microbiome vary with age, location, and a putative bacterial pathogen.

Microbial colonization plays a direct role in host health. Understanding the ecology of the resident microbial community for a given host species is thus an important step for detecting population vulnerabilities like disease. However, the idea of integrating microbiome research into conservation is still relatively new, and wild birds have received less attention in this field than mammals or domesticated animals. Here we examine the composition and function of the gut microbiome of the endangered Galapagos penguin (Spheniscus mendiculus) with the goals of characterizing the normal microbial community and resistome, identifying likely pathogens, and testing hypotheses of structuring forces for this community based on demographics, location, and infection status. We collected fecal samples from wild penguins in 2018 and performed 16S rRNA gene sequencing and whole genome sequencing (WGS) on extracted DNA. 16S sequencing revealed that the bacterial phyla Fusobacteria, Epsilonbacteraeota, Firmicutes, and Proteobacteria dominate the community. Functional pathways were computed from WGS data, showing genetic functional potential primarily focused on metabolism-amino acid metabolism, carbohydrate metabolism, and energy metabolism are the most well-represented functional groups. WGS samples were each screened for antimicrobial resistance, characterizing a resistome made up of nine antibiotic resistance genes. Samples were screened for potential enteric pathogens using virulence factors as indicators; Clostridium perfringens was revealed as a likely pathogen. Overall, three factors appear to be shaping the alpha and beta diversity of the microbial community: penguin developmental stage, sampling location, and C. perfringens. We found that juvenile penguins have significantly lower alpha diversity than adults based on three metrics, as well as significantly different beta diversity. Location effects are minimal, but one site has significantly lower Shannon diversity than the other primary sites. Finally, when samples were grouped by C. perfringens virulence factors, we found dramatic changes in beta diversity based on operational taxonomic units, protein families, and functional pathways. This study provides a baseline microbiome for an endangered species, implicates both penguin age and the presence of a potential bacterial pathogen as primary factors associated with microbial community variance, and reveals widespread antibiotic resistance genes across the population.

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 other's 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.

High-resolution, long-term characterization of bacterial motility using optical tweezers.

We present a single-cell motility assay, which allows the quantification of bacterial swimming in a well-controlled environment, for durations of up to an hour and with a temporal resolution greater than the flagellar rotation rates of approximately 100 Hz. The assay is based on an instrument combining optical tweezers, light and fluorescence microscopy, and a microfluidic chamber. Using this device we characterized the long-term statistics of the run-tumble time series in individual Escherichia coli cells. We also quantified higher-order features of bacterial swimming, such as changes in velocity and reversals of swimming direction.

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.

Genetic Context Significantly Influences the Maintenance and Evolution of Degenerate Pathways.

Understanding the evolution of novel physiological traits is highly relevant for expanding the characterization and manipulation of biological systems. Acquisition of new traits can be achieved through horizontal gene transfer (HGT). Here, we investigate drivers that promote or deter the maintenance of HGT-driven degeneracy, occurring when processes accomplish identical functions through nonidentical components. Subsequent evolution can optimize newly acquired functions; for example, beneficial alleles identified in an engineered Methylorubrum extorquens strain allowed it to utilize a "Foreign" formaldehyde oxidation pathway substituted for its Native pathway for methylotrophic growth. We examined the fitness consequences of interactions between these alleles when they were combined with the Native pathway or both (Dual) pathways. Unlike the Foreign pathway context where they evolved, these alleles were often neutral or deleterious when moved into these alternative genetic backgrounds. However, there were instances where combinations of multiple alleles resulted in higher fitness outcomes than individual allelic substitutions could provide. Importantly, the genetic context accompanying these allelic substitutions significantly altered the fitness landscape, shifting local fitness peaks and restricting the set of accessible evolutionary trajectories. These findings highlight how genetic context can negatively impact the probability of maintaining native and HGT-introduced functions together, making it difficult for degeneracy to evolve. However, in cases where the cost of maintaining degeneracy was mitigated by adding evolved alleles impacting the function of these pathways, we observed rare opportunities for pathway coevolution to occur. Together, our results highlight the importance of genetic context and resulting epistasis in retaining or losing HGT-acquired degenerate functions.

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.

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.

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.

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.

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.

Species interactions differ in their genetic robustness.

Conflict and cooperation between bacterial species drive the composition and function of microbial communities. Stability of these emergent properties will be influenced by the degree to which species' interactions are robust to genetic perturbations. We use genome-scale metabolic modeling to computationally analyze the impact of genetic changes when Escherichia coli and Salmonella enterica compete, or cooperate. We systematically knocked out in silico each reaction in the metabolic network of E. coli to construct all 2583 mutant stoichiometric models. Then, using a recently developed multi-scale computational framework, we simulated the growth of each mutant E. coli in the presence of S. enterica. The type of interaction between species was set by modulating the initial metabolites present in the environment. We found that the community was most robust to genetic perturbations when the organisms were cooperating. Species ratios were more stable in the cooperative community, and community biomass had equal variance in the two contexts. Additionally, the number of mutations that have a substantial effect is lower when the species cooperate than when they are competing. In contrast, when mutations were added to the S. enterica network the system was more robust when the bacteria were competing. These results highlight the utility of connecting metabolic mechanisms and studies of ecological stability. Cooperation and conflict alter the connection between genetic changes and properties that emerge at higher levels of biological organization.

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

BACKGROUND: Due 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. FINDINGS: We 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. CONCLUSIONS: We 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.

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 >10(5) 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.