The Entner-Doudoroff pathway empowers Pseudomonas putida KT2440 with a high tolerance to oxidative stress

Improving 5-(hydroxymethyl)furfural (HMF) tolerance of Pseudomonas taiwanensis VLB120 by automated adaptive laboratory evolution (ALE)

The aldehyde 5-(hydroxymethyl)furfural (HMF) is of great importance for a circular bioeconomy. It is a renewable platform chemical that can be converted into a range of useful compounds to replace petroleum-based products such as the green plastic monomer 2,5-furandicarboxylic acid (FDCA). However, it also exhibits microbial toxicity for example hindering the efficient biotechnological valorization of lignocellulosic hydrolysates. Thus, there is an urgent need for tolerance-improved organisms applicable to whole-cell biocatalysis. Here, we engineer an oxidation-deficient derivative of the naturally robust and emerging biotechnological workhorse P. taiwanensis VLB120 by robotics-assisted adaptive laboratory evolution (ALE). The deletion of HMF-oxidizing enzymes enabled for the first time evolution under constant selection pressure by the aldehyde, yielding strains with consistently improved growth characteristics in presence of the toxicant. Genome sequencing of evolved clones revealed loss-of function mutations in the LysR-type transcriptional regulator-encoding mexT preventing expression of the associated efflux pump mexEF-oprN. This knowledge allowed reverse engineering of strains with enhanced aldehyde tolerance, even in a background of active or overexpressed HMF oxidation machinery, demonstrating a synergistic effect of two distinct tolerance mechanisms.

Chemotaxis-mediated degradation of PAHs and heterocyclic PAHs under low-temperature stress by Pseudomonas fluorescens S01: Insights into the mechanisms of biodegradation and cold adaptation

As wellknown persistent contaminants, polycyclic aromatic hydrocarbons (PAHs) and heterocyclic polyaromatic hydrocarbons (Heterocyclic PAHs)'s fates in cryogenic environments are remains uncertain. Herein, strain S01 was identified as Pseudomonas fluorescens, a novel bacterium tolerant to low temperature and capable of degrading PAHs and heterocyclic PAHs. Strain S01 exhibited growth at 5-40 ℃ and degradation rate of mixed PAHs and heterocyclic PAHs reached 52% under low-temperature. Through comprehensive metabolomic, genomic, and transcriptomic analyses, we reconstructed the biodegradation pathway for PAHs and heterocyclic PAHs in S01 while investigating its response to low temperature. Further experiments involving deletion and replacement of methyl-accepting chemotaxis protein (MCP) confirmed its crucial role in enabling strain S01's adaptation to dual stress of low temperature and pollutants. Additionally, our analysis revealed that MCP was upregulated under cold stress which enhanced strain S01's motility capabilities leading to increased biofilm formation. The establishment of biofilm promoted preservation of distinct cellular membrane stability, thereby enhancing energy metabolism. Consequently, this led to heightened efficiency in pollutant degradation and improved cold resistance capabilities. Our findings provide a comprehensive understanding of the environmental fate of both PAHs and heterocyclic PAHs under low-temperature conditions while also shedding light on cold adaptation mechanism employed by strain S01.

Engineering Pseudomonas putida KT2440 for Dipicolinate Production via the Entner-Doudoroff Pathway

Dipicolinic acid (DPA), a cyclic diacid, has garnered significant interest due to its potential applications in antimicrobial agents, antioxidants, chelating reagents, and polymer precursors. However, its natural bioproduction is limited since DPA is only accumulated in Bacillus and Clostridium species during sporulation. Thus, heterologous production by engineered strains is of paramount importance for developing a sustainable biological route for DPA production. Pseudomonas putida KT2440 has emerged as a promising host for the production of various chemicals thanks to its robustness, metabolic versatility, and genetic tractability. The dominant Entner-Doudoroff (ED) pathway for glucose metabolism in this strain offers an ideal route for DPA production due to the advantage of NADPH generation and the naturally balanced flux between glyceraldehyde-3-phosphate and pyruvate, which are both precursors for DPA synthesis. In this study, DPA production via the ED pathway was in silico designed in P. putida KT2440. The systematically engineered strain produced dipicolinate with a titer of 11.72 g/L from glucose in a 5 L fermentor. This approach not only provides a sustainable green route for DPA production but also expands our understanding of the metabolic potential of the ED pathway in P. putida KT2440.

The etherase system of Novosphingobium sp. MBES04 functions as a sensor of lignin fragments through phenylpropanone production to induce specific transcriptional responses

The MBES04 strain of Novosphingobium accumulates phenylpropanone monomers as end-products of the etherase system, which specifically and reductively cleaves the β-O-4 ether bond (a major bond in lignin molecules). However, it does not utilise phenylpropanone monomers as an energy source. Here, we studied the response to the lignin-related perturbation to clarify the physiological significance of its etherase system. Transcriptome analysis revealed two gene clusters, each consisting of four tandemly linked genes, specifically induced by a lignin preparation extracted from hardwood (Eucalyptus globulus) and a β-O-4-type lignin model biaryl compound, but not by vanillin. The most strongly induced gene was a 2,4'-dihydroxyacetophenone dioxygenase-like protein, which leads to energy production through oxidative degradation. The other cluster was related to multidrug resistance. The former cluster was transcriptionally regulated by a common promoter, where a phenylpropanone monomer acted as one of the effectors responsible for gene induction. These results indicate that the physiological significance of the etherase system of the strain lies in its function as a sensor for lignin fragments. This may be a survival strategy to detect nutrients and gain tolerance to recalcitrant toxic compounds, while the strain preferentially utilises easily degradable aromatic compounds with lower energy demands for catabolism.

Characterization of the Entner-Doudoroff pathway in Pseudomonas aeruginosa catheter-associated urinary tract infections

Pseudomonas aeruginosa is an opportunistic nosocomial pathogen responsible for a subset of catheter-associated urinary tract infections (CAUTI). In a murine model of P. aeruginosa CAUTI, we previously demonstrated that urea within urine suppresses quorum sensing and induces the Entner-Doudoroff (E-D) pathway. The E-D pathway consists of the genes zwf, pgl, edd, and eda. Zwf and Pgl convert glucose-6-phosphate into 6-phosphogluconate. Edd hydrolyzes 6-phosphogluconate to 2-keto-3-deoxy-6-phosphogluconate (KDPG). Finally, Eda cleaves KDPG to glyceraldehyde-3-phosphate and pyruvate, which enters the citric acid cycle. Here, we generated in-frame E-D mutants in the strain PA14 and assessed their growth phenotypes on chemically defined and complex media. These E-D mutants have a growth defect when grown on glucose or gluconate as the sole carbon source, which is similar to results previously reported for PAO1 mutants lacking E-D genes. RNA-sequencing following short exposure to urine revealed minimal gene regulation differences compared to the wild type. In a murine CAUTI model, virulence testing of E-D mutants revealed that two mutants lacking zwf and pgl showed minor fitness defects. Infection with the ∆pgl strain exhibited a 20% increase in host survival, and the ∆zwf strain displayed decreased colonization of the catheter and kidneys. Consequently, our findings suggest that the E-D pathway in P. aeruginosa is dispensable in this model of CAUTI. IMPORTANCE Prior studies have shown that the Entner-Doudoroff pathway is up-regulated when Pseudomonas aeruginosa is grown in urine. Pseudomonads use the Entner-Doudoroff (E-D) pathway to metabolize glucose instead of glycolysis, which led us to ask whether this pathway is required for urinary tract infection. Here, single-deletion mutants of each gene in the pathway were tested for growth on chemically defined media with single-carbon sources as well as complex media. The effect of each mutant on global gene expression in laboratory media and urine was characterized. The virulence of these mutants in a murine model of catheter-associated urinary tract infection revealed that these mutants had similar levels of colonization indicating that glucose is not the primary carbon source utilized in the urinary tract.

Genome-scale metabolic network model and phenome of solvent-tolerant Pseudomonas putida S12

BACKGROUND

Pseudomonas putida S12 is a gram-negative bacterium renowned for its high tolerance to organic solvents and metabolic versatility, making it attractive for various applications, including bioremediation and the production of aromatic compounds, bioplastics, biofuels, and value-added compounds. However, a metabolic model of S12 has yet to be developed.

RESULTS

In this study, we present a comprehensive and highly curated genome-scale metabolic network model of S12 (iSH1474), containing 1,474 genes, 1,436 unique metabolites, and 2,938 metabolic reactions. The model was constructed by leveraging existing metabolic models and conducting comparative analyses of genomes and phenomes. Approximately 2,000 different phenotypes were measured for S12 and its closely related KT2440 strain under various nutritional and environmental conditions. These phenotypic data, combined with the reported experimental data, were used to refine and validate the reconstruction. Model predictions quantitatively agreed well with in vivo flux measurements and the batch cultivation of S12, which demonstrated that iSH1474 accurately represents the metabolic capabilities of S12. Furthermore, the model was simulated to investigate the maximum theoretical metabolic capacity of S12 growing on toxic organic solvents.

CONCLUSIONS

iSH1474 represents a significant advancement in our understanding of the cellular metabolism of P. putida S12. The combined results of metabolic simulation and comparative genome and phenome analyses identified the genetic and metabolic determinants of the characteristic phenotypes of S12. This study could accelerate the development of this versatile organism as an efficient cell factory for various biotechnological applications.

Optimizing Hexose Utilization Pathways of Cupriavidus necator for Improving Growth and L-Alanine Production under Heterotrophic and Autotrophic Conditions

