Competition between growth and shear stress drives intermittency in preferential flow paths in porous medium biofilms

Yielding behaviour of chemically treated Pseudomonas fluorescens biofilms

The mechanics of biofilms are intrinsically shaped by their physicochemical environment. By understanding the influence of the extracellular matrix composition, pH and elevated levels of cationic species on the biofilm rheology, novel living materials with tuned properties can be formulated. In this study, we examine the role of a chaotropic agent (urea), two divalent cations and distilled deionized water on the nonlinear viscoelasticity of a model biofilm Pseudomonas fluorescens. The structural breakdown of each biofilm is quantified using tools of non-linear rheology. Our findings reveal that urea induced a softening response, and displayed strain overshoots comparable to distilled deionized water, without altering the microstructural packing fraction and macroscale morphology. The absorption of divalent ferrous and calcium cations into the biofilm matrix resulted in stiffening and a reduction in normalized elastic energy dissipation, accompanied by macroscale morphological wrinkling and moderate increases in the packing fraction. Notably, ferrous ions induced a predominance of rate dependent yielding, whereas the calcium ions resulted in equal contribution from both rate and strain dependent yielding and structural breakdown of the biofilms. Together, these results indicate that strain rate increasingly becomes an important factor controlling biofilm fluidity with cation-induced biofilm stiffening. The finding can help inform effective biofilm removal protocols and in development of bio-inks for additive manufacturing of biofilm derived materials.

Soft matter physics of the ground beneath our feet

The soft part of the Earth's surface - the ground beneath our feet - constitutes the basis for life and natural resources, yet a general physical understanding of the ground is still lacking. In this critical time of climate change, cross-pollination of scientific approaches is urgently needed to better understand the behavior of our planet's surface. The major topics in current research in this area cross different disciplines, spanning geosciences, and various aspects of engineering, material sciences, physics, chemistry, and biology. Among these, soft matter physics has emerged as a fundamental nexus connecting and underpinning many research questions. This perspective article is a multi-voice effort to bring together different views and approaches, questions and insights, from researchers that work in this emerging area, the soft matter physics of the ground beneath our feet. In particular, we identify four major challenges concerned with the dynamics in and of the ground: (I) modeling from the grain scale, (II) near-criticality, (III) bridging scales, and (IV) life. For each challenge, we present a selection of topics by individual authors, providing specific context, recent advances, and open questions. Through this, we seek to provide an overview of the opportunities for the broad Soft Matter community to contribute to the fundamental understanding of the physics of the ground, strive towards a common language, and encourage new collaborations across the broad spectrum of scientists interested in the matter of the Earth's surface.

In Vivo Microrheology Reveals Local Elastic and Plastic Responses Inside 3D Bacterial Biofilms

Bacterial biofilms are highly abundant 3D living materials capable of performing complex biomechanical and biochemical functions, including programmable growth, self-repair, filtration, and bioproduction. Methods to measure internal mechanical properties of biofilms in vivo with spatial resolution on the cellular scale have been lacking. Here, thousands of cells are tracked inside living 3D biofilms of the bacterium Vibrio cholerae during and after the application of shear stress, for a wide range of stress amplitudes, periods, and biofilm sizes, which revealed anisotropic elastic and plastic responses of both cell displacements and cell reorientations. Using cellular tracking to infer parameters of a general mechanical model, spatially-resolved measurements of the elastic modulus inside the biofilm are obtained, which correlate with the spatial distribution of the polysaccharides within the biofilm matrix. The noninvasive microrheology and force-inference approach introduced here provides a general framework for studying mechanical properties with high spatial resolution in living materials.

The mechanism of biofilm detachment in porous medium under flow field

Biofilms are communities formed by bacteria adhering to surfaces, widely present in porous medium, and their growth can lead to clogging. Our experiment finds that under certain flow conditions, biofilms detach in pores and form a dynamically changing flow path. We define detachment that occurs far from the boundary of the flow path (with a distance greater than 200 μm) as internal detachment and detachment that occurs at the boundary of the flow path as external detachment. To understand the mechanism of biofilm detachment, we study the detachment behaviors of the Bacillus subtilis biofilm in a porous medium in a microfluidic device, where Bacillus subtilis strain is triple fluorescent labeled, which can represent three main phenotypes during the biofilm formation: motile cells, matrix-producing cells, and spores. We find that slow small-scale internal detachment occurs in regions with very few motile cells and matrix-producing cells, and bacterial movement in these areas is disordered. The increase in the number of matrix-producing cells induces clogging, and after clogging, the rapid detachment of the bulk internal biofilm occurs due to the increased pressure difference at the inlet and outlet. When both internal and external detachments occur simultaneously, the number of matrix-producing cells in the internal detachment area is 2.5 times that in the external detachment area. The results indicate that biofilm detachment occurs in areas with fewer matrix-producing cells, as matrix-producing cells can help resist detachment by secreting extracellular polymeric substances (EPSs).

