Search NASA⌕ Search

SEARCH · Search NASA

Results for “nutrient network”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Small Signaling Peptides in Sorghum bicolor : Integrating Phylogeny and Gene Expression to Characterize Roles in Stem Development

Small signaling peptides (SSPs) are important regulators of plant growth, development, and responses to biotic and abiotic stress, yet their role in the C4 grass Sorghum bicolor is largely uncharacterized. To help fill this knowledge gap, 219 sorghum genes that encode SSPs were identified based on SSP sequences previously identified in Arabidopsis thaliana, Zea mays, Oryza sativa, Triticum aestivum , and Brachypodium distachyon . The 219 sorghum SSP-encoding genes were assigned to 19 gene families, analyzed for the presence of motifs, and aligned with genes that encode SSPs in other plants using phylogenetic analysis. Sorghum genes in 12 of the 19 SSP gene families had not been previously characterized. Expression of the 219 SSP-encoding genes in sorghum organs, during stem development, and in stem tissues and cell types revealed distinct spatial, temporal, and developmental patterns of expression. Genes associated with the SbCEP and SbRGF families were preferentially expressed in roots, whereas SbEPF genes were expressed in stem epidermal and pith parenchyma cells and panicles. The expression of genes during bioenergy sorghum stem growth and development was investigated because stems account for ~80% of harvested biomass and serve as conduits for water and nutrient transport between leaves and roots. During stem development, 28 SSP genes in several families ( CLE, EPF, CEP, GASS, PSY, ES, PSK, CAPE, POE ) were expressed at higher levels in zones of cell proliferation. For example, the TDIF homologs SbCLE41 and SbCLE42 were expressed at high levels in nascent stem nodes where they may regulate vascular bundle cambial activity and cell differentiation. A different set of 15 genes in the CIF, POE, CAPE, PSY, CEP, RALF , and CLE families were expressed at higher levels in zones of stem tissue differentiation highlighted by elevated expression of five SbRALFR s in the stem nodal plexus. Cell type–specific expression of many sorghum genes that encode SSPs was observed in fully elongated internodes indicating gene expression is regulated with high spatial resolution. Overall, the results provide a foundation of information for analysis of SSP function in sorghum that can be integrated with knowledge of sorghum gene regulatory networks to modulate traits important for production of sorghum crops.

bioenergy sorghum↗
DOE OSTI2021

All in the timing: epigenetic control of greening

During germination, the emerging seedling relies on a limited supply of energy and nutrients stored within itself and the seed. This often occurs underground in darkness where seedlings undergo skotomorphogenesis, characterized by hypocotyl elongation, reduced root growth, and small closed cotyledons forming an apical hook to protect the shoot apical meristem as the seedling grows though the soil towards sunlight. Such long and pale dark-grown seedlings are commonly referred to as etiolated (from the French étiolier (to make pale)). However, a seedling must eventually find sunlight and make an irreversible shift to photoautotrophic growth (photomorphogenesis or de-etiolation). This leads to a dramatic change in physiology, including reduced hypocotyl elongation, increased root growth, expansion and unfolding of the cotyledons, and the biogenesis of photosynthetic plastids (i.e. chloroplasts). Furthermore, this shift is accompanied by a massive change in gene expression, with about one-third of genes being reprogrammed. A large amount of the transcriptional response is regulated by the plant photoreceptors (phytochromes and cryptochromes) that activate transcriptional networks. Furthermore, this involves the rapid degradation of PHYTOCHROME INTERACTING FACTORS (PIFs), which are basic helix-loop-helix domain-containing transcription factors that repress photomorphogenesis in the dark.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Revealing systematic changes in the transcriptome during the transition from exponential growth to stationary phase