Cupriavidus necator is a versatile microbial chassis to produce high-value products. Blocking the poly-β-hydroxybutyrate synthesis pathway (encoded by the phaC1AB1 operon) can effectively enhance the production of C. necator, but usually decreases cell density in the stationary phase. To address this problem, we modified the hexose utilization pathways of C. necator in this study by implementing strategies such as blocking the Entner-Doudoroff pathway, completing the phosphopentose pathway by expressing the gnd gene (encoding 6-phosphogluconate dehydrogenase), and completing the Embden-Meyerhof-Parnas pathway by expressing the pfkA gene (encoding 6-phosphofructokinase). During heterotrophic fermentation, the OD600 of the phaC1AB1-knockout strain increased by 44.8% with pfkA gene expression alone, and by 93.1% with gnd and pfkA genes expressing simultaneously. During autotrophic fermentation, gnd and pfkA genes raised the OD600 of phaC1AB1-knockout strains by 19.4% and 12.0%, respectively. To explore the effect of the pfkA gene on the production of C. necator, an alanine-producing C. necator was constructed by expressing the NADPH-dependent L-alanine dehydrogenase, alanine exporter, and knocking out the phaC1AB1 operon. The alanine-producing strain had maximum alanine titer and yield of 784 mg/L and 11.0%, respectively. And these values were significantly improved to 998 mg/L and 13.4% by expressing the pfkA gene. The results indicate that completing the Embden-Meyerhof-Parnas pathway by expressing the pfkA gene is an effective method to improve the growth and production of C. necator.

Diamide-based screening method for the isolation of improved oxidative stress tolerance phenotypes in Bacillus mutant libraries

IMPORTANCE

During their life cycle, bacteria are exposed to a range of different stresses that need to be managed appropriately in order to ensure their growth and viability. This applies not only to bacteria in their natural habitats but also to bacteria employed in biotechnological production processes. Oxidative stress is one of these stresses that may originate either from bacterial metabolism or external factors. In biotechnological settings, it is of critical importance that production strains are resistant to oxidative stresses. Accordingly, this also applies to the major industrial cell factory Bacillus subtilis. In the present study, we, therefore, developed a screen for B. subtilis strains with enhanced oxidative stress tolerance. The results show that our approach is feasible and time-, space-, and resource-efficient. We, therefore, anticipate that it will enhance the development of more robust industrial production strains with improved robustness under conditions of oxidative stress.

Characterization of the Entner-Douderoff Pathway in Pseudomonas aeruginosa Catheter-associated Urinary Tract Infections

Pseudomonas aeruginosa is an opportunistic nosocomial pathogen responsible for catheter-associated urinary tract infections (CAUTI). In a murine model of P. aeruginosa CAUTI, we previously demonstrated that urea within urine suppresses quorum sensing and induces the Entner-Douderoff (E-D) pathway. The E-D pathway consists of the genes zwf, pgl, edd, and eda. Zwf and Pgl convert glucose-6-phosphate into 6-phosphogluconate. Edd hydrolyzes 6-phosphogluconate to 2-keto-3-deoxy-6-phosphogluconate (KDPG). Finally, Eda cleaves KDPG to glyceraldehyde-3-phosphate and pyruvate, which enters the citric acid cycle. Here, we generated in-frame E-D mutants in strain PA14 and assessed their growth phenotypes on chemically defined media. These E-D mutants have a growth defect when grown on glucose or gluconate as sole carbon source which are similar to results previously reported for PAO1 mutants lacking E-D genes. RNA-sequencing following short exposure to urine revealed minimal gene regulation differences compared to the wild type. In a murine CAUTI model, virulence testing of E-D mutants revealed that two mutants lacking zwf and pgl showed minor fitness defects. Infection with the ∆pgl strain exhibited a 20% increase in host survival, and the ∆zwf strain displayed decreased colonization of the catheter and kidneys. Consequently, our findings suggest that the E-D pathway in P. aeruginosa is dispensable in this model of CAUTI.

Biofilm formation stabilizes metabolism in a Roseobacteraceae bacterium under temperature increase

Ocean warming profoundly impacts microbes in marine environments; yet, how lifestyle (e.g., free living versus biofilm associated) affects the bacterial response to rising temperature is not clear. Here, we compared transcriptional, enzymatic, and physiological responses of free-living and biofilm-associated Leisingera aquaemixtae M597, a member of the Roseobacteraceae family isolated from marine biofilms, to the increase in temperature from 25℃ to 31℃. Complete genome sequencing and metagenomics revealed the prevalence of M597 in global ocean biofilms. Transcriptomics suggested a significant effect on the expression of genes related to carbohydrate metabolism, nitrogen and sulfur metabolism, and phosphorus utilization of free-living M597 cells due to temperature increase, but such drastic alterations were not observed in its biofilms. In the free-living state, the transcription of the key enzyme participating in the Embden-Meyerhof-Parnas pathway was significantly increased due to the increase in temperature, accompanied by a substantial decrease in the Entner-Doudoroff pathway, but transcripts of these glycolytic enzymes in biofilm-forming strains were independent of the temperature variation. The correlation between the growth condition and the shift in glycolytic pathways under temperature change was confirmed by enzymatic activity assays. Furthermore, the rising temperature affected the growth rate and the production of intracellular reactive oxygen species when M597 cells were free living rather than in biofilms. Thus, biofilm formation stabilizes metabolism in M597 when grown under high temperature and this homeostasis is probably related to the glycolytic pathways.IMPORTANCEBiofilm formation is one of the most successful strategies employed by microbes against environmental fluctuations. In this study, using a marine Roseobacteraceae bacterium, we studied how biofilm formation affects the response of marine bacteria to the increase in temperature. This study enhances our understanding of the function of bacterial biofilms and the microbe-environment interactions in the framework of global climate change.

Time-resolved, deuterium-based fluxomics uncovers the hierarchy and dynamics of sugar processing by Pseudomonas putida

Pseudomonas putida, a microbial host widely adopted for metabolic engineering, processes glucose through convergent peripheral pathways that ultimately yield 6-phosphogluconate. The periplasmic gluconate shunt (PGS), composed by glucose and gluconate dehydrogenases, sequentially transforms glucose into gluconate and 2-ketogluconate. Although the secretion of these organic acids by P. putida has been extensively recognized, the mechanism and spatiotemporal regulation of the PGS remained elusive thus far. To address this challenge, we adopted a dynamic 13C- and 2H-metabolic flux analysis strategy, termed D-fluxomics. D-fluxomics demonstrated that the PGS underscores a highly dynamic metabolic architecture in glucose-dependent batch cultures of P. putida, characterized by hierarchical carbon uptake by the PGS throughout the cultivation. Additionally, we show that gluconate and 2-ketogluconate accumulation and consumption can be solely explained as a result of the interplay between growth rate-coupled and decoupled metabolic fluxes. As a consequence, the formation of these acids in the PGS is inversely correlated to the bacterial growth rate-unlike the widely studied overflow metabolism of Escherichia coli and yeast. Our findings, which underline survival strategies of soil bacteria thriving in their natural environments, open new avenues for engineering P. putida towards efficient, sugar-based bioprocesses.

RB-TnSeq identifies genetic targets for improved tolerance of Pseudomonas putida towards compounds relevant to lignin conversion

Lignin-derived mixtures intended for bioconversion commonly contain high concentrations of aromatic acids, aliphatic acids, and salts. The inherent toxicity of these chemicals places a significant bottleneck upon the effective use of microbial systems for the valorization of these mixtures. Pseudomonas putida KT2440 can tolerate stressful quantities of several lignin-related compounds, making this bacterium a promising host for converting these chemicals to valuable bioproducts. Nonetheless, further increasing P. putida tolerance to chemicals in lignin-rich substrates has the potential to improve bioprocess performance. Accordingly, we employed random barcoded transposon insertion sequencing (RB-TnSeq) to reveal genetic determinants in P. putida KT2440 that influence stress outcomes during exposure to representative constituents found in lignin-rich process streams. The fitness information obtained from the RB-TnSeq experiments informed engineering of strains via deletion or constitutive expression of several genes. Namely, ΔgacAS, ΔfleQ, ΔlapAB, ΔttgR::P_tac:ttgABC, P_tac:PP_1150:PP_1152, ΔrelA, and ΔPP_1430 mutants showed growth improvement in the presence of single compounds, and some also exhibited greater tolerance when grown using a complex chemical mixture representative of a lignin-rich chemical stream. Overall, this work demonstrates the successful implementation of a genome-scale screening tool for the identification of genes influencing stress tolerance against notable compounds within lignin-enriched chemical streams, and the genetic targets identified herein offer promising engineering targets for improving feedstock tolerance in lignin valorization strains of P. putida KT2440.

Fitness-Conditional Genes for Soil Adaptation in the Bioaugmentation Agent Pseudomonas veronii 1YdBTEX2

