Global Rebalancing of Cellular Resources by Pleiotropic Point Mutations Illustrates a Multi-scale Mechanism of Adaptive Evolution

The hallmarks of a tradeoff in transcriptomes that balances stress and growth functions

Fast growth phenotypes are achieved through optimal transcriptomic allocation, in which cells must balance tradeoffs in resource allocation between diverse functions. One such balance between stress readiness and unbridled growth in E. coli has been termed the fear versus greed (f/g) tradeoff. Two specific RNA polymerase (RNAP) mutations observed in adaptation to fast growth have been previously shown to affect the f/g tradeoff, suggesting that genetic adaptations may be primed to control f/g resource allocation. Here, we conduct a greatly expanded study of the genetic control of the f/g tradeoff across diverse conditions. We introduced 12 RNA polymerase (RNAP) mutations commonly acquired during adaptive laboratory evolution (ALE) and obtained expression profiles of each. We found that these single RNAP mutation strains resulted in large shifts in the f/g tradeoff primarily in the RpoS regulon and ribosomal genes, likely through modifying RNAP-DNA interactions. Two of these mutations additionally caused condition-specific transcriptional adaptations. While this tradeoff was previously characterized by the RpoS regulon and ribosomal expression, we find that the GAD regulon plays an important role in stress readiness and ppGpp in translation activity, expanding the scope of the tradeoff. A phylogenetic analysis found the greed-related genes of the tradeoff present in numerous bacterial species. The results suggest that the f/g tradeoff represents a general principle of transcriptome allocation in bacteria where small genetic changes can result in large phenotypic adaptations to growth conditions.IMPORTANCETo increase growth, E. coli must raise ribosomal content at the expense of non-growth functions. Previous studies have linked RNAP mutations to this transcriptional shift and increased growth but were focused on only two mutations found in the protein's central region. RNAP mutations, however, commonly occur over a large structural range. To explore RNAP mutations' impact, we have introduced 12 RNAP mutations found in laboratory evolution experiments and obtained expression profiles of each. The mutations nearly universally increased growth rates by adjusting said tradeoff away from non-growth functions. In addition to this shift, a few caused condition-specific adaptations. We explored the prevalence of this tradeoff across phylogeny and found it to be a widespread and conserved trend among bacteria.

Base editor-mediated large-scale screening of functional mutations in bacteria for industrial phenotypes

Base editing, the targeted introduction of point mutations into cellular DNA, holds promise for improving genome-scale functional genome screening to single-nucleotide resolution. Current efforts in prokaryotes, however, remain confined to loss-of-function screens using the premature stop codons-mediated gene inactivation library, which falls far short of fully releasing the potential of base editors. Here, we developed a base editor-mediated functional single nucleotide variant screening pipeline in Escherichia coli. We constructed a library with 31,123 sgRNAs targeting 462 stress response-related genes in E. coli, and screened for adaptive mutations under isobutanol and furfural selective conditions. Guided by the screening results, we successfully identified several known and novel functional mutations. Our pipeline might be expanded to the optimization of other phenotypes or the strain engineering in other microorganisms.

Serial adaptive laboratory evolution enhances mixed carbon metabolic capacity of Escherichia coli

Microbes have inherent capacities for utilizing various carbon sources, however they often exhibit sub-par fitness due to low metabolic efficiency. To test whether a bacterial strain can optimally utilize multiple carbon sources, Escherichia coli was serially evolved in L-lactate and glycerol. This yielded two end-point strains that evolved first in L-lactate then in glycerol, and vice versa. The end-point strains displayed a universal growth advantage on single and a mixture of adaptive carbon sources, enabled by a concerted action of carbon source-specialists and generalist mutants. The combination of just four variants of glpK, ppsA, ydcI, and rph-pyrE, accounted for more than 80% of end-point strain fitness. In addition, machine learning analysis revealed a coordinated activity of transcriptional regulators imparting condition-specific regulation of gene expression. The effectiveness of the serial adaptive laboratory evolution (ALE) scheme in bioproduction applications was assessed under single and mixed-carbon culture conditions, in which serial ALE strain exhibited superior productivity of acetoin compared to ancestral strains. Together, systems-level analysis elucidated the molecular basis of serial evolution, which hold potential utility in bioproduction applications.

Convergent genetic adaptation of Escherichia coli in minimal media leads to pleiotropic divergence

Adaptation in an environment can either be beneficial, neutral or disadvantageous in another. To test the genetic basis of pleiotropic behaviour, we evolved six lines of E. coli independently in environments where glucose and galactose were the sole carbon sources, for 300 generations. All six lines in each environment exhibit convergent adaptation in the environment in which they were evolved. However, pleiotropic behaviour was observed in several environmental contexts, including other carbon environments. Genome sequencing reveals that mutations in global regulators rpoB and rpoC cause this pleiotropy. We report three new alleles of the rpoB gene, and one new allele of the rpoC gene. The novel rpoB alleles confer resistance to Rifampicin, and alter motility. Our results show how single nucleotide changes in the process of adaptation in minimal media can lead to wide-scale pleiotropy, resulting in changes in traits that are not under direct selection.

Specific codons control cellular resources and fitness

As cellular engineering progresses from simply overexpressing proteins to imparting complex phenotypes through multigene expression, judicious appropriation of cellular resources is essential. Since codon use is degenerate and biased, codons may control cellular resources at a translational level. We investigate how partitioning transfer RNA (tRNA) resources by incorporating dissimilar codon usage can drastically alter interdependence of expression level and burden on the host. By isolating the effect of individual codons' use during translation elongation while eliminating confounding factors, we show that codon choice can trans-regulate fitness of the host and expression of other heterologous or native genes. We correlate specific codon usage patterns with host fitness and derive a coding scheme for multigene expression called the Codon Health Index (CHI, χ). This empirically derived coding scheme (χ) enables the design of multigene expression systems that avoid catastrophic cellular burden and is robust across several proteins and conditions.

Autotrophic growth of Escherichia coli is achieved by a small number of genetic changes

Synthetic autotrophy is a promising avenue to sustainable bioproduction from CO2. Here, we use iterative laboratory evolution to generate several distinct autotrophic strains. Utilising this genetic diversity, we identify that just three mutations are sufficient for Escherichia coli to grow autotrophically, when introduced alongside non-native energy (formate dehydrogenase) and carbon-fixing (RuBisCO, phosphoribulokinase, carbonic anhydrase) modules. The mutated genes are involved in glycolysis (pgi), central-carbon regulation (crp), and RNA transcription (rpoB). The pgi mutation reduces the enzyme's activity, thereby stabilising the carbon-fixing cycle by capping a major branching flux. For the other two mutations, we observe down-regulation of several metabolic pathways and increased expression of native genes associated with the carbon-fixing module (rpiB) and the energy module (fdoGH), as well as an increased ratio of NADH/NAD+ - the cycle's electron-donor. This study demonstrates the malleability of metabolism and its capacity to switch trophic modes using only a small number of genetic changes and could facilitate transforming other heterotrophic organisms into autotrophs.

A multi-scale expression and regulation knowledge base for Escherichia coli

Transcriptomic data is accumulating rapidly; thus, scalable methods for extracting knowledge from this data are critical. Here, we assembled a top-down expression and regulation knowledge base for Escherichia coli. The expression component is a 1035-sample, high-quality RNA-seq compendium consisting of data generated in our lab using a single experimental protocol. The compendium contains diverse growth conditions, including: 9 media; 39 supplements, including antibiotics; 42 heterologous proteins; and 76 gene knockouts. Using this resource, we elucidated global expression patterns. We used machine learning to extract 201 modules that account for 86% of known regulatory interactions, creating the regulatory component. With these modules, we identified two novel regulons and quantified systems-level regulatory responses. We also integrated 1675 curated, publicly-available transcriptomes into the resource. We demonstrated workflows for analyzing new data against this knowledge base via deconstruction of regulation during aerobic transition. This resource illuminates the E. coli transcriptome at scale and provides a blueprint for top-down transcriptomic analysis of non-model organisms.

