Biofilm development and the dynamics of preferential flow paths in porous media

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).

Identification and extraction of cementation patterns in sand modified by MICP: New insights at the pore scale

Microbially induced calcium carbonate precipitation (MICP) is an environmentally friendly technology that improves soil permeability resistance through biocementation. In this study, 2D microscopic analysis and 3D volume reconstruction were performed on river sand after 24 cycles of bio-treatment based on stacked images and computed tomography (CT) scanning data, respectively, to extract biocementation patterns between particles. Based on the mutual validation findings of the two techniques, three patterns in the biocemented sand were identified as G-C-G, G-C, and G-G. Specifically, 2D microscopic analysis showed that G-C-G featured multi-particle encapsulation and bridging, with a pore filling ratio of 81.2%; G-C was characterized by locally coated particle layers, with a pore filling ratio of 19.7%; and the G-G was marked by sporadic filling of interparticle pores, with a pore filling ratio of 11.7%. G-C-G had the best cementation effect and permeability resistance (effective sealing rate of 68.5%), whereas G-C (effective sealing rate of 2.4%) had a relatively minor contribution to pore-filling and flow sealing. 3D volume reconstruction showed that G-C-G had the highest pore filling rate, followed by G-G and G-C. The average filling ratios of area and volume for G-C-G were 83.979% and 77.257%, respectively; for G-G 20.360% and 23.600%; and for G-C 11.545% and 11.250%. The analysis of the representative element volume (REV) was conducted, and the feasibility and reliability of the micro-scale pattern extraction results were confirmed to guide the analysis of macro-scale characteristics. The exploration of the effectiveness of cementation patterns in fluid sealing provides valuable insights into effective biocementation at the pore scale of porous media, which may inspire future research.

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.

The porous structure induced heterogeneous and localized failure of the biofilm in microfluidic channels

Understanding the mechanism of biofilm distribution and detachment is very important to effectively improve water treatment and prevent blockage in porous media. The existing research is more related to the local biofilm evolving around one or few microposts and the lack of the integral biofilm evolution in a micropost array for a longer growth period. This study combines microfluidic experiments and mathematical simulations to study the distribution and detachment of biofilm in porous media. Microfluidic chips with an array of microposts with different sizes are designed to simulate the physical pore structure of soil. The research shows that the initial formation and distribution of biofilm are influenced by bacterial transport velocity gradients within the pore space. Bacteria prefer to aggregate areas with smaller microposts, leading to the development of biofilm in those regions. Consequently, impermeable blockage structures form in this area. By analyzing experimental images of biofilm structures at the later stages, as well as coupling fluid flow and porous medium, and the finite element simulation, we find that the biofilm detachment is correlated with the morphology and permeability (kb) (from 10-15 to 10-9 m2) of the biofilm. The simulations show that there are two modes of biofilm detachment, such as internal detachment and external erosion.

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.

Fluid transport in porous media based on differences in filter media morphology and biofilm growth in bioreactors

To elucidate the effect of pore structure on bioclogging and seepage flow in bioreactors, we used X-ray computed tomography (X-CT) to investigate the changes in seepage flow of porous media in zeolite, gravel and ceramsite bioreactors with biofilm growth by injecting a non-ionic contrast medium iohexol. Based on the X-CT images using a ball-and-stick model, the highest average pore radius (R‾) and the average pore throat radius (r‾) in the ceramsite column were found under the initial conditions, which facilitated its permeability. The pore and throat of the gravel column were small and homogeneous relatively. Biofilm growth decreased the pore and pore throat in the columns. The total throat area of zeolite, gravel and ceramsite columns declined by 74%, 73% and 79% respectively. The zeolite column had the highest average pore throat, which contributed to its maximum conductivity subsequently after biofilm growth. Further, the fractal dimensions of the pore structure increased with biofilm growth, especially in the zeolite and ceramsite columns. The heterogeneity of the porous media was reinforced by the biofilm growth in the zeolite and ceramsite columns due to their higher heterogeneity initially. We also observed that an increase in heterogeneity of porous medium amplified the preferential flow and flow heterogeneities, especially in the zeolite and ceramsite columns.

Bioclogging alleviation for constructed wetland based on the interaction among biofilm growth and hydrodynamics