Strain inoculation (bioaugmentation) is a potentially useful technology to provide microbiomes with new functionalities. However, there is limited understanding of the genetic factors contributing to successful establishment of inoculants. This work aimed to characterize the genes implicated in proliferation of the monoaromatic compound-degrading Pseudomonas veronii 1YdBTEX2 in nonsterile polluted soils. We generated two independent mutant libraries by random minitransposon-delivered marker insertion followed by deep sequencing (Tn-seq) with a total of 5.0 × 10⁵ unique insertions. Libraries were grown in multiple successive cycles for up to 50 generations either in batch liquid medium or in two types of soil microcosms with different resident microbial content (sand or silt) in the presence of toluene. Analysis of gene insertion abundances at different time points (passed generations of metapopulation growth), in comparison to proportions at start and to in silico generated randomized insertion distributions, allowed to define ~800 essential genes common to both libraries and ~2,700 genes with conditional fitness effects in either liquid or soil (195 of which resulted in fitness gain). Conditional fitness genes largely overlapped among all growth conditions but affected approximately twice as many functions in liquid than in soil. This indicates soil to be a more promiscuous environment for mutant growth, probably because of additional nutrient availability. Commonly depleted genes covered a wide range of biological functions and metabolic pathways, such as inorganic ion transport, fatty acid metabolism, amino acid biosynthesis, or nucleotide and cofactor metabolism. Only sparse gene sets were uncovered whose insertion caused fitness decrease exclusive for soils, which were different between silt and sand. Despite detectable higher resident bacteria and potential protist predatory counts in silt, we were, therefore, unable to detect any immediately obvious candidate genes affecting P. veronii biological competitiveness. In contrast to liquid growth conditions, mutants inactivating flagella biosynthesis and motility consistently gained strong fitness advantage in soils and displayed higher growth rates than wild type. In conclusion, although many gene functions were found to be important for growth in soils, most of these are not specific as they affect growth in liquid minimal medium more in general. This indicates that P. veronii does not need major metabolic reprogramming for proliferation in soil with accessible carbon and generally favorable growth conditions. IMPORTANCE Restoring damaged microbiomes is still a formidable challenge. Classical widely adopted approaches consist of augmenting communities with pure or mixed cultures in the hope that these display their intended selected properties under in situ conditions. Ecological theory, however, dictates that introduction of a nonresident microbe is unlikely to lead to its successful proliferation in a foreign system such as a soil microbiome. In an effort to study this systematically, we used random transposon insertion scanning to identify genes and possibly, metabolic subsystems, that are crucial for growth and survival of a bacterial inoculant (Pseudomonas veronii) for targeted degradation of monoaromatic compounds in contaminated nonsterile soils. Our results indicate that although many gene functions are important for proliferation in soil, they are general factors for growth and not exclusive for soil. In other words, P. veronii is a generalist that is not a priori hindered by the soil for its proliferation and would make a good bioaugmentation candidate.

Production of selenium nanoparticles occurs through an interconnected pathway of sulphur metabolism and oxidative stress response in Pseudomonas putida KT2440

The soil bacterium Pseudomonas putida KT2440 has been shown to produce selenium nanoparticles aerobically from selenite; however, the molecular actors involved in this process are unknown. Here, through a combination of genetic and analytical techniques, we report the first insights into selenite metabolism in this bacterium. Our results suggest that the reduction of selenite occurs through an interconnected metabolic network involving central metabolic reactions, sulphur metabolism, and the response to oxidative stress. Genes such as sucA, D2HGDH and PP_3148 revealed that the 2-ketoglutarate and glutamate metabolism is important to convert selenite into selenium. On the other hand, mutations affecting the activity of the sulphite reductase decreased the bacteria's ability to transform selenite. Other genes related to sulphur metabolism (ssuEF, sfnCE, sqrR, sqr and pdo2) and stress response (gqr, lsfA, ahpCF and sadI) were also identified as involved in selenite transformation. Interestingly, suppression of genes sqrR, sqr and pdo2 resulted in the production of selenium nanoparticles at a higher rate than the wild-type strain, which is of biotechnological interest. The data provided in this study brings us closer to understanding the metabolism of selenium in bacteria and offers new targets for the development of biotechnological tools for the production of selenium nanoparticles.

Glucose-6-Phosphate Dehydrogenase, ZwfA, a Dual Cofactor-Specific Isozyme Is Predominantly Involved in the Glucose Metabolism of Pseudomonas bharatica CSV86T

Glucose-6-phosphate dehydrogenase (Zwf) is an important enzyme in glucose metabolism via the Entner-Doudoroff pathway and the first enzyme in the oxidative pentose-phosphate pathway. It generates NAD(P)H during the conversion of glucose-6-phosphate (G6P) to 6-phosphogluconolactone, thus aiding in anabolic processes, energy yield, and oxidative stress responses. Pseudomonas bharatica CSV86T preferentially utilized aromatic compounds over glucose and exhibited a significantly lower growth rate on glucose (0.24 h-1) with a prolonged lag phase (~10 h). In strain CSV86T, glucose was metabolized via the intracellular phosphorylative route only because it lacked an oxidative (gluconate and 2-ketogluconate) route. The genome harbored three genes zwfA, zwfB, and zwfC encoding three Zwf isozymes. The present study aimed to understand gene arrangement, gene expression profiling, and molecular and kinetic properties of the purified enzymes to unveil their physiological significance in the strain CSV86T. The zwfA was found to be a part of the zwfA-pgl-eda operon, which was proximal to other glucose transport and metabolic clusters. The zwfB was found to be arranged as a gnd-zwfB operon, while zwfC was present independently. Among the three, zwfA was transcribed maximally, and the purified ZwfA displayed the highest catalytic efficiency, cooperativity with respect to G6P, and dual cofactor specificity. Isozymes ZwfB and ZwfC were NADP+-preferring and NADP+-specific, respectively. Among other functionally characterized Zwfs, ZwfA from strain CSV86T displayed poor catalytic efficiency and the further absence of oxidative routes of glucose metabolism reflected its lower growth rate on glucose compared to P. putida KT2440 and could be probable reasons for the unique carbon source utilization hierarchy. IMPORTANCE Pseudomonas bharatica CSV86T metabolizes glucose exclusively via the intracellular phosphorylative Entner-Doudoroff pathway leading the entire glucose flux through Zwf as the strain lacks oxidative routes. This may lead to limiting the concentration of downstream metabolic intermediates. The strain CSV86T possesses three isoforms of glucose-6-phosphate dehydrogenase, ZwfA, ZwfB, and ZwfC. The expression profile and kinetic properties of purified enzymes will help to understand glucose metabolism. Isozyme ZwfA dominated in terms of expression and displayed cooperativity with dual cofactor specificity. ZwfB preferred NADP+, and ZwfC was NADP+ specific, which may aid in redox cofactor balance. Such beneficial metabolic flexibility facilitated the regulation of metabolic pathways giving survival/fitness advantages in dynamic environments. Additionally, multiple genes allowed the distribution of function among these isoforms where the primary function was allocated to one of the isoforms.

Dynamic and single cell characterization of a CRISPR-interference toolset in Pseudomonas putida KT2440 for β-ketoadipate production from p-coumarate

Pseudomonas putida KT2440 is a well-studied bacterium for the conversion of lignin-derived aromatic compounds to bioproducts. The development of advanced genetic tools in P. putida has reduced the turnaround time for hypothesis testing and enabled the construction of strains capable of producing various products of interest. Here, we evaluate an inducible CRISPR-interference (CRISPRi) toolset on fluorescent, essential, and metabolic targets. Nuclease-deficient Cas9 (dCas9) expressed with the arabinose (8K)-inducible promoter was shown to be tightly regulated across various media conditions and when targeting essential genes. In addition to bulk growth data, single cell time lapse microscopy was conducted, which revealed intrinsic heterogeneity in knockdown rate within an isoclonal population. The dynamics of knockdown were studied across genomic targets in exponentially-growing cells, revealing a universal 1.75 ± 0.38 h quiescent phase after induction where 1.5 ± 0.35 doublings occur before a phenotypic response is observed. To demonstrate application of this CRISPRi toolset, β-ketoadipate, a monomer for performance-advantaged nylon, was produced at a 4.39 ± 0.5 g/L and yield of 0.76 ± 0.10 mol/mol from p-coumarate, a hydroxycinnamic acid that can be derived from grasses. These cultivation metrics were achieved by using the higher strength IPTG (1K)-inducible promoter to knockdown the pcaIJ operon in the βKA pathway during early exponential phase. This allowed the majority of the carbon to be shunted into the desired product while eliminating the need for a supplemental carbon and energy source to support growth and maintenance.

•

Developed an inducible dCas9-based CRISPR interference toolset in Pseudomonas putida KT2440.

•

Characterized single-cell dynamics of fluorescent and essential gene knockdown.

•

Applied the toolset for glucose-free production of β-ketoadipate from p-coumarate.

•

Produced β-ketoadipate at titer of 4.39 ± 0.5 g/L and 0.76 ± 0.10 mol/mol yield.

Intake of slow-digesting carbohydrates is related to changes in the microbiome and its functional pathways in growing rats with obesity induced by diet

INTRODUCTION

The main cause of insulin resistance in childhood is obesity, which contributes to future comorbidities as in adults. Although high-calorie diets and lack of exercise contribute to metabolic disease development, food quality rather than the quantity of macronutrients is more important than food density. The purpose of the present study was to examine the effects of changing the quality of carbohydrates from rapidly to slowly digestible carbohydrates on the composition of the gut microbiota and the profiles of the functional pathways in growing rats with obesity due to a high-fat diet (HFD).

METHODS

During the course of 4 weeks, rats growing on an HFD-containing carbohydrates with different digestive rates were fed either HFD-containing carbohydrates with a rapid digestion rate (OBE group) or HFD-containing carbohydrates with a slow digestion rate (OBE-ISR group). A non-obese group (NOB) was included as a reference, and rats were fed on a rodent standard diet (AIN93G). An analysis of gut microbiota was conducted using 16S rRNA-based metagenomics; a linear mixed-effects model (LMM) was used to determine changes in abundance between baseline and 4 weeks of treatment, and functional pathways were identified. Gut microbiota composition at bacterial diversity and relative abundance, at phylum and genus levels, and functional profiles were analyzed by integrating the Integrated Microbial Genomes (IMG) database.

RESULTS