Deep mutational scanning reveals the molecular determinants of RNA polymerase-mediated adaptation and tradeoffs

RNA polymerase (RNAP) is emblematic of complex biological systems that control multiple traits involving trade-offs such as growth versus maintenance. Laboratory evolution has revealed that mutations in RNAP subunits, including RpoB, are frequently selected. However, we lack a systems view of how mutations alter the RNAP molecular functions to promote adaptation. We, therefore, measured the fitness of thousands of mutations within a region of rpoB under multiple conditions and genetic backgrounds, to find that adaptive mutations cluster in two modules. Mutations in one module favor growth over maintenance through a partial loss of an interaction associated with faster elongation. Mutations in the other favor maintenance over growth through a destabilized RNAP-DNA complex. The two molecular handles capture the versatile RNAP-mediated adaptations. Combining both interaction losses simultaneously improved maintenance and growth, challenging the idea that growth-maintenance tradeoff resorts only from limited resources, and revealing how compensatory evolution operates within RNAP.

Microbial engineering for shikimate biosynthesis

Shikimate, a precursor to the antiviral drug oseltamivir (Tamiflu®), can influence aromatic metabolites and finds extensive use in antimicrobial, antitumor, and cardiovascular applications. Consequently, various strategies have been developed for chemical synthesis and plant extraction to enhance shikimate biosynthesis, potentially impacting environmental conditions, economic sustainability, and separation and purification processes. Microbial engineering has been developed as an environmentally friendly approach for shikimate biosynthesis. In this review, we provide a comprehensive summary of microbial strategies for shikimate biosynthesis. These strategies primarily include chassis construction, biochemical optimization, pathway remodelling, and global regulation. Furthermore, we discuss future perspectives on shikimate biosynthesis and emphasize the importance of utilizing advanced metabolic engineering tools to regulate microbial networks for constructing robust microbial cell factories.

Trade-offs between the instantaneous growth rate and long-term fitness: Consequences for microbial physiology and predictive computational models

Microbial systems biology has made enormous advances in relating microbial physiology to the underlying biochemistry and molecular biology. By meticulously studying model microorganisms, in particular Escherichia coli and Saccharomyces cerevisiae, increasingly comprehensive computational models predict metabolic fluxes, protein expression, and growth. The modeling rationale is that cells are constrained by a limited pool of resources that they allocate optimally to maximize fitness. As a consequence, the expression of particular proteins is at the expense of others, causing trade-offs between cellular objectives such as instantaneous growth, stress tolerance, and capacity to adapt to new environments. While current computational models are remarkably predictive for E. coli and S. cerevisiae when grown in laboratory environments, this may not hold for other growth conditions and other microorganisms. In this contribution, we therefore discuss the relationship between the instantaneous growth rate, limited resources, and long-term fitness. We discuss uses and limitations of current computational models, in particular for rapidly changing and adverse environments, and propose to classify microbial growth strategies based on Grimes's CSR framework.

Development of heat-shock resistance in Legionella pneumophila modeled by experimental evolution

Because it can grow in buildings with complex hot water distribution systems (HWDS), healthcare facilities recognize the waterborne bacterium Legionella pneumophila as a major nosocomial infection threat and often try to clear the systems with a pasteurization process known as superheat-and-flush. After this treatment, many facilities find that the contaminating populations slowly recover, suggesting the possibility of in situ evolution favoring increased survival in high-temperature conditions. To mimic this process in a controlled environment, an adaptive laboratory evolution model was used to select a wild-type strain of L. pneumophila for survival to transient exposures to temperatures characteristic of routine hot water use or failed pasteurization processes in HWDS. Over their evolution, these populations became insensitive to exposure to 55°C and developed the ability to survive short exposures to 59°C heat shock. Heat-adapted lineages maintained a higher expression of heat-shock genes during low-temperature incubation in freshwater, suggesting a pre-adaptation to heat stress. Although there were distinct mutation profiles in each of the heat-adapted lineages, each acquired multiple mutations in the DnaJ/DnaK/ClpB disaggregase complex, as well as mutations in chaperone htpG and protease clpX. These mutations were specific to heat-shock survival and were not seen in control lineages included in the experimental model without exposure to heat shock. This study supports in situ observations of adaptation to heat stress and demonstrates the potential of L. pneumophila to develop resistance to control measures. IMPORTANCE As a bacterium that thrives in warm water ecosystems, Legionella pneumophila is a key factor motivating regulations on hot water systems. Two major measures to control Legionella are high circulating temperatures intended to curtail growth and the use of superheat-and-flush pasteurization processes to eliminate established populations. Facilities often suffer recolonization of their hot water systems; hospitals are particularly at risk due to the severe nosocomial pneumoniae caused by Legionella. To understand these long-term survivors, we have used an adaptive laboratory evolution model to replicate this process. We find major differences between the mutational profiles of heat-adapted and heat-naïve L. pneumophila populations including mutations in major heat-shock genes like chaperones and proteases. This model demonstrates that well-validated treatment protocols are needed to clear contaminated systems and-in an analog to antibiotic resistance-the importance of complete eradication of the resident population to prevent selection for more persistent bacteria.

Laboratory evolution, transcriptomics, and modeling reveal mechanisms of paraquat tolerance

Relationships between the genome, transcriptome, and metabolome underlie all evolved phenotypes. However, it has proved difficult to elucidate these relationships because of the high number of variables measured. A recently developed data analytic method for characterizing the transcriptome can simplify interpretation by grouping genes into independently modulated sets (iModulons). Here, we demonstrate how iModulons reveal deep understanding of the effects of causal mutations and metabolic rewiring. We use adaptive laboratory evolution to generate E. coli strains that tolerate high levels of the redox cycling compound paraquat, which produces reactive oxygen species (ROS). We combine resequencing, iModulons, and metabolic models to elucidate six interacting stress-tolerance mechanisms: (1) modification of transport, (2) activation of ROS stress responses, (3) use of ROS-sensitive iron regulation, (4) motility, (5) broad transcriptional reallocation toward growth, and (6) metabolic rewiring to decrease NADH production. This work thus demonstrates the power of iModulon knowledge mapping for evolution analysis.

Pervasive genotype-by-environment interactions shape the fitness effects of antibiotic resistance mutations

The fitness effects of antibiotic resistance mutations are a major driver of resistance evolution. While the nutrient environment affects bacterial fitness, experimental studies of resistance typically measure fitness of mutants in a single environment only. We explored how the nutrient environment affected the fitness effects of rifampicin-resistant rpoB mutations in Escherichia coli under several conditions critical for the emergence and spread of resistance-the presence of primary or secondary antibiotic, or the absence of any antibiotic. Pervasive genotype-by-environment (GxE) interactions determined fitness in all experimental conditions, with rank order of fitness in the presence and absence of antibiotics being strongly dependent on the nutrient environment. GxE interactions also affected the magnitude and direction of collateral effects of secondary antibiotics, in some cases so drastically that a mutant that was highly sensitive in one nutrient environment exhibited cross-resistance to the same antibiotic in another. It is likely that the mutant-specific impact of rpoB mutations on the global transcriptome underpins the observed GxE interactions. The pervasive, mutant-specific GxE interactions highlight the importance of doing what is rarely done when studying the evolution and spread of resistance in experimental and clinical work: assessing fitness of antibiotic-resistant mutants across a range of relevant environments.

Reprogramming cellular metabolism to increase the efficiency of microbial cell factories