Bioclogging is the most crucial operation problem of the constructed wetlands, which reduce its removal efficiency and life span. A strategy through properly increasing hydraulic loading is proposed in this study to alleviate the bioclogging for CWs. The two-dimensional porous media flow cell (2D PMFC) test indicated that a quadratic correlation was found between local biofilms growth rate and the near-wall Reynolds number (r > 0.765, p < 0.05). The biofilm growth rate declined with the flowrate when Re exceeded about 6.0. It was also found that the higher flowrate (6 mL/min) lead to the homogeneous biofilm and velocity distribution in the PMFC. The column test indicated that the highest hydraulic loading (9.2 cm/h) produced the smallest decrease in hydraulic conductivity, which was 80 times more than that of low hydraulic load (3.0 cm/h) at the end (40 days) of experiment. Moreover, the relatively homogenized distribution of biofilm was found along the column with the highest hydraulic loading, which confirmed that the proper increase in hydraulic loading can alleviate bioclogging.

Biomass-specific rates as key performance indicators: A nitrogen balancing method for biofilm-based electrochemical conversion

Microbial electrochemical technologies (METs) employ microorganisms utilizing solid-state electrodes as either electron sink or electron source, such as in microbial electrosynthesis (MES). METs reaction rate is traditionally normalized to the electrode dimensions or to the electrolyte volume, but should also be normalized to biomass amount present in the system at any given time. In biofilm-based systems, a major challenge is to determine the biomass amount in a non-destructive manner, especially in systems operated in continuous mode and using 3D electrodes. We developed a simple method using a nitrogen balance and optical density to determine the amount of microorganisms in biofilm and in suspension at any given time. For four MES reactors converting CO₂ to carboxylates, >99% of the biomass was present as biofilm after 69 days of reactor operation. After a lag phase, the biomass-specific growth rate had increased to 0.12-0.16 days^-1. After 100 days of operation, growth became insignificant. Biomass-specific production rates of carboxylates varied between 0.08-0.37 mol_C mol_X ^-1d^-1. Using biomass-specific rates, one can more effectively assess the performance of MES, identify its limitations, and compare it to other fermentation technologies.

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

Bacteria in porous media, such as soils, aquifers, and filters, often form surface-attached communities known as biofilms. Biofilms are affected by fluid flow through the porous medium, for example, for nutrient supply, and they, in turn, affect the flow. A striking example of this interplay is the strong intermittency in flow that can occur when biofilms nearly clog the porous medium. Intermittency manifests itself as the rapid opening and slow closing of individual preferential flow paths (PFPs) through the biofilm-porous medium structure, leading to continual spatiotemporal rearrangement. The drastic changes to the flow and mass transport induced by intermittency can affect the functioning and efficiency of natural and industrial systems. Yet, the mechanistic origin of intermittency remains unexplained. Here, we show that the mechanism driving PFP intermittency is the competition between microbial growth and shear stress. We combined microfluidic experiments quantifying Bacillus subtilis biofilm formation and behavior in synthetic porous media for different pore sizes and flow rates with a mathematical model accounting for flow through the biofilm and biofilm poroelasticity to reveal the underlying mechanisms. We show that the closing of PFPs is driven by microbial growth, controlled by nutrient mass flow. Opposing this, we find that the opening of PFPs is driven by flow-induced shear stress, which increases as a PFP becomes narrower due to microbial growth, causing biofilm compression and rupture. Our results demonstrate that microbial growth and its competition with shear stresses can lead to strong temporal variability in flow and transport conditions in bioclogged porous media.

Time-lagged interspecies interactions prevail during biofilm development in moving bed biofilm reactor

Clarifying the essential succession dynamics of interspecies interactions during biofilm development is crucial for the regulation and application of biofilm-based processes. In this study, regular and time-series phylogenetic molecular ecological networks (pMENs) were constructed to investigate ordinary and time-lagged interspecies interactions during biofilm development in a moving bed biofilm reactor (MBBR). Positive interactions dominated both regular (89.78%) and time-series (77.04%) ecological networks, suggesting that extensive cooperative behaviors facilitated biofilm development. The pronounced directional interactions (72.52%) in the time-series network further indicated that time-lagged interspecies interactions prevailed in the biofilm development process. Specifically, the proportion of directional negative interactions was higher than that of positive interactions, implying that interspecific competition preferred to be time-lagged. The time-series network revealed that module hubs exhibited extensive time-lagged positive interactions with their neighbors, and most of them exhibited altruistic behaviors. Keystone species possessing more positive interactions were positively correlated with biofilm biomass, NO₃ ^- -N concentrations, and the removal efficiencies of NH₄ ⁺ -N and COD. However, keystone species and peripherals that were negatively targeted by their neighbors showed positive correlations with the concentrations of NO₂ ^- -N, polysaccharides, and proteins in the soluble microbial products. The data highlight that the time-series network can provide directional microbial interactions along with the biofilm development process, which would help to predict the tendency of community shifts and propose efficient strategies for the regulation of biofilm-based processes. This article is protected by copyright. All rights reserved.

