Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “microbial 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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Reparameterizing Litter Decomposition Using a Simplified Monte Carlo Method Improves Litter Decay Simulated by a Microbial Model and Alters Bioenergy Soil Carbon Estimates

Litter decomposition determines soil organic matter (SOM) formation and plant-available nutrient cycles. Therefore, accurate model representation of litter decomposition is critical to improving soil carbon (C) projections of bioenergy feedstocks. Soil C models that simulate microbial physiology (i.e., microbial models) are new to bioenergy agriculture, and their parameterization is often based on small datasets or manual calibration to reach benchmarks. Here, we reparameterized litter decomposition in a microbial soil C model (CORPSE - Carbon, Organisms, Rhizosphere, and Protection in the Soil Environment) using the continental-scale Long-term Inter-site Decomposition Experiment Team (LIDET) dataset which documents decomposition across a range of litter qualities over a decade. Here, we conducted a simplified Monte Carlo simulation that constrained parameter values to reduce computational costs. The LIDET-derived parameters improved modeled C and nitrogen (N) remaining, decomposition rates, and litter mean residence times as compared to Baseline parameters. We applied the LIDET litter decomposition parameters to a microbial bioenergy model (Fixation and Uptake of Nitrogen – Bioenergy Carbon, Rhizosphere, Organisms, and Protection) to examine soil C estimates generated by Baseline and LIDET parameters. LIDET parameters increased estimated soil C in bioenergy feedstocks, with even greater increases under elevated plant inputs (i.e., by increasing residue, N fertilization). This was due to the integrated effects of plant litter quantity, quality, and agricultural practices (tillage, fertilization). Collectively, we developed a simple framework for using large-scale datasets to inform the parameterization of microbial models that impacts projections of soil C for bioenergy feedstocks.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Microbial Models for Simulating Soil Carbon Dynamics: A Review

Abstract Soils store the largest amount of carbon (C) in the biosphere, and the C pool in soil is critical to the global C balance. Numerous microbial models have been developed over the last few decades to represent microbial processes that regulate the responses of soil organic carbon (SOC) to climate change. However, the representation of microbial processes varies, and how microbial processes are incorporated into SOC models has not been well explored. Here, we reviewed 71 microbial models to characterize the microbial processes incorporated into SOC models and analyzed variations in mechanistic complexity. We revealed that (a) four processes (microbial‐mediated decomposition, mineral interaction, microbial necromass recycling, and active and dormant microbial dynamics) are commonly incorporated in microbial models, (b) ∼48% of models simulate only one microbial process (i.e., microbial‐mediated decomposition) and 35% of models simulate two microbial processes: for example, microbial‐mediated decomposition and mineral interaction, (c) more than 80% microbial models use nonlinear equations, such as forward Michaelis‐Menten kinetics, to represent SOC decomposition, (d) the concept of persistence of SOC due to its intrinsic properties has been replaced by organo‐mineral interaction (∼39% of microbial models) that protects SOC from decomposition, and (e) various temperature and moisture modifiers and pH effects have been used to explain the environmental effect on microbial processes. In the future, to realistically incorporate microbial processes into Earth System Models, it is imperative to identify experimental evidence on rate limitation processes and firmly ground model structure on the field and laboratory data.

Environmental Sciences & Ecology↗
DOE OSTI2026

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↗
DOE OSTI2024

Modeling microbial carbon fluxes and stocks in global soils from 1901 to 2016

Abstract. Soil microbes play a crucial role in the carbon (C) cycle; however, they have been overlooked in predicting the terrestrial C cycle. We applied a microbial-explicit Earth system model – the Community Land Model-Microbe (CLM-Microbe) – to investigate the dynamics of soil microbes during 1901 to 2016. The CLM-Microbe model was able to reproduce the variations of gross (GPP) and net (NPP) primary productivity, heterotrophic (HR) and soil (SR) respiration, microbial (MBC) biomass C in fungi (FBC) and bacteria (BBC) in the top 30 cm and 1 m, and dissolved (DOC) and soil organic C (SOC) in the top 30 cm and 1 m during 1901–2016. During the study period, simulated C variables increased by approximately 12 PgC yr−1 for HR, 25 PgC yr−1 for SR, 1.0 PgC for FBC and 0.4 PgC for BBC in 0–30 cm, and 1.2 PgC for FBC and 0.7 PgC for BBC in 0–1 m. Increases in microbial C fluxes and pools were widely found, particularly at high latitudes and in equatorial regions, but we also observed their decreases in some grids. Overall, the area-weighted averages of HR, SR, FBC, and BBC in the top 1 m were significantly correlated with those of soil moisture and soil temperature in the top 1 m. These results suggested that microbial C fluxes and pools were jointly governed by vegetation C input and soil temperature and moisture. Our simulations revealed the spatial and temporal patterns of microbial C fluxes and pools in response to environmental change, laying the foundation for an improved understanding of soil microbial roles in the global terrestrial C cycle.

