Microscale advection governs microbial growth and oxygen consumption in macroporous aggregates.

Most microbial life on Earth is found in localized microenvironments that collectively exert a crucial role in maintaining ecosystem health and influencing global biogeochemical cycles. In many habitats such as biofilms in aquatic systems, bacterial flocs in activated sludge, periphyton mats, or particles sinking in the ocean, these microenvironments experience sporadic or continuous flow. Depending on their microscale structure, pores and channels through the microenvironments permit localized flow that shifts the relative importance of diffusive and advective mass transport. How this flow alters nutrient supply, facilitates waste removal, drives the emergence of different microbial niches, and impacts the overall function of the microenvironments remains unclear. Here, we quantify how pores through microenvironments that permit flow can elevate nutrient supply to the resident bacterial community using a microfluidic experimental system and gain further insights from coupled population-based and computational fluid dynamics simulations. We find that the microscale structure determines the relative contribution of advection vs diffusion, and even a modest flow through a pore in the range of 10 µm s-1 can increase the carrying capacity of a microenvironment by 10%. Recognizing the fundamental role that microbial hotspots play in the Earth system, developing frameworks that predict how their heterogeneous morphology and potential interstitial flows change microbial function and collectively alter global scale fluxes is critical.IMPORTANCEMicrobial life is a key driver of global biogeochemical cycles. Similar to the distribution of humans on Earth, they are often not homogeneously distributed in nature but occur in dense clusters that resemble microbial cities. Within and around these clusters, diffusion is often assumed as the sole mass-transfer process that dictates nutrient supply and waste removal. In many natural and engineered systems such as biofilms in aquatic environments, aggregates in bioremediation, or flocs in wastewater treatment plants, these clusters are exposed to flow that elevates mass transfer, a process that is often overlooked. In this study, we show that advective fluxes can increase the local growth of bacteria in a single microenvironment by up to 50% and shape their metabolism by disrupting localized anoxia or supplying nutrients at different rates. Collectively, advection-enhanced mass transport may thus regulate important biogeochemical transformations in both natural and engineered environments.

Uncultivated DPANN archaea are ubiquitous inhabitants of global oxygen-deficient zones with diverse metabolic potential.

Archaea belonging to the DPANN (Diapherotrites, Parvarchaeota, Aenigmarchaeota, Nanoarchaeota, and Nanohaloarchaeota) superphylum have been found in an expanding number of environments and perform a variety of biogeochemical roles, including contributing to carbon, sulfur, and nitrogen cycling. Generally characterized by ultrasmall cell sizes and reduced genomes, DPANN archaea may form mutualistic, commensal, or parasitic interactions with various archaeal and bacterial hosts, influencing the ecology and functioning of microbial communities. While DPANN archaea reportedly comprise a sizeable fraction of the archaeal community within marine oxygen-deficient zone (ODZ) water columns, little is known about their metabolic capabilities in these ecosystems. We report 33 novel metagenome-assembled genomes (MAGs) belonging to the DPANN phyla Nanoarchaeota, Pacearchaeota, Woesearchaeota, Undinarchaeota, Iainarchaeota, and SpSt-1190 from pelagic ODZs in the Eastern Tropical North Pacific and the Arabian Sea. We find these archaea to be permanent, stable residents of all three major ODZs only within anoxic depths, comprising up to 1% of the total microbial community and up to 25%-50% of archaea as estimated from read mapping to MAGs. ODZ DPANN appear to be capable of diverse metabolic functions, including fermentation, organic carbon scavenging, and the cycling of sulfur, hydrogen, and methane. Within a majority of ODZ DPANN, we identify a gene homologous to nitrous oxide reductase. Modeling analyses indicate the feasibility of a nitrous oxide reduction metabolism for host-attached symbionts, and the small genome sizes and reduced metabolic capabilities of most DPANN MAGs suggest host-associated lifestyles within ODZs. IMPORTANCE: Archaea from the DPANN (Diapherotrites, Parvarchaeota, Aenigmarchaeota, Nanoarchaeota, and Nanohaloarchaeota) superphylum have diverse metabolic capabilities and participate in multiple biogeochemical cycles. While metagenomics and enrichments have revealed that many DPANN are characterized by ultrasmall genomes, few biosynthetic genes, and episymbiotic lifestyles, much remains unknown about their biology. We report 33 new DPANN metagenome-assembled genomes originating from the three global marine oxygen-deficient zones (ODZs), the first from these regions. We survey DPANN abundance and distribution within the ODZ water column, investigate their biosynthetic capabilities, and report potential roles in the cycling of organic carbon, methane, and nitrogen. We test the hypothesis that nitrous oxide reductases found within several ODZ DPANN genomes may enable ultrasmall episymbionts to serve as nitrous oxide consumers when attached to a host nitrous oxide producer. Our results indicate DPANN archaea as ubiquitous residents within the anoxic core of ODZs with the potential to produce or consume key compounds.