Recent studies are increasingly focusing on advanced biotechnological tools, self-adjusting smart microorganisms, and artificial intelligent networks, to engineer microorganisms with various functions. Microbial cell factories are a vital platform for improving the bioproduction of medicines, biofuels, and biomaterials from renewable carbon sources. However, these processes are significantly affected by cellular metabolism, and boosting the efficiency of microbial cell factories remains a challenge. In this review, we present a strategy for reprogramming cellular metabolism to enhance the efficiency of microbial cell factories for chemical biosynthesis, which improves our understanding of microbial physiology and metabolic control. Current methods are mainly focused on synthetic pathways, metabolic resources, and cell performance. This review highlights the potential biotechnological strategy to reprogram cellular metabolism and provide novel guidance for designing more intelligent industrial microbes with broader applications in this growing field.

The hallmarks of a tradeoff in transcriptomes that balances stress and growth functions

Fit phenotypes are achieved through optimal transcriptomic allocation. Here, we performed a high-resolution, multi-scale study of the transcriptomic tradeoff between two key fitness phenotypes, stress response (fear) and growth (greed), in Escherichia coli. We introduced twelve RNA polymerase (RNAP) mutations commonly acquired during adaptive laboratory evolution (ALE) and found that single mutations resulted in large shifts in the fear vs. greed tradeoff, likely through destabilizing the rpoB-rpoC interface. RpoS and GAD regulons drive the fear response while ribosomal proteins and the ppGpp regulon underlie greed. Growth rate selection pressure during ALE results in endpoint strains that often have RNAP mutations, with synergistic mutations reflective of particular conditions. A phylogenetic analysis found the tradeoff in numerous bacteria species. The results suggest that the fear vs. greed tradeoff represents a general principle of transcriptome allocation in bacteria where small genetic changes can result in large phenotypic adaptations to growth conditions.

Inference of transcriptome signatures of Escherichia coli in long-term stationary phase

"Non-growing" is a dominant life form of microorganisms in nature, where available nutrients and resources are limited. In laboratory culture systems, Escherichia coli can survive for years under starvation, denoted as long-term stationary phase, where a small fraction of cells manages to survive by recycling resources released from nonviable cells. Although the physiology by which viable cells in long-term stationary phase adapt to prolonged starvation is of great interest, their genome-wide response has not been fully understood. In this study, we analyzed transcriptional profiles of cells exposed to the supernatant of 30-day long-term stationary phase culture and found that their transcriptome profiles displayed several similar responses to those of cells in the 16-h short-term stationary phase. Nevertheless, our results revealed that cells in long-term stationary phase supernatant exhibit higher expressions of stress-response genes such as phage shock proteins (psp), and lower expressions of growth-related genes such as ribosomal proteins than those in the short-term stationary phase. We confirmed that the mutant lacking the psp operon showed lower survival and growth rate in the long-term stationary phase culture. This study identified transcriptional responses for stress-resistant physiology in the long-term stationary phase environment.

Metabolic reprogramming and biosensor-assisted mutagenesis screening for high-level production of L-arginine in Escherichia coli

L-arginine is a value-added amino acid with promising applications in the pharmaceutical and nutraceutical industries. Further unleashing the potential of microbial cell factories to make L-arginine production more competitive remains challenging due to the sophisticated intracellular interaction networks and the insufficient knowledge of global metabolic regulation. Here, we combined multilevel rational metabolic engineering with biosensor-assisted mutagenesis screening to exploit the L-arginine production potential of Escherichia coli. First, multiple metabolic pathways were systematically reprogrammed to redirect the metabolic flux into L-arginine synthesis, including the L-arginine biosynthesis, TCA cycle, and L-arginine export. Specifically, a toggle switch responding to special cellular physiological conditions was designed to dynamically control the expression of sucA and pull more carbon flux from the TCA cycle toward L-arginine biosynthesis. Subsequently, a biosensor-assisted high-throughput screening platform was designed and applied to further exploit the L-arginine production potential. The best-engineered ARG28 strain produced 132 g/L L-arginine in a 5-L bioreactor with a yield of 0.51 g/g glucose and productivity of 2.75 g/(L ⋅ h), which were the highest values reported so far. Through whole genome sequencing and reverse engineering, Frc frameshift mutant, PqiB A78P mutant, and RpoB P564T mutant were revealed for enhancing the L-arginine biosynthesis. Our study exhibited the power of coupling rational metabolic reprogramming and biosensor-assisted mutagenesis screening to unleash the cellular potential for value-added metabolite production.

Laboratory evolution reveals general and specific tolerance mechanisms for commodity chemicals

Although strain tolerance to high product concentrations is a barrier to the economically viable biomanufacturing of industrial chemicals, chemical tolerance mechanisms are often unknown. To reveal tolerance mechanisms, an automated platform was utilized to evolve Escherichia coli to grow optimally in the presence of 11 industrial chemicals (1,2-propanediol, 2,3-butanediol, glutarate, adipate, putrescine, hexamethylenediamine, butanol, isobutyrate, coumarate, octanoate, hexanoate), reaching tolerance at concentrations 60%-400% higher than initial toxic levels. Sequencing genomes of 223 isolates from 89 populations, reverse engineering, and cross-compound tolerance profiling were employed to uncover tolerance mechanisms. We show that: 1) cells are tolerized via frequent mutation of membrane transporters or cell wall-associated proteins (e.g., ProV, KgtP, SapB, NagA, NagC, MreB), transcription and translation machineries (e.g., RpoA, RpoB, RpoC, RpsA, RpsG, NusA, Rho), stress signaling proteins (e.g., RelA, SspA, SpoT, YobF), and for certain chemicals, regulators and enzymes in metabolism (e.g., MetJ, NadR, GudD, PurT); 2) osmotic stress plays a significant role in tolerance when chemical concentrations exceed a general threshold and mutated genes frequently overlap with those enabling chemical tolerance in membrane transporters and cell wall-associated proteins; 3) tolerization to a specific chemical generally improves tolerance to structurally similar compounds whereas a tradeoff can occur on dissimilar chemicals, and 4) using pre-tolerized starting isolates can hugely enhance the subsequent production of chemicals when a production pathway is inserted in many, but not all, evolved tolerized host strains, underpinning the need for evolving multiple parallel populations. Taken as a whole, this study provides a comprehensive genotype-phenotype map based on identified mutations and growth phenotypes for 223 chemical tolerant isolates.

A Systems Biology Approach To Disentangle the Direct and Indirect Effects of Global Transcription Factors on Gene Expression in Escherichia coli

Delineating the pleiotropic effects of global transcriptional factors (TFs) is critical for understanding the system-wide regulatory response in a particular environment. Currently, with the availability of genome-wide TF binding and gene expression data for Escherichia coli, several gene targets can be assigned to the global TFs, albeit inconsistently. Here, using a systematic integrated approach with emphasis on metabolism, we characterized and quantified the direct effects as well as the growth rate-mediated indirect effects of global TFs using deletion mutants of FNR, ArcA, and IHF regulators (focal TFs) under glucose fermentative conditions. This categorization enabled us to disentangle the dense connections seen within the transcriptional regulatory network (TRN) and determine the exact nature of focal TF-driven epistatic interactions with other global and pathway-specific local regulators (iTFs). We extended our analysis to combinatorial deletions of these focal TFs to determine their cross talk effects as well as conserved patterns of regulatory interactions. Moreover, we predicted with high confidence several novel metabolite-iTF interactions using inferred iTF activity changes arising from the allosteric effects of the intracellular metabolites perturbed as a result of the absence of focal TFs. Further, using compendium level computational analyses, we revealed not only the coexpressed genes regulated by these focal TFs but also the coordination of the direct and indirect target expression in the context of the economy of intracellular metabolites. Overall, this study leverages the fundamentals of TF-driven regulation, which could serve as a better template for deciphering mechanisms underlying complex phenotypes. IMPORTANCE Understanding the pleiotropic effects of global TFs on gene expression and their relevance underlying a specific response in a particular environment has been challenging. Here, we distinguish the TF-driven direct effects and growth rate-mediated indirect effects on gene expression using single- and double-deletion mutants of FNR, ArcA, and IHF regulators under anaerobic glucose fermentation. Such dissection assists us in unraveling the precise nature of interactions existing between the focal TF(s) and several other TFs, including those altered by allosteric effects of intracellular metabolites. We were able to recapitulate the previously known metabolite-TF interactions and predict novel interactions with high confidence. Furthermore, we determined that the direct and indirect gene expression have a strong connection with each other when analyzed using the coexpressed- or coregulated-gene approach. Deciphering such regulatory patterns explicitly from the metabolism point of view would be valuable in understanding other unpredicted complex regulation existing in nature.