Agent-based modeling of stress anisotropy driven nematic ordering in growing biofilms

Living active collectives have evolved with remarkable self-patterning capabilities to adapt to the physical and biological constraints crucial for their growth and survival. However, the intricate process by which complex multicellular patterns emerge from a single founder cell remains elusive. In this study, we utilize an agent-based model, validated through single-cell microscopy imaging, to track the three-dimensional (3D) morphodynamics of cells within growing bacterial biofilms encased by agarose gels. The confined growth conditions give rise to a spatiotemporally heterogeneous stress landscape within the biofilm. In the core of the biofilm, where high hydrostatic and low shear stresses prevail, cell packing appears disordered. In contrast, near the gel-cell interface, a state of high shear stress and low hydrostatic stress emerges, driving nematic ordering, albeit with a time delay inherent to shear stress relaxation. Strikingly, we observe a robust spatiotemporal correlation between stress anisotropy and nematic ordering within these confined biofilms. This correlation suggests a mechanism whereby stress anisotropy plays a pivotal role in governing the spatial organization of cells. The reciprocity between stress anisotropy and cell ordering in confined biofilms opens new avenues for innovative 3D mechanically guided patterning techniques for living active collectives, which hold significant promise for a wide array of environmental and biomedical applications.

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.

Hemodynamic regulation allows stable growth of microvascular networks

How do vessels find optimal radii? Capillaries are known to adapt their radii to maintain the shear stress of blood flow at the vessel wall at a set point, yet models of adaptation purely based on average shear stress have not been able to produce complex loopy networks that resemble real microvascular systems. For narrow vessels where red blood cells travel in a single file, the shear stress on vessel endothelium peaks sharply when a red blood cell passes through. We show that stable shear-stress-based adaptation is possible if vessel shear stress set points are cued to the stress peaks. Model networks that respond to peak stresses alone can quantitatively reproduce the observed zebrafish trunk microcirculation, including its adaptive trajectory when hematocrit changes or parts of the network are amputated. Our work reveals the potential for mechanotransduction alone to generate stable hydraulically tuned microvascular networks.

Microfluidic approaches in microbial ecology

Microbial life is at the heart of many diverse environments and regulates most natural processes, from the functioning of animal organs to the cycling of global carbon. Yet, the study of microbial ecology is often limited by challenges in visualizing microbial processes and replicating the environmental conditions under which they unfold. Microfluidics operates at the characteristic scale at which microorganisms live and perform their functions, thus allowing for the observation and quantification of behaviors such as growth, motility, and responses to external cues, often with greater detail than classical techniques. By enabling a high degree of control in space and time of environmental conditions such as nutrient gradients, pH levels, and fluid flow patterns, microfluidics further provides the opportunity to study microbial processes in conditions that mimic the natural settings harboring microbial life. In this review, we describe how recent applications of microfluidic systems to microbial ecology have enriched our understanding of microbial life and microbial communities. We highlight discoveries enabled by microfluidic approaches ranging from single-cell behaviors to the functioning of multi-cellular communities, and we indicate potential future opportunities to use microfluidics to further advance our understanding of microbial processes and their implications.

A versatile micromodel technology to explore biofilm development in porous media flows

Bacterial biofilms that grow in porous media are critical to ecosystem processes and applications ranging from soil bioremediation to bioreactors for treating wastewater or producing value-added products. However, understanding and engineering the complex phenomena that drive the development of biofilms in such systems remains a challenge. Here we present a novel micromodel technology to explore bacterial biofilm development in porous media flows. The technology consists of a set of modules that can be combined as required for any given experiment and conveniently tuned for specific requirements. The core module is a 3D-printed micromodel where biofilm is grown into a perfusable porous substrate. High-precision additive manufacturing, in particular stereolithography, is used to fabricate porous scaffolds with precisely controlled architectures integrating flow channels with diameters down to several hundreds of micrometers. The system is instrumented with: ultraviolet-C light-emitting diodes; on-line measurements of oxygen consumption and pressure drop across the porous medium; camera and spectrophotometric cells for the detection of biofilm detachment events at the outlet. We demonstrate how this technology can be used to study the development of Pseudomonas aeruginosa biofilm for several days within a network of flow channels. We find complex dynamics whereby oxygen consumption reaches a steady-state but not the pressure drop, which instead features a permanent regime with large fluctuations. We further use X-ray computed microtomography to image the spatial distribution of biofilms and computational fluid dynamics to link biofilm development with local flow properties. By combining the advantages of additive manufacturing for the creation of reproducible 3D porous microarchitectures with the flow control and instrumentation accuracy of microfluidics, our system provides a platform to study the dynamics of biofilm development in 3D porous media and to rapidly test new concepts in process engineering.