The groups showed comparable gut microbiota at baseline. At the end of the treatment, animals from the ISR group exhibited differences at the phylum levels by decreasing the diversity of Fisher's index and Firmicutes (newly named as Bacillota), and increasing the Pielou's evenness and Bacteroidetes (newly named as Bacteroidota); at the genus level by increasing Alistipes, Bifidobacterium, Bacteroides, Butyricimonas, Lachnoclostridium, Flavonifractor, Ruminiclostridium 5, and Faecalibaculum and decreasing Muribaculum, Blautia, and Ruminiclostridium 9. Remarkably, relative abundances of genera Tyzzerella and Angelakisella were higher in the OBE group compared to NOB and OBE-ISR groups. In addition, some microbiota carbohydrate metabolism pathways such as glycolysis, glucuronic acid degradation, pentose phosphate pathway, methanogenesis, and fatty acid biosynthesis exhibited increased activity in the OBE-ISR group after the treatment. Higher levels of acetate and propionate were found in the feces of the ISR group compared with the NOB and OBE groups.

CONCLUSION

The results of this study demonstrate that replacing rapidly digestible carbohydrates with slowly digestible carbohydrates within an HFD improve the composition of the gut microbiota. Consequently, metabolic disturbances associated with obesity may be prevented.

Dissolved iron released from nanoscale zero-valent iron (nZVI) activates the defense system in bacterium Pseudomonas putida, leading to high tolerance to oxidative stress

Nanoscale zero-valent iron (nZVI) has increasingly been applied to remediate aquifers polluted by organochlorines or heavy metals. As a result, bacteria in the vicinity of remediate action can be stressed by surplus iron released from nZVI. However, the understanding of the iron stress defense pathways during this process is currently incomplete. Therefore, we aimed to elucidate the physiological and transcriptomic response of the bacterium, Pseudomonas putida NCTC 10936, to 100 mg/L of nZVI and 44.5 µg/L of dissolved iron obtained from nZVI suspension. Cell viability was neither affected by nZVI nor dissolved iron, although the dissolved iron caused stress that altered the cell physiology and caused the generation of smaller cells, whereas cells were elongated in the presence of nZVI. Transcriptomic analysis confirmed the observed stronger physiological effect caused by dissolved iron (in total 3839 differentially expressed genes [DEGs]) than by nZVI (945 DEGs). Dissolved iron (but not nZVI) activated genes involved in oxidative stress-related pathways, antioxidant activity, carbohydrate and energy metabolism, but downregulated genes associated with flagellar assembly proteins and two-component systems involved in sensing external stimuli. As a result, bacteria very effectively faced oxidative insults and cell viability was not affected.

Comparative Genome Analysis of a Novel Alkaliphilic Actinobacterial Species Nesterenkonia haasae

In the present study, a comparative genome analysis of the novel alkaliphilic actinobacterial Nesterenkonia haasae with other members of the genus Nesterenkonia was performed. The genome size of Nesterenkonia members ranged from 2,188,008 to 3,676,111 bp. N. haasae and Nesterenkonia members of the present study encode the essential glycolysis and pentose phosphate pathway genes. In addition, some Nesterenkonia members encode the crucial genes for Entner-Doudoroff pathways. Some Nesterenkonia members possess the genes responsible for sulfate/thiosulfate transport system permease protein/ ATP-binding protein and conversion of sulfate to sulfite. Nesterenkonia members also encode the genes for assimilatory nitrate reduction, nitrite reductase, and the urea cycle. All Nesterenkonia members have the genes to overcome environmental stress and produce secondary metabolites. The present study helps to understand N. haasae and Nesterenkonia members' environmental adaptation and niches specificity based on their specific metabolic properties. Further, based on genome analysis, we propose reclassifying Nesterenkonia jeotgali as a later heterotypic synonym of Nesterenkonia sandarakina.

A synthetic C2 auxotroph of Pseudomonas putida for evolutionary engineering of alternative sugar catabolic routes

Acetyl-coenzyme A (AcCoA) is a metabolic hub in virtually all living cells, serving as both a key precursor of essential biomass components and a metabolic sink for catabolic pathways for a large variety of substrates. Owing to this dual role, tight growth-production coupling schemes can be implemented around the AcCoA node. Building on this concept, a synthetic C2 auxotrophy was implemented in the platform bacterium Pseudomonas putida through an in silico-informed engineering approach. A growth-coupling strategy, driven by AcCoA demand, allowed for direct selection of an alternative sugar assimilation route-the phosphoketolase (PKT) shunt from bifidobacteria. Adaptive laboratory evolution forced the synthetic P. putida auxotroph to rewire its metabolic network to restore C2 prototrophy via the PKT shunt. Large-scale structural chromosome rearrangements were identified as possible mechanisms for adjusting the network-wide proteome profile, resulting in improved PKT-dependent growth phenotypes. ¹³C-based metabolic flux analysis revealed an even split between the native Entner-Doudoroff pathway and the synthetic PKT bypass for glucose processing, leading to enhanced carbon conservation. These results demonstrate that the P. putida metabolism can be radically rewired to incorporate a synthetic C2 metabolism, creating novel network connectivities and highlighting the importance of unconventional engineering strategies to support efficient microbial production.

TCA cycle enhancement and uptake of monomeric substrates support growth of marine Roseobacter at low temperature

Members of the marine Roseobacter group are ubiquitous in global oceans, but their cold-adaptive strategies have barely been studied. Here, as represented by Loktanella salsilacus strains enriched in polar regions, we firstly characterized the metabolic features of a cold-adapted Roseobacter by multi-omics, enzyme activities, and carbon utilization procedures. Unlike in most cold-adapted microorganisms, the TCA cycle is enhanced by accumulating more enzyme molecules, whereas genes for thiosulfate oxidation, sulfate reduction, nitrate reduction, and urea metabolism are all expressed at lower abundance when L. salsilacus was growing at 5 °C in comparison with higher temperatures. Moreover, a carbon-source competition experiment has evidenced the preferential use of glucose rather than sucrose at low temperature. This selective utilization is likely to be controlled by the carbon source uptake and transformation steps, which also reflects an economic calculation balancing energy production and functional plasticity. These findings provide a mechanistic understanding of how a Roseobacter member and possibly others as well counteract polar constraints.

The metabolic adaptation of Loktanella salsilacus strains to cold involves an increase of enzymes involved in the TCA cycle and preferential use of glucose rather than sucrose at low temperature, providing insights into how Roseobacter adapts in polar regions.

Transcriptomic analysis of the response of Photobacterium phosphoreum and Photobacterium carnosum to co-contaminants on chicken meat

This study investigated the impact of Brochothrix (B.) thermosphacta and Pseudomonas (Ps.) fragi on the transcriptomes of Photobacterium (P.) phosphoreum and P. carnosum on chicken meat under modified atmosphere (MA) and air atmosphere (AA). P. phosphoreum TMW2.2103 responded to MA with a reduced transcript number related to cell division and an enhanced number related to oxidative stress. Concomitantly, the analysis revealed upregulation of fermentation and downregulation of respiration. It predicts enhanced substrate competition in presence of co-contaminants/MA. In contrast, the strain upregulated the respiration in AA, supposably due to improved substrate accessibility in this situation. For P. carnosum TMW2.2149 the respiration was downregulated, and the pyruvate metabolism upregulated under MA. MA/co-contaminant resulted in multiple upregulated metabolic routes. Conversely, AA/co-contaminant resulted only in minor regulations, showing inability to cope with fast growing competitors. Observations reveal different strategies of photobacteria to react to co-contaminants on meat.

Differential Genetic Strategies of Burkholderia vietnamiensis and Paraburkholderia kururiensis for Root Colonization of Oryza sativa subsp. japonica and O. sativa subsp. indica , as Revealed by Transposon Mutagenesis Sequencing

Burkholderiaceae are frequent and abundant colonizers of the rice rhizosphere and interesting candidates to investigate for growth promotion. Species of Paraburkholderia have repeatedly been described to stimulate plant growth.

Comparative Genomics Reveals Genetic Diversity and Metabolic Potentials of the Genus Qipengyuania and Suggests Fifteen Novel Species

Members of the genus Qipengyuania are heterotrophic bacteria frequently isolated from marine environments with great application potential in areas such as carotenoid production. However, the genomic diversity, metabolic function, and adaption of this genus remain largely unclear. Here, 16 isolates related to the genus Qipengyuania were recovered from coastal samples and their genomes were sequenced. The phylogenetic inference of these isolates and reference type strains of this genus indicated that the 16S rRNA gene was insufficient to distinguish them at the species level; instead, the phylogenomic reconstruction could provide the reliable phylogenetic relationships and confirm 15 new well-supported branches, representing 15 putative novel genospecies corroborated by the digital DNA-DNA hybridization and average nucleotide identity analyses. Comparative genomics revealed that the genus Qipengyuania had an open pangenome and possessed multiple conserved genes and pathways related to metabolic functions and environmental adaptation, despite the presence of divergent genomic features and specific metabolic potential. Genetic analysis and pigment detection showed that the members of this genus were identified as carotenoid producers, while some proved to be potentially aerobic anoxygenic photoheterotrophs. Collectively, the first insight into the genetic diversity and metabolic potentials of the genus Qipengyuania will contribute to better understanding of the speciation and adaptive evolution in natural environments.

IMPORTANCE The deciphering of the phylogenetic diversity and metabolic features of the abundant bacterial taxa is critical for exploring their ecological importance and application potential. Qipengyuania is a genus of frequently isolated heterotrophic microorganisms with great industrial application potential. Numerous strains related to the genus Qipengyuania have been isolated from diverse environments, but their genomic diversity and metabolic functions remain unclear. Our study revealed a high degree of genetic diversity, metabolic versatility, and environmental adaptation of the genus Qipengyuania using comparative genomics. Fifteen novel species of this genus have been established using a polyphasic taxonomic approach, expanding the number of described species to almost double. This study provided an overall view of the genus Qipengyuania at the genomic level and will enable us to better uncover its ecological roles and evolutionary history.

