Mechanism Across Scales: A Holistic Modeling Framework Integrating Laboratory and Field Studies for Microbial Ecology

Co-exposure to polyethylene fiber and Salmonella enterica serovar Typhimurium alters microbiome and metabolome of in vitro chicken cecal mesocosms

Humans and animals encounter a summation of exposures during their lifetime (the exposome). In recent years, the scope of the exposome has begun to include microplastics. Microplastics (MPs) have increasingly been found in locations, including in animal gastrointestinal tracts, where there could be an interaction with Salmonella enterica serovar Typhimurium, one of the commonly isolated serovars from processed chicken. However, there is limited knowledge on how gut microbiomes are affected by microplastics and if an effect would be exacerbated by the presence of a pathogen. In this study, we aimed to determine if acute exposure to microplastics in vitro altered the gut microbiome membership and activity. The microbiota response to a 24 h co-exposure to Salmonella enterica serovar Typhimurium and/or low-density polyethylene (PE) microplastics in an in vitro broiler cecal model was determined using 16S rRNA amplicon sequencing (Illumina) and untargeted metabolomics. Community sequencing results indicated that PE fiber with and without S. Typhimurium yielded a lower Firmicutes/Bacteroides ratio compared with other treatment groups, which is associated with poor gut health, and overall had greater changes to the cecal microbial community composition. However, changes in the total metabolome were primarily driven by the presence of S. Typhimurium. Additionally, the co-exposure to PE fiber and S. Typhimurium caused greater cecal microbial community and metabolome changes than either exposure alone. Our results indicate that polymer shape is an important factor in effects resulting from exposure. It also demonstrates that microplastic-pathogen interactions cause metabolic alterations to the chicken cecal microbiome in an in vitro chicken cecal mesocosm.

IMPORTANCE

Researching the exposome, a summation of exposure to one's lifespan, will aid in determining the environmental factors that contribute to disease states. There is an emerging concern that microplastic-pathogen interactions in the gastrointestinal tract of broiler chickens may lead to an increase in Salmonella infection across flocks and eventually increased incidence of human salmonellosis cases. In this research article, we elucidated the effects of acute co-exposure to polyethylene microplastics and Salmonella enterica serovar Typhimurium on the ceca microbial community in vitro. Salmonella presence caused strong shifts in the cecal metabolome but not the microbiome. The inverse was true for polyethylene fiber. Polyethylene powder had almost no effect. The co-exposure had worse effects than either alone. This demonstrates that exposure effects to the gut microbial community are contaminant-specific. When combined, the interactions between exposures exacerbate changes to the gut environment, necessitating future experiments studying low-dose chronic exposure effects with in vivo model systems.

Microbiome modeling: a beginner's guide

Microbiomes, comprised of diverse microbial species and viruses, play pivotal roles in human health, environmental processes, and biotechnological applications and interact with each other, their environment, and hosts via ecological interactions. Our understanding of microbiomes is still limited and hampered by their complexity. A concept improving this understanding is systems biology, which focuses on the holistic description of biological systems utilizing experimental and computational methods. An important set of such experimental methods are metaomics methods which analyze microbiomes and output lists of molecular features. These lists of data are integrated, interpreted, and compiled into computational microbiome models, to predict, optimize, and control microbiome behavior. There exists a gap in understanding between microbiologists and modelers/bioinformaticians, stemming from a lack of interdisciplinary knowledge. This knowledge gap hinders the establishment of computational models in microbiome analysis. This review aims to bridge this gap and is tailored for microbiologists, researchers new to microbiome modeling, and bioinformaticians. To achieve this goal, it provides an interdisciplinary overview of microbiome modeling, starting with fundamental knowledge of microbiomes, metaomics methods, common modeling formalisms, and how models facilitate microbiome control. It concludes with guidelines and repositories for modeling. Each section provides entry-level information, example applications, and important references, serving as a valuable resource for comprehending and navigating the complex landscape of microbiome research and modeling.

Reproducible responses of geochemical and microbial successional patterns in the subsurface to carbon source amendment

