2 Citations   Views   Downloads

Numerical investigation of microbial quorum sensing under various flow conditions

View article
Microbiology

Introduction

Microorganisms preferentially reside on solid surfaces, which often leads to a closer proximity of neighboring cells than when in a planktonic form (Donné & Dewilde, 2015). At elevated cell densities, microorganisms need to efficiently coordinate the expression of energetically expensive phenotypes, such as biofilm development, exoenzyme production and microbial dispersal. Efficiency is achieved by producing and detecting relatively cheap signaling molecules which regulate the phenotype expression only when a sufficient signal concentration has been reached (Miller & Bassler, 2001). This microbial decision-making process called “quorum sensing (QS)” was originally understood as a cell-to-cell communication to identify conspecific population density and accomplish cooperative behaviors (Fuqua, Winans & Greenberg, 1994). However, a number of studies have indicated that QS is not necessarily a social trait (Redfield, 2002; Carnes et al., 2010) and depends not only on the population but also on the spatial distribution of microbial cells (Alberghini et al., 2009; Gao et al., 2016). These observations led to an alternative QS concept in which QS depends strictly on the local concentration of signaling molecules (Hense et al., 2007; Hense & Schuster, 2015). This suggests that, to understand QS processes, an integrative approach is required analyzing a multitude of factors including microbial density (Fuqua, Winans & Greenberg, 1994), production and decay kinetics (Lee et al., 2002; Fekete et al., 2010), and transport of signaling molecules through advection and diffusion (Redfield, 2002), as well as the spatial distribution of microorganisms (Alberghini et al., 2009). Thus, spatial constraints and responses may be as important as other biological considerations for the evolution and maintenance of QS. This idea is known to be true in biofilms where cooperative strategies are able to evolve if cooperators are spatially aggregated (Xavier & Foster, 2007).

Individual microbial cells synthesize and release signaling molecules at a basal rate. At low population densities, the concentration of signaling molecules remains low as it degrades both biotically and abiotically (Lee et al., 2002; Yates et al., 2002). At a sufficiently high microbial population density, however, the extracellular concentration of signaling molecules reaches a threshold concentration that activates gene and phenotypes expression (Hense & Schuster, 2015). When QS regulates the production of costly public goods, this balances production cost and the overall benefit (Pai, Tanouchi & You, 2012; Heilmann, Krishna & Kerr, 2015; Schluter et al., 2016), while under nutrient limited conditions, QS can regulate microbial dispersal (Solano, Echeverz & Lasa, 2014; Boyle et al., 2015), improving chances of survival. QS induction also often upregulates genes controlling production of signaling molecules resulting in enhanced signal production (Ward et al., 2001; Fekete et al., 2010; Pérez-Velázquez et al., 2015). Such autoinduction has been thought to confer evolutionary stability and fitness advantages (Brandman et al., 2005; Mitrophanov, Hadley & Groisman, 2010; Gao & Stock, 2018), but its effects on neighboring microbial aggregates and evolutionary benefits in a spatial context have not been fully understood.

QS induction is affected by mass transport characteristics controlling the spatial distribution of signaling molecules. In a confined space, even a single microbial cell can be QS induced if the signaling molecules accumulate to sufficiently high concentration (Carnes et al., 2010). However, higher population densities are required for QS induction in a large open space because the signaling molecules are diluted due to diffusive loss to the surrounding medium (Alberghini et al., 2009; Trovato et al., 2014). Advection may dilute the signaling molecules more effectively than diffusion and repress QS induction. Experimental observations have shown that fast advective flows increase the amount of biomass required for QS induction (Kirisits et al., 2007) and repress QS dependent gene expression (Meyer et al., 2012). Under slower flow conditions, bacteria trapped in a 3D permeable flow cell show more QS dependent gene expressions (Connell et al., 2010). QS induction can be promoted if strong advection is decoupled by heterogeneous pore geometry (e.g., dead-end pores), allowing signaling molecules to accumulate (Kim et al., 2016; Ribbe & Maier, 2016).