ABSTRACT The composition of bacterial transcriptomes is determined by the transcriptional regulatory network (TRN). The TRN regulates the transition from one physiological state to another. Here, we use independent component analysis to monitor the composition of the transcriptome during the transition from the exponential growth phase to the stationary phase. With Escherichia coli K-12 MG1655 as a model strain, we trigger the transition using carbon, nitrogen, and sulfur starvation. We find that (i) the transition to the stationary phase accompanies common transcriptome changes, including increased stringent responses and reduced production of cellular building blocks and energy regardless of the limiting element; (ii) condition-specific changes are strongly associated with transcriptional regulators ( e.g. , Crp, NtrC, CysB, Cbl) responsible for metabolizing the limiting element; and (iii) the shortage of each limiting element differentially affects the production of amino acids and extracellular polymers. This study demonstrates how the combination of genome-scale datasets and new data analytics reveals the fundamental characteristics of a key transition in the life cycle of bacteria. IMPORTANCE Nutrient limitations are critical environmental perturbations in bacterial physiology. Despite its importance, a detailed understanding of how bacterial transcriptomes are adjusted has been limited. By utilizing independent component analysis (ICA) to decompose transcriptome data, this study reveals key regulatory events that enable bacteria to adapt to nutrient limitations. The findings not only highlight common responses, such as the stringent response, but also condition-specific regulatory shifts associated with carbon, nitrogen, and sulfur starvation. The insights gained from this work advance our knowledge of bacterial physiology, gene regulation, and metabolic adaptation.

Lim, Hyun Gyu (ORCID:0000000204692388)↗
DOE OSTI2023

Maximum respiration rates in hyporheic zone sediments are primarily constrained by organic carbon concentration and secondarily by organic matter chemistry

Abstract. River corridors are fundamental components of the Earth system, and their biogeochemistry can be heavily influenced by processes in subsurface zones immediately below the riverbed, referred to as the hyporheic zone. Within the hyporheic zone, organic matter (OM) fuels microbial respiration, and OM chemistry heavily influences aerobic and anaerobic biogeochemical processes. The link between OM chemistry and respiration has been hypothesized to be mediated by OM molecular diversity, whereby respiration is predicted to decrease with increasing diversity. Here we test the specific prediction that aerobic respiration rates will decrease with increases in the number of unique organic molecules (i.e., OM molecular richness, as a measure of diversity). We use publicly available data across the United States from crowdsourced samples taken by the Worldwide Hydrobiogeochemical Observation Network for Dynamic River Systems (WHONDRS) consortium. Our continental-scale analyses rejected the hypothesis of a direct limitation of respiration by OM molecular richness. In turn, we found that organic carbon (OC) concentration imposes a primary constraint over hyporheic zone respiration, with additional potential influences of OM richness. We specifically observed respiration rates to decrease nonlinearly with the ratio of OM richness to OC concentration. This relationship took the form of a constraint space with respiration rates in most systems falling below the constraint boundary. A similar, but slightly weaker, constraint boundary was observed when relating respiration rate to the inverse of OC concentration. These results indicate that maximum respiration rates may be governed primarily by OC concentration, with secondary influences from OM richness. Our results also show that other variables often suppress respiration rates below the maximum associated with the richness-to-concentration ratio. An important focus of future research will identify physical (e.g., sediment grain size), chemical (e.g., nutrient concentrations), and/or biological (e.g., microbial biomass) factors that suppress hyporheic zone respiration below the constraint boundaries observed here.

58 GEOSCIENCES↗
DOE OSTI2025

Hydrological connectivity: a review and emerging strategies for integrating measurement, modeling, and management