Systematic evaluation of PDA/PAM/MAH-modified basalt fibre as biofilm carrier for wastewater treatment

In this study, three kinds of modified basalt fibre (MBF) were used as biofilm carrier to treat wastewater, mainly for the removal of organic matter and nutrient pollutants, which was evaluated the feasibility of modification by pollutants removal performance. Polydopamine modified basalt fibre (PAD-BF) via the surface coating method were obtained. Polyacrylamide modified basalt fibre (PAM-BF) and maleic anhydride-modified basalt fibre (MAH-BF) via the surface grafting method were prepared. The surface physicochemical properties, biomass attachment capacity and pollutants removal efficiency of MBF were systematically investigated. Electron microscope scanning (SEM) revealed that the surface roughness of BF was obviously improved by modification. Besides, fourier transform infrared (FTIR) suggested that the MBF had more surface-active functional groups. The results of sludge immobilization tests showed that PDA/PAM/MAH-BF had higher bio-affinity than ordinary BF with 1.5∼2.3 times on immobilization ratio of microorganisms (IRM). Furthermore, the performances of PDA/PAM/MAH-BF as biofilm carrier for pollutants were significantly higher than that of ordinary BF group. Among them, the highest removal efficiency of COD in PAD-BF biofilm reactor was 95.29 ± 0.99%, while that of BF group was 86.30 ± 3.09%. PAM-BF group had the best removal effect of nutrients with the removal efficiency of 90.83 ± 7.69% for TP and 91.25 ± 6.43% for TN, respectively, while the removal rate of BF group was only about 70%. The improvement of dissolved oxygen (DO) in the MBF reactors was consistent with the enhancement of contaminant removal. Therefore, PDA/PAM/MAH-BF can be used as promising biological carrier fillers in wastewater treatment engineering.

Contributions of biofilm-induced flow heterogeneities to solute retention and anomalous transport features in porous media

Microbial biofilms are ubiquitous within porous media and the dynamics of their growth influence surface and subsurface flow patterns which impacts the physical properties of porous media and large-scale transport of solutes. A two-dimensional pore-scale numerical model was used to evaluate the impact of biofilm-induced flow heterogeneities on conservative transport. Our study integrates experimental biofilm images of Paenibacillus 300A strain in a microfluidic device packed with cylindrical grains in a hexagonal distribution, with mathematical modeling. Biofilm is represented as a synthetic porous structure with locally varying physical properties that honors the impact of biofilm on the porous medium. We find that biofilm plays a major role in shaping the observed conservative transport dynamics by enhancing anomalous characteristics. More specifically, when biofilm is present, the pore structure in our geometry becomes more spatially correlated. We observe intermittent behavior in the Lagrangian velocities that switches between fast transport periods and long trapping events. Our results suggest that intermittency enhances solute spreading in breakthrough curves which exhibit extreme anomalous slope at intermediate times and very marked late solute arrival due to solute retention. The efficiency of solute retention by the biofilm is controlled by a transport regime which can extend the tailing in the concentration breakthrough curves. These results indicate that solute retention by the biofilm exerts a strong control on conservative solute transport at pore-scale, a role that to date has not received enough attention.

Ecological Linkages between a Biofilm Ecosystem and Reactor Performance: The Specificity of Biofilm Development Phases

In biofilm-based engineered ecosystems, the reactor performance was closely linked to interspecies interactions within a biofilm ecosystem, whereas the ecological processes underpinning such linkage were still unenlightened. Herein, the principles of community succession and assembly were integrated to capture the ecological laws of biofilm development by molecular ecological networks and assembly model analysis based on the 16S rRNA sequencing analysis and metagenomics in a well-controlled moving bed biofilm reactor. At the initial colonization phase (days 0-2, driven by initial colonizers), interspecific cooperation (74.18%) facilitated initial biofilm formation, whereas some pioneers, and keystone species disappeared at later phases. At the accumulation phase (days 3-30, rapid biofilm development), interspecific cooperation (81.41 ± 5.07%) contributed to rapid biofilm development and keystone species were mainly involved in quorum sensing or positively correlated with extracellular polymeric substance production. At the maturation phase (days 31-106, a well-adapted quasi-equilibrium state), increased interspecific competition (32.74 ± 4.77%) and higher small-world property facilitated the rapid information transportation and pollutant treatment, and keystone species were positively correlated with the removal of COD and NH₄⁺-N. Homogenizing dispersal diminished the contemporary community dissimilarities, while turnover but rather nestedness governed the temporal variations in the biofilm succession period. This study highlighted the specificity of ecological processes at distinct biofilm development phases, which would advance our understanding on the development-to-function linkages in biofilm-based treatment processes.

