Evolution and engineering of pathways for aromatic O-demethylation in Pseudomonas putida KT2440.

Biological conversion of lignin from biomass offers a promising strategy for sustainable production of fuels and chemicals. However, aromatic compounds derived from lignin commonly contain methoxy groups, and O-demethylation of these substrates is often a rate-limiting reaction that influences catabolic efficiency. Several enzyme families catalyze aromatic O-demethylation, but they are rarely compared in vivo to determine an optimal biocatalytic strategy. Here, two pathways for aromatic O-demethylation were compared in Pseudomonas putida KT2440. The native Rieske non-heme iron monooxygenase (VanAB) and, separately, a heterologous tetrahydrofolate-dependent demethylase (LigM) were constitutively expressed in P. putida, and the strains were optimized via adaptive laboratory evolution (ALE) with vanillate as a model substrate. All evolved strains displayed improved growth phenotypes, with the evolved strains harboring the native VanAB pathway exhibiting growth rates ∼1.8x faster than those harboring the heterologous LigM pathway. Enzyme kinetics and transcriptomics studies investigated the contribution of selected mutations toward enhanced utilization of vanillate. The VanAB-overexpressing strains contained the most impactful mutations, including those in VanB, the reductase for vanillate O-demethylase, PP_3494, a global regulator of vanillate catabolism, and fghA, involved in formaldehyde detoxification. These three mutations were combined into a single strain, which exhibited approximately 5x faster vanillate consumption than the wild-type strain in the first 8 h of cultivation. Overall, this study illuminates the details of vanillate catabolism in the context of two distinct enzymatic mechanisms, yielding a platform strain for efficient O-demethylation of lignin-related aromatic compounds to value-added products.

Selection for Cell Yield Does Not Reduce Overflow Metabolism in Escherichia coli.

Overflow metabolism is ubiquitous in nature, and it is often considered inefficient because it leads to a relatively low biomass yield per consumed carbon. This metabolic strategy has been described as advantageous because it supports high growth rates during nutrient competition. Here, we experimentally evolved bacteria without nutrient competition by repeatedly growing and mixing millions of parallel batch cultures of Escherichia coli. Each culture originated from a water-in-oil emulsion droplet seeded with a single cell. Unexpectedly we found that overflow metabolism (acetate production) did not change. Instead, the numerical cell yield during the consumption of the accumulated acetate increased as a consequence of a reduction in cell size. Our experiments and a mathematical model show that fast growth and overflow metabolism, followed by the consumption of the overflow metabolite, can lead to a higher numerical cell yield and therefore a higher fitness compared with full respiration of the substrate. This provides an evolutionary scenario where overflow metabolism can be favorable even in the absence of nutrient competition.

Revealing oxidative pentose metabolism in new Pseudomonas putida isolates.

The Pseudomonas putida group in the Gammaproteobacteria has been intensively studied for bioremediation and plant growth promotion. Members of this group have recently emerged as promising hosts to convert intermediates derived from plant biomass to biofuels and biochemicals. However, most strains of P. putida cannot metabolize pentose sugars derived from hemicellulose. Here, we describe three isolates that provide a broader view of the pentose sugar catabolism in the P. putida group. One of these isolates clusters with the well-characterized P. alloputida KT2440 (Strain BP6); the second isolate clustered with plant growth-promoting strain P. putida W619 (Strain M2), while the third isolate represents a new species in the group (Strain BP8). Each of these isolates possessed homologous genes for oxidative xylose catabolism (xylDXA) and a potential xylonate transporter. Strain M2 grew on arabinose and had genes for oxidative arabinose catabolism (araDXA). A CRISPR interference (CRISPRi) system was developed for strain M2 and identified conditionally essential genes for xylose growth. A glucose dehydrogenase was found to be responsible for initial oxidation of xylose and arabinose in strain M2. These isolates have illuminated inherent diversity in pentose catabolism in the P. putida group and may provide alternative hosts for biomass conversion.

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.

Kinetic profiling of metabolic specialists demonstrates stability and consistency of in vivo enzyme turnover numbers.

Enzyme turnover numbers (kcats) are essential for a quantitative understanding of cells. Because kcats are traditionally measured in low-throughput assays, they can be inconsistent, labor-intensive to obtain, and can miss in vivo effects. We use a data-driven approach to estimate in vivo kcats using metabolic specialist Escherichia coli strains that resulted from gene knockouts in central metabolism followed by metabolic optimization via laboratory evolution. By combining absolute proteomics with fluxomics data, we find that in vivo kcats are robust against genetic perturbations, suggesting that metabolic adaptation to gene loss is mostly achieved through other mechanisms, like gene-regulatory changes. Combining machine learning and genome-scale metabolic models, we show that the obtained in vivo kcats predict unseen proteomics data with much higher precision than in vitro kcats. The results demonstrate that in vivo kcats can solve the problem of inconsistent and low-coverage parameterizations of genome-scale cellular models.

