Search NASA⌕ Search

SEARCH · Search NASA

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 541 records · Page 30

Modeling temporal and hormonal regulation of plant transcriptional response to wounding

Plants respond to wounding stress by changing gene expression patterns and inducing the production of hormones including jasmonic acid. This wounding transcriptional response activates specialized metabolism pathways such as the glucosinolate pathways in Arabidopsis thaliana. While the regulatory factors and sequences controlling a subset of wound-response genes are known, it remains unclear how wound response is regulated globally. Here, we how these responses are regulated by incorporating putative cis-regulatory elements, known transcription factor binding sites, in vitro DNA affinity purification sequencing, and DNase I hypersensitive sites to predict genes with different wound-response patterns using machine learning. We observed that regulatory sites and regions of open chromatin differed between genes upregulated at early and late wounding time-points as well as between genes induced by jasmonic acid and those not induced. Expanding on what we currently know, we identified cis-elements that improved model predictions of expression clusters over known binding sites. Using a combination of genome editing, in vitro DNA-binding assays, and transient expression assays using native and mutated cis-regulatory elements, we experimentally validated four of the predicted elements, three of which were not previously known to function in wound-response regulation. Our study provides a global model predictive of wound response and identifies new regulatory sequences important for wounding without requiring prior knowledge of the transcriptional regulators.

59 BASIC BIOLOGICAL SCIENCES↗

Identification of potent inhibitors of JUN N-terminal kinases for treatment of endometriosis and associated pain

Endometriosis, defined as the ectopic growth of endometrial tissue outside of the uterine cavity, is an inflammatory and hormone-dependent disease that causes excruciating pelvic pain, infertility, and significantly decreases quality of life in affected patients. The JUN N-terminal kinases (JNKs) are a leading class of nonhormonal therapeutic targets that have been validated in preclinical models of endometriosis and in a Phase 1/2 clinical trial. Despite their therapeutic potential, JNK inhibitors with increased potency and specificity are needed to address the inflammatory pathology of endometriosis and to prevent disease progression. Leveraging a DNA-encoded chemical library collection of ~4 billion compounds, we identified lead inhibitor CDD-2428 and optimized derivatives, CDD-2728 and CDD-3013, with excellent binding affinity to JNK1-3 (K d = 0.12 to 3.7 nM), enhanced selectivity, metabolic stability, and cellular permeability. Crystallographic and biochemical studies confirmed that CDD-3013 exhibited superior kinase selectivity with improved efficacy compared to existing JNK inhibitors. In primary endometriosis cell models, CDD-2728 and CDD-3013 suppressed JNK-dependent inflammatory signaling, dampening pathways linked to pain, invasion, angiogenesis, and macrophage recruitment. In an endometriosis mouse model, both CDD-2728 and CDD-3013 reduced endometriotic lesion size, macrophage infiltration, and cellular proliferation, showing in vivo efficacy. When tested in a lipopolysaccharide-induced hyperalgesia model, CDD-2728 and CDD-3013 decreased markers of induced pain, as measured by changes in a dynamic weight bearing test and Grimace scores. These findings nominate CDD-2728 and CDD-3013 as potent, nonhormonal therapeutic candidates for endometriosis with broad anti-inflammatory and analgesic activity, addressing a critical unmet clinical need.