Inorganic carbon limitation during nitrogen conversions in sponge-bed trickling filters for mainstream treatment of anaerobic effluent

Anaerobic sewage treatment is a proven technology in warm climate regions, and sponge-bed trickling filters (SBTFs) are an important post-treatment technology to remove residual organic carbon and nitrogen. Even though SBTFs can achieve a reasonably good effluent quality, further process optimization is hampered by a lack of mechanistic understanding of the factors influencing nitrogen removal, notably when it comes to mainstream anaerobically treated sewage. In this study, the factors that control the performance of SBTFs following anaerobic (i.e., UASB) reactors for sewage treatment were investigated. A demo-scale SBTF fed with anaerobically pre-treated sewage was monitored for 300 days, showing a median nitrification efficiency of 79% and a median total nitrogen removal efficiency of 26%. Heterotrophic denitrification was limited by the low organic carbon content of the anaerobic effluent. It was demonstrated that nitrification was impaired by a lack of inorganic carbon rather than by alkalinity limitation. To properly describe inorganic carbon limitation in models, bicarbonate was added as a state variable and sigmoidal kinetics were applied. The resulting model was able to capture the overall long-term experimental behaviour. There was no nitrite accumulation, which indicated that nitrite oxidizing bacteria were little or less affected by the inorganic carbon limitation. Overall, this study indicated the vital role of influent characteristics and operating conditions concerning nitrogen conversions in SBTFs treating anaerobic effluent, thus facilitating further process optimization.

A General Model for Biofilm-Driven Microbial Electrosynthesis of Carboxylates From CO2

Up to now, computational modeling of microbial electrosynthesis (MES) has been underexplored, but is necessary to achieve breakthrough understanding of the process-limiting steps. Here, a general framework for modeling microbial kinetics in a MES reactor is presented. A thermodynamic approach is used to link microbial metabolism to the electrochemical reduction of an intracellular mediator, allowing to predict cellular growth and current consumption. The model accounts for CO₂ reduction to acetate, and further elongation to n-butyrate and n-caproate. Simulation results were compared with experimental data obtained from different sources and proved the model is able to successfully describe microbial kinetics (growth, chain elongation, and product inhibition) and reactor performance (current density, organics titer). The capacity of the model to simulate different system configurations is also shown. Model results suggest CO₂ dissolved concentration might be limiting existing MES systems, and highlight the importance of the delivery method utilized to supply it. Simulation results also indicate that for biofilm-driven reactors, continuous mode significantly enhances microbial growth and might allow denser biofilms to be formed and higher current densities to be achieved.

Non-Invasive Measurement, Mathematical Simulation and In Situ Detection of Biofilm Evolution in Porous Media: A Review

The development of biofilms and the related changes in porous media in the subsurface cannot be directly observed and evaluated. The primary reason that the mechanism of biofilm clogging in porous media cannot be clearly demonstrated is due to the opacity and structural complexity of three-dimensional pore space. Interest in exploring methods to overcome this limitation has been increasing. In the first part of this review, we introduce the underlying characteristics of biofilm in porous media. Then, we summarize two approaches, non-invasive measurement methods and mathematical simulation strategies, for studying fluid–biofilm–porous medium interaction with spatiotemporal resolution. We also discuss the advantages and limitations of these approaches. Lastly, we provide a perspective on opportunities for in situ monitoring at the field site.

Modeling the impact of evolving biofilms on flow in porous media inside a microfluidic channel

This study integrates microfluidic experiments and mathematical modeling to study the impacts of biofilms on flow in porous media and to explore approaches to simplify modeling permeability with complicated biofilm geometries. E. coli biofilms were grown in a microfluidic channel packed with a single layer of glass beads to reach three biofilm levels: low, intermediate, and high, with biofilm ratios (β_r) of 2.7%, 17.6%, and 55.2%, respectively. Two-dimensional biofilm structures and distributions in the porous medium were modeled by digitizing confocal images and considering broad ranges of biofilm permeability (k_b) (from 10^-15 m² to 10^-7 m²) and biofilm porosity (ε_b) (from 0.2 to 0.8). The overall permeability of the porous medium (k), the flow pathways and the overall/local pressure gradients were found to be highly dependent on β_r and k_b but were moderately impacted by ε_b when the biofilm levels were high and intermediate with k_b>10^-11 m². When biofilm structures are well developed, simplified biofilm geometries, such as uniform coating and symmetric contact filling, can provide reasonable approximations of k.