This review synthesizes methods for measuring, modeling, and managing hydrologic connectivity, offering pathways to improve practices and address environmental challenges (e.g., climate change) and sustainability. As a key driver of water movement and nutrient cycling, hydrologic connectivity influences flood mitigation, water quality regulation, and biodiversity conservation. However, traditional field-based methods (e.g., dye tracing), indirect measurements (e.g., runoff analysis), and remote sensing techniques (e.g., InSAR) often struggle to capture the complexity of catchment-scale interactions. Similarly, modeling approaches—including process-based and percolation theory-based models, graph theory, and entropy-based metrics—face limitations in fully representing these interconnected processes. Both modeling and measurement techniques are constrained by inadequate spatial and temporal coverage, high data demands, computational complexity, and difficulties in representing subsurface connectivity. Subsequently, we critique current management practices that prioritize isolated variables (e.g., streamflow, sediment transport) over system-wide strategies and emphasize the need for adaptive, connectivity-based approaches in water resource planning and restoration. Moving forward, we highlight the importance of interdisciplinary collaboration, technological innovations (e.g., AI-driven modeling, real-time monitoring), and integrated frameworks to improve connectivity measurement, modeling, and adaptive management to restore fragmented hydrologic networks. This integrated approach sets the stage for transformative water resource management, fostering proactive policy development and stakeholder engagement.

Dwivedi, Dipankar↗
DOE OSTI2023

Graphene oxide exposure alters gut microbial community composition and metabolism in an in vitro human model

Graphene oxide (GO) nanomaterials have unique physicochemical properties that make them highly promising for biomedical, environmental, and agricultural applications. Despite the increasing interest and the use of GO, assessments of its nanotoxicity have largely not interrogated its potential impact on the gut microbiome. This study addresses an important knowledge gap by investigating the impact of GO exposure- both at low (25 ppm) and high (250ppm) doses and fed (nutrient rich) and fasted (nutrient deplete) conditions- on the gut microbial community structure and function, using an in vitro human colon bioreactor model. 16S rRNA amplicon sequencing revealed that GO exposure resulted in a restructuring of community composition. 25 ppm GO induced a marked decrease in the Bacteroidota phylum and increased the ratio of Firmicutes to Bacteroidota (F/B). Untargeted metabolomics on the supernatants indicated that 25 ppm GO impaired microbial utilization and metabolism of substrates (amino acids, carbohydrate metabolites) and reduced production of beneficial microbial metabolites such as 5-hydroxyindole-3-acetic acid and GABA. Exposure to 250 ppm GO resulted in community composition and metabolome profiles that were very similar to the controls that lacked both GO and digestive enzymes, suggesting that high concentrations of GO may interact with digestive enzymes to form protein coronas, causing their depletion in the gut environment. Differential abundance analyses revealed that 3 genera from the phylum Bacteroidota (Bacteroides, Dysgonomonas, and Parabacteroides) were more abundant after 250 ppm GO exposure, irrespective of feed state. Integrative correlation network analysis indicated that the phylum Bacteroidota showed strong positive correlations to multiple microbial metabolites including GABA and 3-indoleacetic acid, are much larger number of correlations compared to other phyla. These results show that GO exposure has a significant impact on gut microbial community composition and metabolism and different mechanisms are at play for low and high GO concentrations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Continental-scale patterns of extracellular enzyme activity in the subsoil: an overlooked reservoir of microbial activity

Abstract Chemical stabilization of microbial-derived products such as extracellular enzymes (EE) onto mineral surfaces has gained attention as a possibly important mechanism leading to the persistence of soil organic carbon (SOC). While the controls on EE activities and their stabilization in the surface soil are reasonably well-understood, how these activities change with soil depth and possibly diverge from those at the soil surface due to distinct physical, chemical, and biotic conditions remains unclear. We assessed EE activity to a depth of 1 m (10 cm increments) in 19 soil profiles across the Critical Zone Observatory Network, which represents a wide range of climates, soil orders, and vegetation types. For all EEs, activities per mass of soil correlated positively with microbial biomass (MB) and SOC, and all three of these variables decreased logarithmically with depth ( p < 0.05). Across all sites, over half of the potential EE activities per mass soil consistently occurred below 20 cm for all measured EEs. Activities per unit MB or SOC were substantially higher at depth (soils below 20 cm accounted for 80% of whole-profile EE activity), suggesting an accumulation of stabilized (i.e. mineral sorbed) EEs in subsoil horizons. The pronounced enzyme stabilization in subsurface horizons was corroborated by mixed-effects models that showed a significant, positive relationship between clay concentration and MB-normalized EE activities in the subsoil. Furthermore, the negative relationships between soil C, N, and P and C-, N-, and P-acquiring EEs found in the surface soil decoupled below 20 cm, which could have also been caused by EE stabilization. This finding suggests that EEs may not reflect soil nutrient availabilities deeper in the soil profile. Taken together, our results suggest that deeper soil horizons hold a significant reservoir of EEs, and that the controls of subsoil EEs differ from their surface soil counterparts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Dominant Controls on Preferential Flow and Their Implications for Future Soil Water Fluxes