Conversion of Escherichia coli into Mixotrophic CO2 Assimilation with Malate and Hydrogen Based on Recombinant Expression of 2-Oxoglutarate:Ferredoxin Oxidoreductase Using Adaptive Laboratory Evolution

We report the mixotrophic growth of Escherichia coli based on recombinant 2-oxoglutarate:ferredoxin oxidoreductase (OGOR) to assimilate CO₂ using malate as an auxiliary carbon source and hydrogen as an energy source. We employ a long-term (~184 days) two-stage adaptive evolution to convert heterotrophic E. coli into mixotrophic E. coli. In the first stage of evolution with serine, diauxic growth emerges as a prominent feature. At the end of the second stage of evolution with malate, the strain exhibits mixotrophy with CO₂ as an essential substrate for growth. We expect this work will open new possibilities in the utilization of OGOR for microbial CO₂ assimilation and future hydrogen-based electro-microbial conversion.

Experimental Evolution Reveals Unifying Systems-Level Adaptations but Diversity in Driving Genotypes

Genotype-fitness maps of evolution have been well characterized for biological components, such as RNA and proteins, but remain less clear for systems-level properties, such as those of metabolic and transcriptional regulatory networks. Here, we take multi-omics measurements of 6 different E. coli strains throughout adaptive laboratory evolution (ALE) to maximal growth fitness. The results show the following: (i) convergence in most overall phenotypic measures across all strains, with the notable exception of divergence in NADPH production mechanisms; (ii) conserved transcriptomic adaptations, describing increased expression of growth promoting genes but decreased expression of stress response and structural components; (iii) four groups of regulatory trade-offs underlying the adjustment of transcriptome composition; and (iv) correlates that link causal mutations to systems-level adaptations, including mutation-pathway flux correlates and mutation-transcriptome composition correlates. We thus show that fitness landscapes for ALE can be described with two layers of causation: one based on system-level properties (continuous variables) and the other based on mutations (discrete variables). IMPORTANCE Understanding the mechanisms of microbial adaptation will help combat the evolution of drug-resistant microbes and enable predictive genome design. Although experimental evolution allows us to identify the causal mutations underlying microbial adaptation, it remains unclear how causal mutations enable increased fitness and is often explained in terms of individual components (i.e., enzyme rate) as opposed to biological systems (i.e., pathways). Here, we find that causal mutations in E. coli are linked to systems-level changes in NADPH balance and expression of stress response genes. These systems-level adaptation patterns are conserved across diverse E. coli strains and thus identify cofactor balance and proteome reallocation as dominant constraints governing microbial adaptation.

Coordination of CcpA and CodY Regulators in Staphylococcus aureus USA300 Strains

The complex cross talk between metabolism and gene regulatory networks makes it difficult to untangle individual constituents and study their precise roles and interactions. To address this issue, we modularized the transcriptional regulatory network (TRN) of the Staphylococcus aureus USA300 strain by applying independent component analysis (ICA) to 385 RNA sequencing samples. We then combined the modular TRN model with a metabolic model to study the regulation of carbon and amino acid metabolism. Our analysis showed that regulation of central carbon metabolism by CcpA and amino acid biosynthesis by CodY are closely coordinated. In general, S. aureus increases the expression of CodY-regulated genes in the presence of preferred carbon sources such as glucose. This transcriptional coordination was corroborated by metabolic model simulations that also showed increased amino acid biosynthesis in the presence of glucose. Further, we found that CodY and CcpA cooperatively regulate the expression of ribosome hibernation-promoting factor, thus linking metabolic cues with translation. In line with this hypothesis, expression of CodY-regulated genes is tightly correlated with expression of genes encoding ribosomal proteins. Together, we propose a coarse-grained model where expression of S. aureus genes encoding enzymes that control carbon flux and nitrogen flux through the system is coregulated with expression of translation machinery to modularly control protein synthesis. While this work focuses on three key regulators, the full TRN model we present contains 76 total independently modulated sets of genes, each with the potential to uncover other complex regulatory structures and interactions. IMPORTANCE Staphylococcus aureus is a versatile pathogen with an expanding antibiotic resistance profile. The biology underlying its clinical success emerges from an interplay of many systems such as metabolism and gene regulatory networks. This work brings together models for these two systems to establish fundamental principles governing the regulation of S. aureus central metabolism and protein synthesis. Studies of these fundamental biological principles are often confined to model organisms such as Escherichia coli. However, expanding these models to pathogens can provide a framework from which complex and clinically important phenotypes such as virulence and antibiotic resistance can be better understood. Additionally, the expanded gene regulatory network model presented here can deconvolute the biology underlying other important phenotypes in this pathogen.

Mitigating Host Burden of Genetic Circuits by Engineering Autonegatively Regulated Parts and Improving Functional Prediction

Mitigating unintended interferences between circuits and host cells is key to realize applications of synthetic regulatory systems both for bacteria and mammalian cells. Here, we demonstrated that growth burden and circuit dysregulation occurred in a concentration-dependent manner for specific transcription factors (CymR*/CymR) in E.coli, and direct negative feedback modules were able to control the concentration of CymR*/CymR, mitigate growth burden, and restore circuit functions. A quantitative design scheme was developed for circuits embedded with autorepression modules. Four key parameters were theoretically identified to determine the performance of autoregulated switches and were experimentally modified by fine-tuning promoter architectures and cooperativity. Using this strategy, we synthesized a number of switches and demonstrated its improvement of product titers and host growth controlling the complex deoxyviolacein biosynthesis pathway. Furthermore, we restored functions of a dysregulated multilayer NOR gate by integrating autorepression modules. Our work provides a blueprint for engineering host-adaptable synthetic systems.

Laboratory evolution of synthetic electron transport system variants reveals a larger metabolic respiratory system and its plasticity

The bacterial respiratory electron transport system (ETS) is branched to allow condition-specific modulation of energy metabolism. There is a detailed understanding of the structural and biochemical features of respiratory enzymes; however, a holistic examination of the system and its plasticity is lacking. Here we generate four strains of Escherichia coli harboring unbranched ETS that pump 1, 2, 3, or 4 proton(s) per electron and characterized them using a combination of synergistic methods (adaptive laboratory evolution, multi-omic analyses, and computation of proteome allocation). We report that: (a) all four ETS variants evolve to a similar optimized growth rate, and (b) the laboratory evolutions generate specific rewiring of major energy-generating pathways, coupled to the ETS, to optimize ATP production capability. We thus define an Aero-Type System (ATS), which is a generalization of the aerobic bioenergetics and is a metabolic systems biology description of respiration and its inherent plasticity.

The bacterial respiratory electron transport system (ETS) is branched to allow condition-specific modulation of energy metabolism. Here the authors examine the systems level properties of aerobic electron transport system using adaptive laboratory evolution and multi-omics analyses.

Regulatory perturbations of ribosome allocation in bacteria reshape the growth proteome with a trade-off in adaptation capacity