Trait-specific dispersal of bacteria in heterogeneous porous environments: from pore to porous medium scale

The dispersal of organisms controls the structure and dynamics of populations and communities, and can regulate ecosystem functioning. Predicting dispersal patterns across scales is important to understand microbial life in heterogeneous porous environments such as soils and sediments. We developed a multi-scale approach, combining experiments with microfluidic devices and time-lapse microscopy to track individual bacterial trajectories and measure the overall breakthrough curves and bacterial deposition profiles: we, then, linked the two scales with a novel stochastic model. We show that motile cells of Pseudomonas putida disperse more efficiently than non-motile mutants through a designed heterogeneous porous system. Motile cells can evade flow-imposed trajectories, enabling them to explore larger pore areas than non-motile cells. While transported cells exhibited a rotation in response to hydrodynamic shear, motile cells were less susceptible to the torque, maintaining their body oriented towards the flow direction and thus changing the population velocity distribution with a significant impact on the overall transport properties. We also found, in a separate set of experiments, that if the suspension flows through a porous system already colonized by a biofilm, P. putida cells are channelled into preferential flow paths and the cell attachment rate is increased. These two effects were more pronounced for non-motile than for motile cells. Our findings suggest that motility coupled with heterogeneous flows can be beneficial to motile bacteria in confined environments as it enables them to actively explore the space for resources or evade regions with unfavourable conditions. Our study also underlines the benefit of a multi-scale approach to the study of bacterial dispersal in porous systems.

Variations in the microbial community of biofilms under different near-wall hydraulic shear stresses in agricultural irrigation systems

The hydraulic characteristics along agricultural irrigation pipelines directly affect the local near-wall hydraulic shear stress and biofilm accumulation. However, the variations in the microbial community during the process remain unknown. Based on the Couette-Taylor reactor, a device was developed to accurately control the hydraulic shear stress. The results indicated that the near-wall hydraulic shear stresses showed quadratic correlations with microbial contents (represented by phospholipid fatty acids r > 0.77, p < 0.05), and the maximum values were obtained under the shear stresses of 0.20-0.35 Pa. For two types of treated wastewater, the mutual operational taxonomic units among different shear stress treatments showed good consistency (>185). Their corresponding response in the microbial community was represented by the quantitative correlations between the near-wall hydraulic shear stresses and the polymorphism indices (r > 0.82, p < 0.05). Among the microorganisms, Firmicutes at the phylum level were significantly affected by the shear stress and significantly influenced the biofilm accumulation process.

Modelling Microbial Fuel Cells Using Lattice Boltzmann Methods

An accurate modelling of bio-electrochemical processes that govern Microbial Fuel Cells (MFCs) and mapping their behavior according to several parameters will enhance the development of MFC technology and enable their successful implementation in well defined applications. The geometry of the electrodes is among key parameters determining efficiency of MFCs due to the formation of a biofilm of anodophilic bacteria on the anode electrode, which is a decisive factor for the functionality of the device. We simulate the bio-electrochemical processes in an MFC while taking into account the geometry of the electrodes. Namely, lattice Boltzmann methods are used to simulate the fluid dynamics and the advection-diffusion phenomena in the anode compartment. The model is verified on voltage and current outputs of a single MFC derived from laboratory experiments under continuous flow. Conclusions can be obtained from a parametric analysis of the model concerning the design of the geometry of the anode compartment, the positioning and microstructure of the anode electrode, in order to achieve more efficient overall performance of the system. An example of such a parametric analysis is presented here, taking into account the positioning of the electrode in the anode compartment.

Testing Anti-Biofilm Polymeric Surfaces: Where to Start?

Present day awareness of biofilm colonization on polymeric surfaces has prompted the scientific community to develop an ever-increasing number of new materials with anti-biofilm features. However, compared to the large amount of work put into discovering potent biofilm inhibitors, only a small number of papers deal with their validation, a critical step in the translation of research into practical applications. This is due to the lack of standardized testing methods and/or of well-controlled in vivo studies that show biofilm prevention on polymeric surfaces; furthermore, there has been little correlation with the reduced incidence of material deterioration. Here an overview of the most common methods for studying biofilms and for testing the anti-biofilm properties of new surfaces is provided.

