Abstract
The vast majority of bacteria live in sessile biofilms that colonize the channels, pores and crevices of confined structures. Flow in these structures carries the nutrients necessary for growth, but also generates stresses and detachment from surfaces. Conversely, bacteria tend to occupy a large part of the available space and, in so doing, increase resistance to flow and modify transport properties. Although the importance of advective transport and hydrodynamic forces on bacteria is well known, the complex feedback effects that control development in confined geome-tries are much less understood. Here, we study how couplings between flow and bacterial development control the spatio-temporal dynamics of Pseudomonas aeruginosa in microchannel flows. We demonstrate that nutrient limitation drives the longitudinal distribution of biomass, while a competition between growth and flow-induced detachment controls the maximum clogging and the temporal dynamics. We find that successive cycles of sloughing and growth cause persistent fluctuations of the hydraulic resistance and prevent the system from ever reaching a true steady-state. Our results indicate that these self-sustained fluctuations are a signature effect of biofilm development in confined flows and could thus be a key component of the spreading of biofilms in infections, environmental processes and engineering applications. Consistent with the description of other bursting events, such as earthquakes or avalanches, we further show that the dynamics of sloughing can be described as a jump stochastic process with a gamma distribution of interevent times. This stochastic modeling approach opens the way towards a new quantitative approach to the characterization of the apparent randomness and irreproducibility of biofilm experiments in such systems.
eLife assessment
This important study integrates microfluidic experiments and mathematical modeling to investigate how flow dynamics and biofilm growth and detachment influence each other. Using Pseudomonas aeruginosa as a model organism, the study identifies several key effects and stages in biofilm development, albeit with some weaknesses in clearly defining the setup and some of their interpretations. The comparison between experimental results and theoretical models is convincing, providing a robust analysis of the biofilm's behavior under varying flow conditions. The findings will be helpful for researchers working on biofilms and their applications.
Introduction
Fluid flow and transport phenomena control many aspects of the life of bacteria (Krsmanovic et al., 2021), from the motility of cells (Busscher and van der Mei, 2006; Rodesney et al., 2017) to the morphology of biofilm colonies (Wang et al., 2022), and even ecological interactions within populations (Battin et al., 2016). As a testimony to the importance of flow, bacteria have evolved specific strategies adapted to their mechanical environment and the rheology of fluids around them (Dufrêne and Persat, 2020). Bacteria have also evolved mechanisms to detect gradients of nutrients or toxic substances and adapt their movement accordingly (Sampedro et al., 2014). Recent studies further suggest that bacteria have mechanosensing and rheosensing capabilities (Dufrêne and Persat, 2020). For example, Pseudomonas aeruginosa has been found to regulate the fro operon in response to the shear rate (Sanfilippo et al., 2019). Shear also modifies the intracellular levels of cyclic-di-GMP in cells of P. aeruginosa attached to a surface, initiating a sessile phenotype (Rodesney et al., 2017).
Even in dense surface-associated colonies known as biofilms (Costerton et al., 1995; Flemming et al., 2016), where cells are partially isolated from the fluid by a matrix of self-produced extracellular polymeric substances (Flemming and Wingender, 2010), flow still plays a fundamental role (Thomen et al., 2017; Rusconi et al., 2011). The flow of a viscous fluid around the matrix generates forces that can induce detachment or remodel the biofilm (Besemer et al., 2007; Stoodley et al., 1998) – the matrix essentially behaves as a viscoelastic material (Peterson et al., 2015) and thus can flow in response to stress (Gloag et al., 2020). Furthermore, key solutes, such as oxygen or nutrients, are transported in the fluid before they can diffuse in the biofilm and be consumed. This can lead to complex couplings between advective transport by flow, diffusion in the fluid and in the biofilm, and uptake by the cells (Taherzadeh et al., 2012; Picioreanu et al., 2000). It has also been recently demonstrated that flow can generate spatial heterogeneities in the activation of quorumsensing within populations, as a result of autoinducers being washed away by flow in zones at high Péclet number (Kim et al., 2016; Emge et al., 2016; Mukherjee and Bassler, 2019).
A vast majority of bacteria live in biofilms that colonize confined geometries (Conrad and Poling-Skutvik, 2018). In such systems, biofilms can occupy a large portion of the available space and thus severely restrict flow. This tends to generate a two-way coupling between biofilm growth and transport phenomena whereby flow and transport mediate the development of the biofilm, but in return the development of the biofilm also modifies flow and transport (Rittmann, 1993; Taylor and Jaffé, 1990a; Taylor et al., 1990; Taylor and Jaffé, 1990b,c; Vandevivere and Baveye, 1992b,a; Stewart and Fogler, 2001; Telgmann et al., 2004; Drescher et al., 2013). In porous media, for instance, Kurz et al. (2022) showed that a Bacillus subtilis biofilm can clog a large part of the pore space, leaving only few preferential flow channels where a competition between shear-induced detachment and growth drives intermittency and spontaneous pressure fluctuations. Such bioclogging is an important component of systems as diverse as bioreactors in the food industry (Verran, 2002), biofilters for wastewater processing (Aslam et al., 2018), pipe flow for water distribution (Cowle et al., 2014), heat exchangers (García and Trueba, 2020), catheters used in medicine (Bixler and Bhushan, 2012), soil bioremediation (Singh et al., 2006), enhanced oil recovery (Sen, 2008; Kryachko, 2018) or biobarriers (Lennox and Ashe, 2009). Understanding the fundamentals of biofilm growth and feedback mechanisms with flow is therefore a crucial step towards developing better approaches in health and engineering.
Here, we investigate the role of couplings between nutrient transport, growth and detachment on biofilm development in microchannel flows. To do so, we developed a microfluidic setup generating a constant flow rate in a microchannel where a P. aeruginosa PAO1 GFP biofilm develops. This microfluidic system is further combined with timelapse microscopy, cellular microbiology and mathematical modeling to study the interactions between biofilm and flow, in particular the dynamics and spatio-temporal fluctuations.
Results
Our experimental approach is summarized in Fig 1. We inoculated cells of P. aeruginosa PAO1 GFP in a microchannel and then continuously flowed a culture medium at constant flow rate to observe biofilm development. In what follows, we first study the effects of nutrient limitation on the spatial distribution of biofilm. We then focus on regimes where nutrients are in excess and biofilm development is primarily driven by the interactions with flow to detail the role of flowinduced detachment on the maximum clogging reachable and the temporal dynamics.
![](https://prod--epp.elifesciences.org/iiif/2/98292%2Fv1%2Fcontent%2F597893v1_fig1.tif/full/max/0/default.jpg)
Schematic representation of the system and of the two main experimental steps. The first step is P. aeruginosa PAO1 GFP culture and inoculation in the microchannel (PDMS on glass, 100 μm × 100 μm cross-section) using a pressure pump. The bacterial suspension is then left for 3 hours without flow to allow cells to adhere. The second step consists in flowing the culture medium at constant flow rate through the microchannel, while recording pressure fluctuations and imaging biofilm development via photonic microscopy. UVC radiation is used during the second step of the experiment to constrain the biofilm in a specific part of the channel.
Nutrient limitation controls the longitudinal distribution of the biofilm
To assess the impact of nutrient limitation, we performed experiments at flow rates over several orders of magnitude (Q = 2 × 10−2, 2 × 10−1, 2 × 100, and 2 × 10+1μL/min) and therefore different total fluxes of nutrients. Fig 2a shows the time-integrated distributions of the GFP fluorescence in the longitudinal direction for the different flow rates. We found that the active biomass expressing GFP is relatively uniform for flow rates between 0.2 and 20μL/min – the sharp decrease on the left-hand side at the outlet corresponds to the effect of the UVCs. For 0.02 μL/min, however, we observed a maximum value of the fluorescence intensity at the inlet on the right-hand side and a distinct decrease when moving towards the outlet.
![](https://prod--epp.elifesciences.org/iiif/2/98292%2Fv1%2Fcontent%2F597893v1_fig2.tif/full/max/0/default.jpg)
Impact of nutrient limitation on the longitudinal distribution of biofilm. (a): Fluorescence intensity integrated for 72 hours for the flow rates Q = 0.02 μL/min, 0.2 μL/min, 2 μL/min and 20 μL/min with 1× concentrated brain heart infusion (BHI) culture medium, along with Q = 0.02 μL/min with either 0.2×BHI or 1×BHI supplemented with 8 g/L glucose. Nutrient limitation is observed only for Q = 0.02 μL/min and is strongly dependent upon the concentration of BHI components. (b): Simulations for Q = 0.02 μL/min for two values of ϕmax and the corresponding Damköhler numbers. The Damköhler numbers for the case with 0.2×BHI were simply obtained by multiplying those for the case 1×BHI by a factor 5. Each experimental curve in (a) and (b) is averaged over 3 replicates. In the model, we considered that the fluorescence signal is proportional to the product of biomass and concentration. We then integrated this signal in time and normalized it with the maximum value as
Considering that the growth of P. aeruginosa PAO1 GFP is aerobic here – P. aeruginosa is a facultative anaerobe and can perform denitrification in anaerobic environments by anoxic respiration using nitrogenated compounds as final electron acceptors (Arat et al., 2015) but this produces less energy (Bartberger et al., 2002) – we hypothesized that this heterogeneity in biofilm development along the channel stems from a limitation of solute species required for growth, either oxygen or one of the nutrients in the growth medium. For the largest flow rates, advective transport through the channel should be fast compared to consumption, so that even bacteria at the outlet receive sufficient levels of nutrients and oxygen for biofilm development. For the lowest flow rate, however, one of the components introduced at the inlet is rapidly consumed by bacteria and becomes limiting, so that growth decreases with the distance from the inlet.
To better understand this limitation, we repeated experiments at 0.02 μL/min either with BHI diluted 5 times or by supplementing 1×BHI with additional glucose. Results in Fig 2a with the glucose supplementation show a similar longitudinal distribution of biofilm compared to assays without glucose, thus suggesting that this carbon source is not the limiting nutrient – it was confirmed by mass spectrometry that only a small fraction of the glucose was consumed. Results with the five times diluted BHI, however, show a much narrower window of biofilm growth, therefore suggesting that one of the components in the BHI is becoming limiting. Although oxygen probably features gradients within the biofilm (Folsom et al., 2010) and in the longitudinal direction, the inlet concentration of oxygen is identical for the two cases 1×BHI and 0.2×BHI, therefore indicating that oxygen is not the primary component limiting growth. Another factor that may further alleviate oxygen limitations is that PDMS is highly permeable to oxygen, so that there are in fact two sources of oxygen in our system: dissolved oxygen in the culture medium and oxygen coming from/through the PDMS.
To explore these different hypotheses and quantify the characteristic times of reactive transport, we simulated the development of the biofilm inside the channels, taking into account the couplings between biofilm growth and nutrient transport. The limiting nutrient is treated as a solute being transported by advection/diffusion and consumed by bacteria. We considered that mass transport is much faster than bacterial growth and division, so that the problem is quasi-steady for solute transport (Picioreanu et al., 2000). The limiting nutrient was thus modeled as
where
![](https://prod--epp.elifesciences.org/iiif/2/98292%2Fv1%2Fcontent%2F597893v1_tbl1.tif/full/max/0/default.jpg)
Summary of various quantities in the experiments and models.
σ0 is the shear stress in empty channel.
![](https://prod--epp.elifesciences.org/iiif/2/98292%2Fv1%2Fcontent%2F597893v1_tbl2.tif/full/max/0/default.jpg)
Values obtained from direct fitting of the experimental data
For the biomass growth, we used a model capturing the coupled effects of growth and flowinduced removal,
where
A number of works express the detachment rate as proportional to the square root of the fluid shear stress at the biofilm fluid interface (Coyte et al., 2017). Here this could be written as
Assuming that ℳ ∈ ]0, 1], the model for the biomass can be written in a simpler form as
with ϕmax the maximum volume fraction of biofilm, defined as ℳ = 1 − ϕmax, ϕmax ∈ [0, 1[ and ϕ ∈ [0, ϕmax]. The only non-trivial steady-state,
reflects an equilibrium between the amount of biomass created and the amount of biomass detached. The trivial steady-state ϕequ = 0 is stable when detachment exceeds growth, i.e. when C⋆ is lower than
Fig 2(b) compares results of the model with those of the experimental fluorescence for the flow rate 0.02 μL/min and K⋆ = 0.1. To roughly assess the sensitivity of our simulation to uncertainties in the value of ϕmax, we further considered two extreme values, ϕmax ≃ 0.6 and ϕmax ≃ 0.9. The model shows good agreement with experiments for DaBHI×1 ≃ 1200 in the case ϕmax ≃ 0.6 and DaBHI×1 ≃ 800 in the case ϕmax ≃ 0.9. The corresponding characteristic reaction time for nutrient uptake, τreac = C0/α, ranges from ≃ 83s to ≃ 125s. We further obtain an excellent correspondence between the model and the experimental data for the 5 times dilution of the BHI without any fitting, simply by multiplying the Damköhler number by a factor 5 – recall that, by definition,
With these estimations of characteristic times for uptake, we can also evaluate whether any form of radial limitation of nutrient is expected by comparing times for transverse diffusion and uptake. The characteristic time for diffusion in the radial direction is
The complete picture for nutrient transport is therefore as follows. In the longitudinal direction, the distribution of biomass is driven by a competition between advective transport and reaction. This can be quantified using the ratio of the Péclet and Damköhler numbers, Pe/Da. When advective transport is faster than uptake, Pe/Da > 1, the distribution of the biofilm is uniform. On the other hand, when Pe/Da < 1, uptake is faster than advective transport so that growth is nutrient limited. Table 1 summarizes the values of Pe/Da for the different cases and shows that the regime Pe/Da < 1 is reached only for the flow rate 0.02 μL/min. In the transverse direction, there is no nutrient limitation and the amount of biomass is determined by an equilibrium between growth and flow-induced detachment. The details of this mechanisms, and its impact on the temporal dynamics, is the core of the following sections.
Hydrodynamic stresses affect all stages of development
We now focus on the flow rates Q = 0.2, 2 and 20 μL/min that featured no strong nutrient limitation and thus allowed us to separate the effects of molecular transport from that of hydrodynamic stresses. We describe the temporal dynamics through two classes of measurements: the time evolution of the pressure and timelapse microscopy imaging. The evolution of the pressure was monitored in the inlet reservoir, while imposing a constant flow rate. From this measurement, we reconstructed the evolution in time of the hydraulic resistance of the zone of interest in the channel, ℛ (t), where biofilm developed (see Material and Methods). This allowed us to determine the dynamics of the ratio
![](https://prod--epp.elifesciences.org/iiif/2/98292%2Fv1%2Fcontent%2F597893v1_fig3.tif/full/max/0/default.jpg)
Temporal dynamics of growth and detachment for Q = 0.2 μL/min, 2 μL/min and 20 μL/min. (a): Summary of the two main stages of biofilm development. The figure shows the evolution of the hydraulic resistance (black solid line) in time, calculated from pressure measurements, and microscopy images corresponding to the different phases. It also shows changes in biofilm colonization extracted from either integrated fluorescence intensity (green squares) or from image segmentation (red solid line). (b), (c) and (d): Temporal dynamics of growth and detachment for the different flow rates. (e), (f) and (g): Wavelet scalograms corresponding to (b), (c) and (d).
We found that biofilm development occurs in different phases, in a way that is consistent across flow rates. We propose to decompose the dynamics in two main components: Stage I that corresponds to the initial adhesion, growth and saturation; and stage II that features large fluctuations in the amount of biomass in the channel. Fig 3(a) illustrates our decomposition in two stages on the temporal evolution of
Early stage I is driven by surface adhesion, motility, division and colony formation
Our inoculation process results in the sparse attachment of individual bacterial cells on the boundaries of the channel. On the glass slide, we evaluated the initial density to be about 38800 cells per mm2. Growth of adhered cells started straight away upon flowing the culture medium, with little to no lag time. The apparent macroscopic lag, as visible in Fig 3, does not stem from a lag at the cellular level, but rather from the sensitivity of the pressure measurements. The relative hydraulic resistance,
To characterize the growth time at the cellular level, we calculated an apparent doubling time at the cellular level from microscopy images of cells attached to the glass slide in the very early stage (from 0 to 3.5 hours, see Fig 4). We first segmented images from differential interference contrast microscopy to identify individual bacteria on the surface and then fitted linearly the log of the number of cells as a function of time. Calculated doubling times were measured as about 198, 95 and 117 minutes on average for, respectively, Q = 0.2 μL/min, 2 μL/min and 20 μL/min. As a reference, the doubling time in liquid culture was measured as about 110 minutes (standard deviation of ~ 10 minutes). Although understanding exactly what generates this dependence of the doubling time upon the flow rate is beyond the scope of this paper, it is worth noting that nutrient limitation is not likely to play a role, as we have previously validated the fact that there is no limitation for Q ≥ 0.2 μL/min, even when biofilm has formed. Mechano-and rheo-sensing could be involved, with bacteria directly regulating their response to perceived shear stress or rate (Sanfilippo et al., 2019; Rodesney et al., 2017), as well as heterogeneities in the rate of division – for instance between motile and adhered cells (Rossy et al., 2019) – or in the attachment/detachment dynamics.
![](https://prod--epp.elifesciences.org/iiif/2/98292%2Fv1%2Fcontent%2F597893v1_fig4.tif/full/max/0/default.jpg)
Flow modifies the apparent doubling time of bacteria on surfaces. (a), (b), (c) and (d): Composite brightfield and GFP images of development stages, starting from single cells that form microcolonies and then evolve towards a biofilm. (e) and (f) show, respectively, the average number of cells on the surface as a function of time and the corresponding doubling time for flow rates (Q = 0.2 μL/min, 2 μL/min and 20 μL/min). For (f), statistical differences were examined by unpaired student test with Gaussian distribution of data and equal standard deviations. Error bars indicate standard error of mean (SEM) and symbols denote statistical significance (****: p-value <0.0001, ***: p-value = 0.0002, ns: p-value > 0.05). The doubling time was calculated by a linear fitting of the logarithm of the number of cells. The slope was used to estimate growth rate and doubling time. Cell count was calculated from image segmentation of four positions in two channels to generate 8 measurements by condition (n= 8) for (Q = 0.2 μL/min, 2 μL/min and 20 μL/min).
After this initial phase, colonies expand and produce large quantities of EPS that yield a much sharper increase of the hydraulic resistance. This increase is linked to a combination between the exponential growth of the microorganisms and the production of the EPS, but also to the nonlinear relationship between
Late Stage I reflects an equilibrium between growth and detachment
We observe an inflection of the signal in late Stage I in Fig 3. This inflection corresponds to the end of the rapid growth phase, with flow-induced detachment progressively increasing, until detachment equilibrates with growth. To further understand this effect, let us consider the mass balance of biofilm in the zone of interest in the channel. The inlet flux of bacteria is zero, since our UVC system prevents growth outside the zone of interest. The source of biomass therefore only results from the uptake of nutrients, division of bacterial cells and EPS production. The sink of biomass is due to flow-induced removal, which includes parts of the biofilm that are displaced out of the zone of interest through the flow/remodelling of the biofilm at the outlet, and parts that are washed away by erosion and seeding (Kaplan, 2010; Krsmanovic et al., 2021). In the absence of nutrient limitation, we can estimate the evolution of the biovolume in the channel from the ordinary differential equation
which is a direct simplification of Eq 3 with ϕ describing an average volume fraction in the entire channel, and thus being only a function of time. This model captures a competition between growth and different types of smooth detachment (flow, seeding, erosion).
Stage II features sloughing-induced jump events followed by re-growth
After growth and detachment start to equilibrate at the end of Stage I, sloughing events become particularly important. These are visible in Fig 3 as jumps in the hydraulic resistance (see also the blue stripes on the wavelet scalograms in Fig 3), which correlate with sharp changes in the biofilm colonization of the channel. The jumps correspond to a range of different events, from minor detachment – where a relatively small portion of the biomass is detached – to major sloughing – where a large portion of the biomass is detached. Fig 5 shows example microscopy images of such events.
![](https://prod--epp.elifesciences.org/iiif/2/98292%2Fv1%2Fcontent%2F597893v1_fig5.tif/full/max/0/default.jpg)
Spatio-temporal dynamics of sloughing in Stage II. (a) and (b): composite bright field and GFP image of a 40 hour biofilm during Stage II. Image before (top) and after (bottom) for (a) a minor sloughing event and (b) a major sloughing event. (c) and (d): kymographs showing the fluorescence intensity (averaged in the radial direction y) as a function of both the longitudinal direction (x) and time. Fluorescence intensity values were normalized by the maximum value. Plots on the right-hand side show the corresponding hydraulic resistance as a function of time.
To better visualize the spatio-temporal dynamics and connect the different observations for the pressure and microscopy, we plotted kymographs in Fig 5 (see also SI figures S7, S10 and S13), showing both variations in time and in space, along with the corresponding hydraulic resistance signal. These graphics show clearly the correlation between the pressure signal on the right-hand side and the detachment events on the kymographs. We can readily identify that large drops in pressure correspond to large events with detachment over almost the entire length of the microchannel. We also visualize a range of detachment events corresponding to various sizes of biofilm being detached, along the rapid increase of resistance and biomass following each sloughing event.
Each sloughing event is followed by a rapid growth phase resembling late Stage I, leading to a cycle of successive sloughing and re-growth. In each re-growth phase, we observe a sigmoid-like shape of the hydraulic resistance signal that is similar to Stage I. As in Stage I, the maximum stems from an equilibrium between growth and detachment. ϕmax thus corresponds to a critical value of the stress, σcrit, that generates enough detachment to compensate growth. From Eqs 2-3, this reads
with ℛmax the maximum value of the hydraulic resistance. The maximum value of the hydraulic resistance over all stages of development decreases by orders of magnitudes from 0.2 to
Pressure becomes prominent in late Stage I and Stage II
Here, we ask the question of the nature of the hydrodynamic stress and, in particular, whether the hydrodynamic stress is dominated by shear or pressure. We are particularly interested in understand the nature of σcrit. Fig 6 shows simple 2D flow simulations illustrating the types of stress induced by the flow of a viscous fluid upon the biofilm. We see that the shear stress becomes particularly strong in the bottlenecks and that the pressure difference also builds up, therefore generating both shear and pressure stresses. In the initial stages of bacterial development, when only individual cells and microcolonies are adhered to surfaces, we expect the dominant stress to be shear. However, what happens in later stages for large fractions of biofilm?
![](https://prod--epp.elifesciences.org/iiif/2/98292%2Fv1%2Fcontent%2F597893v1_fig6.tif/full/max/0/default.jpg)
2D simulations of flow around a biofilm (white) in a channel of 100 μm width and 500 μm length using COMSOL multiphysics. (a): Streamlines and magnitude of the velocity (m/s). (b): Shear rate in (s−1). (c): Pressure field in (Pa). The white arrows indicate the flow direction.
A simple approach to quantifying the relative importance of stresses in our system is to consider the case of uniform film growth between the UVC zones (see details in supplementary information). The tangential stress at the solid/biofilm surface is
Maximum clogging depends on the flow rate and the matrix composition
In our experiments, the maximum values of the volume fraction, estimated from hydraulic resistance, are ϕmax (Q = 0.2) = 0.989 for 0.2, ϕmax (Q = 2) = 0.966 for 2 and ϕmax (Q = 20) = 0.901 for 20 (see Fig 7 and Table 1) – remarkably, this maximum is very reproducible across replicates, with standard deviations for volume fractions in the range of 10−3 in our experiments (Fig 7).
![](https://prod--epp.elifesciences.org/iiif/2/98292%2Fv1%2Fcontent%2F597893v1_fig7.tif/full/max/0/default.jpg)
Evolution of the distribution and maximum of the volume fraction for the different flow rates. (a), (b), (c) and (d): Fractions of biofilm in the microchannel, either calculated from hydraulic resistance (black solid line) or estimated from integrated GFP intensity (green dotted line), for the different flow rates. (e): Distribution of biofilm fraction between 24 and 72 h for all flow rates, represented as whisker boxes. (f): Log-log plot of 1 − ϕmax as a function of the flow rate (Q = 0.2 μL/min, 2 μL/min, 20 μL/min and 200 μL/min). The red dotted line simply shows the slope for an evolution with the square root of the flow rate.
The model Eq 5 for Stage I is consistent with this dependence upon the flow rate. By definition of ϕmax in Eqs 2-3 and in Eq 5, we have
We can express this relatively to a reference flow rate as
Upon considering Qref ≡ 0.2 with ϕmax (Qref) = 0.989, we obtain
To further validate this idea that ϕmax decreases with the square root of the flow rate, we performed a single experiment at 200 μL/min. Results in Fig 7 indicate that the behavior is similar to that of other flow rates, but with a significantly lower value of ϕmax ≃ 0.786. We also see that the scaling remains remarkably similar with
Modifications of the composition and mechanical properties of the biofilm, in particular if they affect 𝒳 or sloughing, should have a significant impact on the value of ϕmax. Among P. aeruginosa PAO1 polysaccharides – Pel, Psl and alginate – Pel and Psl are considered the two primary components of the EPS matrix structure (Ryder et al., 2007), controlling the cohesive and adhesive strength of the biofilm. We therefore performed experiments with mutants that cannot produce either Pel or Psl, with the idea that this should modify the value of 𝒳. The Psl mutant is a pslD deficient ΔpslD obtained from P. aeruginosa PAO1 by non-polar allelic exchange (Colvin et al., 2012). The Pel mutant is a pelF deficient strain ΔpelF also obtained by allelic exchange (Colvin et al., 2011). Liquid culture showed that the lag and doubling times were, respectively, 233.36 and 105.5 minutes (see Fig S21). Fig 8 (see also S14, S15, S16 and S17) shows the behavior of mutants compared to the wild type. We found that the volume fraction of biofilm was significantly lower for both mutants and that the Psl mutant was more strongly affected than the Pel mutant. The Psl mutant only weakly attached to the surface and was more prone to detachment – an observation that is consistent with the prominent role of Psl in the mechanics of PAO1 biofilms (Colvin et al., 2011). For the case 2 μL/min, we even observed an almost complete detachment of the biofilm with only few cells remaining attached and making re-growth possible after a sloughing event. At 20 μL/min, we could not detect the presence of biofilm.
![](https://prod--epp.elifesciences.org/iiif/2/98292%2Fv1%2Fcontent%2F597893v1_fig8.tif/full/max/0/default.jpg)
Fraction of biofilm obtained from image segmentation for the Δpsl (red line) and Δpel strains (green line) compared to the wild-type strain (blue line), for (Q = 0.2 μL/min, 2 μL/min).
Sloughing can be understood as a stochastic process
The model presented so far can describe Stage I and the signal in between two sloughing events, but not the catastrophic jumps associated with sloughing. To characterize the latter, we first performed a frequency analysis of the signal, as shown in Fig 9a. This approach did not prove very informative as it essentially shows a power spectrum typical of noise, with a slope roughly smaller than −2. This result, however, motivated the construction of a more physical representation of the process. The basis of this representation is the observation that fluctuations have a very specific signature on the hydraulic resistance: they first feature a sharp decrease due to sudden sloughing, followed by a slower increase due to growth. Since growth and smooth biofilm detachment are already described in the model, this observation suggests that we only need to capture the sudden sloughing to improve the description. The model takes the form of the following stochastic differential equation,
![](https://prod--epp.elifesciences.org/iiif/2/98292%2Fv1%2Fcontent%2F597893v1_fig9.tif/full/max/0/default.jpg)
Stage II detachment events can be described as a jump process. (a) Power spectrum for the different flow rates (Q = 0.2 μL/min, 2 μL/min and 20 μL/min). Slopes are indicative and were calculated in the interval between 0 and 0.25 Hz. (b) Probability density function of the time between two successive jump events, δt, for (Q = 0.2 μL/min, 2 μL/min and 20 μL/min). The histograms are calculated from experiments, while the solid lines are fitted Gamma distributions (δt)a−1/(baΓ(a)) exp (−δt/b). (c) Probability density function of the relative amplitude of jump events,
with ϕt describing the average fraction of biofilm in the microchannel, N the jumps and
The challenging part of this representation is to accurately describe the randomness of the process N – which is likely strongly dependent upon spatial heterogeneities, for instance in the initial attachment of bacteria. We characterized N in our experiments via the distributions of both the times between two successive jumps and the relative amplitudes of the jumps – i.e. the amplitude of each jump divided by the value of the volume fraction just before the jump ϕt−. Fig 9 shows histograms of these distributions extracted from signal processing of the experimental data (see Material and Methods). We see that the distribution of times between successive jumps, δt, can be well approximated by a Gamma distribution – solid lines in Fig 9 for a fit of the experimental data. The relative amplitude of each jump
Discussion
Regimes of nutrients and oxygen transport
Our experimental setup, which implements a novel UVC strategy to confine biofilm development within a specific zone in the channels, allowed us to finely control boundary conditions and study the role of molecular transport on biofilm development. Our results suggest that, in confined channels, the longitudinal distribution of biofilm ensues from a competition between advective transport of solutes along the channel and uptake by the cells. Modeling molecular transport as an advection-diffusion-reaction equation, this effect can be quantified using the ratio of longitudinal Péclet and Damköhler numbers, Pe/Da, which compares advection to reaction. When Pe/Da < 1, biofilm growth is nutrient-limited and localized close to the inlet. When Pe/Da > 1, biofilm development is close to homogeneous. In the radial direction, the limitation can be characterized in a similar way via a radial Damköhler number, Daradial, defined as the ratio of reaction and radial diffusion. In our experiments, we always had Daradial < 1 with only weak gradients and no radial limitation.
Although the competition between advective/diffusive transport and uptake is a universal feature of biofilm development in flow, the solute that becomes limiting depends on the details of the setup. In our experiments, we found that nutrients were restricting growth, while other works in microchannel flows, such as (Thomen et al., 2017), found that the effect of oxygen was dominant. Several factors play an important role here, including the composition of the growth medium, the catabolism of the microorganisms, the kinetics or even the resilience of bacteria to low oxygen concentration. For instance, Thomen et al. (2017) focused on an E. coli K12 derivative. E. coli can respond strongly to micro-aerobiosis (Colón-González et al., 2004), whereas P. aeruginosa is known to be more resilient (Sabra et al., 2002; Alvarez-Ortega and Harwood, 2007). Factors that affect transport characteristic times, such as the flow velocity or the dimensions of the channels, also modify the competition with uptake. The length of the channel modulates the time for advective transport along the channel and therefore the ratio Pe/Da. The width and height control the characteristic time for radial diffusion. In Thomen et al. (2017), for instance, the width of the channel is 1000 μm so that diffusion takes about 100 times longer and Daradial is 100 times larger than for 100 μm channels.
Temporal dynamics and hydrodynamic stresses
With these regimes clearly identified for molecular transport, we could then study the effect of hydrodynamic stresses on biofilm development and, in particular, the role of flow-induced detachment on the spatio-temporal dynamics. We suppressed any form of nutrient limitation, either radial or longitudinal, by focusing on flow rates ≥ 0.2 μL/min with Pe/Da > 1 and Daradial < 1. In so doing, we found that the time response featured two distinct stages: the hydraulic resistance first followed a relatively standard sigmoid-like pattern (Stage I), but then switched to a different regime with large self-sustained fluctuations (Stage II). As discussed in what follows, we argue that this Stage II is a signature effect of biofilm development in confined flow.
Stage I began just after initial inoculation, as soon as the culture medium started flowing in the microchannel. Early Stage I was dominated by cell adhesion, division, surface motility and microcolony formation. The average velocity in the channel at 0.2 μL/min is above 300 μm/s so bulk motility is likely not playing an important role when Q ≥ 0.2 μL/min (Rossy et al., 2019). Contrary to Thomen et al. (2017) – where no direct initiation of the biofilm of E. coli biofilm was possible above a shear stress of 10 mP a – our results show that P. aeruginosa PAO1 strongly adhered to the surface and could initiate biofilm formation for all shear stresses. We even observed that the apparent doubling time within the first few hours of growth was the smallest in the case at 2 P a. Stage I then evolved towards the standard steps of biofilm formation for P. aeruginosa (Rasamiravaka et al., 2015; Krsmanovic et al., 2021), with clonal microcolonies that extended, joined together and grew into mature biofilms. The exponential growth and EPS production yielded a sharp increase of the hydraulic resistance by several orders of magnitude, which is consistent with previous works (Cunningham et al., 1991; Desmond et al., 2022).
Stage II was characterized by large fluctuations in the hydraulic resistance and in the quantity of biomass in the channel. We showed that these fluctuations are the result of successive cycles of growth and sloughing events: parts of the biofilm get detached by the flow upon reaching a critical stress, generating sharp discrete-like events of detachment, followed by a fast re-growth of the biofilm. In both late Stage I and in Stage II, pressure-induced stresses play a crucial role. Most works have focused on shear stress (Krsmanovic et al., 2021; Kurz et al., 2022; Telgmann et al., 2004; Wagner et al., 2010; Paul et al., 2012) because, in many situations, such as biofilm development in large reactors, the effect of pressure can be neglected. However, this is not the case in confined systems, where biofilm can clog a large proportion of the flow channels and the pressure build-up induced by clogging can generate an important force on the biofilm.
Maximum clogging
One aspect that stands out in our experiments, because of its physical meaning and reproducibility across replicates, is the concept of maximum clogging. In Stage I and after each sloughing in Stage II, curves for the hydraulic resistance had a sigmoid-like shape, with an inflection that was a signature of a competition between growth and smooth detachment due to erosion, seeding and biofilm flow outside the zone of interest. We found that the maximum value of the hydraulic resistance across the entire experiment scales with the inverse of the flow rate, which is consistent with the modeling hypothesis that smooth detachment scales with the square root of the hydrodynamic stress (Rittman, 1982; Duddu et al., 2009; Coyte et al., 2017). This concept of maximum clogging can also be expressed as a critical value of the hydrodynamic stress generating enough detachment to compensate growth. In our framework, this equilibrium stress was in the order of 100-200 Pascals for the wild type and did not vary significantly with the flow rate.
To further understand the role of the EPS matrix composition on maximum clogging, we studied the behavior of Δpsl and Δpel mutants. Pel and Psl are the two main exopolysacharides produced by nonmucoid strains of P. aeruginosa. For both mutants, the maximum amount of biofilm in the microchannel was lower than for the wild type, thus indicating that these mutants were more easily detached by the flow. We also found that Psl had a stronger impact than Pel on the spatiotemporal dynamics and, in particular, on maximum clogging. The Δpsl mutant featured major sloughing events for 2 μL/min that led to a near-complete removal of the biofilm. For 20 μL/min, we observed a complete removal. On the other hand, Pel featured a dynamics similar to that of the wild type, but with a reduced maximum clogging.
Although critical stresses are widely used in characterizing flow-induced detachment of biofilms, they have a broad range of different meanings and values. Critical stresses have been used to characterize a form of complete detachment from the solid surface (Jiang et al., 2021). Ohashi and Harada (1996) described an adhesion strength as the stress required to remove the biofilm – upon conceptualizing the problem as a surface to surface bonding, the adhesive strength has also be defined as the work required to detach the biofilm from the surface per surface area (in Joule per m2) (Chen et al., 1998). Ohashi and Harada (1996) found values for the shear strength in the hundreds of Pascals. Lau et al. (2009) measured an adhesive pressure – adhesive force measured by microbead force spectroscopy divided by contact area – for P. aeruginosa PAO1 that was in the tens of Pascals. Körstgens et al. (2001) studied the yield strength of a mucoid P. aeruginosa strain, discussing its role in the mechanical failure of the biofilm. Considering the biofilm as a viscoelastic gel with plastic flow properties, they found a stress at failure close to 1000 P a. Lee et al. (2023) show that bio-aggregates of E. coli at pore throats become fluidized above a critical value of the shear stress with a yield point at roughly 1.8 P a.
To establish a clear link between these values, our critical stress at 100-200 Pascals and the observed dynamics of sloughing in Stage II, future works will need to further refine our understanding of the role of the different variables on flow-induced detachment, including distributions of stresses within biofilms, adhesive strength to surfaces, mechanical properties of the biofilm, and biological regulation through e.g. quorum sensing (Emerenini et al., 21 juil. 2015) and/or rhamnolipids production (Boles et al., 2005).
Self-sustained fluctuations
Large fluctuations of the hydraulic resistance, similar to that observed in our Stage II, were previously reported in a few works on porous media flows. Howell and Atkinson (1976) showed that, in trickling filters, constant influent concentration and operating conditions can yield large fluctuations induced by sloughing events. Stewart and Scott Fogler (2002) found very large oscillations in the pressure signal due to the formation of successive stable or unstable dextran plugs of Leuconostoc mesenteroides biofilm, behaving in a way similar to a yield stress fluid, and allowing flow only through breakthrough channels – see also discussions on the flow of yield stress fluids through porous media in (Talon, 2022). Sharp et al. (2005) showed that “the pore channels are dynamic, changing in sized, number and location with time” in a flat plate reactor colonized by Vibrio fischeri biofilm. Kurz et al. (2022) found successive cycles of growth and shear-induced detachment in the preferential flow paths for Bacillus subtilis biofilm development in a microfluidic system with cylindrical obstacles. Bottero et al. (2013) modeled the coupled effects of flow, clogging and detachment to study the mechanisms that control these self-sustained fluctuations, in particular what leads to the clogging of a preferential flow path and the unclogging of another one – contrary to the development of stable flow paths.
Although there are similarities with our observations, porous media have an inherent degree of complexity in the geometry that, combined with the nonlinear response of biofilms, can lead to a range of complex mechanisms. For instance, some of the works above-mentioned identified changes in the flow paths associated with the fluctuations. Our work shows that, even in a single channel with continuous flow, these fluctuations occur and prevent the system from reaching a true steady-state. This observations leads us to propose that the self-sustained fluctuations observed in Stage II are a signature of biofilm development in confined flows and may occur in the simplest geometries. The underlying idea is that the strong bioclogging in confined geometries generates a pressure build-up that is the main cause of sloughing. Once a large part of the biofilm gets detached, the pressure decreases and the biofilm can re-grow to get a new cycle started. Our work therefore highlights the ubiquity of these fluctuations, which are likely to play a key role in the development and spreading of biofilms in a range of different systems, such as clogging in pipes, catheters (Stickler, 2014) or stents (Guaglianone et al., 2010), with applications ranging from infections, to environmental processes and engineering systems.
Stochastic modeling
Our approach to modeling was constructed in two main parts. We first proposed differential equations to describe nutrient transport coupled with the dynamics of biofilm development in the channel, including growth and smooth detachment. The sharp jumps corresponding to sloughing were then added to this first layer of the model as a jump stochastic process. Although it seems quite natural to characterize the statistics of detachment (Wilson et al., 2004), there are actually very few attempts to treat such problems in the framework of stochastic processes. Howell and Atkinson (1976) modeled sloughing in trickling filters through a conceptualization as connected filter units describing the discrete pieces of packing material in the filter. They introduced randomness by permitting sloughing at multiples of a fixed time interval. Bohn et al. (2007) used an approach combining a logistic growth with random sloughing events through a stochastic differential equation. They described sloughing via a discrete expression of the amplitude of jumps occurring independently at each time step, which allowed them do describe daily fluctuations in light absorbance data for phototrophic biofilms development in a flow-lane incubator. In our approach, the jump process was characterized by two random variables: the interevent time between two successive jumps and the relative amplitudes of the jumps. Through the analysis of experimental data, we reconstructed the probability density functions for these random variables, with respectively Gamma and log-normal distributions.
There have been previous discussions (Lewandowski et al., 2004) on the determinism of biofilm formation, suggesting that sloughing events are intrinsically random events that generate large fluctuations, prevent the system from reaching a steady state and hinder the reproducibility of longterm experiments. Our work confirms that sloughing is integral to biofilm development (Telgmann et al., 2004) in confined flows but shows that, although a true steady state is never reached, the fluctuations can be precisely characterized using stochastic modeling. This approach paves a way forward in terms of reproducibility: even though the state of the biofilm at any given time may not be reproducible, the randomness of the process may very well be.
Our approach to characterizing bursting events in terms of the distribution of the amplitude and interevent time is reminiscent of the description of other physical systems, such as avalanches (Maaß et al., 2015). For earthquakes, for example, a Gamma distribution for interevent times has been found in many different geographic regions (Corral, 2004). The fact that interevent times for biofilm sloughing also seem to follow a Gamma distribution may point towards specific physical mechanisms (Kumar et al., 2020), which could be used in future works to better understand the physics of sloughing. One interesting perspective of this work is also to assess the universality of these distributions across microorganisms and whether we could define classes of bacteria with specific signatures on the stochastic process. With the same idea, we could further evaluate the impact of ecological interactions in multispecies biofilm or the effect of various molecules, such as biocides, on the sloughing dynamics.
Material and Methods
Bacteria and cultures
Experiments were performed using Pseudomonas aeruginosa PAO1 GFP (ATCC 15692GFP) strain, along with PAO1 GFP ΔpslD and ΔpelF mutants obtained from Colvin et al. (2011, 2012) and built by non-polar deletion through allelic replacement of pslD and pelF operon and harbor pMRP9-1 plasmid expressing GFP. Bacteria were subcultured and grown in brain heart infusion (Sigma Aldrich, Saint-Quentin-Fallavier, France). Cultures were prepared from −80°C frozen aliquots spread on tryptic soy agar plate (Sigma Aldrich, Saint-Quentin-Fallavier, France) supplemented with 300 µg/mL of ampicillin (Sigma Aldrich) then incubated at 30°C during 24 h. Liquid cultures were prepared from the second subcluture on tryptic soy agar by spreading a 24 hours single colony diluted in brain heart infusion media supplemented with 300 µg/mL of ampicillin. 1X concentrated BHI media was prepared by dissolving 37 g of commercial powder in demineralized water and autoclaved with a liquid cycle (121° C for 15 minutes). 0.2Xconcentrated BHI media was obtained by dissolving 7.4 g while a third 1X solution was supplemented by 8 g/L of D-Glucose (Sigma Aldrich).
Microfabrication
Microchannels were fabricated using standard soft lithography techniques. Microchannel molds were prepared by depositing 100 µm SUEX sheets on a silicium wafer via photolithography. The negative mold was cleaned by isopropanol and silanized with trichloromethylsilane (Sigma Aldrich). Square cross-section channels had dimensions of 100 µm height by 100 µm width and 20 mm length. The chips were prepared with a 10% wt/wt cross-linking agent in the polydimethylsiloxane solution (Sylgard 184 Silicone Elastomer Kit, Dow Corning). PDMS was cleaned with isopropanol at 80 °C for 30 minutes and plasma-bonded to a clean glass coverslide.
Inoculation and flow experiments
BHI suspensions were adjusted at optical density at OD620nm= 0.2 (108 CFU/mL) and inoculated inside the microchannels from the outlet, up to approximately ¾ of the channel length in order to keep a clean inlet. The system was let at room temperature (25°C) for 3h under static conditions. Flow experiments were then performed at 0.02, 0.2, 2, 20 and 200 μL/min constant flow rates for 72h in the microchannels at room temperature. For the experiments at 0.2, 2, 20 and 200 μL/min, the fluidic system was based on a sterile culture medium reservoir pressurized by a pressure controller (Fluigent FlowEZ) and connected with a flow rate controller (Fluigent Flow unit). The flow rate was maintained constant by using a controller with a feedback loop adjusting the pressure in the liquid reservoir. The reservoir was connected to the chip using Tygon tubing (Saint Gobain Life Sciences Tygon™ ND 100-80) of 0.52 mm internal diameter and 1.52 mm external diameter, along with PEEK tubing (Cytiva Akta pure) with 0.25 mm inner diameter adapters for flow rate controller. The waste container was also pressurized by another independent pressure controller to reduce air bubble formation in the inlet part. For the experiments at 0.02 μL/min, we used an Harvard Phd2000 syringe pump for the flow.
UVC irradiation
The biofilm is constrained in a part of the microchannel using UVC irradiation directly through the PDMS of the microfluidic chips, thus reducing contamination risk and avoiding unwanted progression/growth of P. aeruginosa in the inlet and tubing for several days of experiment (Ramos et al., 2023). This eliminates parasitic consumption of nutrients in various parts of the fluidic system and maintain a controlled boundary condition with a fixed concentration of nutrients at the inlet of our zone of interest.
To this end, the inlet and outlet of the microchannels were exposed to UVC light by a system of UVC LEDs (Ramos et al., 2023). The system consists of a 3D printed part called the guide carrying in its backside a PCB with a LED light source that delivers UV-C light of 1W power. The light beam follows a straight trajectory until a 45° mirror positioned in the front side of the guide, which reflects the light parallel to the PDMS chip and irradiate it. UVC guides were positioned in both sides of the PDMS microchannel and separated by a 1.2 cm distance to keep a central area unexposed. UVC power is measured by a radiometer after mirror reflection and delivers a power of 200 mW /m2equivalent to 2 mJ /cm2. The guide is elevated to fit with the PDMS chip dimensions and contains a barrier located in the front side that blocks the diffusion of light through the PDMS polymer in the horizontal direction.
Mass spectrometry analysis
Mass spectrometry analysis was performed on a Thermo Scientific Q Exactive Focus Orbitrap LC-MS/MS system connected to a LC device Thermo Scientific Vanquish UPLC system with PDA detection. Analytical separations were performed on a 150 x 2.1 mm Thermo Hypersyl Gold C18 column (1.9 μm) using an MeCN/H2O 1% formic acid gradient. Data were captured in full MS scan mode and processed using Chromeleon 7 software.
Imaging for the biofilm experiment
Bacterial development was imaged for a period of 72 hours with a timestep of 30 minutes at 25°C on an inverted microscope (Ti-2E, Nikon) using a digital camera (back-illuminated PCO edge). Time-lapse images were acquired using brightfield and fluorescence microscopy (Sola light source 10% intensity with 30 ms exposure with 500 nm excitation and 513 nm emission combined with (FITC filter). Images were obtained with a focal plan at the glass/liquid interface. These images had dimensions of 30086 x 154 pixels obtained after multi position scanning using automatic Nikon platform and assembled by Nikon NIS software of single images with 0.65 µm/pixel using a 10X magnification Nikon objective (NA = 0.3).
Image analysis for the biofilm distribution in the longitudinal direction
Fluorescence images were loaded as a matrix (30086 x 154) in MATLAB (MathWorks). For each time acquisition, the signal was first integrated in the radial direction to obtain a mean distribution in the longitudinal direction. In plotting curves to analyze the effect of nutrients, all timepoints were then averaged to obtain one dimensional curves of the mean longitudinal profile. For kymographs, one dimensional curves for each time point were stacked together to describe the spatio-temporal dynamics. The intensity values were normalized to the maximum values of each replicate over all times. For each flow and nutrient condition, three biological replicates were performed.
Image analysis for the biofilm segmentation
To estimate changes in channel colonization, GFP images were binarized using a machine learning software (Ilastik) (Berg et al., 2019). Images were pre-treated with imageJ (Schneider et al., 2012) by normalizing all pixel values between 0 and 65535 grayscale levels. In Ilastik, biofilm structures were first differentiated from the empty flow path using pixel-level manual labeling during pixel classification where visible patches of biofilm and empty channel were annotated manually by mouse cursor. Pixel classification workflow employs a Random Forest classifier, known for its generalization properties. Several samples of biofilm and background images were used to train the classifier by annotating pixels with corresponding labels, allowing the algorithm to learn and make predictions in real-time. The chosen features were color, intensity, edges, and texture. The generated probability maps indicating the likelihood of each class at every pixel were used for the object classification and were thresholded at a value of 0.6 with no size filter. Thresholding is a process involved in converting continuous probability maps generated from pixel classification into binary segmentation images by setting a threshold, where pixels above the threshold are classified as belonging to an object. The size, intensity, position and convexity of the biofilm objects was exported in .cvs format and further analysed in matlab. Volumic fraction of biofilm in microchannel was calculated from the sum of the size of segmented objects divided by the interest growth area (non UVC irradiated central part of microchannel).
Initial adhesion
Separate experiments were performed to study the behavior of cells in the initial phases of attachment in order to increase spatial and temporal resolution. Liquid cultures were prepared following the same protocol as described previously (see Bacteria and cultures). Sterile 1× concentrated BHI culture medium supplemented by 300 µg/mL of ampicillin was flowed under constant flow rate (Q = 0.2 μL/min, 2 μL/min and 20 μL/min) after 3 hours under static conditions at 26° C.
Images were obtained with a 40× Nikon objective (NA = 0.95) using a differential interference contrast (DIC) brightfield with 2 minutes per frame during 3.5 Hours. Replicates were obtained by imaging four positions per channel and each condition was performed in two distinct channels forming two distinct biological replicates to obtain n = 8 replicates by condition (Q = 0.2 μL/min, 2 μL/min and 20 μL/min). The images were further segmented using Ilastik with the same parameters as described before and single cells were selected as objects with an area higher than 20 μm2 to avoid counting dust particles and artefacts that could be considered as distinct objects by Ilastik. The doubling time was calculated in the window between 0 and 3.5 hours by a a linear fitting of the logarithm of the number of cells. The slope was used to estimate growth rate and doubling time. Cell count was calculated from image segmentation of four positions in two channels to generate 8 replicates by condition (n= 8) for (Q = 0.2 μL/min, 2 μL/min and 20 μL/min).
Wavelet analysis
Time series of the pressure data were investigated with a wavelet analysis to identify temporal variations of spectral power (Torrence and Compo, 1998). Wavelet analysis was carried out using a Morlet wavelet, the product of a sinusoidal wave and a Gaussian envelope, with a frequency parameter of 6 and scale width of 300. We then applied a continuous wavelet transform (CWT) results which gives wave power coefficients dependent on the scale or period and the time, as well as a cone of influence (COI), where edge effects become important. The time series were zeropadded to reduce the edge errors. The cone of influence was removed from the scalograms in Fig 3.
Data analysis for the calculation of the hydraulic resistance and volume fraction
Recorded pressure fluctuations in the reservoir were converted to hydraulic resistance in the 10 mm zone between UVC LEDs where biofilm develops. We write the difference between the two pressure reservoirs (one at the inlet and one at the outlet) as ΔPres and express it as
with Q the imposed flow rate,∑i Ri the sum of the hydraulic resistance of the different part of the hydraulic network and R (t) the hydraulic resistance of the the 10 mm zone between UVC LEDs. To estimate R(t), we evaluated ∑i Ri from the mean value of ΔPres, which we write
which yields
with R (t = 0) in the square cross-section of size h0 = 100 μm that can be estimated as
with η the viscosity of the culture medium at 25°C, L = 10 mm the length of the channel and h the width/height of the channel. Since we considered
with b the thickness of the biofilm layer. The volume fraction of biofilm is
Frequency analysis
Welch power spectral analysis (Welch, 1967) of the pressure signal was carried out by dividing the time signal into smaller segments, calculating their periodogram and finally averaging accross the frequencies, resulting in a power spectral density (PSD) estimate. This method allows for PSD estimate that is less noisy than usual periodograms. Here, we used the Matlab tool pwelch to estimate the PSD using hanning window with an overlap of 50%. A linear fit was then applied on the interval in the low frequency range.
Construction of the probability density functions for jumps and fits
Volume fraction data were processed using a homemade Matlab code. Each dataset was first sub-sampled by keeping only 1 in 150 points. The subsampled signal was then differentiated and only negative values corresponding to detachment events were conserved. Jumps were then identified through the selection of local maxima. For each flow rate, all the values for the times between two successive jumps and the relative amplitude of the jumps – amplitude of the jump relative to the value of the volume fraction just before the jump – for the different replicates were aggregated to construct the distributions.
Gamma distributions were then fitted to the experimental data for the times between successive jumps and lognormal distributions were used for the amplitude of the jumps. Values for the different parameters are summarized in Table.2.
Numerical simulation of the stochastic process
Here we describe how we simulated the evolution of the volume fraction ϕ, solving
The simulation was done using a homemade Matlab code. For each case, based on the previously calculated distributions, we first generated two sets of random numbers corresponding to the times between successive jumps and to the relative amplitude of the jumps – therefore allowing us to completely determine the jump process, N. We then simply solved the ordinary differential equation
between jumps using the Matlab solve ode45. Upon reaching a jump, the simulation was stopped and the jump implemented, before proceeding to treating the next interval.
Comsol flow simulations
Flow simulation was performed using a finite element approach in Comsol v6.1. The channel was treated as a rectangle of 850 μm length by 100 μm width. Biofilm, as shown in white in Fig 6, was added from a segmentation of real experimental images. We then solved incompressible Stokes flow with no-slip/no-penetration boundary conditions on the solid and on the biofilm surface. The inlet condition was an imposed velocity, while the outlet was an imposed pressure.
Acknowledgements
This work is part of a project that has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement no 803074). This work was partially supported by the LAAS-CNRS micro and nanotechnologies platform, a member of the French Renatech network. We thank the entire BEBOP team for daily interactions and assistance; Julien Lefort and Emmanuel Libert for their technical support; Terence Desclaux for the introduction to wavelet analysis; Matthew R. Parsek and Maria Zori for providing the Pel and Psl mutants; and Laure Latapie for her support with mass spectrometry.
References
- Responses of Pseudomonas Aeruginosa to Low Oxygen Indicate That Growth in the Cystic Fibrosis Lung Is by Aerobic RespirationMolecular microbiology 65:153–165https://doi.org/10.1111/j.1365-2958.2007.05772.x
- A Network Biology Approach to Denitrification in Pseudomonas AeruginosaPLoS ONE 10https://doi.org/10.1371/journal.pone.0118235
- Recent Developments in Biofouling Control in Membrane Bioreactors for Domestic Wastewater TreatmentSeparation and Purification Technology 206:297–315https://doi.org/10.1016/j.seppur.2018.06.004
- The Reduction Potential of Nitric Oxide (NO) and Its Importance to NO BiochemistryProceedings of the National Academy of Sciences 99:10958–10963https://doi.org/10.1073/pnas.162095599
- The Ecology and Biogeochemistry of Stream BiofilmsNature Reviews Microbiology 14:251–263https://doi.org/10.1038/nrmicro.2016.15
- Ilastik: Interactive Machine Learning for (Bio)Image AnalysisNature Methods 16:1226–1232https://doi.org/10.1038/s41592-019-0582-9
- Biophysical Con-trols on Community Succession in Stream BiofilmsApplied and Environmental Microbiology 73:4966–4974https://doi.org/10.1128/AEM.00588-07
- Review Article: Biofouling: Lessons from NaturePhilosophical transactions Series A, Mathematical, physical, and engineering sciences 370:2381–417https://doi.org/10.1098/rsta.2011.0502
- Stochastic modeling for characterisation of biofilm development with discrete detachment events (sloughing)Water Science and Technology: A Journal of the International Association on Water Pollution Research 55:257–264https://doi.org/10.2166/wst.2007.266
- Rhamnolipids mediate detachment of Pseudomonas aeruginosa from biofilmsMolecular Microbiology 57:1210–1223https://doi.org/10.1111/j.1365-2958.2005.04743.x
- Biofilm Development and the Dynamics of Preferential Flow Paths in Porous MediaBiofouling 29:1069–1086https://doi.org/10.1080/08927014.2013.828284
- Theoretical Microfluidics
- Microbial Adhesion in Flow Displacement SystemsClinical Microbiology Reviews 19:127–141https://doi.org/10.1128/CMR.19.1.127-141.2006
- Dynamics of Biofilm Processes: MethodsWater Research 16:1207–1216https://doi.org/10.1016/0043-1354(82)90139-7
- Direct measurement of the adhesive strength of biofilms in pipes by micromanipulationBiotechnology Techniques 12:875–880https://doi.org/10.1023/A:1008805326385
- Anaerobic Growth Does Not Support Biofilm Formation in Escherichia Coli K-12Research in Microbiology 155:514–521https://doi.org/10.1016/j.resmic.2004.03.004
- The Pel Polysaccharide Can Serve a Structural and Protective Role in the Biofilm Matrix of Pseudomonas AeruginosaPLoS Pathogens 7https://doi.org/10.1371/journal.ppat.1001264
- The Pel and Psl Polysaccharides Provide Pseudomonas Aeruginosa Structural Redundancy within the Biofilm MatrixEnvironmental microbiology 14https://doi.org/10.1111/j.1462-2920.2011.02657.x
- Confined Flow: Consequences and Implications for Bacteria and BiofilmsAnnual Review of Chemical and Biomolecular Engineering 9:175–200https://doi.org/10.1146/annurev-chembioeng-060817-084006
- Long-Term Clustering, Scaling, and Universality in the Temporal Occurrence of EarthquakesPhysical Review Letters 92https://doi.org/10.1103/PhysRevLett.92.108501
- Microbial BiofilmsAnnual Review of Microbiology 49:711–745https://doi.org/10.1146/annurev.mi.49.100195.003431
- Biofilm Development in Water Distribution and Drainage Systems: Dynamics and Implications for Hydraulic EfficiencyEnvironmental Technology Reviews 3:31–47https://doi.org/10.1080/09593330.2014.923517
- Microbial Competition in Porous Environments Can Select against Rapid Biofilm GrowthProceedings of the National Academy of Sciences 114:E161–E170https://doi.org/10.1073/pnas.1525228113
- Influence of biofilm accumulation on porous media hydrodynamicsEnvironmental Science & Technology :1305–1311https://doi.org/10.1021/es00019a013
- Controlling the hydraulic resistance of membrane biofilms by engineering biofilm physical structureWater Research 210https://doi.org/10.1016/j.watres.2021.118031
- Biofilm Streamers Cause Catastrophic Disruption of Flow with Consequences for Environmental and Medical SystemsProceedings of the National Academy of Sciences 110:4345–4350https://doi.org/10.1073/pnas.1300321110
- A Two-dimensional Continuum Model of Biofilm Growth Incorporating Fluid Flow and Shear Stress Based DetachmentBiotechnology and Bioengineering 103:92–104https://doi.org/10.1002/bit.22233
- Mechanomicrobiology: How Bacteria Sense and Respond to ForcesNature Reviews Microbiology 18:227–240https://doi.org/10.1038/s41579-019-0314-2
- A Mathematical Model of Quorum Sensing Induced Biofilm DetachmentPLOS ONE 10https://doi.org/10.1371/journal.pone.0132385
- Resilience of Bacterial Quorum Sensing against Fluid FlowScientific Reports 6https://doi.org/10.1038/srep33115
- The Biofilm MatrixNature Reviews Microbiology 8:623–633https://doi.org/10.1038/nrmicro2415
- Biofilms: An Emergent Form of Bacterial LifeNature Reviews Microbiology 14:563–575https://doi.org/10.1038/nrmicro.2016.94
- Physiology of Pseudomonas Aeruginosa in Biofilms as Revealed by Transcriptome AnalysisBMC Microbiology 10https://doi.org/10.1186/1471-2180-10-294
- Fouling in Heat ExchangersInverse Heat Conduction and Heat Exchangers https://doi.org/10.5772/intechopen.88079
- Biofilm Mechanics: Implications in Infection and SurvivalBiofilm 2https://doi.org/10.1016/j.bioflm.2019.100017
- Microbial Biofilms Associated with Biliary Stent CloggingFEMS Immunology & Medical Microbiology 59:410–420https://doi.org/10.1111/j.1574-695X.2010.00686.x
- Sloughing of microbial film in trickling filtersWater Research 10:307–315https://doi.org/10.1016/0043-1354(76)90172-X
- Searching for the Secret of Stickiness: How Biofilms Adhere to SurfacesFrontiers in Microbiology 12https://doi.org/10.3389/fmicb.2021.686793
- Biofilm Dispersal: Mechanisms, Clinical Implications, and Potential Therapeutic UsesJournal of Dental Research 89:205–218https://doi.org/10.1177/0022034509359403
- Local and Global Consequences of Flow on Bacterial Quorum SensingNature Microbiology 1https://doi.org/10.1038/nmicrobiol.2015.5
- Uniaxial compression measurement device for investigation of the mechanical stability of biofilmsJournal of Microbiological Methods 46:9–17https://doi.org/10.1016/s0167-7012(01)00248-2
- Hydrodynamics and Surface Properties Influence Biofilm ProliferationAdvances in Colloid and Interface Science 288https://doi.org/10.1016/j.cis.2020.102336
- Novel Approaches to Microbial Enhancement of Oil RecoveryJournal of Biotechnology 266:118–123https://doi.org/10.1016/j.jbiotec.2017.12.019
- On Interevent Time Distri-butions of Avalanche DynamicsScientific Reports 10https://doi.org/10.1038/s41598-019-56764-6
- Competition between Growth and Shear Stress Drives Intermittency in Preferential Flow Paths in Porous Medium BiofilmsProceedings of the National Academy of Sciences 119https://doi.org/10.1073/pnas.2122202119
- Absolute Quantitation of Bacterial Biofilm Adhesion and Viscoelasticity by Microbead Force SpectroscopyBiophysical Journal 96:2935–2948https://doi.org/10.1016/j.bpj.2008.12.3943
- Rapid Formation of Bioaggregates and Morphology Transition to Biofilm Streamers Induced by Pore-Throat FlowsProceedings of the National Academy of Sciences 120https://doi.org/10.1073/pnas.2204466120
- Biofilms as BiobarriersThe American Biology Teacher 71:20–26https://doi.org/10.2307/27669358
- Reproducibility of biofilm processes and the meaning of steady state in biofilm reactorsWater Science and Technology: A Journal of the International Association on Water Pollution Research 49:359–364
- Slip Statistics of Dislocation Avalanches under Different Loading ModesPhysical Review E 91https://doi.org/10.1103/PhysRevE.91.042403
- Bacterial Quorum Sensing in Complex and Dynamically Changing EnvironmentsNature Reviews Microbiology 17:371–382https://doi.org/10.1038/s41579-019-0186-5
- A novel concept for evaluation of biofilm adhesion strength by applying tensile force and shear forceWater Science and Technology 34:201–211https://doi.org/10.1016/0273-1223(96)00647-6
- Effect of Shear Stress and Growth Conditions on Detachment and Physical Properties of BiofilmsWater Research 46:5499–5508https://doi.org/10.1016/j.watres.2012.07.029
- Viscoelasticity of Biofilms and Their Recalcitrance to Mechanical and Chemical ChallengesFEMS microbiology reviews 39:234–245https://doi.org/10.1093/femsre/fuu008
- Effect of Diffusive and Convective Substrate Transport on Biofilm Structure Formation: A Two-Dimensional Modeling StudyBiotechnology and Bioengineering 69:504–515https://doi.org/10.1002/1097-0290(20000905)69:5<504::aid-bit5>3.0.co;2-s
- Ultraviolet Control of Bacterial Biofilms in Microfluidic ChipsBiomicrofluidics 17https://doi.org/10.1063/5.0135722
- The Formation of Biofilms by Pseudomonas Aeruginosa : A Review of the Natural and Synthetic Compounds Interfering with Control MechanismsBioMed Research International 2015:1–17https://doi.org/10.1155/2015/759348
- The Effect of Shear Stress on Biofilm Loss RateBiotechnology and Bioengineering 24:501–506https://doi.org/10.1002/bit.260240219
- The significance of biofilms in porous mediaWater Resources Research 29:2195–2202https://doi.org/10.1029/93WR00611
- Mechanosensing of Shear by Pseudomonas Aeruginosa Leads to Increased Levels of the Cyclic-Di-GMP Signal Initiating Biofilm DevelopmentProceedings of the National Academy of Sciences 114:5906–5911https://doi.org/10.1073/pnas.1703255114
- Cellular Advective-Diffusion Drives the Emergence of Bacterial Surface Colonization Patterns and HeterogeneityNature Communications 10https://doi.org/10.1038/s41467-019-10469-6
- Secondary Flow as a Mechanism for the Formation of Biofilm StreamersBiophysical Journal 100:1392–1399https://doi.org/10.1016/j.bpj.2011.01.065
- Role of Polysaccharides in Pseudomonas Aeruginosa Biofilm DevelopmentCurrent Opinion in Microbiology 10:644–648https://doi.org/10.1016/j.mib.2007.09.010
- Physiological Responses of Pseudomonas Aeruginosa PAO1 to Oxidative Stress in Controlled Microaerobic and Aerobic CulturesMicrobiology (Reading, England) 148:3195–3202https://doi.org/10.1099/00221287-148-10-3195
- Pseudomonas ChemotaxisFEMS Microbiology Reviews https://doi.org/10.1111/1574-6976.12081
- Microfluidic-Based Transcriptomics Reveal Force-Independent Bacterial RheosensingNature Microbiology 4:1274–1281https://doi.org/10.1038/s41564-019-0455-0
- NIH Image to ImageJ: 25 Years of Image AnalysisNature methods 9:671–675
- Biotechnology in Petroleum Recovery: The Microbial EORProgress in Energy and Combustion Science 34:714–724https://doi.org/10.1016/j.pecs.2008.05.001
- Visualization and Characterization of Dynamic Patterns of Flow, Growth and Activity of Biofilms Growing in Porous MediaWater Science and Technology 52:85–90https://doi.org/10.2166/wst.2005.0185
- Biofilms: Implications in BioremediationTrends in Microbiology 14:389–397https://doi.org/10.1016/j.tim.2006.07.001
- Biomass Plug Development and Propagation in Porous MediaBiotechnology and Bioengineering 72:353–363https://doi.org/10.1002/1097-0290(20010205)72:3<353::AID-BIT13>3.0.CO;2-U
- Pore-scale Investigation of Biomass Plug Development and Propagation in Porous MediaBiotechnology and Bioengineering 77:577–588https://doi.org/10.1002/bit.10044
- Clinical Complications of Urinary Catheters Caused by Crystalline Biofilms: Something Needs to Be DoneJournal of Internal Medicine 276:120–129https://doi.org/10.1111/joim.12220
- Influence of Hydrodynamics and Nutrients on Biofilm StructureJournal of Applied Microbiology 85:19–28https://doi.org/10.1111/j.1365-2672.1998.tb05279.x
- Mass Transfer Enhancement in Moving Biofilm StructuresBiophysical Journal 102:1483–1492https://doi.org/10.1016/j.bpj.2012.02.033
- On the Determination of a Generalized Darcy Equation for Yield Stress Fluid in Porous MediaScience Talks 3https://doi.org/10.1016/j.sctalk.2022.100042
- Biofilm growth and the related changes in the physical properties of a porous medium: 1. Experimental investigationWater Resources Research 26:2153–2159https://doi.org/10.1029/WR026i009p02153
- Biofilm growth and the related changes in the physical properties of a porous medium: 3Dispersivity and model verification. Water Resources Research 26:2171–2180https://doi.org/10.1029/WR026i009p02171
- Substrate and biomass transport in a porous mediumWater Resources Research 26:2181–2194https://doi.org/10.1029/WR026i009p02181
- Biofilm growth and the related changes in the physical properties of a porous medium: 2Permeability. Water Resources Research 26:2161–2169https://doi.org/10.1029/WR026i009p02161
- Influence of Growth History on Sloughing and Erosion from BiofilmsWater Research 38:3671–3684https://doi.org/10.1016/j.watres.2004.05.020
- Bacterial Biofilm under Flow: First a Physical Struggle to Stay, Then a Matter of BreathingPLoS ONE 12https://doi.org/10.1371/jour-nal.pone.0175197
- A Practical Guide to Wavelet AnalysisBulletin of the American Meteorological Society 79:61–78https://doi.org/10.1175/1520-0477(1998)079<0061:APGTWA>2.0.CO;2
- Effect of bacterial extracellular polymers on the saturated hydraulic conductivity of sand columnsApplied and Environmental Microbiology 58:1690–1698https://doi.org/10.1128/aem.58.5.1690-1698.1992
- Saturated Hydraulic Conductivity Reduction Caused by Aerobic Bacteria in Sand ColumnsSoil Science Society of America Journal 56:1–13https://doi.org/10.2136/ss-saj1992.03615995005600010001x
- Biofouling in Food Processing: Biofilm or Biotransfer PotentialFood and Bioproducts Processing 80:292–298https://doi.org/10.1205/096030802321154808
- Online Assessment of Biofilm Development, Sloughing and Forced Detachment in Tube Reactor by Means of Magnetic Resonance MicroscopyBiotechnology and Bioengineering 107:172–181https://doi.org/10.1002/bit.22784
- Dynamic Changes in Biofilm Structures under Dynamic Flow ConditionsApplied and Environmental Microbiology 88:e01072–22https://doi.org/10.1128/aem.01072-22
- The Use of Fast Fourier Transform for the Estimation of Power Spectra: A Method Based on Time Averaging over Short, Modified PeriodogramsIEEE Transactions on Audio and Electroacoustics 15:70–73https://doi.org/10.1109/TAU.1967.1161901
- Statistical Quantification of Detachment Rates and Size Distributions of Cell Clumps from Wild-Type (PAO1) and Cell Signaling Mutant (JP1) Pseudomonas Aeruginosa BiofilmsApplied and Environmental Microbiology 70:5847–5852https://doi.org/10.1128/AEM.70.10.5847-5852.2004
Article and author information
Version history
- Sent for peer review:
- Preprint posted:
- Reviewed Preprint version 1:
Copyright
© 2024, Benbelkacem et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.