Elucidating aromatic acid tolerance at low pH in Saccharomyces cerevisiae using adaptive laboratory evolution.

Toxicity from the external presence or internal production of compounds can reduce the growth and viability of microbial cell factories and compromise productivity. Aromatic compounds are generally toxic for microorganisms, which makes their production in microbial hosts challenging. Here we use adaptive laboratory evolution to generate Saccharomyces cerevisiae mutants tolerant to two aromatic acids, coumaric acid and ferulic acid. The evolution experiments were performed at low pH (3.5) to reproduce conditions typical of industrial processes. Mutant strains tolerant to levels of aromatic acids near the solubility limit were then analyzed by whole genome sequencing, which revealed prevalent point mutations in a transcriptional activator (Aro80) that is responsible for regulating the use of aromatic amino acids as the nitrogen source. Among the genes regulated by Aro80, ESBP6 was found to be responsible for increasing tolerance to aromatic acids by exporting them out of the cell. Further examination of the native function of Esbp6 revealed that this transporter can excrete fusel acids (byproducts of aromatic amino acid catabolism) and this role is shared with at least one additional transporter native to S. cerevisiae (Pdr12). Besides conferring tolerance to aromatic acids, ESBP6 overexpression was also shown to significantly improve the secretion in coumaric acid production strains. Overall, we showed that regulating the activity of transporters is a major mechanism to improve tolerance to aromatic acids. These findings can be used to modulate the intracellular concentration of aromatic compounds to optimize the excretion of such products while keeping precursor molecules inside the cell.

Membrane transporter identification and modulation via adaptive laboratory evolution.

Membrane transport proteins are potential targets for medical and biotechnological applications. However, more than 30% of reported membrane transporter families are either poorly characterized or lack adequate functional annotation. Here, adaptive laboratory evolution was leveraged to identify membrane transporters for a set of four amino acids as well as specific mutations that modulate the activities of these transporters. Specifically, Escherichia coli was adaptively evolved under increasing concentrations of L-histidine, L-phenylalanine, L-threonine, and L-methionine separately with multiple replicate evolutions. Evolved populations and isolated clones displayed growth rates comparable to the unstressed ancestral strain at elevated concentrations (four-to six-fold increases) of the targeted amino acids. Whole genome sequencing of the evolved strains revealed a diverse number of key mutations, including SNPs, small deletions, and copy number variants targeting the transporters leuE for histidine, yddG for phenylalanine, yedA for methionine, and brnQ and rhtC for threonine. Reverse engineering of the mutations in the ancestral strain established mutation causality of the specific mutations for the tolerant phenotypes. The functional roles of yedA and brnQ in the transport of methionine and threonine, respectively, are novel assignments and their functional roles were validated using a flow cytometry cellular accumulation assay. To demonstrate how the identified transporters can be leveraged for production, an L-phenylalanine overproduction strain was shown to be a superior producer when the identified yddG exporter was overexpressed. Overall, the results revealed the striking efficiency of laboratory evolution to identify transporters and specific mutational mechanisms to modulate their activities, thereby demonstrating promising applicability in transporter discovery efforts and strain engineering.

Machine-learning from Pseudomonas putida KT2440 transcriptomes reveals its transcriptional regulatory network.

Bacterial gene expression is orchestrated by numerous transcription factors (TFs). Elucidating how gene expression is regulated is fundamental to understanding bacterial physiology and engineering it for practical use. In this study, a machine-learning approach was applied to uncover the genome-scale transcriptional regulatory network (TRN) in Pseudomonas putida KT2440, an important organism for bioproduction. We performed independent component analysis of a compendium of 321 high-quality gene expression profiles, which were previously published or newly generated in this study. We identified 84 groups of independently modulated genes (iModulons) that explain 75.7% of the total variance in the compendium. With these iModulons, we (i) expand our understanding of the regulatory functions of 39 iModulon associated TFs (e.g., HexR, Zur) by systematic comparison with 1993 previously reported TF-gene interactions; (ii) outline transcriptional changes after the transition from the exponential growth to stationary phases; (iii) capture group of genes required for utilizing diverse carbon sources and increased stationary response with slower growth rates; (iv) unveil multiple evolutionary strategies of transcriptome reallocation to achieve fast growth rates; and (v) define an osmotic stimulon, which includes the Type VI secretion system, as coordination of multiple iModulon activity changes. Taken together, this study provides the first quantitative genome-scale TRN for P. putida KT2440 and a basis for a comprehensive understanding of its complex transcriptome changes in a variety of physiological states.