Carbon amendments designed to remediate environmental contamination lead to substantial perturbations when injected into the subsurface. For the remediation of uranium contamination, carbon amendments promote reducing conditions to allow microorganisms to reduce uranium to an insoluble, less mobile state. However, the reproducibility of these amendments and underlying microbial community assembly mechanisms have rarely been investigated in the field. In this study, two injections of emulsified vegetable oil were performed in 2009 and 2017 to immobilize uranium in the groundwater at Oak Ridge, TN, USA. Our objectives were to determine whether and how the injections resulted in similar abiotic and biotic responses and their underlying community assembly mechanisms. Both injections caused similar geochemical and microbial succession. Uranium, nitrate, and sulfate concentrations in the groundwater dropped following the injection, and specific microbial taxa responded at roughly the same time points in both injections, including Geobacter, Desulfovibrio, and members of the phylum Comamonadaceae, all of which are well established in uranium, nitrate, and sulfate reduction. Both injections induced a transition from relatively stochastic to more deterministic assembly of microbial taxonomic and phylogenetic community structures based on 16S rRNA gene analysis. We conclude that geochemical and microbial successions after biostimulation are reproducible, likely owing to the selection of similar phylogenetic groups in response to EVO injection.

Contribution of Microorganisms with the Clade II Nitrous Oxide Reductase to Suppression of Surface Emissions of Nitrous Oxide

The sources and sinks of nitrous oxide, as control emissions to the atmosphere, are generally poorly constrained for most environmental systems. Initial depth-resolved analysis of nitrous oxide flux from observation wells and the proximal surface within a nitrate contaminated aquifer system revealed high subsurface production but little escape from the surface. To better understand the environmental controls of production and emission at this site, we used a combination of isotopic, geochemical, and molecular analyses to show that chemodenitrification and bacterial denitrification are major sources of nitrous oxide in this subsurface, where low DO, low pH, and high nitrate are correlated with significant nitrous oxide production. Depth-resolved metagenomes showed that consumption of nitrous oxide near the surface was correlated with an enrichment of Clade II nitrous oxide reducers, consistent with a growing appreciation of their importance in controlling release of nitrous oxide to the atmosphere. Our work also provides evidence for the reduction of nitrous oxide at a pH of 4, well below the generally accepted limit of pH 5.

Environmental stress mediates groundwater microbial community assembly

Community assembly describes how different ecological processes shape microbial community composition and structure. How environmental factors impact community assembly remains elusive. Here we sampled microbial communities and >200 biogeochemical variables in groundwater at the Oak Ridge Field Research Center, a former nuclear waste disposal site, and developed a theoretical framework to conceptualize the relationships between community assembly processes and environmental stresses. We found that stochastic assembly processes were critical (>60% on average) in shaping community structure, but their relative importance decreased as stress increased. Dispersal limitation and 'drift' related to random birth and death had negative correlations with stresses, whereas the selection processes leading to dissimilar communities increased with stresses, primarily related to pH, cobalt and molybdenum. Assembly mechanisms also varied greatly among different phylogenetic groups. Our findings highlight the importance of microbial dispersal limitation and environmental heterogeneity in ecosystem restoration and management.

Molecular mechanisms and environmental adaptations of flagellar loss and biofilm growth of Rhodanobacter under environmental stress

Biofilms aid bacterial adhesion to surfaces via direct and indirect mechanisms, and formation of biofilms is considered as an important strategy for adaptation and survival in suboptimal environmental conditions. However, the molecular underpinnings of biofilm formation in subsurface sediment/groundwater ecosystems where microorganisms often experience fluctuations in nutrient input, pH, and nitrate or metal concentrations are underexplored. We examined biofilm formation under different nutrient, pH, metal, and nitrate regimens of 16 Rhodanobacter strains isolated from subsurface groundwater wells spanning diverse levels of pH (3.5 to 5) and nitrates (13.7 to 146 mM). Eight Rhodanobacter strains demonstrated significant biofilm growth under low pH, suggesting adaptations for survival and growth at low pH. Biofilms were intensified under aluminum stress, particularly in strains possessing fewer genetic traits associated with biofilm formation, findings warranting further investigation. Through random barcode transposon-site sequencing (RB-TnSeq), proteomics, use of specific mutants, and transmission electron microscopy analysis, we discovered flagellar loss under aluminum stress, indicating a potential relationship between motility, metal tolerance, and biofilm growth. Comparative genomic analyses revealed the absence of flagella and chemotaxis genes and the presence of a putative type VI secretion system in the highly biofilm-forming strain FW021-MT20. In this study we identified genetic determinants associated with biofilm growth under metal stress in a predominant environmental genus, Rhodanobacter, and identified traits aiding survival and adaptation to contaminated subsurface environments.

Mixed waste contamination selects for a mobile genetic element population enriched in multiple heavy metal resistance genes