Madasu, Chandrashekhar [Department of Pathology an↗

Carbon source-induced changes in the transcriptional landscape of Novacetimonas hansenii

Novacetimonas hansenii is a model organism for bacterial cellulose (BC) production, a biopolymer with broad industrial applications. To better understand its carbon source-dependent metabolism, we conducted the first time-resolved transcriptomic analysis of N. hansenii ATCC 23769 grown on five carbon sources (glucose, fructose, mannitol, galactose, and sucrose) across five time points (24–120h). Glucose, fructose, and mannitol supported growth and efficient BC production, whereas galactose and sucrose resulted in poor yields, likely due to limited transport and enzymatic capacity. Solid-state NMR revealed that cellulose crystallinity and polymorphism were carbon source-dependent, with reduced bcsD expression correlating with lower crystallinity. Transcriptomic profiling showed distinct metabolic responses. Cultures containing glucose and mannitol activated core metabolic pathways, fructose induced a delayed adaptive shift, and galactose and sucrose triggered stress-related gene expression. These results highlight the metabolic plasticity of N. hansenii and offer insights for enhancing BC production through carbon source optomization.

Bacterial cellulose↗

Formation of a constructed microbial community in a nutrient-rich environment indicates bacterial interspecific competition

ABSTRACT Understanding the organizational principles of microbial communities is essential for interpreting ecosystem stability. Previous studies have investigated the formation of bacterial communities under nutrient-poor conditions or obligate relationships to observe cooperative interactions among different species. How microorganisms form stabilized communities in nutrient-rich environments, without obligate metabolic interdependency for growth, is still not fully disclosed. In this study, three bacterial strains isolated from the Populus deltoides rhizosphere were co-cultured in complex medium, and their growth behavior was tracked. These strains co-exist in mixed culture over serial transfer for multiple growth-dilution cycles. Competition is proposed as an emergent interaction relationship among the three bacteria based on their significantly decreased growth levels. The effects of different initial inoculum ratios, up to three orders of magnitude, on community structure were investigated, and the final compositions of the mixed communities with various starting composition indicate that community structure is not dependent on the initial inoculum ratio. Furthermore, the competitive relationships within the community were not altered by different initial inoculum ratios. The community structure was simulated by generalized Lotka-Volterra and dynamic flux balance analysis to provide mechanistic predictions into emergence of community structure under a nutrient-rich environment. Metaproteomic analyses provide support for the metabolite exchanges predicted by computational modeling and for highly altered physiologies when microbes are grown in co-culture. These findings broaden our understanding of bacterial community dynamics and metabolic diversity in higher-order interactions and could be significant in the management of rhizospheric bacterial communities. IMPORTANCE Bacteria naturally co-exist in multispecies consortia, and the ability to engineer such systems can be useful in biotechnology. Despite this, few studies have been performed to understand how bacteria form a stable community and interact with each other under nutrient-rich conditions. In this study, we investigated the effects of initial inoculum ratios on bacterial community structure using a complex medium and found that the initial inoculum ratio has no significant impact on resultant community structure or on interaction patterns between community members. The microbial population profiles were simulated using computational tools in order to understand intermicrobial relationships and to identify potential metabolic exchanges that occur during stabilization of the bacterial community. Studying microbial community assembly processes is essential for understanding fundamental ecological principles in microbial ecosystems and can be critical in predicting microbial community structure and function.

59 BASIC BIOLOGICAL SCIENCES↗

Hyperaccumulation of Gadolinium by Methylorubrum extorquens AM1 Reveals Impacts of Lanthanides on Cellular Processes Beyond Methylotrophy

Lanthanides (Ln) are a new group of life metals, and many questions remain regarding how they are acquired and used in biology. Methylotrophic bacteria can acquire, transport, biomineralize, and use Ln as part of a cofactor complex with pyrroloquinoline quinone (PQQ) in alcohol dehydrogenases. For most methylotrophic bacteria use is restricted to the light Ln, which range from lanthanum to samarium (atomic numbers 57–62). Understanding how the cell differentiates between light and heavy Ln, and the impacts of these metals on the metabolic network, will advance the field of Ln biochemistry and give insights into enzyme catalysis, stress homeostasis, and metal biomineralization and compartmentalization. We report robust methanol growth with the heavy Ln gadolinium by a genetic variant of the model methylotrophic bacterium Methylorubrum extorquens AM1, named evo -HLn, for “ evo lved for H eavy L antha n ides.” A non-synonymous single nucleotide polymorphism in a cytosolic hybrid histidine kinase/response regulator allowed for sweeping transcriptional alterations to heavy metal stress response, methanol oxidation, and central metabolism. Increased expression of genes for Ln acquisition and uptake, production of the Ln-chelating lanthanophore, PQQ biosynthesis, and phosphate transport and metabolism resulted in gadolinium hyperaccumulation of 36-fold with a trade-off for light Ln accumulation. Gadolinium was hyperaccumulated in an enlarged acidocalcisome-like compartment. This is the first evidence of a bacterial intracellular Ln-containing compartment that we name the “lanthasome.” Carotenoid and toblerol biosynthesis were also upregulated. Due to its unique capabilities, evo -HLn can be used to further magnetic resonance imaging (MRI) and bioremediation technologies. In this regard, we show that gadolinium hyperaccumulation was sufficient to produce MRI contrast in whole cells, and that evo -HLn was able to readily acquire the metal from the MRI contrast agent gadopentetic acid. Finally, hyperaccumulation of gadolinium, differential uptake of light and heavy Ln, increased PQQ levels, and phosphate transport provide new insights into strategies for Ln recovery.

59 BASIC BIOLOGICAL SCIENCES↗

In Vivo Thermodynamic Analysis of Glycolysis in Clostridium thermocellum and Thermoanaerobacterium saccharolyticum Using 13 C and 2 H Tracers

Clostridium thermocellum and Thermoanaerobacterium saccharolyticum are thermophilic anaerobic bacteria with complementary metabolic capabilities that utilize distinct glycolytic pathways for the conversion of cellulosic sugars to biofuels. We integrated quantitative metabolomics with 2 H and 13 C metabolic flux analysis to investigate the in vivo reversibility and thermodynamics of the central metabolic networks of these two microbes. We found that the glycolytic pathway in C. thermocellum operates remarkably close to thermodynamic equilibrium, with an overall drop in Gibbs free energy 5-fold lower than that of T. saccharolyticum or anaerobically grown Escherichia coli . The limited thermodynamic driving force of glycolysis in C. thermocellum could be attributed in large part to the small free energy of the phosphofructokinase reaction producing fructose bisphosphate. The ethanol fermentation pathway was also substantially more reversible in C. thermocellummthan in T. saccharolyticum. These observations help explain the comparatively low ethanol titers of C. thermocellum and suggest engineering interventions that can be used to increase its ethanol productivity and glycolytic rate. In addition to thermodynamic analysis, we used our isotope tracer data to reconstruct the T. saccharolyticum central metabolic network, revealing exclusive use of the Embden-Meyerhof-Parnas (EMP) pathway for glycolysis, a bifurcated tricarboxylic acid (TCA) cycle, and a sedoheptulose bisphosphate bypass active within the pentose phosphate pathway. Thermodynamics constitutes a key determinant of flux and enzyme efficiency in metabolic networks. Here, we provide new insights into the divergent thermodynamics of the glycolytic pathways of C. thermocellum and T. saccharolyticum, two industrially relevant thermophilic bacteria whose metabolism still is not well understood. We report that while the glycolytic pathway in T. saccharolyticum is as thermodynamically favorable as that found in model organisms, such as E. coli or Saccharomyces cerevisiae , the glycolytic pathway of C. thermocellum operates near equilibrium. The use of a near-equilibrium glycolytic pathway, with potentially increased ATP yield, by this cellulolytic microbe may represent an evolutionary adaptation to growth on cellulose, but it has the drawback of being highly susceptible to product feedback inhibition. The results of this study will facilitate future engineering of high-performance strains capable of transforming cellulosic biomass to biofuels at high yields and titers.

59 BASIC BIOLOGICAL SCIENCES↗

Radiation exposure induces cross-species temporal metabolic changes that are mitigated in mice by amifostine

Exposure to acute, damaging radiation may occur through a variety of events from cancer therapy and industrial accidents to terrorist attacks and military actions. Our understanding of how to protect individuals and mitigate the effects of radiation injury or Acute Radiation Syndrome (ARS) is still limited. There are only a few Food and Drug Administration-approved therapies for ARS; whereas, amifostine is limited to treating low dose (0.7–6 Gy) radiation poisoning arising from cancer radiotherapy. An early intervention is critical to treat ARS, which necessitates identifying diagnostic biomarkers to quickly characterize radiation exposure. Towards this end, a multiplatform metabolomics study was performed to comprehensively characterize the temporal changes in metabolite levels from mice and non-human primate serum samples following γ-irradiation. The metabolomic signature of amifostine was also evaluated in mice as a model for radioprotection. The NMR and mass spectrometry metabolomics analysis identified 23 dysregulated pathways resulting from the radiation exposure. These metabolomic alterations exhibited distinct trajectories within glucose metabolism, phospholipid biosynthesis, and nucleotide metabolism. A return to baseline levels with amifostine treatment occurred for these pathways within a week of radiation exposure. Together, our data suggests a unique physiological change that is independent of radiation dose or species. Furthermore, a metabolic signature of radioprotection was observed through the use of amifostine prophylaxis of ARS.

59 BASIC BIOLOGICAL SCIENCES↗

Rubisco regulation in response to altered carbon status in the cyanobacterium Synechococcus elongatus PCC 7942

Abstract Photosynthetic organisms possess a variety of mechanisms to achieve balance between absorbed light (source) and the capacity to metabolically utilize or dissipate this energy (sink). While regulatory processes that detect changes in metabolic status/balance are relatively well studied in plants, analogous pathways remain poorly characterized in photosynthetic microbes. Here, we explored systemic changes that result from alterations in carbon availability in the model cyanobacterium Synechococcus elongatus PCC 7942 by taking advantage of an engineered strain where influx/efflux of a central carbon metabolite, sucrose, can be regulated experimentally. We observed that induction of a high-flux sucrose export pathway leads to depletion of internal carbon storage pools (glycogen) and concurrent increases in estimates of photosynthetic activity. Further, a proteome-wide analysis and fluorescence reporter-based analysis revealed that upregulated factors following the activation of the metabolic sink are concentrated on ribulose-1,5-bisphosphate carboxylase-oxygenase (Rubisco) and auxiliary modules involved in Rubisco maturation. Carboxysome number and Rubisco activity also increased following engagement of sucrose secretion. Conversely, reversing the flux of sucrose by feeding exogenous sucrose through the heterologous transporter resulted in increased glycogen pools, decreased Rubisco abundance, and carboxysome reorganization. Our data suggest that Rubisco activity and organization are key variables connected to regulatory pathways involved in metabolic balancing in cyanobacteria.

Plant Sciences↗

Revisiting diffusion-based moisture functions: why do they fail?

Microbial activities in soil are controlled by the balance between physical and biological mechanisms that are currently implied by soil moisture-respiration relationships, but not explicitly represented. Diffusion based moisture functions focus on physical controls, while substrate kinetics based models primarily address biological drivers. How to harmonize these two types of models for robust predictions of soil respiration and other processes has not been addressed. Here we derived a generalizable analytical solution of diffusion-limited nonlinear microbial uptake kinetics, in which microbial metabolic rates were governed by both physical and biological drivers via substrate diffusion and microbial affinity. We proposed a soil trait concept to integrate both physical and biological perspectives into one unified function, and tested it using both laboratory and field based observations from a variety of soils. This approach provides a tractable and scalable approach that incorporates soil heterogeneity and functional complexity in macroscopic model descriptions.

59 BASIC BIOLOGICAL SCIENCES↗

Thioester synthesis by a designed nickel enzyme models prebiotic energy conversion

The formation of carbon–carbon bonds from prebiotic precursors such as carbon dioxide represents the foundation of all primordial life processes. In extant organisms, this reaction is carried out by the carbon monoxide dehydrogenase (CODH)/acetyl coenzyme A synthase (ACS) enzyme, which performs the cornerstone reaction in the ancient Wood–Ljungdahl metabolic pathway to synthesize the key biological metabolite, acetyl-CoA. Despite its significance, a fundamental understanding of this transformation is lacking, hampering efforts to harness analogous chemistry. To address these knowledge gaps, we have designed an artificial metalloenzyme within the azurin protein scaffold as a structural, functional, and mechanistic model of ACS. We demonstrate the intermediacy of the Ni I species and requirement for ordered substrate binding in the bioorganometallic carbon–carbon bond-forming reaction from the one-carbon ACS substrates. The electronic and geometric structures of the nickel-acetyl intermediate have been characterized using time-resolved optical, electron paramagnetic resonance, and X-ray absorption spectroscopy in conjunction with quantum chemical calculations. Moreover, we demonstrate that the nickel-acetyl species is chemically competent for selective acyl transfer upon thiol addition to biosynthesize an activated thioester. Drawing an analogy to the native enzyme, a mechanism for thioester generation by this ACS model has been proposed. The fundamental insight into the enzymatic process provided by this rudimentary ACS model has implications for the evolution of primitive ACS-like proteins. Ultimately, these findings offer strategies for development of highly active catalysts for sustainable generation of liquid fuels from one-carbon substrates, with potential for broad applications across diverse fields ranging from energy storage to environmental remediation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Daylength variation affects growth, photosynthesis, leaf metabolism, partitioning, and metabolic fluxes

Abstract Daylength, a seasonal and latitudinal variable, exerts a substantial impact on plant growth. However, the relationship between daylength and growth is nonproportional, suggesting the existence of adaptive mechanisms. Thus, our study aimed to comprehensively investigate the adaptive strategies employed by plants in response to daylength variation. We grew false flax (Camelina sativa) plants, a model oilseed crop, under long-day (LD) and short-day (SD) conditions and used growth measurements, gas exchange measurements, and isotopic labeling techniques, including 13C, 14C, and 2H2O, to determine responses to different daylengths. Our findings revealed that daylength influences various growth parameters, photosynthetic physiology, carbon partitioning, metabolic fluxes, and metabolite levels. SD plants employed diverse mechanisms to compensate for reduced CO2 fixation in the shorter photoperiod. These mechanisms included enhanced photosynthetic rates and reduced respiration in the light (RL), leading to increased shoot investment. Additionally, SD plants exhibited reduced rates of the glucose 6-phosphate (G6P) shunt and greater partitioning of sugars into starch, thereby sustaining carbon availability during the longer night. Isotopic labeling results further demonstrated substantial alterations in the partitioning of amino acids and TCA cycle intermediates between rapidly and slowly turning over pools. Overall, the results point to multiple developmental, physiological, and metabolic ways in which plants adapt to different daylengths to maintain growth.

59 BASIC BIOLOGICAL SCIENCES↗

Reshaping the 2‐Pyrone Synthase Active Site for Chemoselective Biosynthesis of Polyketides

Abstract Engineering enzymes with novel reactivity and applying them in metabolic pathways to produce valuable products are quite challenging due to the intrinsic complexity of metabolic networks and the need for high in vivo catalytic efficiency. Triacetic acid lactone (TAL), naturally generated by 2‐pyrone synthase (2PS), is a platform molecule that can be produced via microbial fermentation and further converted into value‐added products. However, these conversions require extra synthetic steps under harsh conditions. We herein report a biocatalytic system for direct generation of TAL derivatives under mild conditions with controlled chemoselectivity by rationally engineering the 2PS active site and then rewiring the biocatalytic pathway in the metabolic network of E. coli to produce high‐value products, such as kavalactone precursors, with yields up to 17 mg/L culture. Computer modeling indicates sterics and hydrogen‐bond interactions play key roles in tuning the selectivity, efficiency and yield.

Zhou, Yu↗

Reshaping the 2‐Pyrone Synthase Active Site for Chemoselective Biosynthesis of Polyketides

Abstract Engineering enzymes with novel reactivity and applying them in metabolic pathways to produce valuable products are quite challenging due to the intrinsic complexity of metabolic networks and the need for high in vivo catalytic efficiency. Triacetic acid lactone (TAL), naturally generated by 2‐pyrone synthase (2PS), is a platform molecule that can be produced via microbial fermentation and further converted into value‐added products. However, these conversions require extra synthetic steps under harsh conditions. We herein report a biocatalytic system for direct generation of TAL derivatives under mild conditions with controlled chemoselectivity by rationally engineering the 2PS active site and then rewiring the biocatalytic pathway in the metabolic network of E. coli to produce high‐value products, such as kavalactone precursors, with yields up to 17 mg/L culture. Computer modeling indicates sterics and hydrogen‐bond interactions play key roles in tuning the selectivity, efficiency and yield.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Establishment of a genome editing tool using CRISPR-Cas9 ribonucleoprotein complexes in the non-model plant pathogen Sphaerulina musiva

CRISPR-Cas9 is a versatile genome editing system widely used since 2013 to introduce site-specific modifications into the genomes of model and non-model species. This technology is used in various applications, from gene knock-outs, knock-ins, and over-expressions to more precise changes, such as the introduction of nucleotides at a targeted locus. CRISPR-Cas9 has been demonstrated to be easy to establish in new species and highly efficient and specific compared to previous gene editing strategies such as Zinc finger nucleases and transcription activator-like effector nucleases. Grand challenges for emerging CRISPR-Cas9 tools in filamentous fungi are developing efficient transformation methods for non-model organisms. In this paper, we have leveraged the establishment of CRISPR-Cas9 genome editing tool that relies on Cas9/sgRNA ribonucleoprotein complexes (RNPs) in the model species Trichoderma reesei and developed the first protocol to efficiently transform the non-model species, Sphaerulina musiva. This fungal pathogen constitutes a real threat to the genus Populus, a foundational bioenergy crop used for biofuel production. Herein, we highlight the general considerations to design sgRNAs and their computational validation. We also describe the use of isolated protoplasts to deliver the CRISPR-Cas9 RNP components in both species and the screening for targeted genome editing events. The development of engineering tools in S. musiva can be used for studying genes involved in diverse processes such as secondary metabolism, establishment, and pathogenicity, among many others, but also for developing genetic mitigation approaches. The approach described here provides guidance for potential development of transformation systems in other non-model spore-bearing ascomycetes.

59 BASIC BIOLOGICAL SCIENCES↗

Insight into industrial hemp ( Cannabis sativa L.) root exudation composition in a simulated soil environment: a rhizosphere-on-a-chip study

Microfluidic technologies provide a reduced complexity and soil-free environment to study plant-soil interactions at the microscale. Traditionally used for model plants such as Arabidopsis thaliana, this study represents the first application of a rhizosphere-on-a-chip (RhizoChip) to investigate root exudation in industrial hemp (Cannabis sativa L.), an agronomic crop with growing economic importance. By incorporating soil-like minerals (kaolinite, potassium feldspar, and biotite), the RhizoChip addresses limitations of previous research. Hemp seedlings grown in mineral-containing chips exhibited significant root growth, emphasizing the critical role of minerals in root development. Using untargeted metabolomics, 170 compounds were identified, including organic acids, amino acids, and secondary metabolites, with distinct profiles across conditions. Metabolic pathway analysis revealed activity in amino acid metabolism, the citric acid cycle, and secondary metabolite biosynthesis. In conclusion, this study highlights the RhizoChip's potential for long-term studies of root exudates in non-model crops and offers insights into rhizosphere processes with implications for sustainable agriculture.

59 BASIC BIOLOGICAL SCIENCES↗

Data and scripts associated with a manuscript modeling microbial regulation of priming effects

This data package is associated with the publication “Modeling Microbial Regulatory Feedback in Organic Matter Decomposition Identifies Copiotrophic Traits as Key Drivers of Positive Priming” published as a preprint on BioRXiv by Ahamed et al. (2026); https://doi.org/10.1101/2024.08.11.607483. The package contains MATLAB scripts and saved simulation outputs used to implement a cybernetic model of microbial regulation during complex organic matter (OM) decomposition governing priming effects. It includes models of (i) single microbial functional groups (copiotrophic or oligotrophic degraders) and (ii) binary consortia composed of degraders and non-degraders with contrasting or common growth traits. Simulation results were generated using Monte Carlo analyses, with randomized key model parameters across a range of environmental mixing fractions of complex and labile OM. The dataset was created to provide a transparent and reusable computational framework for systematically exploring how microbial growth traits, metabolic regulation, and community composition influence OM decomposition dynamics and priming effects. For details on how to navigate data packages generated by this project, see https://data.ess-dive.lbl.gov/portals/PNNLRiverCorridorSFA/About. In addition to a readme, this data package also includes a file-level metadata (FLMD) file that describes each file and a data dictionary (DD) that describes the variable definitions. This package includes: (1) annotated MATLAB code implementing the system of ordinary differential equations and cybernetic control laws; (2) saved output files containing data (e.g., biomass, substrates, enzyme levels, priming metrics); and (3) scripts for processing saved outputs and regenerating figures. Specifically, the data package contains three main MATLAB scripts: runPrimingModel.m, runPlotData.m, and runPlotSuppFigS1.m, along with this readme and supporting documentation. Users should begin with runPrimingModel.m, which contains the annotated code implementing the system of ordinary differential equations and cybernetic control laws. This script runs the Monte Carlo simulations of microbial OM decomposition and allows users to modify microbial trait definitions, adjust parameter distributions, or define new community configurations. Simulation outputs are automatically saved as .mat files in the folder named SavedData, which stores all pre-generated results included in this package. The second script, runPlotData.m, reads files from the SavedData folder and processes them to regenerate the figures presented in the manuscript. The third script, runPlotSuppFigS1.m, specifically generates Figure S1 in the Supplementary Material of the manuscript. The package also includes the aforementioned files in non-proprietary .txt format. If users intend to use them, they should first save the files in their respective .m or .mat formats prior to execution in MATLAB.

Biomass concentration↗

Spatial Study 2022: Water Column, Sediment, and Total Ecosystem Respiration Rates across the Yakima River Basin, Washington, USA (v2)

This dataset supports a broader study examining the drivers of spatial variability in sediment respiration rates in the Yakima River Basin and is associated with the manuscript “Sediment-associated processes account for most of the spatial variation in ecosystem respiration in the Yakima River basin” submitted to Nature Communications Earth & Environment (Garayburu-Caruso et al., in review). The dataset provides ecosystem metabolism estimates generated from streamMetabolizer (Appling et al.; 2018) using data collected during the same five-week period at 48 sites within multiple rivers throughout the Yakima River Basin in Washington, USA. Additionally, it includes the scripts used for the analysis and producing the figures in the manuscript. The contents include streamMetabolizer inputs and outputs and additional relevant data needed to generate the main manuscript results. The data included are: total ecosystem respiration, water respiration, calculated sediment-associated respiration, gross primary production outputs from the river corridor model for the Yakima River Basin, median grain size (d50), depth, dissolved oxygen, water temperature, pressure, and annual oxygen consumption. The associated GitHub repository can be found at https://github.com/river-corridors-sfa/SSS_metabolism. Samples collected during this study were labeled as “Second Spatial Study” or “SSS.” Raw time series sensor data, total suspended solids, and depth data from SSS were published at https://data.ess-dive.lbl.gov/datasets/doi:10.15485/1969566. A subset of data from the SSS samples were published in the contiguous United States (CONUS)-Scale Model-Sample (CM) study data package available at https://data.ess-dive.lbl.gov/view/doi:10.15485/1923689 that presents data from across the CONUS. They include dissolved organic carbon (DOC, measured as non-purgeable organic carbon, NPOC), total nitrogen (TN), grain size, aerobic sediment respiration, dissolved oxygen (DO), and temperature. Parent IDs and Site IDs are consistent between the SSS and CM data packages, and they can be mapped directly so data across packages can be used together. Field metadata for the samples in this da This dataset is comprised of one main data folder with four subfolders. The main data folder contains of (1) file-level metadata; (2) data dictionary; (3) total/water column/sediment respiration; (4) gross primary production (GPP); (5) median grain size (d50); and (6) annual oxygen consumption. The “Figures” subfolder contains the figures used in the paper and all intermediate files (including geospatial files). The “Published_Data” contains a readme directing the user to download the public data to reproduce analyses and figures. The “Scripts” folder contains all scripts used in the analyses that were not part of running StreamMetabolizer. Lastly, the “Stream_Metabolizer” folder contains all files associated with running StreamMetabolizer including (1) model input files, (2) model output files, (3) processing scripts, (4) histogram plots of the outputs, and (5) an R project. All files are .csv, .pdf, .R, .Rmd, .Rproj, .html, .png, .txt, .qgz, .cpg, .dbf, .prj, .shp, .shp.ea.iso.xml, .shp.iso.xml, .shx, .sbn. ta package can be found at either link. We acknowledge the Yakama Nation as owners and caretakers of the lands where we collected these data. We thank the Confederated Tribes and Bands of the Yakama Nation Tribal Council and Yakama Nation Fisheries for working with us to facilitate sample collection and optimization of data usage according to their values and worldview.

54 ENVIRONMENTAL SCIENCES↗

Bacterial community dynamics as a result of growth-yield trade-off and multispecies metabolic interactions toward understanding the gut biofilm niche

Abstract Bacterial communities are ubiquitous, found in natural ecosystems, such as soil, and within living organisms, like the human microbiome. The dynamics of these communities in diverse environments depend on factors such as spatial features of the microbial niche, biochemical kinetics, and interactions among bacteria. Moreover, in many systems, bacterial communities are influenced by multiple physical mechanisms, such as mass transport and detachment forces. One example is gut mucosal communities, where dense, closely packed communities develop under the concurrent influence of nutrient transport from the lumen and fluid-mediated detachment of bacteria. In this study, we model a mucosal niche through a coupled agent-based and finite-volume modeling approach. This methodology enables us to model bacterial interactions affected by nutrient release from various sources while adjusting individual bacterial kinetics. We explored how the dispersion and abundance of bacteria are influenced by biochemical kinetics in different types of metabolic interactions, with a particular focus on the trade-off between growth rate and yield. Our findings demonstrate that in competitive scenarios, higher growth rates result in a larger share of the niche space. In contrast, growth yield plays a critical role in neutralism, commensalism, and mutualism interactions. When bacteria are introduced sequentially, they cause distinct spatiotemporal effects, such as deeper niche colonization in commensalism and mutualism scenarios driven by species intermixing effects, which are enhanced by high growth yields. Moreover, sub-ecosystem interactions dictate the dynamics of three-species communities, sometimes yielding unexpected outcomes. Competitive, fast-growing bacteria demonstrate robust colonization abilities, yet they face challenges in displacing established mutualistic systems. Bacteria that develop a cooperative relationship with existing species typically obtain niche residence, regardless of their growth rates, although higher growth yields significantly enhance their abundance. Our results underscore the importance of bacterial niche dynamics in shaping community properties and succession, highlighting a new approach to manipulating microbial systems.

Microbiology↗