He, Liyuan↗
DOE OSTI2025

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↗
DOE OSTI2024

Scarcity of fixed carbon transfer in a model microbial phototroph–heterotroph interaction

Although the green alga Chlamydomonas reinhardtii has long served as a reference organism, few studies have interrogated its role as a primary producer in microbial interactions. Here, we quantitatively investigated C. reinhardtii’s capacity to support a heterotrophic microbe using the established coculture system with Mesorhizobium japonicum, a vitamin B 12 -producing α-proteobacterium. Using stable isotope probing and nanoscale secondary ion mass spectrometry (nanoSIMS), we tracked the flow of photosynthetic fixed carbon and consequent bacterial biomass synthesis under continuous and diurnal light with single-cell resolution. We found that more 13 C fixed by the alga was taken up by bacterial cells under continuous light, invalidating the hypothesis that the alga’s fermentative degradation of starch reserves during the night would boost M. japonicum heterotrophy. 15 NH 4 assimilation rates and changes in cell size revealed that M. japonicum cells reduced new biomass synthesis in coculture with the alga but continued to divide—a hallmark of nutrient limitation often referred to as reductive division. Despite this sign of starvation, the bacterium still synthesized vitamin B 12 and supported the growth of a B 12 -dependent C. reinhardtii mutant. Finally, we showed that bacterial proliferation could be supported solely by the algal lysis that occurred in coculture, highlighting the role of necromass in carbon cycling. Collectively, these results reveal the scarcity of fixed carbon in this microbial trophic relationship (particularly under environmentally relevant light regimes), demonstrate B 12 exchange even during bacterial starvation, and underscore the importance of quantitative approaches for assessing metabolic coupling in algal–bacterial interactions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Kinetic modeling of microbial growth, enzyme activity, and gene deletions: An integrated model of β-glucosidase function in Cellvibrio japonicus

Understanding the complex growth and metabolic dynamics in microorganisms requires advanced kinetic models containing both metabolic reactions and enzymatic regulation to predict phenotypic behaviors under different conditions and perturbations. Most current kinetic models lack gene expression dynamics and are separately calibrated to distinct media, which consequently makes them unable to account for genetic perturbations or multiple substrates. This challenge limits our ability to gain a comprehensive understanding of microbial processes towards advanced metabolic optimizations that are desired for many biotechnology applications. Here, we present an integrated computational and experimental approach for the development and optimization of mechanistic kinetic models for microbial growth and metabolic and enzymatic dynamics. Our approach integrates growth dynamics, gene expression, protein secretion, and gene-deletion phenotypes. We applied this methodology to build a dynamic model of the growth kinetics in batch culture of the bacterium Cellvibrio japonicus grown using either cellobiose or glucose media. The model parameters were inferred from an experimental data set using an evolutionary computation method. The resulting model was able to explain the growth dynamics of C. japonicus using either cellobiose or glucose media and was also able to accurately predict the metabolite concentrations in the wild-type strain as well as in β-glucosidase gene deletion mutant strains. We validated the model by correctly predicting the non-diauxic growth and metabolite consumptions of the wild-type strain in a mixed medium containing both cellobiose and glucose, made further predictions of mutant strains growth phenotypes when using cellobiose and glucose media, and demonstrated the utility of the model for designing industrially-useful strains. Importantly, the model is able to explain the role of the different β-glucosidases and their behavior under genetic perturbations. Here, this integrated approach can be extended to other metabolic pathways to produce mechanistic models for the comprehensive understanding of enzymatic functions in multiple substrates.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

2020 Multiscale Microbial Dynamics Modeling Course