Bacteria regulate their cellular resource allocation to enable fast growth-adaptation to a variety of environmental niches. We studied the ribosomal allocation, growth, and expression profiles of two sets of fast-growing mutants of Escherichia coli K-12 MG1655. Mutants with only three of the seven copies of ribosomal RNA operons grew faster than the wild-type strain in minimal media and show similar phenotype to previously studied fast-growing rpoB mutants. Comparing these two different regulatory perturbations (rRNA promoters or rpoB mutations), we show how they reshape the proteome for growth with a concomitant fitness cost. The fast-growing mutants shared downregulation of hedging functions and upregulated growth functions. They showed longer diauxic shifts and reduced activity of gluconeogenic promoters during glucose-acetate shifts, suggesting reduced availability of the RNA polymerase for expressing hedging proteome. These results show that the regulation of ribosomal allocation underlies the growth/hedging phenotypes obtained from laboratory evolution experiments.

•

Mutants with only three ribosomal operons grow faster than wild-type in minimal medium

•

Faster growth of mutants is achieved by increased ribosome content

•

Fast-growing mutants display reduced hedging expression and adaptation trade-offs

Global pleiotropic effects in adaptively evolved Escherichia coli lacking CRP reveal molecular mechanisms that define the growth physiology

Evolution facilitates emergence of fitter phenotypes by efficient allocation of cellular resources in conjunction with beneficial mutations. However, system-wide pleiotropic effects that redress the perturbations to the apex node of the transcriptional regulatory networks remain unclear. Here, we elucidate that absence of global transcriptional regulator CRP in Escherichia coli results in alterations in key metabolic pathways under glucose respiratory conditions, favouring stress- or hedging-related functions over growth-enhancing functions. Further, we disentangle the growth-mediated effects from the CRP regulation-specific effects on these metabolic pathways. We quantitatively illustrate that the loss of CRP perturbs proteome efficiency, as evident from metabolic as well as ribosomal proteome fractions, that corroborated with intracellular metabolite profiles. To address how E. coli copes with such systemic defect, we evolved Δcrp mutant in the presence of glucose. Besides acquiring mutations in the promoter of glucose transporter ptsG, the evolved populations recovered the metabolic pathways to their pre-perturbed state coupled with metabolite re-adjustments, which altogether enabled increased growth. By contrast to Δcrp mutant, the evolved strains remodelled their proteome efficiency towards biomass synthesis, albeit at the expense of carbon efficiency. Overall, we comprehensively illustrate the genetic and metabolic basis of pleiotropic effects, fundamental for understanding the growth physiology.

Environmental dependence of competitive fitness in rifampin-resistant rpoB mutants of Bacillus subtilis

RNA polymerase (RNAP) is a highly conserved macromolecular machine that contributes to the flow of genetic information from genotype to phenotype. In Bacillus subtilis , mutations in the rpoB gene encoding the β-subunit of RNAP have been shown to alter a number of global phenotypes including growth, utilization of unusual nutrient sources, sporulation, germination, and production of secondary metabolites. In addition, the spectrum of mutations in rpoB leading to rifampin resistance (Rif R ) can change dramatically depending upon the environment to which B. subtilis cells or spores are exposed. Rif R rpoB mutations have historically been associated with slower growth and reduced fitness; however, these assessments of fitness were conducted on limited collections of mutants in rich laboratory media that poorly reflect natural environments typically inhabited by B. subtilis . Using a novel, deep-sequencing approach in addition to traditional measurements of growth rate, lag time, and pairwise competitions, we demonstrated the competitive advantage of specific rpoB alleles differs depending on the growth environment in which they are determined.

IMPORTANCE

Microbial resistance to antibiotics is a growing threat to public health across the world. Historically, resistance to antibiotics has been associated with reduced fitness. A growing body of evidence indicates that resistance to rifampin, a frontline antibiotic used to treat mycobacterial and biofilm-associated infections, may increase fitness given an appropriate environment even in the absence of the selective antibiotic. Here we experimentally confirm this phenomenon by directly comparing the fitness of multiple rifampin-resistant mutants of Bacillus subtilis in rich LB medium and an asparagine minimal medium. Our research demonstrates that the fitness cost of rifampin resistance can vary greatly depending upon the environment. This has important implications for understanding how microbes develop antimicrobial resistance in the absence of antibiotic selection.

Reconstructing the transcription regulatory network to optimize resource allocation for robust biosynthesis

An ideal microbial cell factory (MCF) should deliver maximal resources to production, which conflicts with the microbe's native growth-oriented resource allocation strategy and can therefore lead to early termination of the high-yield period. Reallocating resources from growth to production has become a critical factor in constructing robust MCFs. Instead of strengthening specific biosynthetic pathways, emerging endeavors are focused on rearranging the gene regulatory network to fundamentally reprogram the resource allocation pattern. Combining this idea with transcriptional regulation within the hierarchical regulatory network, this review discusses recent engineering strategies targeting the transcription machinery, module networks, regulatory edges, and bottom network layer. This global view will help to construct a production-oriented phenotype that fully harnesses the potential of MCFs.

Dissecting the Fitness Costs of Complex Mutations

The fitness cost of complex pleiotropic mutations is generally difficult to assess. On the one hand, it is necessary to identify which molecular properties are directly altered by the mutation. On the other, this alteration modifies the activity of many genetic targets with uncertain consequences. Here, we examine the possibility of addressing these challenges by identifying unique predictors of these costs. To this aim, we consider mutations in the RNA polymerase (RNAP) in Escherichia coli as a model of complex mutations. Changes in RNAP modify the global program of transcriptional regulation, with many consequences. Among others is the difficulty to decouple the direct effect of the mutation from the response of the whole system to such mutation. A problem that we solve quantitatively with data of a set of constitutive genes, those on which the global program acts most directly. We provide a statistical framework that incorporates the direct effects and other molecular variables linked to this program as predictors, which leads to the identification that some genes are more suitable to determine costs than others. Therefore, we not only identified which molecular properties best anticipate fitness, but we also present the paradoxical result that, despite pleiotropy, specific genes serve as more solid predictors. These results have connotations for the understanding of the architecture of robustness in biological systems.

Machine Learning of Bacterial Transcriptomes Reveals Responses Underlying Differential Antibiotic Susceptibility

In vitro antibiotic susceptibility testing often fails to accurately predict in vivo drug efficacies, in part due to differences in the molecular composition between standardized bacteriologic media and physiological environments within the body. Here, we investigate the interrelationship between antibiotic susceptibility and medium composition in Escherichia coli K-12 MG1655 as contextualized through machine learning of transcriptomics data. Application of independent component analysis, a signal separation algorithm, shows that complex phenotypic changes induced by environmental conditions or antibiotic treatment are directly traced to the action of a few key transcriptional regulators, including RpoS, Fur, and Fnr. Integrating machine learning results with biochemical knowledge of transcription factor activation reveals medium-dependent shifts in respiration and iron availability that drive differential antibiotic susceptibility. By extension, the data generation and data analytics workflow used here can interrogate the regulatory state of a pathogen under any measured condition and can be applied to any strain or organism for which sufficient transcriptomics data are available.

IMPORTANCE Antibiotic resistance is an imminent threat to global health. Patient treatment regimens are often selected based on results from standardized antibiotic susceptibility testing (AST) in the clinical microbiology lab, but these in vitro tests frequently misclassify drug effectiveness due to their poor resemblance to actual host conditions. Prior attempts to understand the combined effects of drugs and media on antibiotic efficacy have focused on physiological measurements but have not linked treatment outcomes to transcriptional responses on a systems level. Here, application of machine learning to transcriptomics data identified medium-dependent responses in key regulators of bacterial iron uptake and respiratory activity. The analytical workflow presented here is scalable to additional organisms and conditions and could be used to improve clinical AST by identifying the key regulatory factors dictating antibiotic susceptibility.

Restoration of fitness lost due to dysregulation of the pyruvate dehydrogenase complex is triggered by ribosomal binding site modifications

