A biophysical threshold for biofilm formation
Figures
![](https://iiif.elifesciences.org/lax/76380%2Felife-76380-fig1-v2.tif/full/617,/0/default.jpg)
Competition between motility-mediated dispersal and autoinducer-mediated biofilm formation.
(A) Schematic of chemotactic dispersal: planktonic bacteria (green) consume nutrient (purple) and establish a local gradient that they, in turn, direct their motion in response to. (B) Schematic of positive quorum sensing-controlled biofilm formation: accumulation of produced autoinducer (red) above a threshold concentration causes cells to transition to the biofilm state (blue). (C) Results of an example simulation of Equations 1–4 showing the dynamics of the nutrient, planktonic cells, autoinducer, and biofilm cells from top to bottom, quantified by the normalized concentrations , , , , respectively; , , and represent the initial nutrient concentration, initial bacterial concentration, and autoinducer threshold for biofilm formation, respectively. The position coordinate is represented by the normalized position , where is the width of the initial cellular inoculum. Different shades indicate different time points as listed. The inoculum initially centered about the origin consumes nutrient (purple), establishing a gradient that drives outward dispersal by chemotaxis (outward moving green curves); the cells also produce autoinducer (red) concomitantly. At , sufficient autoinducer has been produced to trigger biofilm formation at the origin; at even longer times (), nutrient depletion limits autoinducer production at this position. However, accumulation of autoinducer by the dispersing planktonic cells triggers partial biofilm formation at as well. This competition between dispersal and biofilm formation leads to a final biofilm fraction of at the final time of . An animated form of this figure is shown in Video 1. The values of the simulation parameters are given in Supplementary file 2.
![](https://iiif.elifesciences.org/lax/76380%2Felife-76380-fig1-figsupp1-v2.tif/full/617,/0/default.jpg)
The exact nature of the temporal dynamics of the arrest in motility while transitioning to the biofilm state does not appreciably influence our model results.
Results of the same simulation as in Figure 1C, but with both motility parameters and smoothly transitioning to zero in time for cells that are transitioning to the biofilm state—instead of these parameters transitioning step-wise to zero as in the simulations presented in the main text. To implement this change, we monitor the amount of time elapsed after planktonic cells are exposed to sufficient autoinducer to drive the transition to the biofilm state (i.e., after the moment when ), and decrease the motility parameters linearly to reach zero after a total time has elapsed: the diffusivity and chemotactic coefficients are given by and , respectively. When, both motility parameters are zero. The results, shown above, are visually identical to those in Figure 1C.
-
Figure 1—figure supplement 1—source data 1
- https://cdn.elifesciences.org/articles/76380/elife-76380-fig1-figsupp1-data1-v2.xlsx
![](https://iiif.elifesciences.org/lax/76380%2Felife-76380-fig1-figsupp2-v2.tif/full/617,/0/default.jpg)
The details of the initial inoculum shape do not appreciably influence our model results.
-
Figure 1—figure supplement 2—source data 1
- https://cdn.elifesciences.org/articles/76380/elife-76380-fig1-figsupp2-data1-v2.xlsx
![](https://iiif.elifesciences.org/lax/76380%2Felife-76380-fig2-v2.tif/full/617,/0/default.jpg)
Faster nutrient consumption limits autoinducer production, leading to complete dispersal.
Retion as in Figure 1C, but for planktonic cells with faster nutrient consumption (larger ). Panels and colors show the same quantities as in Figure 1C. The inoculum initially centered about the origin consumes nutrient (purple), establishing a gradient that drives outward dispersal by chemotaxis (outward moving green curves); the cells also produce autoinducer (red) concomitantly. However, nutrient is depleted at this position more rapidly, limiting autoinducer production; as a result, the population continues to disperse in the planktonic state and the final biofilm fraction is . An animated form of this figure is shown in Video 2. The values of the simulation parameters are given in Supplementary file 2.
![](https://iiif.elifesciences.org/lax/76380%2Felife-76380-fig2-figsupp1-v2.tif/full/617,/0/default.jpg)
Slower nutrient consumption allows greater autoinducer production, leading to more biofilm formation.
Results of the same simulation as in Figure 1C, but for planktonic cells with slower nutrient consumption (smaller ). Panels and colors show the same quantities as in Figure 1C. The inoculum initially centered about the origin slowly consumes nutrient (purple), establishing a slight gradient that allows partial planktonic dispersal (green curves moving outward); the cells also produce autoinducer (red) concomitantly. Because nutrient is consumed slowly, autoinducer production is not limited, resulting in partial biofilm formation (blue). Autoinducer has sufficiently accumulated above the threshold after , which causes a population of biofilm cells to form at the origin (). After , the biofilm population continues to grow, and additionally, autoinducer concentration exceeds the threshold concentration at . Thus, we see a second population of biofilm cells form, centered at.. The slower nutrient consumption results in a greater final biofilm fraction than in Figure 1C—here,.. An animated form of this figure is shown in Video 3. The values of the simulation parameters are given in Supplementary file 2.
-
Figure 2—figure supplement 1—source data 1
- https://cdn.elifesciences.org/articles/76380/elife-76380-fig2-figsupp1-data1-v2.xlsx
![](https://iiif.elifesciences.org/lax/76380%2Felife-76380-fig3-v2.tif/full/617,/0/default.jpg)
Enhanced motility enables cells to disperse before sufficient autoinducer accumulates, leading to complete dispersal.
(A) Results of the same simulation as in Figure 1C, but for faster-moving planktonic cells (larger and ). Panels and colors show the same quantities as in Figure 1C. The inoculum initially centered about the origin consumes nutrient (purple), establishing a gradient that drives outward dispersal by chemotaxis (outward moving green curves); the cells also produce autoinducer (red) concomitantly. More rapid dispersal enables the planktonic cells to ‘outrun’ the growing autoinducer plume, as shown by the extended and magnified view in (B). As a result, the population continues to disperse in the planktonic state and the final biofilm fraction is . An animated form of this figure is shown in Video 4. The values of the simulation parameters are given in Supplementary file 2.
![](https://iiif.elifesciences.org/lax/76380%2Felife-76380-fig3-figsupp1-v2.tif/full/617,/0/default.jpg)
Diminished motility enables autoinducer to accumulate, resulting in increased biofilm formation.
Results of the same simulation as in Figure 1C, but for slower-moving planktonic cells (smaller and ). Panels and colors show the same quantities as in Figure 1C. The inoculum initially centered about the origin consumes nutrient (purple), establishing a slight gradient—however, because the motility parameters are diminished, the planktonic population (green) remains around the origin. The planktonic cells produce autoinducer (red) concomitantly, and after , the autoinducer concentration exceeds the threshold concentration. Thus, some of the planktonic cells transition to biofilm cells, centered at the origin. Both the biofilm cells and planktonic cells continue to grow, produce autoinducer, and consume nutrient; the planktonic cells do not disperse due to their diminished motility, resulting in a larger fraction of biofilm cells () than in Figure 1C. An animated form of this figure is shown in Video 5. The values of the simulation parameters are given in .
-
Figure 3—figure supplement 1—source data 1
- https://cdn.elifesciences.org/articles/76380/elife-76380-fig3-figsupp1-data1-v2.xlsx
![](https://iiif.elifesciences.org/lax/76380%2Felife-76380-fig4-v2.tif/full/617,/0/default.jpg)
The two states of complete dispersal by planktonic cells (green) and complete formation of a biofilm (blue) can be universally described by three dimensionless parameters.
(A) State diagram showing the fraction of biofilm formed, , at the final time () for different values of the nutrient availability and cellular dispersal parameters, and , respectively. The state diagram summarizes the results of 10,983 simulations of Equations 1–4 exploring the full range of parameter values describing different bacterial species/strains and different environmental conditions (Supplementary file 1). Each point represents the mean value of obtained from multiple simulations with different parameter values, but with similar and (identical within each bin defined by the spacing between points). (B) represents the same data, but each point represents the standard deviation of the values of obtained from the same simulations. Despite the vastly differing conditions explored in each simulation, they cluster into the two states of planktonic dispersal (green) and biofilm formation (blue) when parameterized by and . The boundary between the two states can be described by the relation , as shown by the black line; this relation combines the transition between the two states that occurs at both and . Away from this boundary, all simulations for the same and collapse to have the same biofilm fraction , as shown by the points in (B) and examples (i)–(iii) and (v) in (C)—confirming the universality of our parameterization. Near the boundary, we observe some slight differences between simulations, as shown in (B) and examples (iv) and (vi) in (C). The values of the simulation parameters for the examples in (C) are given in Source Data file 1. The data in (A–C) correspond to a fixed value of the third dimensionless parameter , which describes the case of biofilm cells that produce autoinducer rapidly; repeating these simulations for the opposite case of slow autoinducer production by biofilm cells () yields the state diagram shown in (D), but for 14,351 simulation runs; again, (E) shows the standard deviation of the corresponding values of . As shown by (D–E), while the transition between the two states (black line) is unaffected by the change in , the transition to complete biofilm formation is more gradual. Together, the three parameters , , and provide a full description of the onset and extent of biofilm formation across vastly different conditions.
-
Figure 4—source data 1
- https://cdn.elifesciences.org/articles/76380/elife-76380-fig4-data1-v2.xlsx
-
Figure 4—source data 2
- https://cdn.elifesciences.org/articles/76380/elife-76380-fig4-data2-v2.xlsx
-
Figure 4—source data 3
- https://cdn.elifesciences.org/articles/76380/elife-76380-fig4-data3-v2.xlsx
![](https://iiif.elifesciences.org/lax/76380%2Felife-76380-fig4-figsupp1-v2.tif/full/617,/0/default.jpg)
In the case of ‘protected’ nutrient-independent autoinducer production, the transition from planktonic to biofilm states occurs at, independent of .
The model described in the main text considers autoinducer secretion to be nutrient-dependent; here, we consider the alternate case of ‘protected’ nutrient-independent production at the maximal rate per cell. In this case, we do not expect the transition from the planktonic to biofilm state to depend on nutrient availability, as quantified by . Instead, the analysis presented in the main text suggests that the dispersal parameter, , alone should be sufficient to capture this transition. To test this expectation, we perform the same numerical simulations as presented in the main text, but without the Michaelis-Menten function multiplying the autoinducer production term in Equation 4. The resulting state diagrams showing the final () fraction of biofilm formed, , and the standard deviation of the values, as a function of and are shown in (A–B). These diagrams summarize 8709 separate simulations. The results confirm our expectation: in contrast to the results shown in Figure 4, the transition does not appreciably depend on , but instead occurs near (black line). Thus, the theoretical framework developed in this work is more general and can be applied to bacteria with nutrient-dependent or nutrient-independent autoinducer production.
-
Figure 4—figure supplement 1—source data 1
- https://cdn.elifesciences.org/articles/76380/elife-76380-fig4-figsupp1-data1-v2.xlsx
![](https://iiif.elifesciences.org/lax/76380%2Felife-76380-fig4-figsupp2-v2.tif/full/617,/0/default.jpg)
Simulation results are not appreciably influenced by our choice of discretization.
To assess the sensitivity of our results to numerical discretization, we repeat two representative simulations with (A) , corresponding to, and (B) , corresponding to, with varying spatial and temporal step sizes, and, , respectively. The spatial step size is changed by a factor and the corresponding temporal step size is changed by a factor, accordingly, as indicated by the lower and upper horizontal axes, respectively. As shown in the figure, the final biofilm fraction obtained from the simulations is not strongly sensitive to the choice of numerical discretization for spatial and temporal steps smaller than those used for the simulations presented in the main text (indicated by the red stars). Thus, our choice of discretization is sufficiently finely resolved such that the results in the numerical simulations are not appreciably influenced by discretization.
-
Figure 4—figure supplement 2—source data 1
- https://cdn.elifesciences.org/articles/76380/elife-76380-fig4-figsupp2-data1-v2.xlsx
Videos
Animated form of Figure 1C: Results of an example simulation of Equations 1–4 showing the dynamics of the nutrient, planktonic cells, autoinducer, and biofilm cells from top to bottom, quantified by the normalized concentrations , , , and , respectively; , , and represent the initial nutrient concentration, initial bacterial concentration, and autoinducer threshold for biofilm formation, respectively.
The position coordinate is represented by the normalized position , where is the width of the initial cellular inoculum. The inoculum initially centered about the origin consumes nutrient (purple), establishing a gradient that drives outward dispersal by chemotaxis (outward moving green curves); the cells also produce autoinducer (red) concomitantly. At , sufficient autoinducer has been produced to trigger biofilm formation at the origin; at even longer times (), nutrient depletion limits autoinducer production at this position. However, accumulation of autoinducer by the dispersing planktonic cells triggers partial biofilm formation at as well. This competition between dispersal and biofilm formation leads to a final biofilm fraction of at the final time of . The values of the simulation parameters are given in Supplementary file 2. The video displays the profiles every , to retain a manageable file size; however, the temporal step size in the actual simulations is .
Animated form of Figure 2: Results of the same simulation as in Video 1, but for planktonic cells with faster nutrient consumption (larger ).
Panels and colors show the same quantities as in Video 1. The inoculum initially centered about the origin consumes nutrient (purple), establishing a gradient that drives outward dispersal by chemotaxis (outward moving green curves); the cells also produce autoinducer (red) concomitantly. However, nutrient is depleted at this position more rapidly, limiting autoinducer production; as a result, the population continues to disperse in the planktonic state and the final biofilm fraction is . The values of the simulation parameters are given in Supplementary file 2. The video displays the profiles every , to retain a manageable file size; however, the temporal step size in the actual simulations is .
Animated form of Figure 2—figure supplement 1: Results of the same simulation as in Video 1, but for planktonic cells with slower nutrient consumption (smaller ).
Panels and colors show the same quantities as in Video 1. The inoculum initially centered about the origin slowly consumes nutrient (purple), establishing a slight gradient that allows partial planktonic dispersal (green curves moving outward); the cells also produce autoinducer (red) concomitantly. Because nutrient is consumed slowly, autoinducer production is not limited, resulting in partial biofilm formation (blue). Autoinducer has sufficiently accumulated above the threshold after , which causes a population of biofilm cells to form at the origin (). After , the biofilm population continues to grow, and additionally, autoinducer concentration exceeds the threshold concentration at . Thus, we see a second population of biofilm cells form, centered at . The slower nutrient consumption results in a greater final biofilm fraction than in Video 1—here, . The values of the simulation parameters are given in Supplementary file 2. The video displays the profiles every , to retain a manageable file size; however, the temporal step size in the actual simulations is .
Animated form of Figure 3: Results of the same simulation as in Video 1, but for faster-moving planktonic cells (larger and ).
Panels and colors show the same quantities as in Video 1. The inoculum initially centered about the origin consumes nutrient (purple), establishing a gradient that drives outward dispersal by chemotaxis (outward moving green curves); the cells also produce autoinducer (red) concomitantly. More rapid dispersal enables the planktonic cells to ‘outrun’ the growing autoinducer plume, as shown by the extended and magnified view in (B). As a result, the population continues to disperse in the planktonic state and the final biofilm fraction is . The values of the simulation parameters are given in Supplementary file 2. The video displays the profiles every , to retain a manageable file size; however, the temporal step size in the actual simulations is .
Animated form of Figure 3—figure supplement 1: Results of the same simulation as in Video 1, but for slower-moving planktonic cells (smaller and ).
Panels and colors show the same quantities as in Video 1. The inoculum initially centered about the origin consumes nutrient (purple), establishing a slight gradient—however, because the motility parameters are diminished, the planktonic population (green) remains around the origin. The planktonic cells produce autoinducer (red) concomitantly, and after , the autoinducer concentration exceeds the threshold concentration. Thus, some of the planktonic cells transition to biofilm cells, centered at the origin. Both the biofilm cells and planktonic cells continue to grow, produce autoinducer, and consume nutrient; the planktonic cells do not disperse due to their diminished motility, resulting in a larger fraction of biofilm cells () than in Video 1. The values of the simulation parameters are given in . The video displays the profiles every , to retain a manageable file size; however, the temporal step size in the actual simulations is .
Additional files
-
Transparent reporting form
- https://cdn.elifesciences.org/articles/76380/elife-76380-transrepform1-v2.pdf
-
Source data 1
Data for Figure 4C: Excel book summarizing the inputs used for each of the ten simulations shown in Figure 4C (i–vi).
Different sheets within the book correspond to the individual bar graphs—Sheet “Bar Graph (i)” corresponds to all ten runs shown in (i), Sheet “Bar Graph (ii)” corresponds to all ten runs shown in (ii), and this follows for the remaining sheets. The first column identifies the input, its variable, and its units. The following columns correspond to each of the ten runs shown in the bar graph.
- https://cdn.elifesciences.org/articles/76380/elife-76380-data1-v2.xlsx
-
Supplementary file 1
Ranges of the values of the parameters explored in our model with corresponding references.
- https://cdn.elifesciences.org/articles/76380/elife-76380-supp1-v2.zip
-
Supplementary file 2
Specific parameter values for the simulations shown in Figures 1—3, Figure 2—figure supplement 1 and Figure 3—figure supplement 1.
- https://cdn.elifesciences.org/articles/76380/elife-76380-supp2-v2.zip