Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “metabolic models”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

Rubisco supplies pyruvate for the 2- C -methyl-D-erythritol-4-phosphate pathway

RIBULOSE-1,5-BISPHOSPHATE CARBOXYLASE/OXYGENASE (Rubisco) produces pyruvate in the chloroplast through β-elimination of the aci-carbanion intermediate. Here we show that this side reaction supplies pyruvate for isoprenoid, fatty acid and branched-chain amino acid biosynthesis in photosynthetically active tissue. 13 C labelling studies of intact Arabidopsis plants demonstrate that the total carbon commitment to pyruvate is too large for phosphoenolpyruvate to serve as a precursor. Low oxygen stimulates Rubisco carboxylase activity and increases pyruvate production and flux through the 2-C-methyl-d-erythritol-4-phosphate (MEP) pathway, which supplies the precursors for plastidic isoprenoid biosynthesis. Metabolome analysis of mutants defective in phosphoenolpyruvate or pyruvate import and biochemical characterization of isolated chloroplasts further support Rubisco as the main source of pyruvate in chloroplasts. Seedlings incorporated exogenous, 13 C-labelled pyruvate into MEP pathway intermediates, while adult plants did not, underscoring the developmental transition in pyruvate sourcing. Rubisco β-elimination leading to pyruvate constituted 0.7% of the product profile in in vitro assays, which translates to 2% of the total carbon leaving the Calvin–Benson–Bassham cycle. These insights solve the “pyruvate paradox”, improve the fit of metabolic models for central metabolism and connect the MEP pathway directly to carbon assimilation.

Plant physiology↗

Engineering the hyperthermophilic archaeon Pyrococcus furiosus for 1-propanol production

ABSTRACT Society relies heavily on chemicals traditionally produced through the refinement of fossil fuels. The conversion of renewable biomass to value-added chemicals by microbes, particularly hyperthermophiles (T opt ≥80°C), offers a renewable alternative to this traditional approach. Herein, we describe the engineering of the hyperthermophilic archaeon Pyrococcus furiosus , which grows optimally (T opt ) at 100°C, for the conversion of sugar to 1-propanol. This was accomplished by constructing a hybrid metabolic pathway consisting of two native and seven heterologously produced enzymes to convert acetyl-CoA from carbohydrate metabolism to 1-propanol. A total of eleven foreign genes from two other organisms were utilized, one from the thermophilic bacterium Thermoanaerobacter sp. strain X514 and 10 from the thermoacidophilic archaeon Metallosphaera sedula, both of which grow optimally near 70°C. The recombinant P. furiosus strain produced 1-propanol at similar concentrations (up to ~1 mM) when incubated at 75°C to activate the gene products of Thermoanaerobacter sp. strain X514 and M. sedula and by initially incubating at 95°C for P. furiosus growth and then subsequently returning to 75°C to promote 1-propanol formation. Note that 1-propanol was not produced if the culture was grown only at 95°C. This work has the potential for future optimization through harnessing the genome-scale metabolic model of P. furiosus that was used herein to identify engineering targets to increase 1-propanol titer. IMPORTANCE As petroleum reserves become increasingly strained, the development of renewable alternatives to traditional chemical synthesis becomes more important. In this work, a high-temperature biological system for sugar to 1-propanol conversion was demonstrated by metabolic engineering of the hyperthermophilic archaeon Pyrococcus furiosus (T opt 100°C). The engineered strain produced 1-propanol by temperature shifting from 75°C to 95°C and then back to 75°C to accommodate the temperature ranges for native and foreign proteins associated with 1-propanol biosynthesis. Genome-scale metabolic modeling informed the carbon and reductant flux in the system, identified potential factors limiting 1-propanol production, and revealed potential optimization targets.

O'Quinn, Hailey C. (ORCID:0009000274865755)↗

Conjugation-based genome engineering enables rapid prototyping and bioproduction in non-model bacteria