Uncultivated DPANN archaea are ubiquitous inhabitants of global oxygen deficient zones with diverse metabolic potential.

Archaea belonging to the DPANN superphylum have been found within an expanding number of environments and perform a variety of biogeochemical roles, including contributing to carbon, sulfur, and nitrogen cycling. Generally characterized by ultrasmall cell sizes and reduced genomes, DPANN archaea may form mutualistic, commensal, or parasitic interactions with various archaeal and bacterial hosts, influencing the ecology and functioning of microbial communities. While DPANN archaea reportedly comprise 15-26% of the archaeal community within marine oxygen deficient zone (ODZ) water columns, little is known about their metabolic capabilities in these ecosystems. We report 33 novel metagenome-assembled genomes belonging to DPANN phyla Nanoarchaeota, Pacearchaeota, Woesarchaeota, Undinarchaeota, Iainarchaeota, and SpSt-1190 from pelagic ODZs in the Eastern Tropical North Pacific and Arabian Sea. We find these archaea to be permanent, stable residents of all 3 major ODZs only within anoxic depths, comprising up to 1% of the total microbial community and up to 25-50% of archaea. ODZ DPANN appear capable of diverse metabolic functions, including fermentation, organic carbon scavenging, and the cycling of sulfur, hydrogen, and methane. Within a majority of ODZ DPANN, we identify a gene homologous to nitrous oxide reductase. Modeling analyses indicate the feasibility of a nitrous oxide reduction metabolism for host-attached symbionts, and the small genome sizes and reduced metabolic capabilities of most DPANN MAGs suggest host-associated lifestyles within ODZs.

Evaporation-induced hydrodynamics control plasmid transfer during surface-associated microbial growth.

Droplet evaporation is a general process in unsaturated environments that results in micro-scale hydrodynamic flows which in turn determine the spatial distributions of microbial cells across surfaces. These spatial distributions can have significant effects on the development and functioning of surface-associated microbial communities, with consequences for important processes such as the spread of plasmids. Here, we experimentally quantified how evaporation-induced hydrodynamic processes modulate the initial deposition patterns of microbial cells (via the coffee ring effect and Marangoni convection) and how these patterns control the spread of an antibiotic resistance-encoding plasmid during surface-associated growth. We found that plasmid spread is a function of the initial density of cells deposited along the droplet periphery, which is a manifestation of the coffee ring effect. Using an individual-based model, we systematically linked how the different initial cell deposition patterns caused by the relative strengths of the coffee ring effect and Marangoni convection determine the extent of plasmid transfer during surface-associated growth. Our study demonstrates that evaporation-induced hydrodynamic processes that are common in nature can alter crucial ecological properties of surface-associated microbial communities and control the proliferation of plasmids, with consequences on the spread of antibiotic resistance and other plasmid-encoded traits.

Porous marine snow differentially benefits chemotactic, motile, and nonmotile bacteria.

Particulate organic carbon settling through the marine water column is a key process that regulates the global climate by sequestering atmospheric carbon. The initial colonization of marine particles by heterotrophic bacteria represents the first step in recycling this carbon back to inorganic constituents-setting the magnitude of vertical carbon transport to the abyss. Here, we demonstrate experimentally using millifluidic devices that, although bacterial motility is essential for effective colonization of a particle leaking organic nutrients into the water column, chemotaxis specifically benefits at intermediate and higher settling velocities to navigate the particle boundary layer during the brief window of opportunity provided by a passing particle. We develop an individual-based model that simulates the encounter and attachment of bacterial cells with leaking marine particles to systematically evaluate the role of different parameters associated with bacterial run-and-tumble motility. We further use this model to explore the role of particle microstructure on the colonization efficiency of bacteria with different motility traits. We find that the porous microstructure facilitates additional colonization by chemotactic and motile bacteria, and fundamentally alters the way nonmotile cells interact with particles due to streamlines intersecting with the particle surface.

The media composition as a crucial element in high-throughput metabolic network reconstruction.