A study of the flexibility of the carbon catabolic pathways of extremophilic P. aeruginosa san ai exposed to benzoate versus glucose as sole carbon sources by multi omics analytical platform

Polyextremophilic, hydrocarbonoclastic Pseudomonas aeruginosa san ai can survive under extreme environmental challenges in the presence of a variety of pollutants such as organic solvents and hydrocarbons, particularly aromatics, heavy metals, and high pH. To date, the metabolic plasticity of the extremophilic P. aeruginosa, has not been sufficiently studied in regard to the effect of changing carbon sources. Therefore, the present study explores the carbon metabolic pathways of polyextremophilic P. aeruginosa san ai grown on sodium benzoate versus glucose and its potential for aromatic degradation. P. aeruginosa san ai removed/metabolised nearly 430 mg/L of benzoate for 48 h, demonstrating a high capacity for aromatic degradation. Comparative functional proteomics, targeted metabolomics and genomics analytical approaches were employed to study the carbon metabolism of the P. aeruginosa san ai. Functional proteomic study of selected enzymes participating in the β-ketoadipate and the Entner-Doudoroff pathways revealed a metabolic reconfiguration induced by benzoate compared to glucose. Metabolome analysis implied the existence of both catechol and protocatechuate branches of the β-ketoadipate pathway. Enzymatic study of benzoate grown cultures confirmed the activity of the ortho- catechol branch of the β-ketoadipate pathway. Even high concentrations of benzoate did not show increased stress protein synthesis, testifying to its extremophilic nature capable of surviving in harsh conditions. This ability of Pseudomonas aeruginosa san ai to efficiently degrade benzoate can provide a wide range of use of this strain in environmental and agricultural application.

Dual Effect: High NADH Levels Contribute to Efflux-Mediated Antibiotic Resistance but Drive Lethality Mediated by Reactive Oxygen Species

In light of the antibiotic crisis, emerging strategies to sensitize bacteria to available antibiotics should be explored. Several studies on the mechanisms of killing suggest that bactericidal antibiotic activity is enforced through the generation of reactive oxygen species (ROS-lethality hypothesis). Here, we artificially manipulated the redox homeostasis of the model opportunistic pathogen Pseudomonas aeruginosa using specific enzymes that catalyze either the formation or oxidation of NADH. Increased NADH levels led to the activation of antibiotic efflux pumps and high levels of antibiotic resistance. However, higher NADH levels also resulted in increased intracellular ROS and amplified antibiotic killing. Our results demonstrate that growth inhibition and killing activity are mediated via different mechanisms. Furthermore, the profound changes in bioenergetics produced low-virulence phenotypes characterized by reduced interbacterial signaling controlled pathogenicity traits. Our results pave the way for a more effective infection resolution and add an antivirulence strategy to maximize chances to combat devastating P. aeruginosa infections while reducing the overall use of antibiotics. IMPORTANCE The emergence of antibiotic resistance has become one of the major threats to public health. A better understanding of antimicrobial killing mechanisms promises to uncover new ways to resensitize bacteria to commonly used antibiotics. In this context, there is increasing evidence that the metabolic status of the cell plays a fundamental role in reactive oxygen species (ROS)-mediated cell death. In this work, we artificially manipulated the redox balance in Pseudomonas aeruginosa by the expression of two orthologous enzymes. We found that the increase of intracellular NADH concentrations leads to higher antibiotic resistance but also generates a burst in the production of ROS that amplified antimicrobial killing. Our work suggests that the combination of bactericidal antibiotics with agents that disturb the cellular redox homeostasis could significantly enhance antibiotic killing via sensitization of pathogens to currently available antibiotics.

Synthetic Control of Metabolic States in Pseudomonas putida by Tuning Polyhydroxyalkanoate Cycle

Polyhydroxyalkanoates (PHAs) are polyesters produced by numerous microorganisms for energy and carbon storage. Simultaneous synthesis and degradation of PHA drives a dynamic cycle linked to the central carbon metabolism, which modulates numerous and diverse bacterial processes, such as stress endurance, pathogenesis, and persistence. Here, we analyze the role of the PHA cycle in conferring robustness to the model bacterium P. putida KT2440. To assess the effect of this cycle in the cell, we began by constructing a PHA depolymerase (PhaZ) mutant strain that had its PHA cycle blocked. We then restored the flux through the cycle in the context of an engineered library of P. putida strains harboring differential levels of PhaZ. High-throughput phenotyping analyses of this collection of strains revealed significant changes in response to PHA cycle performance impacting cell number and size, PHA accumulation, and production of extracellular (R)-hydroxyalkanoic acids. To understand the metabolic changes at the system level due to PHA turnover, we contextualized these physiological data using the genome-scale metabolic model iJN1411. Model-based predictions suggest successive metabolic steady states during the growth curve and an important carbon flux rerouting driven by the activity of the PHA cycle. Overall, we demonstrate that modulating the activity of the PHA cycle gives us control over the carbon metabolism of P. putida, which in turn will give us the ability to tailor cellular mechanisms driving stress tolerance, e.g., defenses against oxidative stress, and any potential biotechnological applications. IMPORTANCE Despite large research efforts devoted to understanding the flexible metabolism of Pseudomonas beyond the role of key regulatory players, the metabolic basis powering the dynamic control of its biological fitness under disturbance conditions remains largely unknown. Among other metabolic hubs, the so-called PHA cycle, involving simultaneous synthesis and degradation of PHAs, is emerging as a pivotal metabolic trait powering metabolic robustness and resilience in this bacterial group. Here, we provide evidence suggesting that metabolic states in Pseudomonas can be anticipated, controlled, and engineered by tailoring the flux through the PHA cycle. Overall, our study suggests that the PHA cycle is a promising metabolic target toward achieving control over bacterial metabolic robustness. This is likely to open up a broad range of applications in areas as diverse as pathogenesis and biotechnology.

The OxyR and SoxR transcriptional regulators are involved in a broad oxidative stress response in Paraburkholderia xenovorans LB400

Background

Aerobic metabolism generates reactive oxygen species that may cause critical harm to the cell. The aim of this study is the characterization of the stress responses in the model aromatic-degrading bacterium Paraburkholderia xenovorans LB400 to the oxidizing agents paraquat and H₂O₂.

Methods

Antioxidant genes were identified by bioinformatic methods in the genome of P. xenovorans LB400, and the phylogeny of its OxyR and SoxR transcriptional regulators were studied. Functionality of the transcriptional regulators from strain LB400 was assessed by complementation with LB400 SoxR of null mutant P. aeruginosa ΔsoxR, and the construction of P. xenovorans pIZoxyR that overexpresses OxyR. The effects of oxidizing agents on P. xenovorans were studied measuring bacterial susceptibility, survival and ROS formation after exposure to paraquat and H₂O₂. The effects of these oxidants on gene expression (qRT-PCR) and the proteome (LC–MS/MS) were quantified.

Results

P. xenovorans LB400 possesses a wide repertoire of genes for the antioxidant defense including the oxyR, ahpC, ahpF, kat, trxB, dpsA and gorA genes, whose orthologous genes are regulated by the transcriptional regulator OxyR in E. coli. The LB400 genome also harbors the soxR, fumC, acnA, sodB, fpr and fldX genes, whose orthologous genes are regulated by the transcriptional regulator SoxR in E. coli. The functionality of the LB400 soxR gene was confirmed by complementation of null mutant P. aeruginosa ΔsoxR. Growth, susceptibility, and ROS formation assays revealed that LB400 cells were more susceptible to paraquat than H₂O₂. Transcriptional analyses indicated the upregulation of the oxyR, ahpC1, katE and ohrB genes in LB400 cells after exposure to H₂O₂, whereas the oxyR, fumC, ahpC1, sodB1 and ohrB genes were induced in presence of paraquat. Proteome analysis revealed that paraquat induced the oxidative stress response proteins AhpCF and DpsA, the universal stress protein UspA and the RNA chaperone CspA. Both oxidizing agents induced the Ohr protein, which is involved in organic peroxide resistance. Notably, the overexpression of the LB400 oxyR gene in P. xenovorans significantly decreased the ROS formation and the susceptibility to paraquat, suggesting a broad OxyR-regulated antioxidant response.

Conclusions

This study showed that P. xenovorans LB400 possess a broad range oxidative stress response, which explain the high resistance of this strain to the oxidizing compounds paraquat and H₂O₂.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40659-022-00373-7.

Methylotrophs and Hydrocarbon-Degrading Bacteria Are Key Players in the Microbial Community of an Abandoned Century-Old Oil Exploration Well

In this work, we studied the microbial community and the physicochemical conditions prevailing in an exploratory oil well, abandoned a century ago, located in the Cahuita National Park (Costa Rica). According to our analysis, Cahuita well is characterized by a continuous efflux of methane and the presence of a mixture of hydrocarbons including phenanthrene/anthracene, fluoranthene, pyrene, dibenzothiophene, tricyclic terpanes, pyrene, sesquiterpenes, sterane, and n-alkanes. Based on the analysis of 16S rRNA gene amplicons, we detected a significant abundance of methylotrophic bacteria such as Methylobacillus (6.3-26.0% of total reads) and Methylococcus (4.1-30.6%) and the presence of common genera associated with hydrocarbon degradation, such as Comamonas (0.8-4.6%), Hydrogenophaga (1.5-3.3%) Rhodobacter (1.0-4.9%), and Flavobacterium (1.1-6.5%). The importance of C1 metabolism in this niche was confirmed by amplifying the methane monooxygenase (MMO)-encoding gene (pmo) from environmental DNA and the isolation of two strains closely related to Methylorubrum rhodesianum and Paracoccus communis with the ability to growth using methanol and formate as sole carbon source respectively. In addition, we were able to isolated 20 bacterial strains from the genera Pseudomonas, Acinetobacter, and Microbacterium which showed the capability to grow using the hydrocarbons detected in the oil well as sole carbon source. This work describes the physicochemical properties and microbiota of an environment exposed to hydrocarbons for 100 years, and it not only represents a contribution to the understanding of microbial communities in environments with permanently high concentrations of these compounds but also has biotechnological implications for bioremediation of petroleum-polluted sites.