Abstract Soil water flow, particularly preferential flow (PF), is a critical control on hydrological and biogeochemical processes, including groundwater recharge, contaminant transport, and carbon cycling. However, it remains challenging to predict PF occurrence across large environmental gradients. Here, we developed a deep learning (DL) model to estimate event‐scale soil water flow velocity and the probability of PF occurrence using high‐frequency soil moisture and precipitation data from 33 sites across the National Ecological Observatory Network. The model demonstrated high skill in predicting the binary occurrence of PF (91% F1‐score; 85% accuracy) but the performance was limited in predicting soil water velocity ( R 2 = 0.31). We found that precipitation characteristics (duration, volume, and intensity) were the most important predictors for soil water velocity. Among the non‐precipitation event variables, sand content showed relatively high predictive skill, though differences among non‐event climate variables were generally modest. Lower sand content was associated with increased predicted soil water velocity, a finding that highlights the role of soil structure in producing more non‐uniform flow, which contrasts with traditional uniform flow models. Projecting a reduced DL model under both moderate and high‐emissions future climate scenarios (2060–2099 Representative Concentration Pathways 4.5 and 8.5), we found ∼7.3% increase under RCP4.5 and ∼15% under RCP8.5 of soil water velocities compared to the historical simulation, while modeled likelihood of PF changed little. These findings suggest climate change is not making PF more frequent, but it is making existing PF pathways more efficient with important consequences for associated nutrient and contaminant transport under climate change. Plain Language Summary Water movement in soil is critical for water quality. While often modeled as a uniform flow process, in reality water moves rapidly through cracks and burrows in what is called “preferential flow” (PF), which limits natural filtration and can transport pollutants. We developed a deep learning model, trained on data from 33 U.S. sites, to predict when and how fast this PF occurs based on precipitation, soil, and climate data. The model showed that precipitation characteristics (duration, intensity, volume) were the most important predictors of PF. Lower soil sand content/higher clay content was associated with faster water flow, likely due to clay soils forming aggregates and cracks that water moves through rather than infiltrating uniformly. Further analyses based on climate projections suggest that the speed at which PF occurs will become more rapid under future climate scenarios compared to historical simulation. This highlights the need to represent PF in soil water models when assessing future water quality. Key Points The effect of precipitation peak intensity on soil water velocities declined with increasing precipitation intensity Antecedent soil moisture failed to predict preferential flow (PF), contrasting the high predictive power of sand content Climate predictions suggest that soil water velocities through PF paths will increase ∼15% by 2099

Li, Bonan↗
DOE OSTI2022

Chemotaxonomic patterns in intracellular metabolites of marine microbial plankton

Most biological diversity on Earth is contained within microbial communities. In the ocean, these communities dominate processes related to carbon fixation and nutrient recycling. Yet, specific factors that determine community composition and metabolic activity are difficult to resolve in complex microbial populations, complicating predictions of microbial processes in a changing ocean. Microbial metabolism generates small organic molecules that reflect both the biochemical and physiological diversity as well as the taxonomic specificity of these biological processes. These small molecules serve as the conduit for taxon-specific signaling and exchange. Here, we use liquid chromatography-mass spectrometry (LC-MS)-based metabolomics to taxonomically categorize 111 metabolites that include small molecules in central and secondary metabolism across 42 taxa representing numerically dominant and metabolically important lineages of microbial autotrophs and heterotrophs. Patterns in metabolite presence-absence broadly reflected taxonomic lineages. A subset of metabolites that includes osmolytes, sulfur-containing metabolites, sugars, and amino acid derivatives provided chemotaxonomic information among phytoplankton taxa. A variety of phytohormones and signaling molecules were predominantly found in the heterotrophic bacteria and archaea, expanding knowledge of metabolites implicated in modulating interactions between microbes. This chemotaxonomic inventory of marine microbial metabolites is a key step in deciphering metabolic networks that influence ocean biogeochemical cycles.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Spatial turnover of soil viral populations and genotypes overlain by cohesive responses to moisture in grasslands