Pyruvate dehydrogenase complex (PDC) functions as the main determinant of the respiro-fermentative balance because it converts pyruvate to acetyl-coenzyme A (CoA), which then enters the TCA (tricarboxylic acid cycle). PDC is repressed by the pyruvate dehydrogenase complex regulator (PdhR) in Escherichia coli. The deletion of the pdhR gene compromises fitness in aerobic environments. We evolve the E. coli pdhR deletion strain to examine its achievable growth rate and the underlying adaptive strategies. We find that (1) optimal proteome allocation to PDC is critical in achieving optimal growth rate; (2) expression of PDC in evolved strains is reduced through mutations in the Shine-Dalgarno sequence; (3) rewiring of the TCA flux and increased reactive oxygen species (ROS) defense occur in the evolved strains; and (4) the evolved strains adapt to an efficient biomass yield. Together, these results show how adaptation can find alternative regulatory mechanisms for a key cellular process if the primary regulatory mode fails.

Global Transcriptional Regulators Fine-Tune the Translational and Metabolic Efficiency for Optimal Growth of Escherichia coli

Global transcriptional regulators coordinate complex genetic interactions that bestow better adaptability for an organism against external and internal perturbations. These transcriptional regulators are known to control an enormous array of genes with diverse functionalities. However, regulator-driven molecular mechanisms that underpin precisely tuned translational and metabolic processes conducive for rapid exponential growth remain obscure. Here, we comprehensively reveal the fundamental role of global transcriptional regulators FNR, ArcA, and IHF in sustaining translational and metabolic efficiency under glucose fermentative conditions in Escherichia coli By integrating high-throughput gene expression profiles and absolute intracellular metabolite concentrations, we illustrate that these regulators are crucial in maintaining nitrogen homeostasis, govern expression of otherwise unnecessary or hedging genes, and exert tight control on metabolic bottleneck steps. Furthermore, we characterize changes in expression and activity profiles of other coregulators associated with these dysregulated metabolic pathways, determining the regulatory interactions within the transcriptional regulatory network. Such systematic findings emphasize their importance in optimizing the proteome allocation toward metabolic enzymes as well as ribosomes, facilitating condition-specific phenotypic outcomes. Consequentially, we reveal that disruption of this inherent trade-off between ribosome and metabolic proteome economy due to the loss of regulators resulted in lowered growth rates. Moreover, our findings reinforce that the accumulations of intracellular metabolites in the event of proteome repartitions negatively affects the glucose uptake. Overall, by extending the three-partition proteome allocation theory concordant with multi-omics measurements, we elucidate the physiological consequences of loss of global regulators on central carbon metabolism restraining the organism to attain maximal biomass synthesis.IMPORTANCE Cellular proteome allocation in response to environmental or internal perturbations governs their eventual phenotypic outcome. This entails striking an effective balance between amino acid biosynthesis by metabolic proteins and its consumption by ribosomes. However, the global transcriptional regulator-driven molecular mechanisms that underpin their coordination remains unexplored. Here, we emphasize that global transcriptional regulators, known to control expression of a myriad of genes, are fundamental for precisely tuning the translational and metabolic efficiencies that define the growth optimality. Towards this, we systematically characterized the single deletion effect of FNR, ArcA, and IHF regulators of Escherichia coli on exponential growth under anaerobic glucose fermentative conditions. Their absence disrupts the stringency of proteome allocation, which manifests as impairment in key metabolic processes and an accumulation of intracellular metabolites. Furthermore, by incorporating an extension to the empirical growth laws, we quantitatively demonstrate the general design principles underlying the existence of these regulators in E. coli.

Independent component analysis recovers consistent regulatory signals from disparate datasets

The availability of bacterial transcriptomes has dramatically increased in recent years. This data deluge could result in detailed inference of underlying regulatory networks, but the diversity of experimental platforms and protocols introduces critical biases that could hinder scalable analysis of existing data. Here, we show that the underlying structure of the E. coli transcriptome, as determined by Independent Component Analysis (ICA), is conserved across multiple independent datasets, including both RNA-seq and microarray datasets. We subsequently combined five transcriptomics datasets into a large compendium containing over 800 expression profiles and discovered that its underlying ICA-based structure was still comparable to that of the individual datasets. With this understanding, we expanded our analysis to over 3,000 E. coli expression profiles and predicted three high-impact regulons that respond to oxidative stress, anaerobiosis, and antibiotic treatment. ICA thus enables deep analysis of disparate data to uncover new insights that were not visible in the individual datasets.

Cells adapt to diverse environments by regulating gene expression. Genome-wide measurements of gene expression levels have exponentially increased in recent years, but successful integration and analysis of these datasets are limited. Recently, we showed that independent component analysis (ICA), a signal deconvolution algorithm, can separate a large bacterial gene expression dataset into groups of co-regulated genes. This previous study focused on data generated by a standardized pipeline and did not address whether ICA extracts the same quantitative co-expression signals across expression profiling platforms. In this study, we show that ICA finds similar co-regulation patterns underlying multiple gene expression datasets and can be used as a tool to integrate and interpret diverse datasets. Using a dataset containing over 3,000 expression profiles, we predicted three new regulons and characterized their activities. Since large, standardized expression datasets only exist for a few bacterial strains, these results broaden the possible applications of this tool to better understand transcriptional regulation across a wide range of microbes.

OxyR Is a Convergent Target for Mutations Acquired during Adaptation to Oxidative Stress-Prone Metabolic States

Oxidative stress is concomitant with aerobic metabolism. Thus, bacterial genomes encode elaborate mechanisms to achieve redox homeostasis. Here we report that the peroxide-sensing transcription factor, oxyR, is a common mutational target using bacterial species belonging to two genera, Escherichia coli and Vibrio natriegens, in separate growth conditions implemented during laboratory evolution. The mutations clustered in the redox active site, dimer interface, and flexible redox loop of the protein. These mutations favor the oxidized conformation of OxyR that results in constitutive expression of the genes it regulates. Independent component analysis of the transcriptome revealed that the constitutive activity of OxyR reduces DNA damage from reactive oxygen species, as inferred from the activity of the SOS response regulator LexA. This adaptation to peroxide stress came at a cost of lower growth, as revealed by calculations of proteome allocation using genome-scale models of metabolism and macromolecular expression. Further, identification of similar sequence changes in natural isolates of E. coli indicates that adaptation to oxidative stress through genetic changes in oxyR can be a common occurrence.

Adaptive laboratory evolution of Pseudomonas putida KT2440 improves p-coumaric and ferulic acid catabolism and tolerance

Pseudomonas putida KT2440 is a promising bacterial chassis for the conversion of lignin-derived aromatic compound mixtures to biofuels and bioproducts. Despite the inherent robustness of this strain, further improvements to aromatic catabolism and toxicity tolerance of P. putida will be required to achieve industrial relevance. Here, tolerance adaptive laboratory evolution (TALE) was employed with increasing concentrations of the hydroxycinnamic acids p-coumaric acid (pCA) and ferulic acid (FA) individually and in combination (pCA ​+ ​FA). The TALE experiments led to evolved P. putida strains with increased tolerance to the targeted acids as compared to wild type. Specifically, a 37 ​h decrease in lag phase in 20 ​g/L pCA and a 2.4-fold increase in growth rate in 30 ​g/L FA was observed. Whole genome sequencing of intermediate and endpoint evolved P. putida populations revealed several expected and non-intuitive genetic targets underlying these aromatic catabolic and toxicity tolerance enhancements. PP_3350 and ttgB were among the most frequently mutated genes, and the beneficial contributions of these mutations were verified via gene knockouts. Deletion of PP_3350, encoding a hypothetical protein, recapitulated improved toxicity tolerance to high concentrations of pCA, but not an improved growth rate in high concentrations of FA. Deletion of ttgB, part of the TtgABC efflux pump, severely inhibited growth in pCA ​+ ​FA TALE-derived strains but did not affect growth in pCA ​+ ​FA in a wild type background, suggesting epistatic interactions. Genes involved in flagellar movement and transcriptional regulation were often mutated in the TALE experiments on multiple substrates, reinforcing ideas of a minimal and deregulated cell as optimal for domesticated growth. Overall, this work demonstrates increased tolerance towards and growth rate at the expense of hydroxycinnamic acids and presents new targets for improving P. putida for microbial lignin valorization.

