Search for: All records

Bacteriophage auxiliary metabolic genes (AMGs) alter host metabolism by hijacking reactions, but previous studies mostly inferred their roles from annotations, ignoring system-wide impacts and phage production. Here we integrate AMGs and phage assembly into a genome-scale metabolic model ofProchloroccocus marinusMED4 infected by P-HM2. We show that 17 directly hijacked reactions substantially affect more than 30% of the reactions in MED4 metabolism, including carbon fixation, photosynthesis, and nucleotide synthesis, distinguishing these AMGs as either phage aligned—shifting feasible reaction velocities in accordance with maximal phage production—or phage antialigned. Pareto optimization reveals that phage-aligned reactions alter phage-host growth trade-offs, while phage-antialigned reactions do not. We experimentally validate our predictions of system-level AMG impacts by measuring the N-dependent effect of P-HM2cp12expression on growth in a model relative of the genetically intractable MED4,Synechococcus elongatus. We also show that AMGs’ indirect impacts are synergistically and antagonistically coupled, providing systems-level insight into AMG perturbations and highlighting how nontrivial cascading effects shape host metabolism. 

Stochastic nature of gene expression leads to the complex formation of the bacterial transcriptome and proteome. In contrast to typical transcriptome studies, we employ a near wild-type, Syn1.0, of the naturally genome-reducedMycoplasmas, and the dramatically further genome-reduced JCVI-syn3A thus avoiding additional contributions from many non-essential cellular functions. To aid in profiling the transcriptional landscape within these bacteria, we present a bioinformatic analysis of the genetic sequence motifs implicated in modulating the stochastic gene expression events, coupled with genome-wide short-read (Illumina) and long-read (Oxford Nanopore Technologies and Pacfic Biosciences) RNA sequencing. The bioinformatic analysis coupled with information from structural studies assigns strengths of the Shine-Dalgarno signatures and identifies both transcription initiation and termination sites, leading to predictions of RNA isoforms in Syn1.0 (and related organisms). The long-read and short-read RNA sequencing characterized the predicted transcriptional activity, and the long-read methods provide direct insight into the RNA isoform complexity within Syn1.0. Comparison of the RNA sequencing results with that of the bioinformatic analysis highlights the inability of bioinformatics alone to capture the results of bacterial transcription without including effects of RNA degradation. This study emphasizes the need for comparative analysis and potential dangers of genome reduction, exemplified through the discovery of altered gene expression patterns of JCVI-syn1.0 and JCVI-syn3A, achieved via the union of our transcriptome study with their proteomics data. Analysis of the transcriptomics data sets through a Jupyter notebook allows any genomic region to be easily examined. Table of Content Image 

Abstract Photosynthetic microorganisms rely on multiple pathways in central carbon metabolism to adapt to fluctuating light and energy availability across diel cycles. Mechanistic insight into the regulatory dynamics of this adaptation requires integrating processes spanning disparate timescales, from rapid redox-dependent post-translational modifications (PTMs) to slower changes in protein expression and metabolic pathway usage. To address this complexity beyond genome-based inference and traditional modeling, we develop a whole-cell four-dimensional (3D + time) model of the marine cyanobacteriumProchlorococcus marinusMED4 that explicitly represents the spatial organization of enzymatic and molecular processes in central carbon metabolism under light perturbation. We employ a perturbation-based research design to experimentally generate time-series, multi-omics measurements that provide molecular descriptors and cryo-ET images as constraints for this dynamic 4D framework. The integration of experiments and modeling across defined light regimes enables quantitative validation of system-level responses and forecasting under distinct light disturbances. We test the hypothesis that light-dependent redox PTMs regulating the structural assembly of a protein megacomplex, the “dark complex,” modulate metabolic flux at a conserved regulatory node of the Calvin–Benson cycle (CBC) in cyanobacteria. Our model shows that subcellular spatial organization buffers rapid light-induced changes in thylakoid reaction rates, which are followed by redox-PTM-mediated sequestration or release of CBC enzymes in the dark complex, ultimately impacting carbon fixation dynamics within carboxysomes. Comparison with an equivalently parameterized well-mixed stochastic model demonstrates that post-translational regulation not only buffers transcriptional noise and diffusion-driven fluctuations but also stabilizes phenotypic outcomes, underscoring the importance of spatial heterogeneity in phenotypic robustness. This ability to probe adaptive, spatiotemporally resolved mechanisms in photosynthetic machinery and central carbon metabolism addresses a critical gap in genotype-to-phenotype inference and expands modeling and design capabilities for understudied or genetically intractable autotrophs such asP. marinusMED4. Significance StatementThis work advances 4D whole-cell modeling by presenting the first spatiotemporal simulation of a photosynthetic autotroph using the Lattice Microbes platform. UsingProchlorococcus marinusMED4, we show that subcellular spatial organization of organelles, diffusion constraints, and redox regulation collectively shape central carbon metabolism across orders of magnitude in space and time. Through a perturbation-based strategy that generates multi-omics data sets over time, we construct and validate a spatially and temporally resolved model of MED4, constrained by high-resolution (10 nm) cryo-electron tomography. Our results highlight the importance of localized biochemical reactions and redox-dependent post-translational modification of enzymes in regulating carbon fixation in a noisy environment under light disturbance. This study establishes a spatiotemporal, whole-cell physiology modeling framework as a transformative tool for uncovering multiscale regulatory responses to environmental gradients.