Viruses shape microbial communities, food web dynamics, and carbon and nutrient cycling in diverse ecosystems. However, little is known about the patterns and drivers of viral community composition, particularly in soil, precluding a predictive understanding of viral impacts on terrestrial habitats. To investigate soil viral community assembly processes, here we analyzed 43 soil viromes from a rainfall manipulation experiment in a Mediterranean grassland in California. We identified 5,315 viral populations (viral operational taxonomic units [vOTUs] with a representative sequence ≥10 kbp) and found that viral community composition exhibited a highly significant distance–decay relationship within the 200-m 2 field site. This pattern was recapitulated by the intrapopulation microheterogeneity trends of prevalent vOTUs (detected in ≥90% of the viromes), which tended to exhibit negative correlations between spatial distance and the genomic similarity of their predominant allelic variants. Although significant spatial structuring was also observed in the bacterial and archaeal communities, the signal was dampened relative to the viromes, suggesting differences in local assembly drivers for viruses and prokaryotes and/or differences in the temporal scales captured by viromes and total DNA. Despite the overwhelming spatial signal, evidence for environmental filtering was revealed in a protein-sharing network analysis, wherein a group of related vOTUs predicted to infect actinobacteria was shown to be significantly enriched in low-moisture samples distributed throughout the field. Overall, our results indicate a highly diverse, dynamic, active, and spatially structured soil virosphere capable of rapid responses to changing environmental conditions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Measurement of Volatile Compounds for Real-Time Analysis of Soil Microbial Metabolic Response to Simulated Snowmelt

Snowmelt dynamics are a significant determinant of microbial metabolism in soil and regulate global biogeochemical cycles of carbon and nutrients by creating seasonal variations in soil redox and nutrient pools. With an increasing concern that climate change accelerates both snowmelt timing and rate, obtaining an accurate characterization of microbial response to snowmelt is important for understanding biogeochemical cycles intertwined with soil. However, observing microbial metabolism and its dynamics non-destructively remains a major challenge for systems such as soil. Microbial volatile compounds (mVCs) emitted from soil represent information-dense signatures and when assayed non-destructively using state-of-the-art instrumentation such as Proton Transfer Reaction-Time of Flight-Mass Spectrometry (PTR-TOF-MS) provide time resolved insights into the metabolism of active microbiomes. In this study, we used PTR-TOF-MS to investigate the metabolic trajectory of microbiomes from a subalpine forest soil, and their response to a simulated wet-up event akin to snowmelt. Using an information theory approach based on the partitioning of mutual information, we identified mVC metabolite pairs with robust interactions, including those that were non-linear and with time lags. The biological context for these mVC interactions was evaluated by projecting the connections onto the Kyoto Encyclopedia of Genes and Genomes (KEGG) network of known metabolic pathways. Simulated snowmelt resulted in a rapid increase in the production of trimethylamine (TMA) suggesting that anaerobic degradation of quaternary amine osmo/cryoprotectants, such as glycine betaine, may be important contributors to this resource pulse. Unique and synergistic connections between intermediates of methylotrophic pathways such as dimethylamine, formaldehyde and methanol were observed upon wet-up and indicate that the initial pulse of TMA was likely transformed into these intermediates by methylotrophs. Increases in ammonia oxidation signatures (transformation of hydroxylamine to nitrite) were observed in parallel, and while the relative role of nitrifiers or methylotrophs cannot be confirmed, the inferred connection to TMA oxidation suggests either a direct or indirect coupling between these processes. Overall, it appears that such mVC time-series from PTR-TOF-MS combined with causal inference represents an attractive approach to non-destructively observe soil microbial metabolism and its response to environmental perturbation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Linking metabolic phenotypes to pathogenic traits among “ Candidatus Liberibacter asiaticus” and its hosts