The signaling molecules transported either via advection or diffusion can induce QS in neighboring cells (Frederick et al., 2010; Pérez-Velázquez, Gölgeli & García-Contreras, 2016). Because the signal concentration decreases with distance from its source, cells should be located close to each other in order to send and receive enough signaling molecules to and from their neighbors (Hense et al., 2007; Matur et al., 2015). The distance between two QS induced microbial cells or aggregates is referred to as the “calling distance” and has been reported to be 5–78 μm between individual cells (Gantner et al., 2006) and ~180 μm between microbial aggregates (Darch et al., 2018). However, the dependance of QS processes on advection and diffusion suggests that transport regimes affect calling distances, highlighting the importance of relative positioning of microorganisms coupled with the mass transport characteristics of a habitat.

Here, we evaluate the effect of combined diffusive and advective transport on QS processes in environmentally relevant conditions using a reactive transport modeling approach. The advection–diffusion–reaction equation was nondimensionalized to capture the characteristic properties of QS systems (i.e., production rates of signaling molecules, cell density, mass transport and spatial distribution of microbial aggregates) and used to formulate empirical expressions describing concentration profiles of signaling molecules under various flow conditions. Using these relationships, we evaluate calling distances and threshold biochemical conditions for QS induction of a single microbial aggregate under various flow conditions. Then, we investigate QS interactions between heterogeneously distributed microbial aggregates. Finally, we demonstrate the importance of autoinduction for coordinated microbial behaviors in advection-dominated environments. This study quantifies the effect of flow velocities, autoinduction, and relative position of microbial aggregates to calling distances in a 2D flow channel.

Materials and Methods