Genome-scale metabolic modeling reveals key features of a minimal gene set.

Mesoplasma florum, a fast-growing near-minimal organism, is a compelling model to explore rational genome designs. Using sequence and structural homology, the set of metabolic functions its genome encodes was identified, allowing the reconstruction of a metabolic network representing ˜ 30% of its protein-coding genes. Growth medium simplification enabled substrate uptake and product secretion rate quantification which, along with experimental biomass composition, were integrated as species-specific constraints to produce the functional iJL208 genome-scale model (GEM) of metabolism. Genome-wide expression and essentiality datasets as well as growth data on various carbohydrates were used to validate and refine iJL208. Discrepancies between model predictions and observations were mechanistically explained using protein structures and network analysis. iJL208 was also used to propose an in silico reduced genome. Comparing this prediction to the minimal cell JCVI-syn3.0 and its parent JCVI-syn1.0 revealed key features of a minimal gene set. iJL208 is a stepping-stone toward model-driven whole-genome engineering.

Engineering Pseudomonas putida for improved utilization of syringyl aromatics.

Lignin is a largely untapped source for the bioproduction of value-added chemicals. Pseudomonas putida KT2440 has emerged as a strong candidate for bioprocessing of lignin feedstocks due to its resistance to several industrial solvents, broad metabolic capabilities, and genetic amenability. Here we demonstrate the engineering of P. putida for the ability to metabolize syringic acid, one of the major products that comes from the breakdown of the syringyl component of lignin. The rational design was first applied for the construction of strain Sy-1 by overexpressing a native vanillate demethylase. Subsequent adaptive laboratory evolution (ALE) led to the generation of mutations that achieved robust growth on syringic acid as a sole carbon source. The best mutant showed a 30% increase in the growth rate over the original engineered strain. Genomic sequencing revealed multiple mutations repeated in separate evolved replicates. Reverse engineering of mutations identified in agmR, gbdR, fleQ, and the intergenic region of gstB and yadG into the parental strain recaptured the improved growth of the evolved strains to varied extent. These findings thus reveal the ability of P. putida to utilize lignin more fully as a feedstock and make it a more economically viable chassis for chemical production.

Coupling S-adenosylmethionine-dependent methylation to growth: Design and uses.

We present a selection design that couples S-adenosylmethionine-dependent methylation to growth. We demonstrate its use in improving the enzyme activities of not only N-type and O-type methyltransferases by 2-fold but also an acetyltransferase of another enzyme category when linked to a methylation pathway in Escherichia coli using adaptive laboratory evolution. We also demonstrate its application for drug discovery using a catechol O-methyltransferase and its inhibitors entacapone and tolcapone. Implementation of this design in Saccharomyces cerevisiae is also demonstrated.

Tracing compartmentalized NADPH metabolism in the cytosol and mitochondria of mammalian cells.

Eukaryotic cells compartmentalize biochemical processes in different organelles, often relying on metabolic cycles to shuttle reducing equivalents across intracellular membranes. NADPH serves as the electron carrier for the maintenance of redox homeostasis and reductive biosynthesis, with separate cytosolic and mitochondrial pools providing reducing power in each respective location. This cellular organization is critical for numerous functions but complicates analysis of metabolic pathways using available methods. Here we develop an approach to resolve NADP(H)-dependent pathways present within both the cytosol and the mitochondria. By tracing hydrogen in compartmentalized reactions that use NADPH as a cofactor, including the production of 2-hydroxyglutarate by mutant isocitrate dehydrogenase enzymes, we can observe metabolic pathway activity in these distinct cellular compartments. Using this system we determine the direction of serine/glycine interconversion within the mitochondria and cytosol, highlighting the ability of this approach to resolve compartmentalized reactions in intact cells.

Adaptive evolution reveals a tradeoff between growth rate and oxidative stress during naphthoquinone-based aerobic respiration.