A pore-scale reactive transport modeling study for quorum sensing-driven biofilm dispersal in heterogeneous porous media

Microorganisms regulate the expression of energetically expensive phenotypes via a collective decision-making mechanism known as quorum sensing (QS). This study investigates the intricate dynamics of biofilm growth and QS-controlled biofilm dispersal in heterogeneous porous media, employing a pore-scale reactive transport modeling approach. Model simulations carried out under various fluid flow conditions and biofilm growth scenarios reveal that QS processes are influenced not only by the biomass density of biofilm colonies but also by a complex interplay between pore architecture, flow velocity, and the rates of biofilm growth and dispersal. This study demonstrates that pore architecture controls the initiation of QS processes and advection gives rise to oscillatory growth of biofilms. Such oscillation is suppressed if biofilm dynamics are in favor of sustaining a sufficiently high signal concentration, such as fast growth or slow dispersal rates. By establishing a mathematical framework, this study contributes to the fundamental understanding of QS-controlled biofilm dynamics in complex environments.

Controlling biofilm transport with porous metamaterials designed with Bayesian learning

Biofilm growth and transport in confined systems frequently occur in natural and engineered systems. Designing customizable engineered porous materials for controllable biofilm transportation properties could significantly improve the rapid utilization of biofilms as engineered living materials for applications in pollution alleviation, material self-healing, energy production, and many more. We combine Bayesian optimization (BO) and individual-based modeling to conduct design optimizations for maximizing different porous materials' (PM) biofilm transportation capability. We first characterize the acquisition function in BO for designing 2-dimensional porous membranes. We use the expected improvement acquisition function for designing lattice metamaterials (LM) and 3-dimensional porous media (3DPM). We find that BO is 92.89% more efficient than the uniform grid search method for LM and 223.04% more efficient for 3DPM. For all three types of structures, the selected characterization simulation tests are in good agreement with the design spaces approximated with Gaussian process regression. All the extracted optimal designs exhibit better biofilm growth and transportability than unconfined space without substrates. Our comparison study shows that PM stimulates biofilm growth by taking up volumetric space and pushing biofilms' upward growth, as evidenced by a 20% increase in bacteria cell numbers in unconfined space compared to porous materials, and 128% more bacteria cells in the target growth region for PM-induced biofilm growth compared with unconfined growth. Our work provides deeper insights into the design of substrates to tune biofilm growth, analyzing the optimization process and characterizing the design space, and understanding biophysical mechanisms governing the growth of biofilms.

Microfluidic investigation of the impacts of flow fluctuations on the development of Pseudomonas putida biofilms

Biofilms play critical roles in wastewater treatment, bioremediation, and medical-device-related infections. Understanding the dynamics of biofilm formation and growth is essential for controlling and exploiting their properties. However, the majority of current studies have focused on the impact of steady flows on biofilm growth, while flow fluctuations are common in natural and engineered systems such as water pipes and blood vessels. Here, we reveal the effects of flow fluctuations on the development of Pseudomonas putida biofilms through systematic microfluidic experiments and the development of a theoretical model. Our experimental results showed that biofilm growth under fluctuating flow conditions followed three phases: lag, exponential, and fluctuation phases. In contrast, biofilm growth under steady-flow conditions followed four phases: lag, exponential, stationary, and decline phases. Furthermore, we demonstrated that low-frequency flow fluctuations promoted biofilm growth, while high-frequency fluctuations inhibited its development. We attributed the contradictory impacts of flow fluctuations on biofilm growth to the adjustment time (T0) needed for biofilm to grow after the shear stress changed from high to low. Furthermore, we developed a theoretical model that explains the observed biofilm growth under fluctuating flow conditions. Our insights into the mechanisms underlying biofilm development under fluctuating flows can inform the design of strategies to control biofilm formation in diverse natural and engineered systems.

A novel model for biofilm initiation in porous media flow