Mobile genetic elements (MGEs) like plasmids, viruses, and transposable elements can provide fitness benefits to their hosts for survival in the presence of environmental stressors. Heavy metal resistance genes (HMRGs) are frequently observed on MGEs, suggesting that MGEs may be an important driver of adaptive evolution in environments contaminated with heavy metals. Here, we report the meta-mobilome of the heavy metal-contaminated regions of the Oak Ridge Reservation subsurface. This meta-mobilome was compared with one derived from samples collected from unimpacted regions of the Oak Ridge Reservation subsurface. We assembled 1615 unique circularized DNA elements that we propose to be MGEs. The circular elements from the highly contaminated subsurface were enriched in HMRG clusters relative to those from the nearby unimpacted regions. Additionally, we found that these HMRGs were associated with Gamma and Betaproteobacteria hosts in the contaminated subsurface and potentially facilitate the persistence and dominance of these taxa in this region. Finally, the HMRGs were associated with conjugative elements, suggesting their potential for future lateral transfer. We demonstrate how our understanding of MGE ecology, evolution, and function can be enhanced through the genomic context provided by completed MGE assemblies.

Microbially mediated climate feedbacks from wetland ecosystems

Wetlands are crucial nodes in the carbon cycle, emitting approximately 20% of global CH4 while also sequestering 20%-30% of all soil carbon. Both greenhouse gas fluxes and carbon storage are driven by microbial communities in wetland soils. However, these key players are often overlooked or overly simplified in current global climate models. Here, we first integrate microbial metabolisms with biological, chemical, and physical processes occurring at scales from individual microbial cells to ecosystems. This conceptual scale-bridging framework guides the development of feedback loops describing how wetland-specific climate impacts (i.e., sea level rise in estuarine wetlands, droughts and floods in inland wetlands) will affect future climate trajectories. These feedback loops highlight knowledge gaps that need to be addressed to develop predictive models of future climates capturing microbial contributions. We propose a roadmap connecting environmental scientific disciplines to address these knowledge gaps and improve the representation of microbial processes in climate models. Together, this paves the way to understand how microbially mediated climate feedbacks from wetlands will impact future climate change.

Challenges in upscaling laboratory studies to ecosystems in soil microbiology research

Soil microbiology has entered into the big data era, but the challenges in bridging laboratory-, field-, and model-based studies of ecosystem functions still remain. Indeed, the limitation of factors in laboratory experiments disregards interactions of a broad range of in situ environmental drivers leading to frequent contradictions between laboratory- and field-based studies, which may consequently mislead model development and projections. Upscaling soil microbiology research from laboratory to ecosystems represents one of the grand challenges facing environmental scientists, but with great potential to inform policymakers toward climate-smart and resource-efficient ecosystems. The upscaling is not only a scale problem, but also requires disentangling functional relationships and processes on each level. We point to three potential reasons for the gaps between laboratory- and field-based studies (i.e., spatiotemporal dynamics, sampling disturbances, and plant-soil-microbial feedbacks), and three key issues of caution when bridging observations and model predictions (i.e., across-scale effect, complex-process coupling, and multi-factor regulation). Field-based studies only cover a limited range of environmental variation that must be supplemented by laboratory and mesocosm manipulative studies when revealing the underlying mechanisms. The knowledge gaps in upscaling soil microbiology from laboratory to ecosystems should motivate interdisciplinary collaboration across experimental, observational, theoretic, and modeling research.

Ecophysiological and genomic analyses of a representative isolate of highly abundant Bacillus cereus strains in contaminated subsurface sediments

Bacillus cereus strain CPT56D-587-MTF (CPTF) was isolated from the highly contaminated Oak Ridge Reservation (ORR) subsurface. This site is contaminated with high levels of nitric acid and multiple heavy metals. Amplicon sequencing of the 16S rRNA genes (V4 region) in sediment from this area revealed an amplicon sequence variant (ASV) with 100% identity to the CPTF 16S rRNA sequence. Notably, this CPTF-matching ASV had the highest relative abundance in this community survey, with a median relative abundance of 3.77% and comprised 20%-40% of reads in some samples. Pangenomic analysis revealed that strain CPTF has expanded genomic content compared to other B. cereus species-largely due to plasmid acquisition and expansion of transposable elements. This suggests that these features are important for rapid adaptation to native environmental stressors. We connected genotype to phenotype in the context of the unique geochemistry of the site. These analyses revealed that certain genes (e.g. nitrate reductase, heavy metal efflux pumps) that allow this strain to successfully occupy the geochemically heterogenous microniches of its native site are characteristic of the B. cereus species while others such as acid tolerance are mobile genetic element associated and are generally unique to strain CPTF.