Candidatus Liberibacter asiaticus (CLas) has been associated with Huanglongbing, a lethal vector-borne disease affecting citrus crops worldwide. While comparative genomics has provided preliminary insights into the metabolic capabilities of this uncultured microorganism, a comprehensive functional characterization is currently lacking. Here, we reconstructed and manually curated genome-scale metabolic models for the six CLas strains A4, FL17, gxpsy, Ishi-1, psy62, and YCPsy, in addition to a model of the closest related culturable microorganism, L. crescens BT-1. Predictions about nutrient requirements and changes in growth phenotypes of CLas were confirmed using in vitro hairy root-based assays, while the L. crescens BT-1 model was validated using cultivation assays. Host-dependent metabolic phenotypes were revealed using expression data obtained from CLas-infected citrus trees and from the CLas-harboring psyllid Diaphorina citri Kuwayama. These results identified conserved and unique metabolic traits, as well as strain-specific interactions between CLas and its hosts, laying the foundation for the development of model-driven Huanglongbing management strategies.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

A continental scale analysis reveals widespread root bimodality

An improved understanding of root vertical distribution is crucial for assessing plant-soil-atmosphere interactions and their influence on the land carbon sink. Here, we analyze a continental-scale dataset of fine roots reaching 2 meters depth, spanning from Alaskan tundra to Puerto Rican forests. Contrary to the expectation that fine root abundance decays exponentially with depth, we found root bimodality at ~20% of 44 sites, with secondary biomass peaks often below 1m. Root bimodality was more likely in areas with low total fine root biomass and was more frequent in shrublands than grasslands. Notably, secondary peaks coincided with high soil nitrogen content at depth. Our analyses suggest that deep soil nutrients tend to be underexploited, while root bimodality offers plants a mechanism to tap into deep soil resources. Our findings add to the growing recognition that deep soil dynamics are systematically overlooked, and calls for more research attention to this deep frontier in the face of global environmental change.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Bacterial hemophilin homologs and their specific type eleven secretor proteins have conserved roles in heme capture and are diversifying as a family

Cellular life relies on enzymes that require metals, which must be acquired from extracellular sources. Bacteria utilize surface and secreted proteins to acquire such valuable nutrients from their environment. These include the cargo proteins of the type eleven secretion system (T11SS), which have been connected to host specificity, metal homeostasis, and nutritional immunity evasion. This Sec-dependent, Gram-negative secretion system is encoded by organisms throughout the phylum Proteobacteria, including human pathogens Neisseria meningitidis, Proteus mirabilis, Acinetobacter baumannii, and Haemophilus influenzae. Experimentally verified T11SS-dependent cargo include transferrin-binding protein B (TbpB), the hemophilin homologs heme receptor protein C (HrpC), hemophilin A (HphA), the immune evasion protein factor-H binding protein (fHbp), and the host symbiosis factor nematode intestinal localization protein C (NilC). Here, we examined the specificity of T11SS systems for their cognate cargo proteins using taxonomically distributed homolog pairs of T11SS and hemophilin cargo and explored the ligand binding ability of those hemophilin cargo homologs. In vivo expression in Escherichia coli of hemophilin homologs revealed that each is secreted in a specific manner by its cognate T11SS protein. Sequence analysis and structural modeling suggest that all hemophilin homologs share an N-terminal ligand-binding domain with the same topology as the ligand-binding domains of the Haemophilus haemolyticus heme binding protein (Hpl) and HphA. We term this signature feature of this group of proteins the hemophilin ligand-binding domain. Network analysis of hemophilin homologs revealed five subclusters and representatives from four of these showed variable heme-binding activities, which, combined with sequence-structure variation, suggests that hemophilins are diversifying in function.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Expanding standards in viromics: in silico evaluation of dsDNA viral genome identification, classification, and auxiliary metabolic gene curation