In recent years, metagenome-assembled genomes (MAGs) have provided glimpses into the intra- and interspecies genetic diversity and interactions that form the bases of complex microbial communities. High-throughput reconstruction of genome-scale metabolic networks (GEMs) from MAGs is a promising avenue to disentangle the myriad trophic interactions stabilizing these communities. However, high-throughput reconstruction of GEMs relies on accurate gap filling of metabolic pathways using automated algorithms. Here, we systematically explore how the composition of the media (specification of the available nutrients and metabolites) during gap filling influences the resulting GEMs concerning predicted auxotrophies for fully sequenced model organisms and environmental isolates. We expand this analysis by using 106 MAGs from the same species with differing quality. We find that although the completeness of MAGs influences the fraction of gap-filled reactions, the composition of the media plays the dominant role in the accurate prediction of auxotrophies that form the basis of myriad community interactions. We propose that constraining the media composition for gap filling through both experimental approaches and computational approaches will increase the reliability of high-throughput reconstruction of genome-scale metabolic models from MAGs and paves the way for culture independent prediction of trophic interactions in complex microbial communities.

Bacterial age distribution in soil - Generational gaps in adjacent hot and cold spots.

Resource patchiness and aqueous phase fragmentation in soil may induce large differences local growth conditions at submillimeter scales. These are translated to vast differences in bacterial age from cells dividing every thirty minutes in close proximity to plant roots to very old cells experiencing negligible growth in adjacent nutrient poor patches. In this study, we link bacterial population demographics with localized soil and hydration conditions to predict emerging generation time distributions and estimate mean bacterial cell ages using mechanistic and heuristic models of bacterial life in soil. Results show heavy-tailed distributions of generation times that resemble a power law for certain conditions, suggesting that we may find bacterial cells of vastly different ages living side by side within small soil volumes. Our results imply that individual bacteria may exist concurrently with all of their ancestors, resulting in an archive of bacterial cells with traits that have been gained (and lost) throughout time-a feature unique to microbial life. This reservoir of bacterial strains and the potential for the reemergence of rare strains with specific functions may be critical for ecosystem stability and function.

Rare and localized events stabilize microbial community composition and patterns of spatial self-organization in a fluctuating environment.

Spatial self-organization is a hallmark of surface-associated microbial communities that is governed by local environmental conditions and further modified by interspecific interactions. Here, we hypothesize that spatial patterns of microbial cell-types can stabilize the composition of cross-feeding microbial communities under fluctuating environmental conditions. We tested this hypothesis by studying the growth and spatial self-organization of microbial co-cultures consisting of two metabolically interacting strains of the bacterium Pseudomonas stutzeri. We inoculated the co-cultures onto agar surfaces and allowed them to expand (i.e. range expansion) while fluctuating environmental conditions that alter the dependency between the two strains. We alternated between anoxic conditions that induce a mutualistic interaction and oxic conditions that induce a competitive interaction. We observed co-occurrence of both strains in rare and highly localized clusters (referred to as "spatial jackpot events") that persist during environmental fluctuations. To resolve the underlying mechanisms for the emergence of spatial jackpot events, we used a mechanistic agent-based mathematical model that resolves growth and dispersal at the scale relevant to individual cells. While co-culture composition varied with the strength of the mutualistic interaction and across environmental fluctuations, the model provides insights into the formation of spatially resolved substrate landscapes with localized niches that support the co-occurrence of the two strains and secure co-culture function. This study highlights that in addition to spatial patterns that emerge in response to environmental fluctuations, localized spatial jackpot events ensure persistence of strains across dynamic conditions.

Spatiotemporal metabolic modeling of bacterial life in complex habitats.

The combination of genome-scale metabolic networks with spatially explicit representation of microbial habitats (spatiotemporal metabolic network modeling) paves the way to predict complex metabolic landscapes to a hitherto unparalleled detail, thus providing new insights into trophic interactions occurring at different scales. Placing detailed bacterial metabolism in realistic physical environment highlights the roles of physical barriers and diffusional bottlenecks on bacterial community interactions, structure and stability. We review recent advances in spatiotemporal metabolic network modeling using a few illustrative examples that highlight the immense potential of these novel approaches to interpret and design metabolic mediated interactions in structures (natural and engineered) environments.

Spatial organization in microbial range expansion emerges from trophic dependencies and successful lineages.

Evidence suggests that bacterial community spatial organization affects their ecological function, yet details of the mechanisms that promote spatial patterns remain difficult to resolve experimentally. In contrast to bacterial communities in liquid cultures, surface-attached range expansion fosters genetic segregation of the growing population with preferential access to nutrients and reduced mechanical restrictions for cells at the expanding periphery. Here we elucidate how localized conditions in cross-feeding bacterial communities shape community spatial organization. We combine experiments with an individual based mathematical model to resolve how trophic dependencies affect localized growth rates and nucleate successful cell lineages. The model tracks individual cell lineages and attributes these with trophic dependencies that promote counterintuitive reproductive advantages and result in lasting influences on the community structure, and potentially, on its functioning. We examine persistence of lucky lineages in structured habitats where expansion is interrupted by physical obstacles to gain insights into patterns in porous domains.