Challenges of biofilm control and utilization: lessons from mathematical modelling

This article reviews modern applications of mathematical descriptions of biofilm formation. The focus is on theoretically obtained results which have implications for areas including the medical sector, food industry and wastewater treatment. Examples are given as to how models have contributed to the overall knowledge on biofilms and how they are used to predict biofilm behaviour. We conclude that the use of mathematical models of biofilms has demonstrated over the years the ability to significantly contribute to the vast field of biofilm research. Among other things, they have been used to test various hypotheses on the nature of interspecies interactions, viability of biofilm treatment methods or forces behind observed biofilm pattern formations. Mathematical models can also play a key role in future biofilm research. Many models nowadays are analysed through computer simulations and continue to improve along with computational capabilities. We predict that models will keep on providing answers to important challenges involving biofilm formation. However, further strengthening of the ties between various disciplines is necessary to fully use the tools of collective knowledge in tackling the biofilm phenomenon.

A bilayer coarse-fine infiltration system minimizes bioclogging: The relevance of depth-dynamics

Bioclogging is a main concern in infiltration systems as it may significantly shorten the service life of these low-technology water treatment methods. In porous media, biofilms grow to clog partially or totally the pore network. Dynamics of biofilm accumulation (e.g., by attachment, detachment, advective transport in depth) and their impact on both surface and deep bioclogging are still not yet fully understood. To address this concern, a 104 day-long outdoor infiltration experiment in sand tanks was performed, using secondary treated wastewater and two grain size distributions (GSDs): a monolayer system filled with fine sand, and a bilayer one composed by a layer of coarse sand placed on top of a layer of fine sand. Biofilm dynamics as a function of GSD and depth were studied through cross-correlations and multivariate statistical analyses using different parameters from biofilm biomass and activity indices, plus hydraulic parameters measured at different depths. Bioclogging (both surface and deep) was found more significant in the monolayer fine system than in the bilayer coarse-fine one, possibly due to an early low-cohesive biofilm formation in the former, driven by lower porosity and lower fluxes; under such conditions biomass is favorably detached from the top layer, transported and accumulated in depth, so that new biomass might colonize the surface. On the other hand, in the bilayer system, fluxes are highest, and the biofilm is still in a growing phase, with low biofilm detachment capability from the top sand layer and high microbial activity in depth, resulting in low bioclogging. Overall, the bilayer coarse-fine system allows infiltrating higher volume of water per unit of surface area than the monolayer fine one, minimizing surface and deep bioclogging, and thus increasing the longevity and efficiency of infiltration systems.

Determination of mechanical properties of biofilms by modelling the deformation measured using optical coherence tomography

The advantage of using non-invasive imaging such as optical coherence tomography (OCT) to asses material properties from deformed biofilm geometries can be compromised by the assumptions made on fluid forces acting on the biofilm. This study developed a method for the determination of elastic properties of biofilms by modelling the biofilm deformation recorded by OCT imaging with poroelastic fluid-structure interaction computations. Two-dimensional biofilm geometries were extracted from OCT scans of non-deformed and deformed structures as a result of hydrodynamic loading. The biofilm geometries were implemented in a model coupling fluid dynamics with elastic solid mechanics and Darcy flow in the biofilm. The simulation results were compared with real deformed geometries and a fitting procedure allowed estimation of the Young's modulus in given flow conditions. The present method considerably improves the estimation of elastic moduli of biofilms grown in mini-fluidic rectangular channels. This superior prediction is based on the relaxation of several simplifying assumptions made in past studies: shear stress is not anymore taken constant over the biofilm surface, total stress including also pressure is accounted for, any biofilm shape can be used in the determinations, and non-linear behavior of mechanical properties can be estimated. Biofilm elastic moduli between 70 and 700 Pa were obtained and biofilm hardening at large applied stress due to increasing flow velocity was quantified. The work performed here opens the way for in-situ determination of other mechanical properties (e.g., viscoelastic properties, relaxation times, plastic yields) and provides data for modelling biofilm deformation and detachment with eventual applications in biofilm control and removal strategies.

Towards standardized mechanical characterization of microbial biofilms: analysis and critical review