Analogous Metabolic Decoupling in Pseudomonas putida and Comamonas testosteroni Implies Energetic Bypass to Facilitate Gluconeogenic Growth

Gluconeogenic carbon metabolism is not well understood, especially within the context of flux partitioning between energy generation and biomass production, despite the importance of gluconeogenic carbon substrates in natural and engineered carbon processing. Here, using multiple omics approaches, we elucidate the metabolic mechanisms that facilitate gluconeogenic fast-growth phenotypes in Pseudomonas putida and Comamonas testosteroni, two Proteobacteria species with distinct metabolic networks. In contrast to the genetic constraint of C. testosteroni, which lacks the enzymes required for both sugar uptake and a complete oxidative pentose phosphate (PP) pathway, sugar metabolism in P. putida is known to generate surplus NADPH by relying on the oxidative PP pathway within its characteristic cyclic connection between the Entner-Doudoroff (ED) and Embden-Meyerhoff-Parnas (EMP) pathways. Remarkably, similar to the genome-based metabolic decoupling in C. testosteroni, our 13C-fluxomics reveals an inactive oxidative PP pathway and disconnected EMP and ED pathways in P. putida during gluconeogenic feeding, thus requiring transhydrogenase reactions to supply NADPH for anabolism in both species by leveraging the high tricarboxylic acid cycle flux during gluconeogenic growth. Furthermore, metabolomics and proteomics analyses of both species during gluconeogenic feeding, relative to glycolytic feeding, demonstrate a 5-fold depletion in phosphorylated metabolites and the absence of or up to a 17-fold decrease in proteins of the PP and ED pathways. Such metabolic remodeling, which is reportedly lacking in Escherichia coli exhibiting a gluconeogenic slow-growth phenotype, may serve to minimize futile carbon cycling while favoring the gluconeogenic metabolic regime in relevant proteobacterial species. IMPORTANCE Glycolytic metabolism of sugars is extensively studied in the Proteobacteria, but gluconeogenic carbon sources (e.g., organic acids, amino acids, aromatics) that feed into the tricarboxylic acid (TCA) cycle are widely reported to produce a fast-growth phenotype, particularly in species with biotechnological relevance. Much remains unknown about the importance of glycolysis-associated pathways in the metabolism of gluconeogenic carbon substrates. Here, we demonstrate that two distinct proteobacterial species, through genetic constraints or metabolic regulation at specific metabolic nodes, bypass the oxidative PP pathway during gluconeogenic growth and avoid unnecessary carbon fluxes by depleting protein investment into connected glycolysis pathways. Both species can leverage instead the high TCA cycle flux during gluconeogenic feeding to meet NADPH demand. Importantly, lack of a complete oxidative pentose phosphate pathway is a widespread metabolic trait in Proteobacteria with a gluconeogenic carbon preference, thus highlighting the important relevance of our findings toward elucidating the metabolic architecture in these bacteria.

Metabolic Model of the Nitrogen-Fixing Obligate Aerobe Azotobacter vinelandii Predicts Its Adaptation to Oxygen Concentration and Metal Availability

There is considerable interest in promoting biological nitrogen fixation (BNF) as a mechanism to reduce the inputs of nitrogenous fertilizers in agriculture, but considerable fundamental knowledge gaps still need to be addressed. BNF is catalyzed by nitrogenase, which requires a large input of energy in the form of ATP and low potential electrons. Diazotrophs that respire aerobically have an advantage in meeting the ATP demands of BNF but face challenges in protecting nitrogenase from inactivation by oxygen. Here, we constructed a genome-scale metabolic model of the nitrogen-fixing bacterium Azotobacter vinelandii, which uses a complex respiratory protection mechanism to consume oxygen at a high rate to keep intracellular conditions microaerobic. Our model accurately predicts growth rate under high oxygen and substrate concentrations, consistent with a large electron flux directed to the respiratory protection mechanism. While a partially decoupled electron transport chain compensates for some of the energy imbalance under high-oxygen conditions, it does not account for all substrate intake, leading to increased maintenance rates. Interestingly, the respiratory protection mechanism is required for accurate predictions even when ammonia is supplemented during growth, suggesting that the respiratory protection mechanism might be a core principle of metabolism and not just used for nitrogenase protection. We have also shown that rearrangement of flux through the electron transport system allows A. vinelandii to adapt to different oxygen concentrations, metal availability, and genetic disruption, which cause an ammonia excretion phenotype. Accurately determining the energy balance in an aerobic nitrogen-fixing metabolic model is required for future engineering approaches.

The potential of Pseudomonas for bioremediation of oxyanions

Non-metal, metal and metalloid oxyanions occur naturally in minerals and rocks of the Earth's crust and are mostly found in low concentrations or confined in specific regions of the planet. However, anthropogenic activities including urban development, mining, agriculture, industrial activities and new technologies have increased the release of oxyanions to the environment, which threatens the sustainability of natural ecosystems, in turn affecting human development. For these reasons, the implementation of new methods that could allow not only the remediation of oxyanion contaminants but also the recovery of valuable elements from oxyanions of the environment is imperative. From this perspective, the use of microorganisms emerges as a strategy complementary to physical, mechanical and chemical methods. In this review, we discuss the opportunities that the Pseudomonas genus offers for the bioremediation of oxyanions, which is derived from its specialized central metabolism and the high number of oxidoreductases present in the genomes of these bacteria. Finally, we review the current knowledge on the transport and metabolism of specific oxyanions in Pseudomonas species. We consider that the Pseudomonas genus is an excellent starting point for the development of biotechnological approaches for the upcycling of oxyanions into added-value metal and metalloid byproducts.

Genome analysis of the metabolically versatile Pseudomonas umsongensis GO16: the genetic basis for PET monomer upcycling into polyhydroxyalkanoates

The throwaway culture related to the single-use materials such as polyethylene terephthalate (PET) has created a major environmental concern. Recycling of PET waste into biodegradable plastic polyhydroxyalkanoate (PHA) creates an opportunity to improve resource efficiency and contribute to a circular economy. We sequenced the genome of Pseudomonas umsongensis GO16 previously shown to convert PET-derived terephthalic acid (TA) into PHA and performed an in-depth genome analysis. GO16 can degrade a range of aromatic substrates in addition to TA, due to the presence of a catabolic plasmid pENK22. The genetic complement required for the degradation of TA via protocatechuate was identified and its functionality was confirmed by transferring the tph operon into Pseudomonas putida KT2440, which is unable to utilize TA naturally. We also identified the genes involved in ethylene glycol (EG) metabolism, the second PET monomer, and validated the capacity of GO16 to use EG as a sole source of carbon and energy. Moreover, GO16 possesses genes for the synthesis of both medium and short chain length PHA and we have demonstrated the capacity of the strain to convert mixed TA and EG into PHA. The metabolic versatility of GO16 highlights the potential of this organism for biotransformations using PET waste as a feedstock.

Oxidative Stress Response in Pseudomonas aeruginosa

Pseudomonas aeruginosa is a Gram-negative environmental and human opportunistic pathogen highly adapted to many different environmental conditions. It can cause a wide range of serious infections, including wounds, lungs, the urinary tract, and systemic infections. The high versatility and pathogenicity of this bacterium is attributed to its genomic complexity, the expression of several virulence factors, and its intrinsic resistance to various antimicrobials. However, to thrive and establish infection, P. aeruginosa must overcome several barriers. One of these barriers is the presence of oxidizing agents (e.g., hydrogen peroxide, superoxide, and hypochlorous acid) produced by the host immune system or that are commonly used as disinfectants in a variety of different environments including hospitals. These agents damage several cellular molecules and can cause cell death. Therefore, bacteria adapt to these harsh conditions by altering gene expression and eliciting several stress responses to survive under oxidative stress. Here, we used PubMed to evaluate the current knowledge on the oxidative stress responses adopted by P. aeruginosa. We will describe the genes that are often differently expressed under oxidative stress conditions, the pathways and proteins employed to sense and respond to oxidative stress, and how these changes in gene expression influence pathogenicity and the virulence of P. aeruginosa. Understanding these responses and changes in gene expression is critical to controlling bacterial pathogenicity and developing new therapeutic agents.

Pseudomonas putida KT2440 endures temporary oxygen limitations

The obligate aerobic nature of Pseudomonas putida, one of the most prominent whole-cell biocatalysts emerging for industrial bioprocesses, questions its ability to be cultivated in large-scale bioreactors, which exhibit zones of low dissolved oxygen tension. P. putida KT2440 was repeatedly subjected to temporary oxygen limitations in scale-down approaches to assess the effect on growth and an exemplary production of rhamnolipids. At those conditions, the growth and production of P. putida KT2440 were decelerated compared to well-aerated reference cultivations, but remarkably, final biomass and rhamnolipid titers were similar. The robust growth behavior was confirmed across different cultivation systems, media compositions, and laboratories, even when P. putida KT2440 was repeatedly exposed to dual carbon and oxygen starvation. Quantification of the nucleotides ATP, ADP, and AMP revealed a decrease of intracellular ATP concentrations with increasing duration of oxygen starvation, which can, however, be restored when re-supplied with oxygen. Only small changes in the proteome were detected when cells encountered oscillations in dissolved oxygen tensions. Concluding, P. putida KT2440 appears to be able to cope with repeated oxygen limitations as they occur in large-scale bioreactors, affirming its outstanding suitability as a whole-cell biocatalyst for industrial-scale bioprocesses.