The 2020 Multiscale Microbial Dynamics course is adapted from the virtual 2020 Mutliscale Microbial Dynamics Summer School that was hosted by Environmental Molecular Sciences Laboratory (EMSL), a U.S. Department of Energy (DOE) science user facility located on the Pacific Northwest National Laboratory (PNNL) campus, in collaboration with the Joint Genome Institute (JGI) and the DOE Systems Biology Knowledgebase (KBase). The course course covers how to incorporate microbial metagenomic and environmental metabolite data from watershed ecosystems into metabolic and community modeling using computational frameworks, such as KBase and PFLOTRAN. The curriculum includes lectures and software and data analysis tutorials. All materials are freely accessible to the community as part of the 2020 Microbial Dynamics Summer School Organization in KBase.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Leveraging observed soil heterotrophic respiration fluxes as a novel constraint on global-scale models

Microbially-explicit models may improve understanding and projections of carbon dynamics in response to future climate change, but their fidelity in simulating global-scale soil heterotrophic respiration (RH), a stringent test for soil biogeochemical models, has never been evaluated. We used statistical global RH products, as well as 7,821 daily site-scale RH measurements, to evaluate the spatio-temporal performance of one first-order decay model (CASA-CNP) and two microbially-explicit biogeochemical models (CORPSE and MIMICS) that were forced by two different climate datasets. CORPSE and MIMICS did not provide any measurable performance improvement; instead, the models were highly sensitive to the meteorological input data used to drive them. Spatial RH variability was generally well simulated except in the northern middle latitudes (~50°N) and arid regions; models captured the seasonal variability of RH well, but showed more divergence in tropic and arctic regions. Our results demonstrate that the next generation of biogeochemical models shows promise, but also needs to be improved for realistic spatio-temporal variability of RH. Finally, we emphasize the importance of net primary production, soil moisture, and soil temperature inputs, and that jointly evaluating soil models for their spatial (global scale) and temporal (site scale) performance provides crucial benchmarks for improving biogeochemical models.

Jian, Jinshi↗
DOE OSTI2021

Soil enzymes as indicators of soil function: A step toward greater realism in microbial ecological modeling

Soil carbon (C) and nitrogen (N) cycles and their complex responses to environmental changes have received increasing attention. However, large uncertainties in model predictions remain, partially due to the lack of explicit representation and parameterization of microbial processes. One great challenge is to effectively integrate rich microbial functional traits into ecosystem modeling for better predictions. Here, using soil enzymes as indicators of soil function, we developed a competitive dynamic enzyme allocation scheme and detailed enzyme-mediated soil inorganic N processes in the Microbial-ENzyme Decomposition (MEND) model. We conducted a rigorous calibration and validation of MEND with diverse soil C-N fluxes, microbial C:N ratios, and functional gene abundances from a 12-year CO 2 × N grassland experiment (BioCON) in Minnesota, USA. In addition to accurately simulating soil CO 2 fluxes and multiple N variables, the model correctly predicted microbial C:N ratios and their negative response to enriched N supply. Model validation further showed that, compared to the changes in simulated enzyme concentrations and decomposition rates, the changes in simulated activities of eight C-N-associated enzymes were better explained by the measured gene abundances in responses to elevated atmospheric CO 2 concentration. In conclusion, our results demonstrated that using enzymes as indicators of soil function and validating model predictions with functional gene abundances in ecosystem modeling can provide a basis for testing hypotheses about microbially mediated biogeochemical processes in response to environmental changes. Further development and applications of the modeling framework presented here will enable microbial ecologists to address ecosystem-level questions beyond empirical observations, toward more predictive understanding, an ultimate goal of microbial ecology.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Exploring the determinants of organic matter bioavailability through substrate-explicit thermodynamic modeling

Microbial decomposition of organic matter (OM) in river corridors is a major driver of nutrient and energy cycles in natural ecosystems. Recent advances in omics technologies enabled high-throughput generation of molecular data that could be used to inform biogeochemical models. With ultrahigh-resolution OM data becoming more readily available, in particular, the substrate-explicit thermodynamic modeling (SXTM) has emerged as a promising approach due to its ability to predict OM degradation and respiration rates from chemical formulae of compounds. This model implicitly assumes that all detected organic compounds are bioavailable, and that aerobic respiration is driven solely by thermodynamics. Despite promising demonstrations in previous studies, these assumptions may not be universally valid because OM degradation is a complex process governed by multiple factors. To identify key drivers of OM respiration, we performed a comprehensive analysis of diverse river systems using Fourier-transform ion cyclotron resonance mass spectrometry OM data and associated respiration measurements collected by the Worldwide Hydrobiogeochemistry Observation Network for Dynamic River Systems (WHONDRS) consortium. In support of our argument, we found that the incorporation of all compounds detected in the samples into the SXTM resulted in a poor correlation between the predicted and measured respiration rates. The data-model consistency was significantly improved by the selective use of a small subset (i.e., only about 5%) of organic compounds identified using an optimization method. Through a subsequent comparative analysis of the subset of compounds (which we presume as bioavailable) against the full set of compounds, we identified three major traits that potentially determine OM bioavailability, including: (1) thermodynamic favorability of aerobic respiration, (2) the number of C atoms contained in compounds, and (2) carbon/nitrogen (C/N) ratio. We found that all three factors serve as “filters” in that the compounds with undesirable properties in any of these traits are strictly excluded from the bioavailable fraction. This work highlights the importance of accounting for the complex interplay among multiple key traits to increase the predictive power of biogeochemical and ecosystem models.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Geomechanical and flow implications with continued bioconversion of coal to methane: Experiments and modeling