A Defined Medium for Cultivation and Exometabolite Profiling of Soil Bacteria

Exometabolomics is an approach to assess how microorganisms alter, or react to their environments through the depletion and production of metabolites. It allows the examination of how soil microbes transform the small molecule metabolites within their environment, which can be used to study resource competition and cross-feeding. This approach is most powerful when used with defined media that enable tracking of all metabolites. However, microbial growth media have traditionally been developed for the isolation and growth of microorganisms but not metabolite utilization profiling through Liquid Chromatography Tandem Mass Spectrometry (LC-MS/MS). Here, we describe the construction of a defined medium, the Northen Lab Defined Medium (NLDM), that not only supports the growth of diverse soil bacteria but also is defined and therefore suited for exometabolomic experiments. Metabolites included in NLDM were selected based on their presence in R2A medium and soil, elemental stoichiometry requirements, as well as knowledge of metabolite usage by different bacteria. We found that NLDM supported the growth of 108 of the 110 phylogenetically diverse (spanning 36 different families) soil bacterial isolates tested and all of its metabolites were trackable through LC–MS/MS analysis. These results demonstrate the viability and utility of the constructed NLDM medium for growing and characterizing diverse microbial isolates and communities.

A Framework for Investigating Rules of Life by Establishing Zones of Influence

The incredible complexity of biological processes across temporal and spatial scales hampers defining common underlying mechanisms driving the patterns of life. However, recent advances in sequencing, big data analysis, machine learning, and molecular dynamics simulation have renewed the hope and urgency of finding potential hidden rules of life. There currently exists no framework to develop such synoptic investigations. Some efforts aim to identify unifying rules of life across hierarchical levels of time, space, and biological organization, but not all phenomena occur across all the levels of these hierarchies. Instead of identifying the same parameters and rules across levels, we posit that each level of a temporal and spatial scale and each level of biological organization has unique parameters and rules that may or may not predict outcomes in neighboring levels. We define this neighborhood, or the set of levels, across which a rule functions as the zone of influence. Here, we introduce the zone of influence framework and explain using three examples: (a) randomness in biology, where we use a Poisson process to describe processes from protein dynamics to DNA mutations to gene expressions, (b) island biogeography, and (c) animal coloration. The zone of influence framework may enable researchers to identify which levels are worth investigating for a particular phenomenon and reframe the narrative of searching for a unifying rule of life to the investigation of how, when, and where various rules of life operate.

Genomic Features and Pervasive Negative Selection in Rhodanobacter Strains Isolated from Nitrate and Heavy Metal Contaminated Aquifer

Rhodanobacter species dominate in the Oak Ridge Reservation (ORR) subsurface environments contaminated with acids, nitrate, metal radionuclides, and other heavy metals. To uncover the genomic features underlying adaptations to these mixed-waste environments and to guide genetic tool development, we sequenced the whole genomes of eight Rhodanobacter strains isolated from the ORR site. The genome sizes ranged from 3.9 to 4.2 Mb harboring 3,695 to 4,035 protein-coding genes and GC contents approximately 67%. Seven strains were classified as R. denitrificans and one strain, FW510-R12, as R. thiooxydans based on full length 16S rRNA sequences. According to gene annotation, the top two Cluster of Orthologous Groups (COGs) with high pan-genome expansion rates (Pan/Core gene ratio) were “replication, recombination and repair” and “defense mechanisms.” The denitrifying genes had high DNA homologies except the predicted protein structure variances in NosZ. In contrast, heavy metal resistance genes were diverse with between 7 to 34% of them were located in genomic islands, and these results suggested origins from horizontal gene transfer. Analysis of the methylation patterns in four strains revealed the unique 5mC methylation motifs. Most orthologs (78%) had ratios of nonsynonymous to synonymous substitutions (dN/dS) less than one when compared to the type strain 2APBS1, suggesting the prevalence of negative selection. Overall, the results provide evidence for the important roles of horizontal gene transfer and negative selection in genomic adaptation at the contaminated field site. The complex restriction-modification system genes and the unique methylation motifs in Rhodanobacter strains suggest the potential recalcitrance to genetic manipulation.