Reduced gravity promotes bacterially mediated anoxic hotspots in unsaturated porous media.

Human endeavours into deep space exploration and the prospects of establishing colonies on nearby planets would invariably involve components of bioregenerative life support for food production, cabin atmosphere renewal, and waste recycling. Growing plants and their microbiomes in porous media under different gravitational fields may present new challenges due to effects of liquid distribution on gaseous exchange with roots and microorganisms. We provide the first direct evidence that capillary driven liquid reconfiguration in porous media under reduced gravity conditions reduces oxygen diffusion pathways and enhances anoxic conditions within bacterial hotspots. Parabolic flight experiments using model porous media inoculated with aerobic and facultative anaerobic bacteria reveal the systematic enhancement of anoxic conditions during the reduced gravity periods in the presence but not in the absence of bacterial activity. The promotion of anoxic conditions under reduced gravity may lead to higher nitrous oxide and methane emissions relative to Earth conditions, on the other hand, anoxic conditions could be beneficial for perchlorate bioremediation of Martian soil. The results highlight changes in soil bacterial microhabitats under reduced gravity and the challenges of managing bioregenerative life support systems in space.

The engineering of spatially linked microbial consortia - potential and perspectives.

Traditional biotechnological applications of microorganisms employ mono-cultivation or co-cultivation in well-mixed vessels disregarding the potential of spatially organized cultures. Metabolic specialization and guided species interactions facilitated through spatial isolation would enable consortia of microbes to accomplish more complex functions than currently possible, for bioproduction as well as biodegradation processes. Here, we review concepts of spatially linked microbial consortia in which spatial arrangement is optimized to increase control and facilitate new species combinations. We highlight that genome-scale metabolic network models can inform the design and tuning of synthetic microbial consortia and suggest that a standardized assembly of such systems allows the combination of 'incompatibles', potentially leading to countless novel applications.

Modeling metabolic networks of individual bacterial agents in heterogeneous and dynamic soil habitats (IndiMeSH).

Natural soil is characterized as a complex habitat with patchy hydrated islands and spatially variable nutrients that is in a constant state of change due to wetting-drying dynamics. Soil microbial activity is often concentrated in sparsely distributed hotspots that contribute disproportionally to macroscopic biogeochemical nutrient cycling and greenhouse gas emissions. The mechanistic representation of such dynamic hotspots requires new modeling approaches capable of representing the interplay between dynamic local conditions and the versatile microbial metabolic adaptations. We have developed IndiMeSH (Individual-based Metabolic network model for Soil Habitats) as a spatially explicit model for the physical and chemical microenvironments of soil, combined with an individual-based representation of bacterial motility and growth using adaptive metabolic networks. The model uses angular pore networks and a physically based description of the aqueous phase as a backbone for nutrient diffusion and bacterial dispersal combined with dynamic flux balance analysis to calculate growth rates depending on local nutrient conditions. To maximize computational efficiency, reduced scale metabolic networks are used for the simulation scenarios and evaluated strategically to the genome scale model. IndiMeSH was compared to a well-established population-based spatiotemporal metabolic network model (COMETS) and to experimental data of bacterial spatial organization in pore networks mimicking soil aggregates. IndiMeSH was then used to strategically study dynamic response of a bacterial community to abrupt environmental perturbations and the influence of habitat geometry and hydration conditions. Results illustrate that IndiMeSH is capable of representing trophic interactions among bacterial species, predicting the spatial organization and segregation of bacterial populations due to oxygen and carbon gradients, and provides insights into dynamic community responses as a consequence of environmental changes. The modular design of IndiMeSH and its implementation are adaptable, allowing it to represent a wide variety of experimental and in silico microbial systems.

Spatial organization of bacterial populations in response to oxygen and carbon counter-gradients in pore networks.

Microbial activity in soil is spatially heterogeneous often forming spatial hotspots that contribute disproportionally to biogeochemical processes. Evidence suggests that bacterial spatial organization contributes to the persistence of anoxic hotspots even in unsaturated soils. Such processes are difficult to observe in situ at the microscale, hence mechanisms and time scales relevant for bacterial spatial organization remain largely qualitative. Here we develop an experimental platform based on glass-etched micrometric pore networks that mimics resource gradients postulated in soil aggregates to observe spatial organization of fluorescently tagged aerobic and facultative anaerobic bacteria. Two initially intermixed bacterial species, Pseudomonas putida and Pseudomonas veronii, segregate into preferential regions promoted by opposing gradients of carbon and oxygen (such persistent coexistence is not possible in well-mixed cultures). The study provides quantitative visualization and modeling of bacterial spatial organization within aggregate-like hotspots, a key step towards developing a mechanistic representation of bacterial community organization in soil pores.