Metaproteomics Reveals Alteration of the Gut Microbiome in Weaned Piglets Due to the Ingestion of the Mycotoxins Deoxynivalenol and Zearalenone

The ingestion of mycotoxins can cause adverse health effects and represents a severe health risk to humans and livestock. Even though several acute and chronic effects have been described, the effect on the gut metaproteome is scarcely known. For that reason, we used metaproteomics to evaluate the effect of the mycotoxins deoxynivalenol (DON) and zearalenone (ZEN) on the gut microbiome of 15 weaned piglets. Animals were fed for 28 days with feed contaminated with different concentrations of DON (DONlow: 870 μg DON/kg feed, DONhigh: 2493 μg DON/kg feed) or ZEN (ZENlow: 679 μg ZEN/kg feed, ZENhigh: 1623 μg ZEN/kg feed). Animals in the control group received uncontaminated feed. The gut metaproteome composition in the high toxin groups shifted compared to the control and low mycotoxin groups, and it was also more similar among high toxin groups. These changes were accompanied by the increase in peptides belonging to Actinobacteria and a decrease in peptides belonging to Firmicutes. Additionally, DONhigh and ZENhigh increased the abundance of proteins associated with the ribosomes and pentose-phosphate pathways, while decreasing glycolysis and other carbohydrate metabolism pathways. Moreover, DONhigh and ZENhigh increased the abundance of the antioxidant enzyme thioredoxin-dependent peroxiredoxin. In summary, the ingestion of DON and ZEN altered the abundance of different proteins associated with microbial metabolism, genetic processing, and oxidative stress response, triggering a disruption in the gut microbiome structure.

Model-guided dynamic control of essential metabolic nodes boosts acetyl-coenzyme A-dependent bioproduction in rewired Pseudomonas putida

Pseudomonas putida is evolutionarily endowed with features relevant for bioproduction, especially under harsh operating conditions. The rich metabolic versatility of this species, however, comes at the price of limited formation of acetyl-coenzyme A (CoA) from sugar substrates. Since acetyl-CoA is a key metabolic precursor for a number of added-value products, in this work we deployed an in silico-guided rewiring program of central carbon metabolism for upgrading P. putida as a host for acetyl-CoA-dependent bioproduction. An updated kinetic model, integrating fluxomics and metabolomics datasets in addition to manually-curated information of enzyme mechanisms, identified targets that would lead to increased acetyl-CoA levels. Based on these predictions, a set of plasmids based on clustered regularly interspaced short palindromic repeats (CRISPR) and dead CRISPR-associated protein 9 (dCas9) was constructed to silence genes by CRISPR interference (CRISPRi). Dynamic reduction of gene expression of two key targets (gltA, encoding citrate synthase, and the essential accA gene, encoding subunit A of the acetyl-CoA carboxylase complex) mediated an 8-fold increase in the acetyl-CoA content of rewired P. putida. Poly(3-hydroxybutyrate) (PHB) was adopted as a proxy of acetyl-CoA availability, and two synthetic pathways were engineered for biopolymer accumulation. By including cell morphology as an extra target for the CRISPRi approach, fully rewired P. putida strains programmed for PHB accumulation had a 5-fold increase in PHB titers in bioreactor cultures using glucose. Thus, the strategy described herein allowed for rationally redirecting metabolic fluxes in P. putida from central metabolism towards product biosynthesis-especially relevant when deletion of essential pathways is not an option.

Molecular cloning and biochemical characterization of a NAD-dependent sorbitol dehydrogenase from cold-adapted Pseudomonas mandelii

Sugar alcohols (polyols) have important roles as nutrients, anti-freezing agents and scavengers of free radicals in cold-adapted bacteria, but the characteristics of polyol dehydrogenases in cold-adapted bacteria remain largely unknown. In this study, based on the observation that a cold-adapted bacterium Pseudomonas mandelii JR-1 predominantly utilized d-sorbitol as its carbon source, among the four polyols examined (d-galactitol, d-mannitol, d-sorbitol and d-xylitol), we cloned and characterized a sorbitol dehydrogenase (SDH, EC 1.1.1.14) belonging to the short-chain dehydrogenase/reductase family from this bacterium (the SDH hereafter referred to as PmSDH). PmSDH contained Asn111, Ser140, Tyr153 and Lys157 as catalytic active site residues and existed as an ∼67-kDa dimer in size-exclusion chromatography. PmSDH converted d-sorbitol to d-fructose using nicotinamide adenine dinucleotide (NAD+) as a cofactor and, vice versa, d-fructose to d-sorbitol using nicotinamide adenine dinucleotide reduced (NADH) as a cofactor. PmSDH maintained its conformational flexibility, secondary and tertiary structures, and thermal stability at 4-25°C. These results indicate that PmSDH, which has a flexible structure and a high catalytic activity at colder temperatures, is well suited to sorbitol utilization in the cold-adapted bacterium P. mandelii JR-1.

The anoxic electrode-driven fructose catabolism of Pseudomonas putida KT2440

Pseudomonas putida (P. putida) is a microorganism of interest for various industrial processes, yet its strictly aerobic nature limits application. Despite previous attempts to adapt P. putida to anoxic conditions via genetic engineering or the use of a bioelectrochemical system (BES), the problem of energy shortage and internal redox imbalance persists. In this work, we aimed to provide the cytoplasmic metabolism with different monosaccharides, other than glucose, and explored the physiological response in P. putida KT2440 during bioelectrochemical cultivation. The periplasmic oxidation cascade was found to be able to oxidize a wide range of aldoses to their corresponding (keto-)aldonates. Unexpectedly, isomerization of the ketose fructose to mannose also enabled oxidation by glucose dehydrogenase, a new pathway uncovered for fructose metabolism in P. putida KT2440 in BES. Besides the isomerization, the remainder of fructose was imported into the cytoplasm and metabolized. This resulted in a higher NADPH/NADP+ ratio, compared to glucose. Comparative proteomics further revealed the upregulation of proteins in the lower central carbon metabolism during the experiment. These findings highlight that the choice of a substrate in BES can target cytosolic and periplasmic oxidation pathways, and that electrode-driven redox balancing can drive these pathways in P. putida under anaerobic conditions.

Micro-aerobic production of isobutanol with engineered Pseudomonas putida

Pseudomonas putida KT2440 is emerging as a promising microbial host for biotechnological industry due to its broad range of substrate affinity and resilience to physicochemical stresses. Its natural tolerance towards aromatics and solvents qualifies this versatile microbe as promising candidate to produce next generation biofuels such as isobutanol. In this study, we scaled-up the production of isobutanol with P. putida from shake flask to fed-batch cultivation in a 30 L bioreactor. The design of a two-stage bioprocess with separated growth and production resulted in 3.35 gisobutanol L-1. Flux analysis revealed that the NADPH expensive formation of isobutanol exceeded the cellular catabolic supply of NADPH finally causing growth retardation. Concomitantly, the cell counteracted to the redox imbalance by increased formation of 2-ketogluconic thereby providing electrons for the respiratory ATP generation. Thus, P. putida partially uncoupled ATP formation from the availability of NADH. The quantitative analysis of intracellular pyridine nucleotides NAD(P)+ and NAD(P)H revealed elevated catabolic and anabolic reducing power during aerobic production of isobutanol. Additionally, the installation of micro-aerobic conditions during production doubled the integral glucose-to-isobutanol conversion yield to 60 mgisobutanol gglucose -1 while preventing undesired carbon loss as 2-ketogluconic acid.

Genome-Wide Metabolic Reconstruction of the Synthesis of Polyhydroxyalkanoates from Sugars and Fatty Acids by Burkholderia Sensu Lato Species

Burkholderia sensu lato (s.l.) species have a versatile metabolism. The aims of this review are the genomic reconstruction of the metabolic pathways involved in the synthesis of polyhydroxyalkanoates (PHAs) by Burkholderia s.l. genera, and the characterization of the PHA synthases and the pha genes organization. The reports of the PHA synthesis from different substrates by Burkholderia s.l. strains were reviewed. Genome-guided metabolic reconstruction involving the conversion of sugars and fatty acids into PHAs by 37 Burkholderia s.l. species was performed. Sugars are metabolized via the Entner-Doudoroff (ED), pentose-phosphate (PP), and lower Embden-Meyerhoff-Parnas (EMP) pathways, which produce reducing power through NAD(P)H synthesis and PHA precursors. Fatty acid substrates are metabolized via β-oxidation and de novo synthesis of fatty acids into PHAs. The analysis of 194 Burkholderia s.l. genomes revealed that all strains have the phaC, phaA, and phaB genes for PHA synthesis, wherein the phaC gene is generally present in ≥2 copies. PHA synthases were classified into four phylogenetic groups belonging to class I II and III PHA synthases and one outlier group. The reconstruction of PHAs synthesis revealed a high level of gene redundancy probably reflecting complex regulatory layers that provide fine tuning according to diverse substrates and physiological conditions.

Complete Genome Sequence and Function Gene Identify of Prometryne-Degrading Strain Pseudomonas sp. DY-1