Developing reliable anti-biofilm strategies or efficient biofilm-based bioprocesses strongly depends on having a clear understanding of the mechanisms underlying biofilm development, and knowledge of the relevant mechanical parameters describing microbial biofilm behavior. Many varied mechanical testing methods are available to assess these parameters. The mechanical properties thus identified can then be used to compare protocols such as antibiotic screening. However, the lack of standardization in both mechanical testing and the associated identification methods for a given microbiological goal remains a blind spot in the biofilm community. The pursuit of standardization is problematic, as biofilms are living structures, i.e., both complex and dynamic. Here, we review the main available methods for characterizing the mechanical properties of biofilms through the lens of the relationship linking experimental testing to the identification of mechanical parameters. We propose guidelines for characterizing biofilms according to microbiological objectives that will help the reader choose an appropriate test and a relevant identification method for measuring any given mechanical parameter. The use of a common methodology for the mechanical characterization of biofilms will enable reliable analysis and comparison of microbiological protocols needed for improvement of engineering process and screening.

Confined Flow: Consequences and Implications for Bacteria and Biofilms

Bacteria overwhelmingly live in geometrically confined habitats that feature small pores or cavities, narrow channels, or nearby interfaces. Fluid flows through these confined habitats are ubiquitous in both natural and artificial environments colonized by bacteria. Moreover, these flows occur on time and length scales comparable to those associated with motility of bacteria and with the formation and growth of biofilms, which are surface-associated communities that house the vast majority of bacteria to protect them from host and environmental stresses. This review describes the emerging understanding of how flow near surfaces and within channels and pores alters physical processes that control how bacteria disperse, attach to surfaces, and form biofilms. This understanding will inform the development and deployment of technologies for drug delivery, water treatment, and antifouling coatings and guide the structuring of bacterial consortia for production of chemicals and pharmaceuticals.

Pore-Scale Hydrodynamics in a Progressively Bioclogged Three-Dimensional Porous Medium: 3-D Particle Tracking Experiments and Stochastic Transport Modeling

Biofilms are ubiquitous bacterial communities that grow in various porous media including soils, trickling, and sand filters. In these environments, they play a central role in services ranging from degradation of pollutants to water purification. Biofilms dynamically change the pore structure of the medium through selective clogging of pores, a process known as bioclogging. This affects how solutes are transported and spread through the porous matrix, but the temporal changes to transport behavior during bioclogging are not well understood. To address this uncertainty, we experimentally study the hydrodynamic changes of a transparent 3-D porous medium as it experiences progressive bioclogging. Statistical analyses of the system's hydrodynamics at four time points of bioclogging (0, 24, 36, and 48 h in the exponential growth phase) reveal exponential increases in both average and variance of the flow velocity, as well as its correlation length. Measurements for spreading, as mean-squared displacements, are found to be non-Fickian and more intensely superdiffusive with progressive bioclogging, indicating the formation of preferential flow pathways and stagnation zones. A gamma distribution describes well the Lagrangian velocity distributions and provides parameters that quantify changes to the flow, which evolves from a parallel pore arrangement under unclogged conditions, toward a more serial arrangement with increasing clogging. Exponentially evolving hydrodynamic metrics agree with an exponential bacterial growth phase and are used to parameterize a correlated continuous time random walk model with a stochastic velocity relaxation. The model accurately reproduces transport observations and can be used to resolve transport behavior at intermediate time points within the exponential growth phase considered.

Flow environment and matrix structure interact to determine spatial competition in Pseudomonas aeruginosa biofilms

Bacteria often live in biofilms, which are microbial communities surrounded by a secreted extracellular matrix. Here, we demonstrate that hydrodynamic flow and matrix organization interact to shape competitive dynamics in Pseudomonas aeruginosa biofilms. Irrespective of initial frequency, in competition with matrix mutants, wild-type cells always increase in relative abundance in planar microfluidic devices under simple flow regimes. By contrast, in microenvironments with complex, irregular flow profiles – which are common in natural environments – wild-type matrix-producing and isogenic non-producing strains can coexist. This result stems from local obstruction of flow by wild-type matrix producers, which generates regions of near-zero shear that allow matrix mutants to locally accumulate. Our findings connect the evolutionary stability of matrix production with the hydrodynamics and spatial structure of the surrounding environment, providing a potential explanation for the variation in biofilm matrix secretion observed among bacteria in natural environments.

DOI:http://dx.doi.org/10.7554/eLife.21855.001