Bacteria often form biofilms in porous environments where an external flow is present, such as soil or filtration systems. To understand the initial stages of biofilm formation, one needs to study the interactions between cells, the fluid and the confining geometries. Here, we present an agent based numerical model for bacteria that takes into account the planktonic stage of motile cells as well as surface attachment and biofilm growth in a lattice Boltzmann fluid. In the planktonic stage we show the importance of the interplay between complex flow and cell motility when determining positions of surface attachment. In the growth stage we show the applicability of our model by investigating how external flow and biofilm stiffness determine qualitative colony morphologies as well as quantitative measurements of, e.g., permeability.

Internal Biofilm Heterogeneities Enhance Solute Mixing and Chemical Reactions in Porous Media

Bacterial biofilms can form in porous media that are of interest in industrial applications ranging from medical implants to biofilters as well as in environmental applications such as in situ groundwater remediation, where they can be critical locations for biogeochemical reactions. The presence of biofilms modifies porous media topology and hydrodynamics by clogging pores and consequently solutes transport and reactions kinetics. The interplay between highly heterogeneous flow fields found in porous media and microbial behavior, including biofilm growth, results in a spatially heterogeneous biofilm distribution in the porous media as well as internal heterogeneity across the thickness of the biofilm. Our study leverages highly resolved three-dimensional X-ray computed microtomography images of bacterial biofilms in a tubular reactor to numerically compute pore-scale fluid flow and solute transport by considering multiple equivalent stochastically generated internal permeability fields for the biofilm. We show that the internal heterogeneous permeability mainly impacts intermediate velocities when compared with homogeneous biofilm permeability. While the equivalent internal permeability fields of the biofilm do not impact fluid-fluid mixing, they significantly control a fast reaction. For biologically driven reactions such as nutrient or contaminant uptake by the biofilm, its internal permeability field controls the efficiency of the process. This study highlights the importance of considering the internal heterogeneity of biofilms to better predict reactivity in industrial and environmental bioclogged porous systems.

Morphogenesis of Biofilms in Porous Media and Control on Hydrodynamics

The functioning of natural and engineered porous media, like soils and filters, depends in many cases on the interplay between biochemical processes and hydrodynamics. In such complex environments, microorganisms often form surface-attached communities known as biofilms. Biofilms can take the shape of clusters, which alter the distribution of fluid flow velocities within the porous medium, subsequently influencing biofilm growth. Despite numerous experimental and numerical efforts, the control of the biofilm clustering process and the resulting heterogeneity in biofilm permeability is not well understood, limiting our predictive abilities for biofilm-porous medium systems. Here, we use a quasi-2D experimental model of a porous medium to characterize biofilm growth dynamics for different pore sizes and flow rates. We present a method to obtain the time-resolved biofilm permeability field from experimental images and use the obtained permeability field to compute the flow field through a numerical model. We observe a biofilm cluster size distribution characterized by a spectrum slope evolving in time between -2 and -1, a fundamental measure that can be used to create spatio-temporal distributions of biofilm clusters for upscaled models. We find a previously undescribed biofilm permeability distribution, which can be used to stochastically generate permeability fields within biofilms. An increase in velocity variance for a decrease in physical heterogeneity shows that the bioclogged porous medium behaves differently than expected from studies on heterogeneity in abiotic porous media.

Rapid formation of bioaggregates and morphology transition to biofilm streamers induced by pore-throat flows

Bioaggregates are condensed porous materials comprising microbes, organic and inorganic matters, and water. They are commonly found in natural and engineered porous media and often cause clogging. Despite their importance, the formation mechanism of bioaggregates in porous media systems is largely unknown. Through microfluidic experiments and direct numerical simulations of fluid flow, we show that the rapid bioaggregation is driven by the interplay of the viscoelastic nature of biomass and hydrodynamic conditions at pore throats. At an early stage, unique flow structures around a pore throat promote the biomass attachment at the throat. Then, the attached biomass fluidizes when the shear stress at the partially clogged pore throat reaches a critical value. After the fluidization, the biomass is displaced and accumulated in the expansion region of throats forming bioaggregates. We further find that such criticality in shear stress triggers morphological changes in bioaggregates from rounded- to streamer-like shapes. This knowledge was used to control the clogging of throats by tuning the flow conditions: When the shear stress at the throat exceeded the critical value, clogging was prevented. The bioaggregation process did not depend on the detailed pore-throat geometry, as we reproduced the same dynamics in various pore-throat geometries. This study demonstrates that pore-throat structures, which are ubiquitous in porous media systems, induce bioaggregation and can lead to abrupt disruptions in flow.