The genus Pseudomonas is widely recognized for its potential for environmental remediation and plant growth promotion. Pseudomonas sp. DY-1 was isolated from the agricultural soil contaminated five years by prometryne, it manifested an outstanding prometryne degradation efficiency and an untapped potential for plant resistance improvement. Thus, it is meaningful to comprehend the genetic background for strain DY-1. The whole genome sequence of this strain revealed a series of environment adaptive and plant beneficial genes which involved in environmental stress response, heavy metal or metalloid resistance, nitrate dissimilatory reduction, riboflavin synthesis, and iron acquisition. Detailed analyses presented the potential of strain DY-1 for degrading various organic compounds via a homogenized pathway or the protocatechuate and catechol branches of the β-ketoadipate pathway. In addition, heterologous expression, and high efficiency liquid chromatography (HPLC) confirmed that prometryne could be oxidized by a Baeyer-Villiger monooxygenase (BVMO) encoded by a gene in the chromosome of strain DY-1. The result of gene knock-out suggested that the sulfate starvation-induced (SSI) genes in this strain might also involve in the process of prometryne degradation. These results would provide the molecular basis for the application of strain DY-1 in various fields and would contribute to the study of prometryne biodegradation mechanism as well.

Reconfiguration of metabolic fluxes in Pseudomonas putida as a response to sub-lethal oxidative stress

As a frequent inhabitant of sites polluted with toxic chemicals, the soil bacterium and plant-root colonizer Pseudomonas putida can tolerate high levels of endogenous and exogenous oxidative stress. Yet, the ultimate reason of such phenotypic property remains largely unknown. To shed light on this question, metabolic network-wide routes for NADPH generation-the metabolic currency that fuels redox-stress quenching mechanisms-were inspected when P. putida KT2440 was challenged with a sub-lethal H₂O₂ dose as a proxy of oxidative conditions. ¹³C-tracer experiments, metabolomics, and flux analysis, together with the assessment of physiological parameters and measurement of enzymatic activities, revealed a substantial flux reconfiguration in oxidative environments. In particular, periplasmic glucose processing was rerouted to cytoplasmic oxidation, and the cyclic operation of the pentose phosphate pathway led to significant NADPH-forming fluxes, exceeding biosynthetic demands by ~50%. The resulting NADPH surplus, in turn, fueled the glutathione system for H₂O₂ reduction. These properties not only account for the tolerance of P. putida to environmental insults-some of which end up in the formation of reactive oxygen species-but they also highlight the value of this bacterial host as a platform for environmental bioremediation and metabolic engineering.

Portable bacterial CRISPR transcriptional activation enables metabolic engineering in Pseudomonas putida

CRISPR-Cas transcriptional programming in bacteria is an emerging tool to regulate gene expression for metabolic pathway engineering. Here we implement CRISPR-Cas transcriptional activation (CRISPRa) in P. putida using a system previously developed in E. coli. We provide a methodology to transfer CRISPRa to a new host by first optimizing expression levels for the CRISPRa system components, and then applying rules for effective CRISPRa based on a systematic characterization of promoter features. Using this optimized system, we regulate biosynthesis in the biopterin and mevalonate pathways. We demonstrate that multiple genes can be activated simultaneously by targeting multiple promoters or by targeting a single promoter in a multi-gene operon. This work will enable new metabolic engineering strategies in P. putida and pave the way for CRISPR-Cas transcriptional programming in other bacterial species.

Refactoring the Conjugation Machinery of Promiscuous Plasmid RP4 into a Device for Conversion of Gram-Negative Isolates to Hfr Strains

Chromosomal exchange and subsequent recombination of the cognate DNA between bacteria was one of the most useful genetic tools (e.g., Hfr strains) for genetic analyses of E. coli before the genomic era. In this paper, yeast assembly has been used to recruit the conjugation machinery of environmentally promiscuous RP4 plasmid into a minimized, synthetic construct that enables transfer of chromosomal segments between donor/recipient strains of P. putida KT2440 and potentially many other Gram-negative bacteria. The synthetic device features [i] a R6K suicidal plasmid backbone, [ii] a mini-Tn5 transposon vector, and [iii] the minimal set of genes necessary for active conjugation (RP4 Tra1 and Tra2 clusters) loaded as cargo in the mini-Tn5 mobile element. Upon insertion of the transposon in different genomic locations, the ability of P. putida-TRANS (transference of RP4-activated nucleotide segments) donor strains to mobilize genomic stretches of DNA into neighboring bacteria was tested. To this end, a P. putida double mutant ΔpyrF (uracil auxotroph) Δedd (unable to grow on glucose) was used as recipient in mating experiments, and the restoration of the pyrF+/edd+ phenotypes allowed for estimation of chromosomal transfer efficiency. Cells with the inserted transposon behaved in a manner similar to Hfr-like strains and were able to transfer up to 23% of their genome at frequencies close to 10-6 exconjugants per recipient cell. The hereby described TRANS device not only expands the molecular toolbox for P. putida, but it also enables a suite of genomic manipulations which were thus far only possible with domesticated laboratory strains and species.

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.

Comparative Proteomics Reveals the Anaerobic Lifestyle of Meat-Spoiling Pseudomonas Species

The ability of certain Pseudomonas (P.) species to grow or persist in anoxic habitats by either denitrification, acetate fermentation, or arginine fermentation has been described in several studies as a special property. Previously, we had isolated strains belonging to the species P. lundensis, P. weihenstephanensis, and P. fragi from anoxic modified atmosphere packaged (MAP) minced beef and further proved their anaerobic growth in vitro on agar plates. This follow-up study investigated the anaerobic growth of two strains per respective species in situ on inoculated chicken breast filet under 100% N₂ modified atmosphere. We were able to prove anaerobic growth of all six strains on chicken breast filet with cell division rates of 0.2–0.8/day. Furthermore, we characterized the anaerobic metabolic lifestyle of these Pseudomonas strains by comparative proteomics, upon their cultivation in meat simulation media, which were constantly gassed with either air or 100% N₂ atmospheres. From these proteomic predictions, and respective complementation by physiological experiments, we conclude that the Pseudomonas strains P. fragi, P. weihenstephanensis, P. lundensis exhibit a similar anaerobic lifestyle and employ arginine fermentation via the arginine deiminase (ADI) pathway to grow anaerobically also on MAP meats. Furthermore, glucose fermentation to ethanol via the ED-pathway is predicted to enable long term survival but no true growth, while respiratory growth with nitrate as alternative electron acceptor or glucose fermentation to acetate could be excluded due to absence of essential genes. The citric acid cycle is partially bypassed by the glyoxylate shunt, functioning as the gluconeogenetic route without production of NADH₂ under carbon limiting conditions as e.g., in packaged meats. Triggered by an altered redox balance, we also detected upregulation of enzymes involved in protein folding as well as disulfide bonds isomerization under anoxic conditions as a counteracting mechanism to reduce protein misfolding. Hence, this study reveals the mechanisms enabling anaerobic grow and persistence of common meat-spoiling Pseudomonas species, and further complements the hitherto limited knowledge of the anaerobic lifestyle of Pseudomonas species in general.

Cofactor Specificity of Glucose-6-Phosphate Dehydrogenase Isozymes in Pseudomonas putida Reveals a General Principle Underlying Glycolytic Strategies in Bacteria

Glucose-6-phosphate dehydrogenase (G6PDH) is widely distributed in nature and catalyzes the first committing step in the oxidative branch of the pentose phosphate (PP) pathway, feeding either the reductive PP or the Entner-Doudoroff pathway. Besides its role in central carbon metabolism, this dehydrogenase provides reduced cofactors, thereby affecting redox balance. Although G6PDH is typically considered to display specificity toward NADP⁺, some variants accept NAD⁺ similarly or even preferentially. Furthermore, the number of G6PDH isozymes encoded in bacterial genomes varies from none to more than four orthologues. On this background, we systematically analyzed the interplay of the three G6PDH isoforms of the soil bacterium Pseudomonas putida KT2440 from genomic, genetic, and biochemical perspectives. P. putida represents an ideal model to tackle this endeavor, as its genome harbors gene orthologues for most dehydrogenases in central carbon metabolism. We show that the three G6PDHs of strain KT2440 have different cofactor specificities and that the isoforms encoded by zwfA and zwfB carry most of the activity, acting as metabolic “gatekeepers” for carbon sources that enter at different nodes of the biochemical network. Moreover, we demonstrate how multiplication of G6PDH isoforms is a widespread strategy in bacteria, correlating with the presence of an incomplete Embden-Meyerhof-Parnas pathway. The abundance of G6PDH isoforms in these species goes hand in hand with low NADP⁺ affinity, at least in one isozyme. We propose that gene duplication and relaxation in cofactor specificity is an evolutionary strategy toward balancing the relative production of NADPH and NADH.

IMPORTANCE Protein families have likely arisen during evolution by gene duplication and divergence followed by neofunctionalization. While this phenomenon is well documented for catabolic activities (typical of environmental bacteria that colonize highly polluted niches), the coexistence of multiple isozymes in central carbon catabolism remains relatively unexplored. We have adopted the metabolically versatile soil bacterium Pseudomonas putida KT2440 as a model to interrogate the physiological and evolutionary significance of coexisting glucose-6-phosphate dehydrogenase (G6PDH) isozymes. Our results show that each of the three G6PDHs in this bacterium display distinct biochemical properties, especially at the level of cofactor preference, impacting bacterial physiology in a carbon source-dependent fashion. Furthermore, the presence of multiple G6PDHs differing in NAD⁺ or NADP⁺ specificity in bacterial species strongly correlates with their predominant metabolic lifestyle. Our findings support the notion that multiplication of genes encoding cofactor-dependent dehydrogenases is a general evolutionary strategy toward achieving redox balance according to the growth conditions.