Viruses influence global patterns of microbial diversity and nutrient cycles. Though viral metagenomics (viromics), specifically targeting dsDNA viruses, has been critical for revealing viral roles across diverse ecosystems, its analyses differ in many ways from those used for microbes. To date, viromics benchmarking has covered read pre-processing, assembly, relative abundance, read mapping thresholds and diversity estimation, but other steps would benefit from benchmarking and standardization. Here we use in silico-generated datasets and an extensive literature survey to evaluate and highlight how dataset composition (i.e., viromes vs bulk metagenomes) and assembly fragmentation impact (i) viral contig identification tool, (ii) virus taxonomic classification, and (iii) identification and curation of auxiliary metabolic genes (AMGs). The in silico benchmarking of five commonly used virus identification tools show that gene-content-based tools consistently performed well for long (≥3 kbp) contigs, while k -mer- and blast-based tools were uniquely able to detect viruses from short (≤3 kbp) contigs. Notably, however, the performance increase of k -mer- and blast-based tools for short contigs was obtained at the cost of increased false positives (sometimes up to ~5% for virome and ~75% bulk samples), particularly when eukaryotic or mobile genetic element sequences were included in the test datasets. Furthermore, for viral classification, variously sized genome fragments were assessed using gene-sharing network analytics to quantify drop-offs in taxonomic assignments, which revealed correct assignations ranging from ~95% (whole genomes) down to ~80% (3 kbp sized genome fragments). A similar trend was also observed for other viral classification tools such as VPF-class, ViPTree and VIRIDIC, suggesting that caution is warranted when classifying short genome fragments and not full genomes. Finally, we highlight how fragmented assemblies can lead to erroneous identification of AMGs and outline a best-practices workflow to curate candidate AMGs in viral genomes assembled from metagenomes. Together, these benchmarking experiments and annotation guidelines should aid researchers seeking to best detect, classify, and characterize the myriad viruses ‘hidden’ in diverse sequence datasets.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Meta Biome: a multiscale model integrating agent-based and metabolic networks to reveal spatial regulation in gut mucosal microbial communities

ABSTRACT Mucosal microbial communities (MMCs) are complex ecosystems near the mucosal layers of the gut essential for maintaining health and modulating disease states. Despite advances in high-throughput omics technologies, current methodologies struggle to capture the dynamic metabolic interactions and spatiotemporal variations within MMCs. In this work, we presentMetaBiome, a multiscale model integrating agent-based modeling (ABM), finite volume methods, and constraint-based models to explore the metabolic interactions within these communities. Integrating ABM allows for the detailed representation of individual microbial agents each governed by rules that dictate cell growth, division, and interactions with their surroundings. Through a layered approach—encompassing microenvironmental conditions, agent information, and metabolic pathways—we simulated different communities to showcase the potential of the model. Using ourin-silicoplatform, we explored the dynamics and spatiotemporal patterns of MMCs in the proximal small intestine and the cecum, simulating the physiological conditions of the two gut regions. Our findings revealed how specific microbes adapt their metabolic processes based on substrate availability and local environmental conditions, shedding light on spatial metabolite regulation and informing targeted therapies for localized gut diseases.MetaBiome provides a detailed representation of microbial agents and their interactions, surpassing the limitations of traditional grid-based systems. This work marks a significant advancement in microbial ecology, as it offers new insights into predicting and analyzing microbial communities. IMPORTANCE Our study presents a novel multiscale model that combines agent-based modeling, finite volume methods, and genome-scale metabolic models to simulate the complex dynamics of mucosal microbial communities in the gut. This integrated approach allows us to capture spatial and temporal variations in microbial interactions and metabolism that are difficult to study experimentally. Key findings from our model include the following: (i) prediction of metabolic cross-feeding and spatial organization in multi-species communities, (ii) insights into how oxygen gradients and nutrient availability shape community composition in different gut regions, and (iii) identification of spatiallyregulated metabolic pathways and enzymes inE. coli. We believe this work represents a significant advance in computational modeling of microbial communities and provides new insights into the spatial regulation of gut microbiome metabolism. The multiscale modeling approach we have developed could be broadly applicable for studying other complex microbial ecosystems.