Evolution fine-tunes biological pathways to achieve a robust cellular physiology. Two and a half billion years ago, rapidly rising levels of oxygen as a byproduct of blooming cyanobacterial photosynthesis resulted in a redox upshift in microbial energetics. The appearance of higher-redox-potential respiratory quinone, ubiquinone (UQ), is believed to be an adaptive response to this environmental transition. However, the majority of bacterial species are still dependent on the ancient respiratory quinone, naphthoquinone (NQ). Gammaproteobacteria can biosynthesize both of these respiratory quinones, where UQ has been associated with aerobic lifestyle and NQ with anaerobic lifestyle. We engineered an obligate NQ-dependent γ-proteobacterium, Escherichia coli ΔubiC, and performed adaptive laboratory evolution to understand the selection against the use of NQ in an oxic environment and also the adaptation required to support the NQ-driven aerobic electron transport chain. A comparative systems-level analysis of pre- and postevolved NQ-dependent strains revealed a clear shift from fermentative to oxidative metabolism enabled by higher periplasmic superoxide defense. This metabolic shift was driven by the concerted activity of 3 transcriptional regulators (PdhR, RpoS, and Fur). Analysis of these findings using a genome-scale model suggested that resource allocation to reactive oxygen species (ROS) mitigation results in lower growth rates. These results provide a direct elucidation of a resource allocation tradeoff between growth rate and ROS mitigation costs associated with NQ usage under oxygen-replete condition.

Cellular responses to reactive oxygen species are predicted from molecular mechanisms.

Catalysis using iron-sulfur clusters and transition metals can be traced back to the last universal common ancestor. The damage to metalloproteins caused by reactive oxygen species (ROS) can prevent cell growth and survival when unmanaged, thus eliciting an essential stress response that is universal and fundamental in biology. Here we develop a computable multiscale description of the ROS stress response in Escherichia coli, called OxidizeME. We use OxidizeME to explain four key responses to oxidative stress: 1) ROS-induced auxotrophy for branched-chain, aromatic, and sulfurous amino acids; 2) nutrient-dependent sensitivity of growth rate to ROS; 3) ROS-specific differential gene expression separate from global growth-associated differential expression; and 4) coordinated expression of iron-sulfur cluster (ISC) and sulfur assimilation (SUF) systems for iron-sulfur cluster biosynthesis. These results show that we can now develop fundamental and quantitative genotype-phenotype relationships for stress responses on a genome-wide basis.

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.

ALEdb 1.0: a database of mutations from adaptive laboratory evolution experimentation.