Bacteria often live together – attached to surfaces like river rocks, water pipes, the lining of the gut and catheters – in communities called biofilms. These groups of bacteria are small-scale ecosystems in which cells cooperate and compete with one another to obtain resources, such as food and space to grow. Within a biofilm, a sticky glue-like substance called the matrix binds the cells to each other and to the surface. Cells that make the matrix typically have an advantage over those that do not because they can better resist the shearing forces experienced when liquid flows over the surface. The matrix also helps cells to capture nutrients from the passing liquid. Nevertheless, not all strains of bacteria make matrix, despite its advantages.

Because of where they can grow, biofilms are fundamentally important in the environment, in industry and in infections. Resolving why some bacteria make matrix while others do not could therefore allow scientists and engineers to re-design the surfaces involved in these settings to discourage harmful biofilms or to encourage beneficial ones.

Nadell, Ricaurte et al. have now used a bacterium called Pseudomonas aeruginosa to explore how the properties of the surface and the flowing liquid affect matrix production among cells in biofilms. P. aeruginosa typically lives in soil and can cause infections in people, especially in hospital patients and people who have weakened immune systems. Nadell, Ricaurte et al. studied normal P. aeruginosa bacteria and a mutant strain that is unable to make matrix. The strains were labeled with fluorescent markers and put into special chambers that simulated different environments. The proportion of each strain was measured after three days of biofilm growth. When biofilms were grown under flowing liquid in simple environments with flat surfaces, matrix producers always outcompeted non-producers. However, the two strains coexisted in more complex and porous environments, like those found in soil.

Nadell, Ricaurte et al. went on to show that the strains could co-exist because the matrix producers made biofilms that created areas within the environment where the liquid flows very slowly or not at all. In these regions, non-producing cells could compete successfully because resistance to shearing forces is less important when flow is weak or absent, and so the non-producing cells were not washed away. The results begin to explain why matrix production among cells in environmental settings is diverse and highlight that the environment is important in the evolution of bacterial biofilms.

DOI:http://dx.doi.org/10.7554/eLife.21855.002

Effect of water temperature on biofouling development in reverse osmosis membrane systems

Understanding the factors that determine the spatial and temporal biofilm development is a key to formulate effective control strategies in reverse osmosis membrane systems for desalination and wastewater reuse. In this study, biofilm development was investigated at different water temperatures (10, 20, and 30 °C) inside a membrane fouling simulator (MFS) flow cell. The MFS studies were done at the same crossflow velocity with the same type of membrane and spacer materials, and the same feed water type and nutrient concentration, differing only in water temperature. Spatially resolved biofilm parameters such as oxygen decrease rate, biovolume, biofilm spatial distribution, thickness and composition were measured using in-situ imaging techniques. Pressure drop (PD) increase in time was used as a benchmark as to when to stop the experiments. Biofilm measurements were performed daily, and experiments were stopped once the average PD increased to 40 mbar/cm. The results of the biofouling study showed that with increasing feed water temperature (i) the biofilm activity developed faster, (ii) the pressure drop increased faster, while (iii) the biofilm thickness decreased. At an average pressure drop increase of 40 mbar/cm over the MFS for the different feed water temperatures, different biofilm activities, structures, and quantities were found, indicating that diagnosis of biofouling of membranes operated at different or varying (seasonal) feed water temperatures may be challenging. Membrane installations with a high temperature feed water are more susceptible to biofouling than installations fed with low temperature feed water.

Filaments in curved streamlines: Rapid formation of Staphylococcus aureus biofilm streamers

Biofilms are surface-associated conglomerates of bacteria that are highly resistant to antibiotics. These bacterial communities can cause chronic infections in humans by colonizing, for example, medical implants, heart valves, or lungs. Staphylococcus aureus, a notorious human pathogen, causes some of the most common biofilm-related infections. Despite the clinical importance of S. aureus biofilms, it remains mostly unknown how physical effects, in particular flow, and surface structure influence biofilm dynamics. Here we use model microfluidic systems to investigate how environmental factors, such as surface geometry, surface chemistry, and fluid flow affect biofilm development in S. aureus. We discovered that S. aureus rapidly forms flow-induced, filamentous biofilm streamers, and furthermore if surfaces are coated with human blood plasma, streamers appear within minutes and clog the channels more rapidly than if the channels are uncoated. To understand how biofilm streamer filaments reorient in flows with curved streamlines to bridge the distances between corners, we developed a mathematical model based on resistive force theory of slender filaments. Understanding physical aspects of biofilm formation in S. aureus may lead to new approaches for interrupting biofilm formation of this pathogen.