•

Pseudomonas putida was evolved in increasing hydroxycinnamic acid concentrations.

•

The lag phase in high p-coumaric acid concentration was reduced.

•

The growth rate and tolerance to high ferulic acid concentration was increased.

•

PP_3350 and ttgB contribute to the observed phenotypes.

Independent component analysis of E. coli's transcriptome reveals the cellular processes that respond to heterologous gene expression

Achieving the predictable expression of heterologous genes in a production host has proven difficult. Each heterologous gene expressed in the same host seems to elicit a different host response governed by unknown mechanisms. Historically, most studies have approached this challenge by manipulating the properties of the heterologous gene through methods like codon optimization. Here we approach this challenge from the host side. We express a set of 45 heterologous genes in the same Escherichia coli strain, using the same expression system and culture conditions. We collect a comprehensive RNAseq set to characterize the host's transcriptional response. Independent Component Analysis of the RNAseq data set reveals independently modulated gene sets (iModulons) that characterize the host response to heterologous gene expression. We relate 55% of variation of the host response to: Fear vs Greed (16.5%), Metal Homeostasis (19.0%), Respiration (6.0%), Protein folding (4.5%), and Amino acid and nucleotide biosynthesis (9.0%). If these responses can be controlled, then the success rate with predicting heterologous gene expression should increase.

A universal trade-off between growth and lag in fluctuating environments

The rate of cell growth is crucial for bacterial fitness and a main driver of proteome allocation^1,2, but it is unclear what ultimately determines growth rates in different environmental conditions. Increasing evidence suggests that other objectives also play key roles^3–7, such as the rate of physiological adaptation to changing environments^8,9. The challenge for cells is that these objectives often cannot be independently optimized, and maximizing one often reduces another. Many such tradeoffs have indeed been hypothesized, based on qualitative correlative studies^8–11. Here we report the occurrence of a tradeoff between steady-state growth rate and physiological adaptability for Escherichia coli, upon abruptly shifting a growing culture from a preferred carbon source to fermentation products such as acetate. Such transitions, common for enteric bacteria, are often accompanied by multi-hour lags before growth resumes. Metabolomic analysis revealed that the long lags resulted from the depletion of key metabolites due to the sudden reversal of central carbon flux imposed by these nutrient shifts. A model of sequential flux limitation not only explained the observed universal tradeoff between growth and adaptability, but also generated quantitative predictions that were validated experimentally. The observed trade-off reflects the opposing enzyme requirements for glycolysis versus gluconeogenesis.

A quantitative method for proteome reallocation using minimal regulatory interventions

Engineering resource allocation in biological systems is an ongoing challenge. Organisms allocate resources for ensuring survival, reducing the productivity of synthetic biology functions. Here we present a new approach for engineering the resource allocation of Escherichia coli by rationally modifying its transcriptional regulatory network. Our method (ReProMin) identifies the minimal set of genetic interventions that maximizes the savings in cell resources. To this end, we categorized transcription factors according to the essentiality of its targets and we used proteomic data to rank them. We designed the combinatorial removal of transcription factors that maximize the release of resources. Our resulting strain containing only three mutations, theoretically releasing 0.5% of its proteome, had higher proteome budget, increased production of an engineered metabolic pathway and showed that the regulatory interventions are highly specific. This approach shows that combining proteomic and regulatory data is an effective way of optimizing strains using conventional molecular methods.

Revealing 29 sets of independently modulated genes in Staphylococcus aureus, their regulators, and role in key physiological response

Staphylococcus aureus infections impose an immense burden on the healthcare system. To establish a successful infection in a hostile host environment, S. aureus must coordinate its gene expression to respond to a wide array of challenges. This balancing act is largely orchestrated by the transcriptional regulatory network. Here, we present a model of 29 independently modulated sets of genes that form the basis for a segment of the transcriptional regulatory network in clinical USA300 strains of S. aureus. Using this model, we demonstrate the concerted role of various cellular systems (e.g., metabolism, virulence, and stress response) underlying key physiological responses, including response during blood infection.

The ability of Staphylococcus aureus to infect many different tissue sites is enabled, in part, by its transcriptional regulatory network (TRN) that coordinates its gene expression to respond to different environments. We elucidated the organization and activity of this TRN by applying independent component analysis to a compendium of 108 RNA-sequencing expression profiles from two S. aureus clinical strains (TCH1516 and LAC). ICA decomposed the S. aureus transcriptome into 29 independently modulated sets of genes (i-modulons) that revealed: 1) High confidence associations between 21 i-modulons and known regulators; 2) an association between an i-modulon and σS, whose regulatory role was previously undefined; 3) the regulatory organization of 65 virulence factors in the form of three i-modulons associated with AgrR, SaeR, and Vim-3; 4) the roles of three key transcription factors (CodY, Fur, and CcpA) in coordinating the metabolic and regulatory networks; and 5) a low-dimensional representation, involving the function of few transcription factors of changes in gene expression between two laboratory media (RPMI, cation adjust Mueller Hinton broth) and two physiological media (blood and serum). This representation of the TRN covers 842 genes representing 76% of the variance in gene expression that provides a quantitative reconstruction of transcriptional modules in S. aureus, and a platform enabling its full elucidation.

Metabolic role of pyrophosphate-linked phosphofructokinase pfk for C1 assimilation in Methylotuvimicrobium alcaliphilum 20Z

BACKGROUND

Methanotrophs is a promising biocatalyst in biotechnological applications with their ability to utilize single carbon (C1) feedstock to produce high-value compounds. Understanding the behavior of biological networks of methanotrophic bacteria in different parameters is vital to systems biology and metabolic engineering. Interestingly, methanotrophic bacteria possess the pyrophosphate-dependent 6-phosphofructokinase (PPi-PFK) instead of the ATP-dependent 6-phosphofructokinase, indicating their potentials to serve as promising model for investigation the role of inorganic pyrophosphate (PPi) and PPi-dependent glycolysis in bacteria. Gene knockout experiments along with global-omics approaches can be used for studying gene functions as well as unraveling regulatory networks that rely on the gene product.

RESULTS

In this study, we performed gene knockout and RNA-seq experiments in Methylotuvimicrobium alcaliphilum 20Z to investigate the functional roles of PPi-PFK in C1 metabolism when cells were grown on methane and methanol, highlighting its metabolic importance in C1 assimilation in M. alcaliphilum 20Z. We further conducted adaptive laboratory evolution (ALE) to investigate regulatory architecture in pfk knockout strain. Whole-genome resequencing and RNA-seq approaches were performed to characterize the genetic and metabolic responses of adaptation to pfk knockout. A number of mutations, as well as gene expression profiles, were identified in pfk ALE strain to overcome insufficient C1 assimilation pathway which limits the growth in the unevolved strain.

CONCLUSIONS

This study first revealed the regulatory roles of PPi-PFK on C1 metabolism and then provided novel insights into mechanism of adaptation to the loss of this major metabolic enzyme as well as an improved basis for future strain design in type I methanotrophs.

Conversion of Escherichia coli to Generate All Biomass Carbon from CO2