Microbiology↗
DOE OSTI2020

A multirate mass transfer model to represent the interaction of multicomponent biogeochemical processes between surface water and hyporheic zones (SWAT-MRMT-R 1.0)

Abstract. Surface water quality along river corridors can be modulated by hyporheic zones (HZs) that are ubiquitous and biogeochemically active. Watershed management practices often ignore the potentially important role of HZs as a natural reactor. To investigate the effect of hydrological exchange and biogeochemical processes on the fate of nutrients in surface water and HZs, a novel model, SWAT-MRMT-R, was developed coupling the Soil and Water Assessment Tool (SWAT) watershed model and the reaction module from a flow and reactive transport code (PFLOTRAN). SWAT-MRMT-R simulates concurrent nonlinear multicomponent biogeochemical reactions in both the channel water and its surrounding HZs, connecting the channel water and HZs through hyporheic exchanges using multirate mass transfer (MRMT) representation. Within the model, HZs are conceptualized as transient storage zones with distinguished exchange rates and residence times. The biogeochemical processes within HZs are different from those in the channel water. Hyporheic exchanges are modeled as multiple first-order mass transfers between the channel water and HZs. As a numerical example, SWAT-MRMT-R is applied to the Hanford Reach of the Columbia River, a large river in the United States, focusing on nitrate dynamics in the channel water. Major nitrate contaminants entering the Hanford Reach include those from the legacy waste, irrigation return flows (irrigation water that is not consumed by crops and runs off as point sources to the stream), and groundwater seepage resulting from irrigated agriculture. A two-step reaction sequence for denitrification and an aerobic respiration reaction is assumed to represent the biogeochemical transformations taking place within the HZs. The spatially variable hyporheic exchange rates and residence times in this example are estimated with the basin-scale Networks with EXchange and Subsurface Storage (NEXSS) model. Our simulation results show that (1), given a residence time distribution, how the exchange fluxes to HZs are approximated when using MRMT can significantly change the amount of nitrate consumption in HZs through denitrification and (2) source locations of nitrate have a different impact on surface water quality due to the spatially variable hyporheic exchanges.

58 GEOSCIENCES↗
DOE OSTI2022

Molecular Regulatory Mechanisms Drive Emergent Pathogenetic Properties of Neisseria gonorrhoeae

Neisseria gonorrhoeae is the causative agent of the sexually transmitted infection (STI) gonorrhea, with an estimated 87 million annual cases worldwide. N. gonorrhoeae predominantly colonizes the male and female genital tract (FGT). In the FGT, N. gonorrhoeae confronts fluctuating levels of nutrients and oxidative and non-oxidative antimicrobial defenses of the immune system, as well as the resident microbiome. One mechanism utilized by N. gonorrhoeae to adapt to this dynamic FGT niche is to modulate gene expression primarily through DNA-binding transcriptional regulators. Here, we describe the major N. gonorrhoeae transcriptional regulators, genes under their control, and how these regulatory processes lead to pathogenic properties of N. gonorrhoeae during natural infection. We also discuss the current knowledge of the structure, function, and diversity of the FGT microbiome and its influence on gonococcal survival and transcriptional responses orchestrated by its DNA-binding regulators. We conclude with recent multi-omics data and modeling tools and their application to FGT microbiome dynamics. Understanding the strategies utilized by N. gonorrhoeae to regulate gene expression and their impact on the emergent characteristics of this pathogen during infection has the potential to identify new effective strategies to both treat and prevent gonorrhea.

59 BASIC BIOLOGICAL SCIENCES↗