We used the Lattice Boltzmann (LB) method to implement a numerical model for the transport of signaling molecules due to diffusion and advection. The LB method is a mesoscopic approach solving the Boltzmann equation across a defined set of particles which recovers the macroscopic Navier–Stokes equation (NSE) and advection–diffusion–reaction equation (ADRE) (Tang et al., 2013; Krüger et al., 2017). First, we obtained the flow field by solving the particle distribution function f: fi(r+ciΔt,t+Δt)=fi(r,t)+Δtτ(fieq(r,t)fi(r,t))where particles fi(r, t) travel in the direction i with the lattice velocity ci (c0 = (0, 0), c1 = (1, 0), c2 = (0, 1), c3 = (−1, 0), c4 = (0, −1), c5 = (1, 1), c6 = (−1, 1), c7 = (−1, −1), c8 = (1, −1)) to a new position r + ciΔt after a time step Δt. The relaxation time (τ) was described by the commonly used Bhatnagar–Gross–Krook collision operator (Bhatnagar, Gross & Krook, 1954) and the D2Q9 lattice with the corresponding equilibrium distribution function: fieq(r,t)=ωiρ(1+ucics2+(uci2)2cs4uu2cs2)where ωi are lattice weights (ω0 = 4/9, ω1–4 = 1/9, ω5–8 = 1/36), cs is a lattice dependent constant (here, cs2 = 1/3), and u is the macroscopic flow velocity. The moments of the discretized mesoscopic particles retrieve the macroscopic density ρ=fi and momentum ρu=cifi. The Chapman-Enskog expansion showed that this LB approach recovers the incompressible NSE with the viscosity ν=cs2(τΔt2) (Krüger et al., 2017). Once the flow field was obtained, we simulated solute transport with a particle distribution function g, using the regularized LB algorithm (RLB) for numerical accuracy (Latt & Chopard, 2006; Latt, 2007) and the D2Q5 lattice for numerical efficiency (Li, Mei & Klausner, 2017): gi(r+ciΔt,t+Δt)=gieq(r,t)+(1Δtτ)ωi2cs4Qi:Πineq+ΩiRXN(r,t)where ci are the lattice velocities (c0 = (0, 0), c1 = (1, 0), c2 = (0, 1), c3 = (−1, 0), c4 = (0, −1), c5 = (1,1)) corresponding to the lattice weights ωi0 = 1/3, ω1–4 = 1/6), and Qi : Πineq is the tensor contraction of the two tensors Qi=ciciTcs2I and Πineq=jciciT(gj(r,t)gjeq(r,t)). The reaction term in the Eq. (3) describes the production of signaling molecules: ΩiRXN(r,t)=Δtωi(1+FH[A^θ^])k^B^where F represents a multiplication factor which was set to either 0 or 10 to reflect the magnitude of autoinduced signal production (Fekete et al., 2010), A^ is a concentration of signaling molecules, θ^ is the QS induction threshold, k^ is the basal production rate constant of signaling molecules, and B^ is the microbial density. QS induction often displays a switch-like behavior (Fujimoto & Sawai, 2013; Heilmann, Krishna & Kerr, 2015; Hense & Schuster, 2015), which is represented in the model by a step function with a higher signal production rate above the threshold concentration of signaling molecule: H[A^θ^]={1,(A^θ^)0,(A^<θ^)

With the imposed flow field from Eq. (1), the LB transport solver (Eq. 3) recovers the following ADRE: A^t^+u^^A^=D^^2A^+(1+FH[A^θ^])k^B^with the molecular diffusivity D^=cs2(τΔt2). Note that we are ignoring the breakdown of signaling molecules (Lee et al., 2002), limiting us to settings where production and transport are the dominant processes.

To describe the characteristic properties of a microbial system across various flow and reaction conditions, Eq. (6) was recast by introducing dimensionless quantities A=A^θ^, t=D^t^l^2, =^l^, B=B^B^θ, u=u^U^, where l^ is a characteristic length (i.e., the width of the flow channel), U^ is a characteristic fluid velocity (here, the average pore fluid velocity), and B^θ is a threshold biomass density required for QS induction, resulting in: At+PeuA=2A+Da

This nondimensionalized ADRE is fully characterized by the Péclet number, expressing the magnitude of advective flow relative to diffusion (Pe=U^l^D^), and the diffusive Damköhler number, comparing reaction to diffusion (Da=kBl^2θ^D^; where k=(1+FH[A^θ^])k^B^θ). A system with high Da, either due to high k’ (i.e., fast signal production), high B (i.e., high microbial density), or low θ^ (i.e., high sensitivity to signaling molecules) - is more likely to be QS induced.

An important property of Eq. (7) is that its solution linearly scales in Da (Lin, Xu & Zhang, 2020). For example, if Da is increased 2-fold at a fixed Pe condition, the concentrations of signaling molecule are doubled. This linearity allows to calculate the concentration distribution of signaling molecules for any Da from a single simulation result with an arbitrary Da at a given Pe. However, this simple approach cannot be applied to the flow conditions because the solution is not linear in Pe. Therefore, multiple numerical simulations were carried out with 24 Pe conditions (Pe ∈ {0.5, 0.6, 0.7, 0.8, 0.9, 1, 1.5, 2, 2.5, 3, 3.5, 4, 4.5, 5, 5.5, 6, 6.5, 7, 7.5, 8, 8.5, 9, 9.5, 10}) while Da was fixed at 5. For the 2D simulations in a straight channel (Fig. 1), the flow field was established by imposing pressures at in- and outlet and no flow conditions at the top and bottom boundaries, resulting in a flow from left to right. Fixed concentration (A|left boundary, x=0 = 0) and no-gradient (A/x|rightboundary,x=4=0) boundary conditions were imposed at the inlet and outlet boundaries, with no-flux at the top and bottom boundaries, respectively. All simulations were run to steady state.

Mathematical investigation of microbial quorum sensing under various flow conditions steady state concentration fields of signaling molecules at three Peclet numbers.

Figure 1: Mathematical investigation of microbial quorum sensing under various flow conditions steady state concentration fields of signaling molecules at three Peclet numbers.

Concentration fields of signal concentration (A) produced by microorganisms located at x = 1 and y = 0 with Da = 5 and (A) Pe = 1, (B) Pe = 5 and (C) Pe = 10, without autoinduction (F = 0). Note the difference in scale on the horizontal axis.

Simulations were conducted for a 2D flow channel of non-dimensional length of 4 and a width of 2, discretized with 2,000 × 1,000 grid elements. The flow field (Eq. 1) was generated by imposing fixed pressures at inlet (x = 0) and outlet (x = 4) with no flow boundaries in both normal and tangential direction at the bottom (y = 0) and top (y = 2) of the domain resulting in parabolic Poiseuille flows. Simulations were carried out under low Mach numbers (Ma = u/cs ≪ 1) to ensure incompressible flow conditions (Krüger et al., 2017).

Results and Discussion

QS processes of a single microbial aggregate

The effect of various flow conditions on the distribution of signaling molecules (A) produced from a single microbial aggregate assumed a source constrained to a single grid cell located at x = 1 was investigated under various Pe conditions (0.5 ≤ Pe ≤ 10) while Da was fixed at 5 (Fig. 1). The environmentally relevant range of Pe was chosen (Battiato et al., 2009, 2011) while Da is arbitrary because of the linearity of Eq. (7) in Da. The QS induction enhancing the signal production rate was not considered.

The signal concentration fields developed under various advective flows show maximum concentrations (Amax = A(x=1)) decreasing with increasing Pe (i.e., faster advective flow): Amax decreased from 1.68 (Pe = 1) to 1.35 (Pe = 5) and 1.21 (Pe = 10). However, Amax of all of the simulations with Da = 5 exceeded 1 (i.e., A^θ^), indicating the potential for QS induction. The threshold Da for QS induction (Daθ), where Amax = 1, can easily be computed using the linearity of the nondimensionalized ADRE in Da (Eq. 7). For example, Daθ at Pe = 1 was calculated by dividing Da = 5 by its corresponding Amax = 1.68 which resulted in Daθ = 2.98. Thus, at Pe = 1, conditions for which Da ≥ 2.98 lead to or exceed the concentration of signaling molecules needed for QS induction. Figure 2 shows the calculated Daθ for each simulated Pe condition.

Threshold Damkohler numbers under a range of Peclet numbers for quorum sensing induction.

Figure 2: Threshold Damkohler numbers under a range of Peclet numbers for quorum sensing induction.

The relationship between the threshold Da for QS induction (Daθ) and Pe. The simulation results (block dots) were fitted using the power regression (red line; Eq. 8).

The regression analysis revealed that the simulated Daθ for QS induction varies as a function of Pe following the power law: Daθ=1.3812Pe0.2626+1.592

The increasing Daθ along with the increasing Pe indicates higher Da (i.e., higher microbial density (B), higher signal production rate constant (k’), or lower QS induction threshold (θ^)) is required for QS induction under higher Pe. This result corresponds to the observed repressed QS induction under the presence of advection (Vaughan, Smith & Chopp, 2010; Meyer et al., 2012; Kim et al., 2016) and matches the pattern of biomass required for QS under varying flow conditions (Kirisits et al., 2007). Equation (8) was further evaluated by applying the experimentally measured QS parameters of Pseudomonas putida (k^ = 2.3 × 10−10 nmol/cell/h, and θ^ = 70 nmol/L (Fekete et al., 2010)) in a flow system where l^ = 1 cm and D^ = 3.0 × 10−10 m2/s (Dilanji et al., 2012). Our results show B^θ of 9.77, 12.2 and 13.5 × 106 cells/mL at Pe = 1, 5 and 10, respectively. If Eq. (8) is extrapolated to diffusion only transport condition (Pe = 0, Daθ = 1.592), B^θ is estimated as 5.23 × 106 cells/mL which largely agree with the experimental observation of 2.69~6.23 × 106 cells/mL where signal concentration starts to show a strong spike (Table S1 in Fekete et al., 2010).

In addition to reducing Amax, advection also influenced the spatial distribution of signaling molecules. We define the “transport distance” (d) as the distance between the point of production (x0) and the point (x1) where the signal concentration reaches a certain value A* (i.e., d = |x0x1|), distinguishing it from the “calling distance” between two QS induced microbial cells or aggregates. If the signal transport occurred only through diffusion, transport distances would be isotropic (Alberghini et al., 2009). However, advection resulted in anisotropic concentration distribution where upstream transport distances (dup) are much shorter than the downstream distances (ddn). Moreover, fast advective flows (i.e., high Pe) reduced overall transport distances which are illustrated in Figs. 1A1C as the shrinking areas covered by contour lines. For example, the (nondimensional) transport distances to the location where A = 0.1 are dup = 0.08 and ddn = 0.62 at Pe = 1 and Da = 5 (Fig. 1A). These values decrease to dup = 0.033 and ddn = 0.27 at Pe = 5 (Fig. 1B) and to dup = 0.023 and ddn = 0.19 at Pe = 10 (Fig. 1C).

Empirical approximation of concentration profiles

Obtaining transport distances for different Pe conditions requires running numerical simulations for each of the corresponding Pe. However, this may be avoided if we can express the concentration profiles as a function of Pe. For this purpose, parametric regression analysis was applied to the numerically obtained concentration profiles along the bottom of the flow channel (Fig. 3).

Constructing empirical relationships between the transport distance of signaling molecules and Peclet numbers.

Figure 3: Constructing empirical relationships between the transport distance of signaling molecules and Peclet numbers.

(A) Simulated (dots) and reconstructed (lines) concentration profiles along the bottom of the flow channel (y = 0.001) at Pe = 1 and Da = Daθ = 2.98. The upstream (x ≤ 1; Aup(x); solid line) and downstream (x > 1; Adn(x); dashed line) concentration profiles were obtained from Eqs. (9) and (10), respectively. (B–D) The coefficients for Aup(x) (aup, bup and cup) and Adn(x) (adn, bdn and cdn) obtained from the parametric regressions of the simulated concentration profiles at each simulated Pe conditions with Eq. (9) (black dots) and Eq. (10) (white dots), respectively. The solid and dashed lines are the exponential (Eqs. 1113) and power fits (Eqs. 1416) of the estimated coefficients as a function of Pe.

Several parametric regression models (linear, power, exponential and polynomial models) were tested to the upstream (Aup(x); 0 ≤ x ≤ 1) and downstream (Adn(x); 1 < x ≤ 4) signal concentration profiles. Among the tested regression models, the exponential (Eq. 9) and power-law models (Eq. 10) provided the best fit for log-transformed upstream and downstream signal concentration profiles, respectively. In the regression analysis of upstream profiles, only the locations where A(x) > 0.001 were used to improve the fitting quality and the signal concentration at x = 1 was fixed as 1. The additional regression analysis was then carried out for the coefficients (a, b and c) obtained from simulated profiles at 24 Pe conditions to construct a relationship between the coefficients and Pe (Figs. 3B3D). The exponential and power-law models provided the best fit for the upstream (Eqs. 1113) and downstream coefficients (Eqs. 1416), respectively: Aup(x)|x1=exp(aup(xbupxcup)) Adn(x)|x>1=exp(adnln(x)bdn+cdn)where Aup and Adn are 0 in the down- and up-stream directions, respectively, and aup=0.376exp(2.5975Pe)+2.7165exp(0.0244Pe)

bup=20.311exp(0.0115Pe)13.38exp(0.6121Pe) cup=7.1289exp(0.0348Pe)+5.9469exp(0.4272Pe) 𝑎dn=8.6156𝑃𝑒0.066813.3056 bdn=0.1051Pe0.2522+0.1082 cdn=7.5322Pe0.0464+8.7195

Equations (9) and (10) can be used as approximations of the concentration profiles along a pore channel without running simulations for various Pe conditions, with the microbial aggregate located at x = 1. Due to the linearity in Da, the concentration profiles at different Da conditions can be calculated simply by multiplying Da/Daθ to Eqs. (9) and (10), so that A(x)=DaDaθ(Aup(x)+Adn(x))

These analytical expressions are applicable not only to QS but also to other chemical processes subject to zero-order production reactions (e.g., Bezemer et al., 2000; Tang et al., 2015). The equations become less accurate at low Pe as under low flow conditions, the estimates from Eq. (17) in a flow channel with a small width (i.e., low l^ and Pe) could underestimate the actual concentration because the confined channel width would push the produced chemical further upstream and downstream.

The effect of QS induced signal production on transport distances

QS often involves autoinduction which substantially increases signal production rates. The effect of autoinduction on transport distances was investigated by using Eq. (17) for the conditions without (F = 0; Da = Daθ) and with (F = 10; Da = 11Daθ) enhanced signal production. The transport distances from a single microbial aggregate under various Pe were then calculated using Eq. (17) for the location x.

Figure 4 shows the transport distances without (Fig. 4A) and with (Figs. 4B and 4C) the enhanced signal production at Pe = 1, 5 and 10. The concentration ratios (0.1 ≤ A/Amax ≤ 0.9) were used instead of absolute concentrations to generalize transport distances for various Da conditions. For example, the transport distance (dA) for A/Amax = 0.5 indicates that A(x0 + dA) = 0.5 if Da = Daθ while A(x0 + dA) = 0.05 when Da = 0.1Daθ. The consequence of the enhanced signal production was the significant increase of dup and ddn. Without the enhanced signal production, dup and ddn for A/Amax = 0.4 at Pe = 1 were estimated as 0.021 and 0.024, respectively (Fig. 4A). These values increased to dup = 0.1 and ddn = 1.28 with the enhanced signal production (Figs. 4B and 4C). The downstream transport distance of 1.28 is translated into 6.4 mm in a flow channel with l^ = 1 cm. This result is much longer than the generally observed ranges of calling distances (Whiteley, Diggle & Greenberg, 2017). However, we emphasize again that the transport distance merely indicates the distance of signaling molecules transported from a source location while the calling distance involves QS induced microbial cells or aggregates.

Transport distances of signaling molecules with and without autoinduction.

Figure 4: Transport distances of signaling molecules with and without autoinduction.

Upstream (dup) and downstream (ddn) transport distances (A) without (F = 0) and (B) with (F = 10) enhanced signal production for the concentration ratios (0.1 ≤ A/Amax ≤ 0.9) at Pe = 1, 5 and 10, and (C) the enlarged barplot of upstream transport distances with F = 10. Note the different scale of the horizontal axes between panels.

QS induction between spatially distributed multiple microbial aggregates

QS processes of multiple aggregates were investigated by constructing the concentration profiles using Eq. (18). Concentration fields of signaling molecules with multiple microbial aggregates can be calculated as the superposition of the concentration profile produced by each individual aggregate: A(x)=i=1nDaiDaθ(Aup(x+di0)+Adn(x+di0))where n is the number of aggregates, di0 is the distance between xi and x0 (di0 = xix0), xi is the location of ith aggregate, x0 is the reference location (x0 = 1), Dai is the Da calculated only with the density of ith microbial aggregate (i.e., microscopic Da), and Aup and Adn are Eqs. (9) and (10), respectively. Here, an example system with macroscopic Da(DaT=Dai) = 3.2Daθ consist of four aggregates (A1–4) located at x1 = 0.4, x2 = 1, x3 = 1.096 and x4 = 1.7 with the evenly distributed microscopic Dai (i.e. Da1 = Da2 = Da3 = Da4 = 0.8Daθ) was tested. In using Eq. (18), the profile was first constructed for Dai = Da* that does not consider autoinduction (F = 0). Then, if there is an aggregate with A(xi) ≥ 1, the profile was recalculated with updated Dai = (1 + F) × Da* until all Dai with A ≥ 1 were updated.

The signal concentration profile produced by four aggregates without the enhanced signal production (F = 0) illustrates the crucial importance of relative positioning of microbial aggregates for QS induction with respect not only to each aggregate but also to the flow direction (Fig. 5A). The microscopic Dai was set such that the maximum concentration produced by a single aggregate was 0.8, as observed at the most upstream location (A1 at x1 = 0.4). But due to transport, the local concentration at A2 reached 0.879, receiving A of 0.048 and 0.031 from A1 and A3, respectively. A3 received slightly less signaling molecules from A1 (A = 0.044) due to the longer distance of A3 than A2 from A1. However, A2 provided much more signaling molecules (A = 0.157) to A3 than was provided by A3 because of advective flows favoring downstream transport of signaling molecules (Figs. 2 and 4). As a consequence, A3 exceeded the QS threshold (A(x3) = 0.044 from A1 + 0.157 from A2 + 0.8 from A3 + 0 from A4 = 1.001 > 1) while the upstream located A2 did not. The QS induction of A3 demonstrates the importance of transport distances. QS induction was achieved because of the upstream aggregates located within the transport distance of 0.696. However, the calling distance would have been estimated as the length of a grid voxel (0.002) because only A3 was QS induced. Therefore, considering only the calling distance could lead to the wrong conclusion that the local Da condition at A3 (i.e., Da3 = 0.8Daθ) is a sufficient condition for QS induction. Although A4 did not reach the QS induction threshold, it received A from all the other aggregates resulting in a concentration (A(x4) = 0.029 + 0.044 + 0.048 + 0.8 = 0.921) that was higher than at A2 despite the longest separation distance from other aggregates.

Quorum sensing amongst multiple microbial aggregates.

Figure 5: Quorum sensing amongst multiple microbial aggregates.

Concentration profile (A) without (F = 0) and (B) with (F = 10) enhanced signal production where four aggregates are located at x1 = 0.4, x2 = 1, x3 = 1.096 and x4 = 1.7. Black dots are the simulated results and the yellow lines represent the profile from Eq. (18).

Accounting for QS induction (F = 10) increased the transport distances and hence induced other aggregates (Fig. 5b). With the same spatial distribution, QS-induced A3 produced signaling molecules much more and faster (i.e. k’ = 11k^ and Da3 = 8.8Daθ) and provided more signaling molecules to A2. As a result, A(x2) exceeded the QS threshold (0.048 + 0.8 + 0.335 + 0 = 1.183). The QS induction of A2 and A3 resulted in the final signal concentrations of A(x2) = 9.183 (= 0.048 + 8.8 + 0.035 + 0) and A(x3) = 10.569 (= 0.044 + 1.725 + 8.8 + 0). While A4 still did not contribute signaling molecules to any of upstream aggregates, enhanced contribution from A2 and A3 QS induced A4, A(x4) = 9.839 (= 0.029 + 0.48 + 0.53 + 8.8). Despite increased transport distances by QS induction, A1 was still too far away from the other aggregates thus the signal concentration at A1 remained unchanged A(x1) = 0.8. As a result of the QS induction of A2-4, DaT had increased from the initial 3.2Daθ (= 0.8Daθ × 4) to 27.2Daθ (= 0.8Daθ + 3 × 11 × 0.8Daθ).

This example illustrates the importance of enhanced signal production on the spatial propagation of QS induction. While only A3 experienced signaling molecule levels that could induce QS when all the aggregates produce signaling molecules at the basal production rate, the enhanced signal production of A3 when considering induced production (F = 10) provided more signaling molecules to its adjacent microbial aggregates and resulted in the QS induction of neighboring aggregates, A2 and A4. It may be counterintuitive that the upstream-located A2 was also QS-induced by the contribution from A3 despite the contracted upstream transport distances under the presence of advective flows. This result shows that the enhanced signal production can overcome the influence of advection and promote QS induction, and provide a way to provoke upstream microbial aggregates, for example, to slow down the substrate consumption to ensure efficient resource utilization in crowded environments (An et al., 2014).

Conclusions and Perspectives

This study has demonstrated that advection and the enhanced signal production can determine the spatial extent of QS induction. Reactive transport simulation results reveal that fast flow conditions dilute signaling molecules and thus higher Daθ (i.e., faster signal production rate, higher microbial density, or lower QS induction threshold concentration) is required for QS induction. Reduced upstream delivery of signaling molecules under advective flow limits propagation of QS; it can be relaxed if autoinduction increases signal production rates. Our study results highlight the importance of relative positioning of microbial aggregates with respect to flow directions and the role of autoinduction to overcome advection for upstream signal propagation.

The simulations focused on the effect of various flow conditions on QS and assumed that microbial aggregates have a negligible impact on flow fields, which is a reasonable approximation for low microbial density conditions. However, it may not hold when large aggregates producing extracellular polymeric substances (EPS) perturb flows substantially. In such a case, estimating signal transport requires fully resolving nonlinear feedback between cell activity and fluid flow (Thullner & Baveye, 2008; Carrel et al., 2018; Jung & Meile, 2019), including diffusion limitation (Stewart, 2003). Finally, accounting for degradation of signaling molecules (Lee et al., 2002; Yates et al., 2002) and increased spreading of signaling molecules in 3D systems than 2D, would result in shorter transport distances than this study.

Although QS mediated gene expression has been understood as evolutionarily beneficial collective behaviors, long transport distances observed in this study suggests that it may not be always true. The transport of signaling molecules, especially in downstream direction, combined with enhanced signal production, suggests that QS induction can be decoupled from microbial density. In the above example (Fig. 5B), any microbial cell located where A > 1 (e.g., A(x = 2.5) = 1.05) would have been QS-induced, independent of the local cell density. This could lead to detrimental impacts on a microbial population, unless there are other counteracting mechanisms such as differential QS induction sensitivity to signal concentration even in within a clonal population (Darch et al., 2018) or biofilm formation modifying local transport characteristics (Emge et al., 2016). Future investigations should explicitly examine the evolutionary consequences of QS strategies in spatially heterogeneous environments under advective–diffusion–reaction dynamics.