Adaptive Laboratory Evolution (ALE) has emerged as an experimental approach to discover causal mutations that confer desired phenotypic functions. ALE not only represents a controllable experimental approach to systematically discover genotype-phenotype relationships, but also allows for the revelation of the series of genetic alterations required to acquire the new phenotype. Numerous ALE studies have been published, providing a strong impetus for developing databases to warehouse experimental evolution information and make it retrievable for large-scale analysis. Here, the first step towards establishing this resource is presented: ALEdb (http://aledb.org). This initial release contains over 11 000 mutations that have been discovered from eleven ALE publications. ALEdb (i) is a web-based platform that comprehensively reports on ALE acquired mutations and their conditions, (ii) reports key mutations using previously established trends, (iii) enables a search-driven workflow to enhance user mutation functional analysis through mutation cross-reference, (iv) allows exporting of mutation query results for custom analysis, (v) includes a bibliome describing the databased experiment publications and (vi) contains experimental evolution mutations from multiple model organisms. Thus, ALEdb is an informative platform which will become increasingly revealing as the number of reported ALE experiments and identified mutations continue to expand.

Dissecting the genetic and metabolic mechanisms of adaptation to the knockout of a major metabolic enzyme in Escherichia coli.

Unraveling the mechanisms of microbial adaptive evolution following genetic or environmental challenges is of fundamental interest in biological science and engineering. When the challenge is the loss of a metabolic enzyme, adaptive responses can also shed significant insight into metabolic robustness, regulation, and areas of kinetic limitation. In this study, whole-genome sequencing and high-resolution 13C-metabolic flux analysis were performed on 10 adaptively evolved pgi knockouts of Escherichia coliPgi catalyzes the first reaction in glycolysis, and its loss results in major physiological and carbon catabolism pathway changes, including an 80% reduction in growth rate. Following adaptive laboratory evolution (ALE), the knockouts increase their growth rate by up to 3.6-fold. Through combined genomic-fluxomic analysis, we characterized the mutations and resulting metabolic fluxes that enabled this fitness recovery. Large increases in pyridine cofactor transhydrogenase flux, correcting imbalanced production of NADPH and NADH, were enabled by direct mutations to the transhydrogenase genes sthA and pntAB The phosphotransferase system component crr was also found to be frequently mutated, which corresponded to elevated flux from pyruvate to phosphoenolpyruvate. The overall energy metabolism was found to be strikingly robust, and what have been previously described as latently activated Entner-Doudoroff and glyoxylate shunt pathways are shown here to represent no real increases in absolute flux relative to the wild type. These results indicate that the dominant mechanism of adaptation was to relieve the rate-limiting steps in cofactor metabolism and substrate uptake and to modulate global transcriptional regulation from stress response to catabolism.

Causal mutations from adaptive laboratory evolution are outlined by multiple scales of genome annotations and condition-specificity.

BACKGROUND: Adaptive Laboratory Evolution (ALE) has emerged as an experimental approach to discover mutations that confer phenotypic functions of interest. However, the task of finding and understanding all beneficial mutations of an ALE experiment remains an open challenge for the field. To provide for better results than traditional methods of ALE mutation analysis, this work applied enrichment methods to mutations described by a multiscale annotation framework and a consolidated set of ALE experiment conditions. A total of 25,321 unique genome annotations from various sources were leveraged to describe multiple scales of mutated features in a set of 35 Escherichia coli based ALE experiments. These experiments totalled 208 independent evolutions and 2641 mutations. Additionally, mutated features were statistically associated across a total of 43 unique experimental conditions to aid in deconvoluting mutation selection pressures. RESULTS: Identifying potentially beneficial, or key, mutations was enhanced by seeking coding and non-coding genome features significantly enriched by mutations across multiple ALE replicates and scales of genome annotations. The median proportion of ALE experiment key mutations increased from 62%, with only small coding and non-coding features, to 71% with larger aggregate features. Understanding key mutations was enhanced by considering the functions of broader annotation types and the significantly associated conditions for key mutated features. The approaches developed here were used to find and characterize novel key mutations in two ALE experiments: one previously unpublished with Escherichia coli grown on glycerol as a carbon source and one previously published with Escherichia coli tolerized to high concentrations of L-serine. CONCLUSIONS: The emergent adaptive strategies represented by sets of ALE mutations became more clear upon observing the aggregation of mutated features across small to large scale genome annotations. The clarification of mutation selection pressures among the many experimental conditions also helped bring these strategies to light. This work demonstrates how multiscale genome annotation frameworks and data-driven methods can help better characterize ALE mutations, and thus help elucidate the genotype-to-phenotype relationship of the studied organism.

High-quality genome-scale metabolic modelling of Pseudomonas putida highlights its broad metabolic capabilities.

Genome-scale reconstructions of metabolism are computational species-specific knowledge bases able to compute systemic metabolic properties. We present a comprehensive and validated reconstruction of the biotechnologically relevant bacterium Pseudomonas putida KT2440 that greatly expands computable predictions of its metabolic states. The reconstruction represents a significant reactome expansion over available reconstructed bacterial metabolic networks. Specifically, iJN1462 (i) incorporates several hundred additional genes and associated reactions resulting in new predictive capabilities, including new nutrients supporting growth; (ii) was validated by in vivo growth screens that included previously untested carbon (48) and nitrogen (41) sources; (iii) yielded gene essentiality predictions showing large accuracy when compared with a knock-out library and Bar-seq data; and (iv) allowed mapping of its network to 82 P. putida sequenced strains revealing functional core that reflect the large metabolic versatility of this species, including aromatic compounds derived from lignin. Thus, this study provides a thoroughly updated metabolic reconstruction and new computable phenotypes for P. putida, which can be leveraged as a first step toward understanding the pan metabolic capabilities of Pseudomonas.

Adaptive laboratory evolution of Escherichia coli under acid stress.

The ability of Escherichia coli to tolerate acid stress is important for its survival and colonization in the human digestive tract. Here, we performed adaptive laboratory evolution of the laboratory strain E. coli K-12 MG1655 at pH 5.5 in glucose minimal medium. After 800 generations, six independent populations under evolution had reached 18.0 % higher growth rates than their starting strain at pH 5.5, while maintaining comparable growth rates to the starting strain at pH 7. We characterized the evolved strains and found that: (1) whole genome sequencing of isolated clones from each evolved population revealed mutations in rpoC appearing in five of six sequenced clones; and (2) gene expression profiles revealed different strategies to mitigate acid stress, which are related to amino acid metabolism and energy production and conversion. Thus, a combination of adaptive laboratory evolution, genome resequencing and expression profiling revealed, on a genome scale, the strategies that E. coli uses to mitigate acid stress.