Abstract Non-model bacteria offer unique metabolic capabilities for sustainable bioproduction, yet their limited genetic accessibility hinders systematic strain development. Here we present conjugation-based serine recombinase-assisted genome engineering (cSAGE), a broad-host-range platform that enables predictable, iterative genomic integration in transformation-resistant bacteria. cSAGE combines conjugative DNA delivery, standardized low-copy vectors, orthogonal recombinases, and modular genetic parts to support rapid pathway assembly and cross-host benchmarking. Using purple nonsulfur bacteria as a testbed, we integrate promoter engineering, multi-payload genome modification, and genome-scale metabolic modeling to empirically evaluate host-dependent pathway performance. Applying this workflow, we identify strain-specific differences in photosynthetic conversion of lignin-derived p -coumarate to the thermoplastic precursor p -vinylphenol. By enabling genome engineering and functional comparison across diverse bacteria using a single plasmid system, cSAGE provides a general framework for non-model strain prototyping and biotransformation discovery.

Guzman, Michael S. [Department of Chemical Enginee↗

Current State, Challenges, and Opportunities in Genome-Scale Resource Allocation Models: A Mathematical Perspective

Stoichiometric genome-scale metabolic models (generally abbreviated GSM, GSMM, or GEM) have had many applications in exploring phenotypes and guiding metabolic engineering interventions. Nevertheless, these models and predictions thereof can become limited as they do not directly account for protein cost, enzyme kinetics, and cell surface or volume proteome limitations. Lack of such mechanistic detail could lead to overly optimistic predictions and engineered strains. Initial efforts to correct these deficiencies were by the application of precursor tools for GSMs, such as flux balance analysis with molecular crowding. In the past decade, several frameworks have been introduced to incorporate proteome-related limitations using a genome-scale stoichiometric model as the reconstruction basis, which herein are called resource allocation models (RAMs). This review provides a broad overview of representative or commonly used existing RAM frameworks. This review discusses increasingly complex models, beginning with stoichiometric models to precursor to RAM frameworks to existing RAM frameworks. RAM frameworks are broadly divided into two categories: coarse-grained and fine-grained, with different strengths and challenges. Discussion includes pinpointing their utility, data needs, highlighting framework strengths and limitations, and appropriateness to various research endeavors, largely through contrasting their mathematical frameworks. Finally, promising future applications of RAMs are discussed.

59 BASIC BIOLOGICAL SCIENCES↗

Developing a media formulation to sustain ex vivo chloroplast function

Chloroplasts are critical organelles in plants and algae responsible for accumulating biomass through photosynthetic carbon fixation and cellular maintenance through metabolism in the cell. Chloroplasts are increasingly appreciated for their role in biomanufacturing, as they can produce many useful molecules, and a deeper understanding of chloroplast regulation and function would provide more insight for the biotechnological applications of these organelles. However, traditional genetic approaches to manipulate chloroplasts are slow, and generation of transgenic organisms to study their function can take weeks to months, significantly delaying the pace of research. To develop chloroplasts themselves as a quicker and more defined platform, we isolated chloroplasts from the green algae, Chlamydomonas reinhardtii, and examined their photosynthetic function after extraction. Combined with a metabolic modeling approach using flux-balance analysis, we identified key metabolic reactions essential to chloroplast function and leveraged this information into reagents that can be used in a “chloroplast media” capable of maintaining chloroplast photosynthetic function over time ex vivo compared to buffer alone. We envision this could serve as a model platform to enable more rapid design-build-test-learn cycles to study and improve chloroplast function in combination with genetic modifications and potentially as a starting point for the bottom-up design of a synthetic organelle-containing cell.

Chlamydomonas reinhardtii↗

Leptothrix ochracea genomes reveal potential for mixotrophic growth on Fe(II) and organic carbon

ABSTRACT Leptothrix ochracea creates distinctive iron-mineralized mats that carpet streams and wetlands. Easily recognized by its iron-mineralized sheaths, L. ochracea was one of the first microorganisms described in the 1800s. Yet it has never been isolated and does not have a complete genome sequence available, so key questions about its physiology remain unresolved. It is debated whether iron oxidation can be used for energy or growth and if L. ochracea is an autotroph, heterotroph, or mixotroph. To address these issues, we sampled L. ochracea -rich mats from three of its typical environments (a stream, wetlands, and a drainage channel) and reconstructed nine high-quality genomes of L. ochracea from metagenomes. These genomes contain iron oxidase genes cyc2 and mtoA, showing that L. ochracea has the potential to conserve energy from iron oxidation. Sox genes confer potential to oxidize sulfur for energy. There are genes for both carbon fixation (RuBisCO) and utilization of sugars and organic acids (acetate, lactate, and formate). In silico stoichiometric metabolic models further demonstrated the potential for growth using sugars and organic acids. Metatranscriptomes showed a high expression of genes for iron oxidation; aerobic respiration; and utilization of lactate, acetate, and sugars, as well as RuBisCO, supporting mixotrophic growth in the environment. In summary, our results suggest that L. ochracea has substantial metabolic flexibility. It is adapted to iron-rich, organic carbon-containing wetland niches, where it can thrive as a mixotrophic iron oxidizer by utilizing both iron oxidation and organics for energy generation and both inorganic and organic carbon for cell and sheath production. IMPORTANCE Winogradsky's observations of L. ochracea led him to propose autotrophic iron oxidation as a new microbial metabolism, following his work on autotrophic sulfur-oxidizers. While much culture-based research has ensued, isolation proved elusive, so most work on L. ochracea has been based in the environment and in microcosms. Meanwhile, the autotrophic Gallionella became the model for freshwater microbial iron oxidation, while heterotrophic and mixotrophic iron oxidation is not well-studied. Ecological studies have shown that Leptothrix overtakes Gallionella when dissolved organic carbon content increases, demonstrating distinct niches. This study presents the first near-complete genomes of L. ochracea , which share some features with autotrophic iron oxidizers, while also incorporating heterotrophic metabolisms. These genome, metabolic modeling, and transcriptome results give us a detailed metabolic picture of how the organism may combine lithoautotrophy with organoheterotrophy to promote Fe oxidation and C cycling and drive many biogeochemical processes resulting from microbial growth and iron oxyhydroxide formation in wetlands.

59 BASIC BIOLOGICAL SCIENCES↗

Charting the state of GEMs in microalgae: progress, challenges, and innovations

Genome-scale metabolic models (GEMs) provide a systems-level framework for understanding and engineering microalgal metabolism. This review explores the evolution of GEMs in microalgae, highlighting advances in light modeling, automation, and multi-omics integration. Special emphasis is placed on Chlamydomonas reinhardtii as a model species. Limitations of current models, particularly for microalgae, are discussed, alongside promising developments in dynamic modeling and machine learning. Together, these innovations chart a path toward more predictive, adaptable GEMs that can accelerate biotechnological applications of microalgae in sustainable production systems.

Plant Sciences↗

Development and Validation of a Process Model and Open-Source Process Simulator for Microalgae-Based Tertiary Phosphorus Recovery

Microalgae-based tertiary wastewater treatment has the potential to meet stringent effluent phosphorus limits, with the added benefit of producing a marketable feedstock. However, the lack of validated mechanistic models and their implementation in process simulators have limited the adoption of this technology. In this study, an updated lumped pathway metabolic model (Phototrophic-Mixotrophic Process Model, PM 2 ), including both photoautotrophic and heterotrophic metabolisms of microalgae, was developed to predict effluent phosphorus concentration and biomass yield in response to dynamic influent and varying environmental conditions. The model was implemented in QSDsan – an open-source, Python-based design and simulation platform – for robust simulation under uncertainty. A global sensitivity analysis was performed to prioritize model parameters for calibration. The model was then calibrated and validated using batch experimental data and 45 days of continuous online monitoring data from a full-scale (568 m 3 ·d -1 ) microalgae-based tertiary wastewater treatment plant (EcoRecover process). In particular, along with dynamic influent composition, temperature and light intensity data with diel variation were provided as model inputs to reflect the microalgal behavior under day-night cycling. Overall, the QSDsan-based microalgae process simulator was able to predict effluent phosphorus within 0.02–0.04 mg-P·L -1 , while also capturing the general trends of state variables according to nutrient availability.

Lumped pathway metabolic model↗

Spatially structured bacterial interactions alter algal carbon flow to bacteria

Phytoplankton account for nearly half of global photosynthetic carbon fixation, and the fate of that carbon is regulated in large part by microbial food web processing. We currently lack a mechanistic understanding of how interactions among heterotrophic bacteria impact the fate of photosynthetically fixed carbon. Here, we used a set of bacterial isolates capable of growing on exudates from the diatom Phaeodactylum tricornutum to investigate how bacteria-bacteria interactions affect the balance between exudate remineralization and incorporation into biomass. With exometabolomics and genome-scale metabolic modeling, we estimated the degree of resource competition between bacterial pairs. In a sequential spent media experiment, we found that pairwise interactions were more beneficial than predicted based on resource competition alone, and 30% exhibited facilitative interactions. To link this to carbon fate, we used single-cell isotope tracing in a custom cultivation system to compare the impact of different "primary" bacterial strains in close proximity to live P. tricornutum on a distal "secondary" strain. We found that a primary strain with a high degree of competition decreased secondary strain carbon drawdown by 51% at the single-cell level, providing a quantitative metric for the "cost" of competition on algal carbon fate. Additionally, a primary strain classified as facilitative based on sequential interactions increased total algal-derived carbon assimilation by 7.6 times, integrated over all members, compared to the competitive primary strain. Our findings suggest that the degree of interaction between bacteria along a spectrum from competitive to facilitative is directly linked to algal carbon drawdown.

genome-scale metabolic model↗

Prediction of non-intuitive metabolic targets with bayesian metabolic control analysis to improve 3-hydroxypropionic acid production in Aspergillus niger

Development of efficient bioconversion processes is limited by the ability to predictably improve metabolic flux. Here we deployed Bayesian Metabolic Control Analysis as a platform to integrate multi-omics data with metabolic modeling and evaluated its ability to predict genetic interventions that improve metabolic flux. Global Metabolomics and proteomics data was collected from 17 Aspergillus niger strains engineered to produce the platform biochemical 3-hydroxypropionic acid from which seven actional genetic interventions were predicted from significant flux control coefficients. Of the suggested genetic interventions, two were present within the intuitively designed strains used for training (malonic semialdehyde dehydrogenase and pyruvate carboxylase) while five predicted targets were present within non-intuitive areas of the metabolic network including 5-formyltetrahydrofolate deformylase and four mitochondrial enzymes, alcohol dehydrogenase, succinyl-CoA ligase, aspartate aminotransferase, and malate dehydrogenase. Six of the targets were validated in the highest performing 3-HP strain used for multi-omics data generation which contained a prior disruption of the highest scoring target malonic semialdehyde dehydrogenase. Predicted directional perturbation of five of the six tested targets significantly improved titer and rate of 3-HP production and two significantly improved yield. The greatest improvements were observed following disruption of the non-intuitive target succinyl-CoA ligase which increased titer by 39% and yield by 29% (to 20.4 g/L 3-HP and 0.31 g 3-HP/g glucose) over the strains used for training. This study demonstrates the utility of Bayesian Metabolic Control Analysis and highlights the ability to predict meaningful genetic targets in unexpected areas of metabolism to improve engineered strains for bioconversion.

3-hydroxypropionic acid↗

A multi-omic characterization of the physiological responses to salt stress in Scenedesmus obliquus UTEX393

Scenedesmus obliquus UTEX393 is a promising microalgal candidate for sustainable biomanufacturing but its limited halotolerance hinders large-scale cultivation in saline environments. To investigate the molecular basis of salt stress responses, we conducted a comprehensive multi-omic analysis integrating genomics, transcriptomics, proteomics, lipidomics, metabolomics, and DNA affinity purification sequencing (DAP-seq). An improved nuclear genome assembly and annotation yielded 19,017 gene models and a 97% BUSCO completeness score, enabling construction of a genome-scale metabolic model. Comparing 15 ppt salinity stress to 5 ppt control, growth and productivity were significantly reduced, accompanied by widespread transcriptomic and proteomic changes. Transcriptomic analysis revealed downregulation of photosynthetic machinery and energy conservation genes, and upregulation of stress-responsive elements such as expansins, flavodoxins, and osmoprotectants. Lipidomic profiling showed accumulation of triacylglycerols (TAGs) and degradation of galactosyl lipids, consistent with a shift toward lipid biosynthesis to mitigate redox imbalance. Depletion of key polar metabolites and branched-chain amino acids suggested a rerouting of central carbon metabolism under stress. DAP-seq identified key transcription factors, including LHY1 and SPL12, that target central metabolic enzymes involved in redox balancing, such as glyceraldehyde-3-phosphate dehydrogenase (GAPDH) and malate dehydrogenase (MDH). These findings establish a regulatory-metabolic framework linking redox stress to lipid accumulation and reveal potential engineering targets to enhance salt tolerance. Overall, the multi-omic analysis supports the “overflow” hypothesis, where impaired photosynthesis results in excess reducing equivalents being diverted into TAG synthesis and highlights transcriptional regulators as candidates for improving algal robustness in brackish environments.

09 BIOMASS FUELS↗

Gene-Metabolite Association Prediction with Interactive Knowledge Transfer Enhanced Graph for Metabolite Production

Identifying gene targets for enhancing metabolite production in metabolic engineering is challenging due to the vast research literature and the approximation in genome-scale metabolic model (GEM) simulations. Here, to address this, we propose the Gene-Metabolite Association Prediction task, which automates gene discovery for given metabolite-gene pairs, accompanied by a benchmark dataset of 2474 metabolites and 1947 genes for Saccharomyces cerevisiae (SC) and Issatchenkia orientalis (IO). This task is complicated by incomplete metabolic graphs and metabolic heterogeneity. We introduce an Interactive Knowledge Transfer mechanism based on Metabolism Graphs (IKT4Meta) to enhance prediction accuracy by integrating cross-metabolism knowledge. Using Pretrained Language Models (PLMs) to generate inter-graph links mitigates heterogeneity issues, while intra-graph links are propagated via these anchors. Gene-metabolite predictions are then performed on the enriched graphs integrating multiple microorganisms’ knowledge. Experiments show that IKT4Meta outperforms baselines by up to 12.3% in link prediction.

59 BASIC BIOLOGICAL SCIENCES↗

Coarse-Grained Modeling Elucidates Differential Metabolism of Saccharomyces cerevisiae under Varied Nutrient Limitations

Microorganisms such as Saccharomyces cerevisiae have a native ability to adapt their metabolism to varying nutrient conditions. Understanding their responses to nutrient limitations is critical for decoding cellular physiology and designing strategies for metabolic engineering. While the influence of carbon availability on yeast metabolism has been extensively studied, the role of nitrogen availability remains relatively underexplored. In this study, we utilized a coarse-grained kinetic model to systematically analyze and compare the effects of carbon and nitrogen limitations on yeast metabolism. Our model successfully revealed the differential metabolic characteristics of S. cerevisiae under carbon- and nitrogen-limited chemostat conditions. It also highlighted the significance of protein activity regulation at varying carbon-to-nitrogen ratios, and elucidated distinct strategies employed to maintain ATP homeostasis. Furthermore, this study provides a computational tool for investigating yeast physiology under nutrient limitations and offers quantitative and mechanistic insights into yeast metabolism.

59 BASIC BIOLOGICAL SCIENCES↗

Moisture Metaphenome Incubation Analysis Results

The Birch effect, a pulse of CO2 release that occurs when dry soil is rewet, is commonly observed, yet the underlying biogeochemistry remains elusive. Using multi-omics data, real-time mass spectrometry and modeling approaches, we investigated the molecular response to rewetting of a soil microbiome exposed to drought for one and two weeks. The microbiome response was evaluated through analysis of transcript, protein, metabolite, and respiration profiles and metabolic modeling using an enhanced version of the Metabolite Expression Metabolic Network Integration for Pathway Identification and Selection (MEMPIS) algorithm (Roy Chowdhury et al, mSystems, 2019).

Lipton, Mary S [Pacific Northwest National Laborat↗

pnnl-predictive-phenomics/SElon-GEM

The software was developed to solve the problem of continuous validation for genome-scale model of S. elongatus. The origin of the work was due to the difficulty of curating genome-scale metabolic models, and memote was invented to solve this problem.

Pino, James [Pacific Northwest National Laboratory↗

A Computational Workflow of Elucidating Viral Impact on Mediating Microbial Response to In-situ Experimental Warming: Bridging microbial modeling to carbon and mineral modeling

Viruses are abundant in soils and shape microbial communities in ways that can potentially influence ecosystem processes, yet their contributions to carbon cycling and mineral transformations remain poorly understood. Here we present a multi-phase framework that links virus-host interactions to soil biogeochemistry by combining ecological simulations, genome- and community-scale metabolic modeling, and statistical and machine-learning analyses. We first calibrated microbial abundance profiles under explicit infection scenarios to capture how viral pressure alters community structure, then explored alternative interaction strategies, including kill-the-winner, piggyback-the-winner, and mixed lytic-lysogenic modes, through forward simulations. These ecological shifts were translated into metabolic consequences using exchange fluxes summarized into biologically meaningful categories, while integrated statistical and machine-learning screens elevated subtle but consistent signals. Application of this framework revealed that viral infections shift the balance between organic and inorganic fluxes, redirecting metabolism from diffuse organic transformations toward inorganic pools such as protons and CO 2 , directly linking viral regulation to respiration and soil carbon balance. The roll-up analysis also isolated perturbations in critical mineral ions, including magnesium, manganese, zinc, and copper, which serve as essential enzymatic cofactors. In piggyback-the-winner scenarios, uptake of these ions was strongly suppressed. Contrasting viral strategies produced distinct community structures and metabolic outcomes, from broad suppression under kill-the-winner dynamics to dramatic redistributions under high-lytic and high-gain lysogenic regimes that collapsed vulnerable microbial populations while promoting opportunists. Together, these results provide a tractable path to trace viral perturbations from host abundance shifts to metabolic flux adjustments and ecosystem-scale processes, offering a practical way to include viruses in earth system models.

54 ENVIRONMENTAL SCIENCES↗

Addressing genome scale design tradeoffs in Pseudomonas putida for bioconversion of an aromatic carbon source

Genome-scale metabolic models (GSMM) are commonly used to identify gene deletion sets that result in growth coupling and pairing product formation with substrate utilization and can improve strain performance beyond levels typically accessible using traditional strain engineering approaches. However, sustainable feedstocks pose a challenge due to incomplete high-resolution metabolic data for non-canonical carbon sources required to curate GSMM and identify implementable designs. Here we address a four-gene deletion design in the Pseudomonas putida KT2440 strain for the lignin-derived non-sugar carbon source, p-coumarate (p-CA), that proved challenging to implement. We examine the performance of the fully implemented design for p-coumarate to glutamine, a useful biomanufacturing intermediate. In this study glutamine is then converted to indigoidine, an alternative sustainable pigment and a model heterologous product that is commonly used to colorimetrically quantify glutamine concentration. Through proteomics, promoter-variation, and growth characterization of a fully implemented gene deletion design, we provide evidence that aromatic catabolism in the completed design is rate-limited by fumarase hydratase (FUM) enzyme activity in the citrate cycle and requires careful optimization of another fumarate hydratase protein (PP_0897) expression to achieve growth and production. A double sensitivity analysis also confirmed a strict requirement for fumarate hydratase activity in the strain where all genes in the growth coupling design have been implemented. Metabolic cross-feeding experiments were used to examine the impact of complete removal of the fumarase hydratase reaction and revealed an unanticipated nutrient requirement, suggesting additional functions for this enzyme. While a complete implementation of the design was achieved, this study highlights the challenge of completely inactivating metabolic reactions encoded by under-characterized proteins, especially in the context of multi-gene edits.

59 BASIC BIOLOGICAL SCIENCES↗

Multi-strain analysis of Pseudomonas putida reveals the metabolic and genetic diversity of the species

Pseudomonas putida is a gram-negative bacterial species increasingly utilized in biotechnology due to its robust growth, ability to degrade aromatic compounds, solvent tolerance, and genetic tractability. In this study, we report a comprehensive multi-strain analysis of 164 P. putida strains based on the reconstruction of a pan-putida metabolic network and the formulation of strain-specific genome-scale metabolic models (GEMs). We performed whole-genome sequencing and hybrid assembly for 40 strains, contributing a ~8% increase to the available genomic data for P. putida . Furthermore, high-throughput phenotypic profiling using the Biolog phenotype microarray system for 24 strains on 190 unique carbon sources, along with 15 aromatic compounds not present on Biolog plates, yielded 4,920 unique strain-phenotype measurements. These data were leveraged to curate GEMs for 24 representative strains, including a refined model for strain KT2440, which comprised 1,480 genes and 2,191 metabolites, achieving a prediction accuracy of 91.2% in carbon utilization. Systematic comparison of genomes and GEMs revealed both conserved core pathways and significant allelic and functional divergence across strains, highlighting strain-specific variation in aromatic degradation. While pathways for protocatechuate and phenylacetate degradation were widely conserved, metabolic capabilities for compounds such as ferulate, phenol, and cresols varied markedly, suggesting adaptation to distinct ecological niches. Alleleome analysis of enzymes, such as PcaI and PcaJ, revealed distinct, functionally similar clades, indicating possible convergent evolution or horizontal gene transfer. These results provide computable resources and informative models for selecting P. putida strains with desired traits for biomanufacturing and bioremediation and offer insights into the evolution and phylogeny of the P. putida species.

aromatics utilization↗