The living world is largely divided into autotrophs that convert CO₂ into biomass and heterotrophs that consume organic compounds. In spite of widespread interest in renewable energy storage and more sustainable food production, the engineering of industrially relevant heterotrophic model organisms to use CO₂ as their sole carbon source has so far remained an outstanding challenge. Here, we report the achievement of this transformation on laboratory timescales. We constructed and evolved Escherichia coli to produce all its biomass carbon from CO₂. Reducing power and energy, but not carbon, are supplied via the one-carbon molecule formate, which can be produced electrochemically. Rubisco and phosphoribulokinase were co-expressed with formate dehydrogenase to enable CO₂ fixation and reduction via the Calvin-Benson-Bassham cycle. Autotrophic growth was achieved following several months of continuous laboratory evolution in a chemostat under intensifying organic carbon limitation and confirmed via isotopic labeling.

Directed Metabolic Pathway Evolution Enables Functional Pterin-Dependent Aromatic-Amino-Acid Hydroxylation in Escherichia coli

Tetrahydrobiopterin-dependent hydroxylation of aromatic amino acids is the first step in the biosynthesis of many neuroactive compounds in humans. A fundamental challenge in building these pathways in Escherichia coli is the provision of the non-native hydroxylase cofactor, tetrahydrobiopterin. To solve this, we designed a genetic selection that relies on the tyrosine synthesis activity of phenylalanine hydroxylase. Using adaptive laboratory evolution, we demonstrate the use of this selection to discover: (1) a minimum set of heterologous enzymes and a host folE (T198I) mutation for achieving this type of hydroxylation chemistry in whole cells, (2) functional complementation of tetrahydrobiopterin by indigenous cofactors, and (3) a tryptophan hydroxylase mutation for improving protein abundance. Thus, the goal of having functional aromatic-amino-acid hydroxylation in E. coli was achieved through directed metabolic pathway evolution.

The Escherichia coli transcriptome mostly consists of independently regulated modules

Underlying cellular responses is a transcriptional regulatory network (TRN) that modulates gene expression. A useful description of the TRN would decompose the transcriptome into targeted effects of individual transcriptional regulators. Here, we apply unsupervised machine learning to a diverse compendium of over 250 high-quality Escherichia coli RNA-seq datasets to identify 92 statistically independent signals that modulate the expression of specific gene sets. We show that 61 of these transcriptomic signals represent the effects of currently characterized transcriptional regulators. Condition-specific activation of signals is validated by exposure of E. coli to new environmental conditions. The resulting decomposition of the transcriptome provides: a mechanistic, systems-level, network-based explanation of responses to environmental and genetic perturbations; a guide to gene and regulator function discovery; and a basis for characterizing transcriptomic differences in multiple strains. Taken together, our results show that signal summation describes the composition of a model prokaryotic transcriptome.

Mechanistic insight into the regulation of transcriptional modules remains scarce. Here, the authors identify statistically independent gene sets by applying independent component analysis to a high-quality E. coli RNA-seq data compendium and find that most gene sets represent the effects of specific transcriptional regulators.

Common principles and best practices for engineering microbiomes

Despite broad scientific interest in harnessing the power of Earth's microbiomes, knowledge gaps hinder their efficient use for addressing urgent societal and environmental challenges. We argue that structuring research and technology developments around a design-build-test-learn (DBTL) cycle will advance microbiome engineering and spur new discoveries of the basic scientific principles governing microbiome function. In this Review, we present key elements of an iterative DBTL cycle for microbiome engineering, focusing on generalizable approaches, including top-down and bottom-up design processes, synthetic and self-assembled construction methods, and emerging tools to analyse microbiome function. These approaches can be used to harness microbiomes for broad applications related to medicine, agriculture, energy and the environment. We also discuss key challenges and opportunities of each approach and synthesize them into best practice guidelines for engineering microbiomes. We anticipate that adoption of a DBTL framework will rapidly advance microbiome-based biotechnologies aimed at improving human and animal health, agriculture and enabling the bioeconomy.

Trade-offs between gene expression, growth and phenotypic diversity in microbial populations

•

Limitations in molecular resources for gene expression influence bacterial physiology.

•

Bacteria optimise trade-offs between resource allocation and growth.

•

Resource allocation plays a role in the emergence of phenotypic heterogeneity.

•

Trade-offs between bet-hedging and growth can be harnessed in biotechnology.

Bacterial cells have a limited number of resources that can be allocated for gene expression. The intracellular competition for these resources has an impact on the cell physiology. Bacteria have evolved mechanisms to optimize resource allocation in a variety of scenarios, showing a trade-off between the resources used to maximise growth (e.g. ribosome synthesis) and the rest of cellular functions. Limitations in gene expression also play a role in generating phenotypic diversity, which is advantageous in fluctuating environments, at the expenses of decreasing growth rates. Our current understanding of these trade-offs can be exploited for biotechnological applications benefiting from the selective manipulation of the allocation of resources.

Energy metabolism controls phenotypes by protein efficiency and allocation

Cells require energy for growth and maintenance and have evolved to have multiple pathways to produce energy in response to varying conditions. A basic question in this context is how cells organize energy metabolism, which is, however, challenging to elucidate due to its complexity, i.e., the energy-producing pathways overlap with each other and even intertwine with biomass formation pathways. Here, we propose a modeling concept that decomposes energy metabolism into biomass formation and ATP-producing pathways. The latter can be further decomposed into a high-yield and a low-yield pathway. This enables independent estimation of protein efficiency for each pathway. With this concept, we modeled energy metabolism for Escherichia coli and Saccharomyces cerevisiae and found that the high-yield pathway shows lower protein efficiency than the low-yield pathway. Taken together with a fixed protein constraint, we predict overflow metabolism in E. coli and the Crabtree effect in S. cerevisiae, meaning that energy metabolism is sufficient to explain the metabolic switches. The static protein constraint is supported by the findings that protein mass of energy metabolism is conserved across conditions based on absolute proteomics data. This also suggests that enzymes may have decreased saturation or activity at low glucose uptake rates. Finally, our analyses point out three ways to improve growth, i.e., increasing protein allocation to energy metabolism, decreasing ATP demand, or increasing activity for key enzymes.

Generation of an E. coli platform strain for improved sucrose utilization using adaptive laboratory evolution

Background

Sucrose is an attractive industrial carbon source due to its abundance and the fact that it can be cheaply generated from sources such as sugarcane. However, only a few characterized Escherichia coli strains are able to metabolize sucrose, and those that can are typically slow growing or pathogenic strains.

Methods

To generate a platform strain capable of efficiently utilizing sucrose with a high growth rate, adaptive laboratory evolution (ALE) was utilized to evolve engineered E. coli K-12 MG1655 strains containing the sucrose utilizing csc genes (cscB, cscK, cscA) alongside the native sucrose consuming E. coli W.

Results

Evolved K-12 clones displayed an increase in growth and sucrose uptake rates of 1.72- and 1.40-fold on sugarcane juice as compared to the original engineered strains, respectively, while E. coli W clones showed a 1.4-fold increase in sucrose uptake rate without a significant increase in growth rate. Whole genome sequencing of evolved clones and populations revealed that two genetic regions were frequently mutated in the K-12 strains; the global transcription regulatory genes rpoB and rpoC, and the metabolic region related to a pyrimidine biosynthetic deficiency in K-12 attributed to pyrE expression. These two mutated regions have been characterized to confer a similar benefit when glucose is the main carbon source, and reverse engineering revealed the same causal advantages on M9 sucrose. Additionally, the most prevalent mutation found in the evolved E. coli W lineages was the inactivation of the cscR gene, the transcriptional repression of sucrose uptake genes.

Conclusion

The generated K-12 and W platform strains, and the specific sets of mutations that enable their phenotypes, are available as valuable tools for sucrose-based industrial bioproduction in the facile E. coli chassis.

Electronic supplementary material

The online version of this article (10.1186/s12934-019-1165-2) contains supplementary material, which is available to authorized users.