Microbial conversion of coal to methane is a promising technology during transition of coal-based energy from conventional coal mining to natural gas recovery from coal. Significant research advances have been made towards engineering ideal microbial communities and nutrients for bio-stimulation of coal. However, actual field applications require geomechanical and flow behavior characterization of coal during the bioconversion process as well as gas production over the life of the created “biogenic gas reservoirs”. This work presents the results of an experimental investigation to estimate gas/methane production using bioconversion of coal and analyses of variations in bulk modulus, strain and permeability. Bulk modulus of coal demonstrated time-dependent behavior with continued bioconversion of solid coal. The phenomenon was modeled using mass balancing in a closed environment and the logistic equation based biogenic gas production with time. The modeled gas production and changes in modulus of coal showed excellent agreement with the corresponding experimental results. Next, numerical simulation of biogenic conversion of coal in a constant stress in-situ condition, replicating the Huff ‘n Puff method under field conditions, was carried out. The results showed that bioconversion can lead to decrease in effective stress, increase in permeability and decrease in the modulus of coal with time. Repeated recharge of coal with nutrients and recovering the produced gas, that is, sequential cycles of Huff n’ Puff, showed further increase in permeability and decrease in coal strength, potentially leading to coal failure in-situ, further increasing the permeability, thus enhancing the prospect of field application of the technology.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2022

How do soil microbes shape ecosystem biogeochemistry in the context of global change?

Earth system models disagree on the future of global soil carbon C, with projections ranging from losses of 70 Pg to gains of 250 Pg by the year 2100 (Todd-Brown et al., 2014). This disagreement is due, in part, to the high degree of uncertainty in the mechanisms that control soil C, the largest dynamic pool of organic C on Earth. It remains unclear how climate change, rising atmospheric CO2 levels, and other environmental changes will combine to affect global soil C. As key engineers of biogeochemical cycling, microbes play a crucial role in soil C responses to these environmental changes. Yet, the path toward representing microbes in soil models remains uncertain. Furthermore, we address four burning questions relevant for modelling microbial ecological, evolutionary, and biogeochemical processes in soils: (1) What is the role of soil microbes in biogeochemical cycling? (2) Does microbial community composition matter for soil carbon cycling? (3) Which ecological and evolutionary processes contribute to functional changes? (4) How should we model microbial eco-evolutionary mechanisms in soil C models?

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Hysteretic temperature sensitivity of wetland CH4 fluxes explained by substrate availability and microbial activity: Model Archive