IMPORTANCE Despite the dominance of Rhodanobacter species in the subsurface of the contaminated Oak Ridge Reservation (ORR) site, very little is known about the mechanisms underlying their adaptions to the various stressors present at ORR. Recently, multiple Rhodanobacter strains have been isolated from the ORR groundwater samples from several wells with varying geochemical properties. Using Illumina, PacBio, and Oxford Nanopore sequencing platforms, we obtained the whole genome sequences of eight Rhodanobacter strains. Comparison of the whole genomes demonstrated the genetic diversity, and analysis of the long nanopore reads revealed the heterogeneity of methylation patterns in strains isolated from the same well. Although all strains contained a complete set of denitrifying genes, the predicted tertiary structures of NosZ differed. The sequence comparison results demonstrate the important roles of horizontal gene transfer and negative selection in adaptation. In addition, these strains may be recalcitrant to genetic manipulation due to the complex restriction-modification systems and methylations.

Microfluidics as an Emerging Platform for Exploring Soil Environmental Processes: A Critical Review

Investigating environmental processes, especially those occurring in soils, calls for innovative and multidisciplinary technologies that can provide insights at the microscale. The heterogeneity, opacity, and dynamics make the soil a "black box" where interactions and processes are elusive. Recently, microfluidics has emerged as a powerful research platform and experimental tool which can create artificial soil micromodels, enabling exploring soil processes on a chip. Micro/nanofabricated microfluidic devices can mimic some of the key features of soil with highly controlled physical and chemical microenvironments at the scale of pores, aggregates, and microbes. The combination of various techniques makes microfluidics an integrated approach for observation, reaction, analysis, and characterization. In this review, we systematically summarize the emerging applications of microfluidic soil platforms, from investigating soil interfacial processes and soil microbial processes to soil analysis and high-throughput screening. We highlight how innovative microfluidic devices are used to provide new insights into soil processes, mechanisms, and effects at the microscale, which contribute to an integrated interrogation of the soil systems across different scales. Critical discussions of the practical limitations of microfluidic soil platforms and perspectives of future research directions are summarized. We envisage that microfluidics will represent the technological advances toward microscopic, controllable, and in situ soil research.

Deciphering Microbial Metal Toxicity Responses via Random Bar Code Transposon Site Sequencing and Activity-Based Metabolomics

To uncover metal toxicity targets and defense mechanisms of the facultative anaerobe Pantoea sp. strain MT58 (MT58), we used a multiomic strategy combining two global techniques, random bar code transposon site sequencing (RB-TnSeq) and activity-based metabolomics. MT58 is a metal-tolerant Oak Ridge Reservation (ORR) environmental isolate that was enriched in the presence of metals at concentrations measured in contaminated groundwater at an ORR nuclear waste site. The effects of three chemically different metals found at elevated concentrations in the ORR contaminated environment were investigated: the cation Al³⁺, the oxyanion CrO₄^2-, and the oxycation UO₂²⁺. Both global techniques were applied using all three metals under both aerobic and anaerobic conditions to elucidate metal interactions mediated through the activity of metabolites and key genes/proteins. These revealed that Al³⁺ binds intracellular arginine, CrO₄^2- enters the cell through sulfate transporters and oxidizes intracellular reduced thiols, and membrane-bound lipopolysaccharides protect the cell from UO₂²⁺ toxicity. In addition, the Tol outer membrane system contributed to the protection of cellular integrity from the toxic effects of all three metals. Likewise, we found evidence of regulation of lipid content in membranes under metal stress. Individually, RB-TnSeq and metabolomics are powerful tools to explore the impact various stresses have on biological systems. Here, we show that together they can be used synergistically to identify the molecular actors and mechanisms of these pertubations to an organism, furthering our understanding of how living systems interact with their environment. IMPORTANCE Studying microbial interactions with their environment can lead to a deeper understanding of biological molecular mechanisms. In this study, two global techniques, RB-TnSeq and activity metabolomics, were successfully used to probe the interactions between a metal-resistant microorganism, Pantoea sp. strain MT58, and metals contaminating a site where the organism can be located. A number of novel metal-microbe interactions were uncovered, including Al³⁺ toxicity targeting arginine synthesis, which could lead to a deeper understanding of the impact Al³⁺ contamination has on microbial communities as well as its impact on higher-level organisms, including plants for whom Al³⁺ contamination is an issue. Using multiomic approaches like the one described here is a way to further our understanding of microbial interactions and their impacts on the environment overall.