This Modeling Archive is in support of an NGEE Arctic publication "Hysteretic temperature sensitivity of wetland CH4 fluxes explained by substrate availability and microbial activity" in the Journal Biogeosciences (https://doi.org/10.5194/bg-17-5849-2020), which includes the model data used in the publication. CH4 emissions from terrestrial systems are posited to increase, which can offset mitigation efforts and accelerate climate change. Yet, the accuracy of modeled CH4 emissions is sensitive to the prescribed CH4 production (or emission) temperature dependencies that are currently uncertain. Here, we use a comprehensive biogeochemistry model (ecosys) to investigate factors modulating CH4 production and emission rates across a permafrost thaw gradient encompassing a partly thawed bog and a fully thawed fen. We find that seasonally varying substrate availability drives lower and higher modeled methanogen biomass and activity, and thereby CH4 production, during the earlier and later periods of the thawed season, respectively. Package follows the Model data archiving guidelines with data in a *.zip file with three subfolders containing *.csv files and raw model output files; a data dictionary table (data_dictionary.csv) and two model output description tables (ecosys_plantspecies_output_notes.csv and ecosys_soil_ouput_notes.csv) to explain the format and meaning of individual output variables; and a user guide as a *.pdf. A detailed model description can be found in the supplement of (Grant, 2013). The ecosys source code is available at https://doi:10.5281/zenodo.3906642. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Targeted rRNA depletion enables efficient mRNA sequencing in diverse bacterial species and complex co-cultures

Bacterial mRNA sequencing is inefficient due to the abundance of ribosomal RNA that is challenging to deplete. While commercial kits target rRNA from common bacterial species, they are frequently inefficient when applied to divergent species, including those from environmental isolates. Similarly, other methods typically employ large probe sets that tile the entire length of rRNAs; however, such approaches are infeasible when applied to many species. Therefore, we present EMBR-seq+, which requires fewer than 10 oligonucleotides per rRNA by combining rRNA blocking primers with RNase H-mediated depletion to achieve rRNA removal efficiencies of up to 99% in diverse bacterial species. Furthermore, in more complex microbial co-cultures between Fibrobacter succinogenes strain UWB7 and anaerobic fungi, EMBR-seq+ depleted both bacterial and fungal rRNA, with a fourfold improvement in bacterial rRNA depletion compared with a commercial kit, thereby demonstrating that the method can be applied to non-model microbial mixtures. Notably, for microbes with unknown rRNA sequences, EMBR-seq+ enables rapid iterations in probe design without requiring to start experiments from total RNA. Finally, efficient depletion of rRNA enabled systematic quantification of the reprogramming of the bacterial transcriptome when cultured in the presence of the anaerobic fungi Anaeromyces robustus or Caecomyces churrovis. We observed that F. succinogenes strain UWB7 downregulated several lignocellulose-degrading carbohydrate-active enzymes in the presence of anaerobic gut fungi, suggesting close interactions between two cellulolytic species that specialize in different aspects of biomass breakdown. Thus, EMBR-seq+ enables efficient, cost-effective, and rapid quantification of the transcriptome to gain insights into non-model microbial systems.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

How microscopic epistasis and clonal interference shape the fitness trajectory in a spin glass model of microbial long-term evolution

The adaptive dynamics of evolving microbial populations takes place on a complex fitness landscape generated by epistatic interactions. The population generically consists of multiple competing strains, a phenomenon known as clonal interference. Microscopic epistasis and clonal interference are central aspects of evolution in microbes, but their combined effects on the functional form of the population’s mean fitness are poorly understood. Here, we develop a computational method that resolves the full microscopic complexity of a simulated evolving population subject to a standard serial dilution protocol. Through extensive numerical experimentation, we find that stronger microscopic epistasis gives rise to fitness trajectories with slower growth independent of the number of competing strains, which we quantify with power-law fits and understand mechanistically via a random walk model that neglects dynamical correlations between genes. We show that increasing the level of clonal interference leads to fitness trajectories with faster growth (in functional form) without microscopic epistasis, but leaves the rate of growth invariant when epistasis is sufficiently strong, indicating that the role of clonal interference depends intimately on the underlying fitness landscape. The simulation package for this work may be found at https://github.com/nmboffi/spin_glass_evodyn .

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

How microscopic epistasis and clonal interference shape the fitness trajectory in a spin glass model of microbial long-term evolution

The adaptive dynamics of evolving microbial populations takes place on a complex fitness landscape generated by epistatic interactions. The population generically consists of multiple competing strains, a phenomenon known as clonal interference. Microscopic epistasis and clonal interference are central aspects of evolution in microbes, but their combined effects on the functional form of the population’s mean fitness are poorly understood. Here, we develop a computational method that resolves the full microscopic complexity of a simulated evolving population subject to a standard serial dilution protocol. Through extensive numerical experimentation, we find that stronger microscopic epistasis gives rise to fitness trajectories with slower growth independent of the number of competing strains, which we quantify with power-law fits and understand mechanistically via a random walk model that neglects dynamical correlations between genes. We show that increasing the level of clonal interference leads to fitness trajectories with faster growth (in functional form) without microscopic epistasis, but leaves the rate of growth invariant when epistasis is sufficiently strong, indicating that the role of clonal interference depends intimately on the underlying fitness landscape. The simulation package for this work may be found at https://github.com/nmboffi/spin_glass_evodyn .

59 BASIC BIOLOGICAL SCIENCES↗