Computational and Mathematical Methods in Medicine
Volume 2015 (2015), Article ID 581829, 10 pages
http://dx.doi.org/10.1155/2015/581829

Research Article

Modeling of the Bacillus subtilis Bacterial Biofilm Growing on an Agar Substrate

Xiaoling Wang,1,2 Guoqing Wang,¹ and Mudong Hao¹

¹School of Mechanical Engineering, University of Science and Technology Beijing, Beijing 100083, China
²School of Engineering and Applied Sciences, Harvard University, Cambridge, MA 02138, USA

Received 13 May 2015; Revised 21 July 2015; Accepted 26 July 2015

Academic Editor: Chung-Min Liao

Copyright © 2015 Xiaoling Wang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Bacterial biofilms are organized communities composed of millions of microorganisms that accumulate on almost any kinds of surfaces. In this paper, a biofilm growth model on an agar substrate is developed based on mass conservation principles, Fick’s first law, and Monod’s kinetic reaction, by considering nutrient diffusion between biofilm and agar substrate. Our results show biofilm growth evolution characteristics such as biofilm thickness, active biomass, and nutrient concentration in the agar substrate. We quantitatively obtain biofilm growth dependence on different parameters. We provide an alternative mathematical method to describe other kinds of biofilm growth such as multiple bacterial species biofilm and also biofilm growth on various complex substrates.

1. Introduction

A biofilm is an assemblage of microbial cells that is irreversibly associated with a surface and enclosed in an extracellular matrix which are secreted by bacteria [1]. Biofilm is ubiquitous in natural and industrial environment, and it has both good and bad effects on us. Biofilm can cause problems in energy losses due to increased fluid frictional resistance of the hull surfaces and increased heat transfer resistance in heat exchange equipment, device damage, food contamination, and medical infections; however, biofilms also have a lot of applications such as biofiltration of wastewater and remediation of contaminated soil and groundwater [2–4].

Scientists did a lot of work on biofilm, such as biofilm morphology, structures, and growth process. They proposed models on biofilm to better understand its growth. These models are established for two types of biofilms: one is growing in liquid environment and the other is growing on solid nutritive substrate. The difference between two types is the way of nutrient supply. The biofilm in liquid environment is fed on nutrient from the sides, top, and substrate; the models for them described their heterogeneous structures and growth characteristics [5–8], due to osmotic pressure or extracellular polymeric substances (EPS) [9–11]. The biofilm growing on solid nutritive substrate is exposed to air and only fed through bottom substrate, which contains ions salts and nutrient for biofilm growth, and models for this kind of biofilm described a fractal colony with various morphologies through the diffusion-limited aggregation process [12, 13], studying the role of osmotic stress in biofilm spreading based on the physics of polymer solutions [14].

There are some key factors deciding whole biofilm growth, for example, the yield coefficient, which is an important parameter in biofilm growth and degradation processes [15–17], the maximum specific growth rate, which indicates the microbial inherent growth characteristic [18], nutrient diffusion between biofilm and solid substrate or liquid, which is influenced by extracellular polysaccharides and cellular layers, and the biofilm density, which affects diffusion coefficient and thus has an influence on biofilm growth characteristic [19, 20]. We will consider the above four main parameters in our biofilm growth modeling; for simplicity, we assume diffusion coefficient between biofilm and substrate is a constant.

The existing models for biofilm growth on an agar substrate did not analyze influences of the above main factors on biofilm growth. In our paper, we establish a mathematical biofilm growth model on an agar plate based on mass conservation principle, Fick’s first law, and Monod’s kinetic reaction. From the numerical analysis, we can obtain dynamic evolution of the biofilm thickness, the volume fraction of active biomass and nutrient concentration in the agar substrate, and parameter effect on biofilm growth. Our model can be used as a tool for describing of other kinds of biofilm growth such as multiple bacterial species biofilm and also biofilm growth on various complex substrates.

2. Materials and Methods

Bacillus subtilis strain NCIB 3610 is used for all experiments. We grow colonies on 1.5 wt% agar gel with minimal media, MSgg, designed to induce biofilm formation: 5 mM potassium phosphate (pH 7); 100 mM MOPS (pH 7); 2 mM MgCl₂; 700 μM CaCl₂; 50 μM MnCl₂; 50 μM FeCl₂; 1 μM ZnCl₂; 2 μM thiamine; 0.5% glycerol; 0.5% glutamate; 50 g/mL tryptophan; 50 g/mL phenylalanine. The agar solution is cooled to 55°C before adding the remaining ingredients. We typically use 100 mm diameter Petri dishes containing 12 mL of media to obtain a 16 mm diameter biofilm after two days. The plates are covered with lids and cooled overnight at room temperature and then spotted within 24 h.

The triple reporter Bacillus subtilis strain was obtained from the Kolter lab. This strain was created using phage transduction to fuse a fluorescent gene to the featured promoter into the bacterial chromosome following the established procedure [21]: Red fluorescent protein mKate2 is linked to the hag promoter, which is responsible for flagellin production in motile cells. Blue fluorescent protein cfp is linked to the tapA promoter, which produces the amyloid proteins in matrix-producing cells. Green fluorescent protein citrus is linked to the sspB promoter, which is encoded in late period in sporulating cells. When a featured promoter for a phenotype is expressed, the relevant fluorescent protein is also produced. There is no spectrum overlap of the three fluorophores, except for some minor autofluorescence with CFP filter set. The autofluorescence is because of the pigment production of the colony after 48 h; so most of our experiments last only 48 h.

We transfer the bacteria to the agar surface by spotting with 0.1 μL of bacterial culture at OD₆₀₀ = 1; before inoculating the plates, we remove the lids and allow the surface to dry for 5–10 minutes. We allow the drop to dry for 5–10 additional minutes with the lid off until the meniscus of the initial drop is no longer visible and the bacteria are left in a “coffee ring” around the perimeter. For two days’ time-lapse movies, we grow biofilm colonies in a Tupperware container stuffed with wet paper towels and sealed around the microscope using Glad Press’n’Seal plastic wrap to prevent evaporation. The temperature of the microscope is maintained at 32°C using heating elements and fans.

3. The Biofilm Growth Model

We assume the biofilm consists of active biomass, inactive biomass, and water, which are shown in Figure 1; here the biomass is treated as a homogeneous continuum [22]. Inactive biomass is related to endogenous decay and the fraction of the active biomass that is not biodegradable, which also indicates the inert biomass [23, 24]. We consider nutrient diffusion between biofilm and agar substrate, which is shown as red arrows in Figure 1. In our experiment and model, the height of agar substrate is 0.5 cm, and the biofilm area is 1 cm². The volume of the agar substrate is 0.5 cm³.

Figure 1: The simplified schematic of biofilm model on the agar substrate. A biofilm is assumed to be made up of active biomass, inactive biomass, and water. The red circle represents active biomass, the yellow rhombus represents inactive biomass, and the green region is water. Nutrient diffuses from agar substrate to biofilm; red arrows indicate the diffusion direction.

The unknown dependent variables determined by the model are given in Table 1 and the other variables and the parameters used to develop this model are given in Table 2. Through our paper, the fundamental units of nutrient mass , biofilm mass , length , and time are identified. Derivation of biofilm model on agar substrate is shown as in Figure 1.

Table 1: Unknown dependent variables and their fundamental units in the biofilm growth model on agar substrate.

Variables
Physical quantity
Units

Biofilm thickness
L

Nutrient concentration in the agar substrate

Average nutrient concentration in the biofilm

Average volume fraction of active biomass in the biofilm
None

Average volume fraction of inactive biomass in the biofilm
None

Table 2: Parameters used in the biofilm growth model and their fundamental units.

Parameters
Physical quantity
Units

Area of the film-agar interface

Diffusion coefficient for the substrate through

the laminar diffusional sublayer

Thickness of laminar diffusional sublayer

Maximum specific growth rate for Monod kinetics

Monod constant

Biomass density

Inactivation coefficient

Yield coefficient

3.1. Mass Balance of Active Biomass

The rate of change of the active biomass depends on the increment rate of active biomass due to cell growth and the inactivation rate of active biomass. The rate of change of the active biomass is given by

The equation can be further reduced to

3.2. Mass Balance of Inactive Biomass

The inactive biomass increases as the active biomass becomes inactive. The rate of change of inactive biomass can be given by

which simplifies to

Adding this to (1) and using the assumption that ( is the volume fraction of water in the biofilm) [26], we get the rate of change of the biofilm thickness:

3.3. Mass Balance of Nutrient in Biofilm

The average rate of change of the nutrient in the biofilm depends on the consumption rate of the nutrient and the diffusion rate of the nutrient into the biofilm. A mathematical expression for the rate of change of the nutrient in the biofilm may be given by

which reduces to

Here, , refers to the nutrient flux through the diffusion layer (), according to Fick’s first law [22].

3.4. Mass Balance of the Nutrient Concentration in the Agar Substrate

The rate of change of the nutrient concentration in the agar substrate is

This reduces to

Equations (2), (5), (7), and (9) and initial conditions shown in Table 3 constitute our model. And this model can be simplified to

Table 3: Parameter values and initial values for the biofilm growth model simulations on agar substrate.

Term
Value
Units

0.0001
mg/cm³

12.96
1/day [27]

0.5
mg/mg

0.25
1/day

0.864
cm²/day [28]

1
cm²

0.8
None

300
mg/cm³

0.5
cm³

0.3
cm

15
mg/cm³ [14]

0.0001
cm [14]

0.00004
mg/cm³

0.15
None

4. Results

4.1. Experimental Results

We make time-lapse movies of the growing biofilm by recording the three fluorescent channels and the transmitted channel. After about 12 h the biofilm becomes visible to the naked eye and continues to grow in a circular fashion such that after 48 h its diameter is 16 mm as shown in Figure 2(a). The images show that three different phenotypes are spatially and temporally organized; details will be described below. To estimate the biofilm thickness from the transmission images, we develop a novel calibration procedure involving cross sections of biofilms. Similar to previous studies [14, 29, 30], we use the Beer-Lambert law to estimate the biofilm thickness from the optical density, OD, through the agar and bacterial colony:

where is the attenuation length. The optical density is defined as

where is intensity of the transmitted light through the substrate and colony and is that of the substrate alone, which is transparent agar. To determine , we compare the transmission (taken from above) with the height obtained from a side view. We reinforce the biofilm for subsequent manipulations by covering it with agar. We cut a thin slab of the biofilm on the agar substrate, whose top view is shown by the transmission image in the top of Figure 2(b). The transmission, , is determined from the ratio of the transmitted light relative to that of the agar substrate alone and averaged along the narrow, that is, transverse, direction. We next obtain the biofilm’s height by flipping the slab onto its side and image its cross section based on the constitutive fluorescent channel as shown in the bottom of Figure 2(b). Finally, we characterize the biofilm’s growth based upon the height profiles obtained from the optical density, as shown in Figures 2(c) and 2(b).

(a)

(b)

(c)

Figure 2: Characterization of growing Bacillus subtilis biofilm. (a) Time-lapse images of a growing Bacillus subtilis biofilm growing on top of nutritive agar medium at 12-hour intervals. Upon excitation, the motile, matrix-producing, and sporulating phenotypes produce different fluorescent colors and are observed using red, blue, and yellow filters, respectively. The transmission images are shown at the bottom and are used to determine the optical density. Green line serves as guide indicating at 48 hours that the extension of the matrix-producing region is comparable to the motile cells and exceeds that of the sporulating cells.

Scale bar is 10 mm. (b) Top panel: top view biofilm cross section; here yellow color region is biofilm and red color region is agar substrate. Bottom panel: side view of red constitutive fluorescent labeled biofilm; blue color is agar substrate. (c) Height profile derived from calibration of transmitted images.

4.2. Simulation Results

Given the parameters and initial values in Table 3, we numerically solve the main equations (10) by using Matlab package “ODE23s.”

4.2.1. Biofilm Thickness Evolution

Biofilm thickness increases from its initial value 0.0001 cm to its steady-state value 0.06 cm, which roughly agrees with our experiment, as shown in Figures 3(a) and 3(c). Biofilm thickness in our experiment firstly increases from its initial 0.025 cm to a maximum and then decreases to zero, as shown in Figure 3(b). By calculating biofilm approximate growth rate through the linear fitting, our results show a good coincidence with our experiment, as shown in inset of Figures 3(a) and 3(c).

(a)

(b)

(c)

Figure 3: The thickness profile of biofilm. (a) The biofilm thickness in 10 days for biofilm growth model on agar substrate. (b) The biofilm thickness change with its radius in our experiment. (c) The biofilm thickness in 2 days in our experiment. In addition, the rates of biofilm growth in our model and experiment are shown in inset of (a) and (c), which shows a good coincidence.

In addition, experiments showed that the biofilm resistance to antimicrobial agents was associated with the rate of cell growth in biofilm and biofilm age. The faster the rate of cell growth, the more rapid the rate of inactivation by ciprofloxacin [31]. Anwar et al. [32] found that 10-day-old chemostat-grown P. aeruginosa biofilms are significantly more resistant to tobramycin and piperacillin than 2-day-old biofilm. The biofilm stages with different growth rate in our model can be estimated, as shown in Figures 3(a) and 3(c); accordingly, we can test biocides more efficiently.

4.2.2. Volume Fraction of Active Biomass Evolution

The volume fraction of active biomass in biofilm increases from initial 0.15 rapidly to maximum 0.19 in about 1 day, which means that the yield rate of active biomass is larger than inactive rate of active biomass in this period. Then volume fraction of active biomass decreases until biofilms are occupied completely by inactive biomass, as shown in Figure 4.

Figure 4: The volume fraction of active biomass in biofilm in 30 days for biofilm growth model on agar substrate.

4.2.3. Nutrient Concentration in the Agar Substrate Evolution

Nutrient concentration in the agar substrate decreases all the time due to nutrient consumption by biofilm continuous growth and ultimately decreases to zero, as shown in Figure 5. At about 2.5 days, the nutrient concentration in the agar substrate decreases to zero and simultaneously the biofilm thickness increases to a steady-state value due to the nutrient depletion, as shown in Figure 3(a). The rate of nutrient consumption in agar substrate increases with the increase of nutrient diffusion coefficient, as shown in inset of Figure 5.

Figure 5: The nutrient concentration in the agar substrate in 10 days for the biofilm growth model. Inset: the nutrient concentration in agar substrate changes with diffusion coefficient   ().

4.2.4. Effect of Influenced Parameters on Biofilm Growth

We find that when yield coefficient and the maximum specific growth rate increase, the biofilm thickness increases accordingly, as shown in Figures 6(a) and 6(b), while yield coefficient has significant effect on biofilm final thickness, as shown in Figure 6(a). It was approved that a larger yield coefficient would enlarge cell growth rate, and increase of the maximum specific growth rate would increase time-rate-of-change of biomass [33]. The maximum specific growth rate can have positive effect on the maximum value of active biomass volume fraction during the first day; at the same time, active biofilm volume fraction decreases more rapidly with higher value of after the first day, as shown in inset of Figure 7(a). The reason is that the higher the is, the more the bacteria multiply, which causes more nutrient consumption and competition between bacteria cells. Active biomass volume fraction also decreases with increase of rate of active biomass inactivation , as shown in Figure 7(b). Similarly, Kluge et al. found that the inactivation rate of active biomass can cause a negative growth rate of active biomass [34].

(a)

(b)

(c)

(d)

Figure 6: The parameter effect on the biofilm thickness. (a) The yield coefficient    describes the ratio of the amount of biomass produced to the amount of substrate consumed [25]. (b) The maximum specific growth rate    describes the maximum proliferation ability of microorganisms. (c) The diffusion coefficient    reflects the nutrient diffusion ability into biofilm and (d) the biofilm biomass density    indicates the biomass in unit volume. Moreover, the biofilm thickness as a function of time for different (b) () and (c) () at a steady state (inset) is shown.

(a)

(b)

Figure 7: The active biomass volume fraction changes with (a) the maximum specific growth rate    in 25 days and (b) inactivation coefficient    in 30 days; the insets show the change of active biomass volume fraction with the two kinds of parameters in 2 days.

The biofilm thickness increases with increase of the diffusion coefficient at first 1.5 days, as shown in Figure 6(c). When the diffusion coefficient increases, the transportation rate of nutrient increases, and biofilm grows rapidly. At the same time, the time needed to nutrient complete depletion decreases, as shown in inset of Figure 5. Gonpot et al. [35] also found that an increase in diffusion coefficient favours growth of the biofilm. We find that when the biofilm reaches the steady state, the thickness is inversely proportional to diffusion coefficient . Because bacteria cannot fully absorb all nutrients from agar substrate when nutrient diffusion coefficient is high, the differences of final biofilm thickness with various diffusion coefficients are slight at 0.002 cm. Biofilm thickness is inversely proportional to biofilm density , as shown in Figure 6(d). Liu found that active biofilm thickness decreases significantly with the increased biofilm density [36].

5. Conclusion and Discussion

We obtain the biofilm thickness and the volume fraction of active biomass in biofilm evolution from our model. The biofilm thickness from our model is comparable with experimental result, which is shown in Figure 3. In fact, the biofilm thickness is inhomogeneous in both azimuthal and radial directions; it is much thicker between center and edge, as shown in Figures 2(b) and 3(b). Bacillus subtilis cells can differentiate into multiple phenotypes with different functions during biofilm formation, the main phenotypes including motile, matrix-producing, and sporulating types. We make time-lapse movies of the growing biofilm by recording the three fluorescent channels and the transmitted channel. We observe that the motile cells stay in a circular region about the inoculation spot, as shown in the second row of Figure 2(a). Matrix-producing cells grow radially outwards from the motile region as shown in the third row of Figure 2(a). Sporulating cells occur in a region inside the ring of the matrix-producing cells and outside the region of motile cells, as shown in the last row of Figure 2(a). There are complicated wrinkle patterns distributed in biofilm, which can transport water and nutrient; reasons for wrinkles structure formation are vertical force induced by cell death and lateral force induced from active cells. Wrinkle patterns also depend on mechanical properties of agar substrate and biofilm itself [37–39].

In our model, we obtain the quantitative relationship between some parameters and biofilm evolution; we can control biofilm growth by adjusting these parameters in future experiment. To get biofilm with different thickness and active biomass volume fraction, firstly we can adjust biomass yield coefficient by changing the amount of oxygen in biofilm [40, 41], biomass decay, and different nutrient concentration [2, 42]. Secondly, we can adjust diffusion coefficient and biofilm density [20, 43]. The biofilm density depends on many factors, such as different types of microorganisms [44], physical forces created by hydrodynamic conditions [45], increasing detachment forces related to particle-particle collisions [46, 47], and EPS [48]. Finally, we can adjust the maximum special growth rate of Monod equation by using different microbial species and initial nutrient concentration in agar substrate [49, 50].

In our model, we think nutrient acquisition of biofilm growth mainly depends on nutrient diffusion; the experiment found that the EPS triggering mechanism might bring in nutrient more efficiently than diffusion [29]. We need to consider this point for complete biofilm growth model.

Biofilm growth is complicated and is affected by many factors, including the specific bacteria strain [51, 52], material surface properties, and environmental parameters such as the pH and temperature [1]. To better understand its growth, we need to consider the influence of biofilm channels on biofilm growth, especially the effect of nutrient convection on biofilm thickness and radial expansion.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The work was supported by Harvard University, the National Natural Science Foundation of China (11272002), and Beijing Education Young Elite Teacher Project (YETP0363).

References

R. M. Donlan, “Biofilms: microbial life on surfaces,” Emerging Infectious Diseases, vol. 8, no. 9, pp. 881–890, 2002.

M. G. Trulear and W. G. Characklis, “Dynamics of biofilm processes,” Journal of the Water Pollution Control Federation, vol. 54, no. 9, pp. 1288–1301, 1982.

W. G. Characklis, M. J. Nimmons, and B. F. Picologlou, “Influence of fouling biofilms on heat transfer,” Heat Transfer Engineering, vol. 3, no. 1, pp. 23–37, 1981.

Q. Wang and T. Zhang, “Review of mathematical models for biofilms,” Solid State Communications, vol. 150, no. 21-22, pp. 1009–1022, 2010.

C. Picioreanu, M. C. M. Van Loosdrecht, and J. J. Heijnen, “Effect of diffusive and convective substrate transport on biofilm structure formation: a two-dimensional modeling study,” Biotechnology and Bioengineering, vol. 69, no. 5, pp. 504–515, 2000.

E. Alpkvist and I. Klapper, “A multidimensional multispecies continuum model for heterogeneous biofilm development,” Bulletin of Mathematical Biology, vol. 69, no. 2, pp. 765–789, 2007.

P. Cumsille, J. A. Asenjo, and C. Conca, “A novel model for biofilm growth and its resolution by using the hybrid immersed interface-level set method,” Computers and Mathematics with Applications, vol. 67, no. 1, pp. 34–51, 2014.

T. Zhang, N. Cogan, and Q. Wang, “Phase-field models for biofilms II. 2-D numerical simulations of biofilm-flow interaction,” Communications in Computational Physics, vol. 4, no. 1, pp. 72–101, 2008.

U. Ahmed and K. Vafai, “Analysis of biofilm growth in the presence of osmotic pressure and temperature effects,” International Journal of Heat and Mass Transfer, vol. 55, no. 21-22, pp. 5709–5721, 2012.

N. G. Cogan and J. P. Keener, “The role of the biofilm matrix in structural development,” Mathematical Medicine and Biology, vol. 21, no. 2, pp. 147–166, 2004.

E. Alpkvist, C. Picioreanu, M. C. M. van Loosdrecht, and A. Heyden, “Three-dimensional biofilm model with individual cells and continuum EPS matrix,” Biotechnology and Bioengineering, vol. 94, no. 5, pp. 961–979, 2006.

H. Fujikawa, “Diversity of the growth patterns of Bacillus subtilis colonies on agar plates,” FEMS Microbiology Ecology, vol. 13, no. 3, pp. 159–167, 1994.

H. Fujikawa and M. Matsushita, “Fractal growth of Bacillus subtilis on agar plates,” Journal of the Physical Society of Japan, vol. 58, no. 11, pp. 3875–3878, 1989.

A. Seminara, T. E. Angelini, J. N. Wilking et al., “Osmotic spreading of Bacillus subtilis biofilms driven by an extracellular matrix,” Proceedings of the National Academy of Sciences of the United States of America, vol. 109, no. 4, pp. 1116–1121, 2012.

M. Henze, “Characterization of wastewater for modelling of activated sludge processes,” Water Science and Technology, vol. 25, no. 6, pp. 1–15, 1992.

G. A. Ekama, P. L. Dold, and G. V. R. Marais, “Procedures for determining influent COD fractions and the maximum specific growth rate of heterotrophs in activated sludge systems,” Water Science & Technology, vol. 18, no. 6, pp. 91–114, 1986.

W. Gujer and M. Henze, “Activated sludge modelling and simulation,” Water Science & Technology, vol. 23, no. 4–6, pp. 1011–1023, 1991.

P. Dalgaard and K. Koutsoumanis, “Comparison of maximum specific growth rates and lag times estimated from absorbance and viable count data by different mathematical models,” Journal of Microbiological Methods, vol. 43, no. 3, pp. 183–196, 2001.

R. S. Renslow, P. D. Majors, J. S. McLean, J. K. Fredrickson, B. Ahmed, and H. Beyenal, “In situ effective diffusion coefficient profiles in live biofilms using pulsed-field gradient nuclear magnetic resonance,” Biotechnology and Bioengineering, vol. 106, no. 6, pp. 928–937, 2010.

L.-S. Fan, R. Leyva-Ramos, K. D. Wisecarver, and B. J. Zehner, “Diffusion of phenol through a biofilm grown on activated carbon particles in draft-tube three-phase fluidized-bed bioreactor,” Biotechnology and Bioengineering, vol. 35, no. 3, pp. 279–286, 1990.

S. Chang and S. N. Cohen, “High frequency transformation of Bacillus subtilis protoplasts by plasmid DNA,” Molecular and General Genetics, vol. 168, no. 1, pp. 111–115, 1979.

V. Gadani, P. Villon, J. Manem, and B. Rittmann, “A new method to solve a non-steady-state multispecies biofilm model,” Bulletin of Mathematical Biology, vol. 55, no. 6, pp. 1039–1061, 1993.

C. S. Laspidou and B. E. Rittmann, “Non-steady state modeling of extracellular polymeric substances, soluble microbial products, and active and inert biomass,” Water Research, vol. 36, no. 8, pp. 1983–1992, 2002.

C. S. Laspidou and B. E. Rittmann, “A unified theory for extracellular polymeric substances, soluble microbial products, and active and inert biomass,” Water Research, vol. 36, no. 11, pp. 2711–2720, 2002.

http://www.lenntech.com/wwtp/coefficienty.htm#ixzz3XTbbgDwu

S. K. Tiwari and K. L. Bowers, “Modeling biofilm growth for porous media applications,” Mathematical and Computer Modelling, vol. 33, no. 1–3, pp. 299–319, 2001.

K. Kovarova, A. J. Zehnder, and T. Egli, “Temperature-dependent growth kinetics of Escherichia coli ML 30 in glucose-limited continuous culture,” Journal of Bacteriology, vol. 178, pp. 4530–4539, 1996.

N. N. Sinha, Visualizing patterns of gene expression in growing Bacillus subtilis biofilms [Ph.D. thesis], Harvard University, Cambridge, Mass, USA, 2013.

W. Zhang, A. Seminara, M. Suaris, M. P. Brenner, D. A. Weitz, and T. E. Angelini, “Nutrient depletion in Bacillus subtilis biofilms triggers matrix production,” New Journal of Physics, vol. 16, Article ID 015028, 2014.

J. Dervaux, J. C. Magniez, and A. Libchaber, “On growth and form of Bacillus subtilis biofilms,” Interface Focus, vol. 4, no. 6, Article ID 20130051, 2014.

I. G. Duguid, E. Evans, M. R. W. Brown, and P. Gilbert, “Growth-rate-independent killing by ciprofloxacin of biofilm-derived Staphylococcus epidermidis; evidence for cell-cycle dependency,” Journal of Antimicrobial Chemotherapy, vol. 30, no. 6, pp. 791–802, 1992.

H. Anwar, J. L. Strap, K. Chen, and J. W. Costerton, “Dynamic interactions of biofilms of mucoid Pseudomonas aeruginosa with tobramycin and piperacillin,” Antimicrobial Agents and Chemotherapy, vol. 36, no. 6, pp. 1208–1214, 1992.

http://www.cs.montana.edu/webworks/projects/stevesbook/contents/chapters/chapter002/section002/black/page001.html

M. Kluge, D. Siegmund, H. Diekmann, and M. Thoma, “A model for penicillin production with and without temperature shift after the growth phase,” Applied Microbiology and Biotechnology, vol. 36, no. 4, pp. 446–451, 1992.

P. Gonpot, R. Smith, and A. Richter, “Diffusion limited biofilm growth,” Modelling and Simulation in Materials Science and Engineering, vol. 8, no. 5, pp. 707–726, 2000.

Y. Liu, “Estimating minimum fixed biomass concentration and active thickness of nitrifying biofilm,” Journal of Environmental Engineering, vol. 123, no. 2, pp. 198–202, 1997.

J. N. Wilking, T. E. Angelini, A. Seminara, M. P. Brenner, and D. A. Weitz, “Biofilms as complex fluids,” MRS Bulletin, vol. 36, no. 5, pp. 385–391, 2011.

M. Asally, M. Kittisopikul, P. Rué et al., “Localized cell death focuses mechanical forces during 3D patterning in a biofilm,” Proceedings of the National Academy of Sciences of the United States of America, vol. 109, no. 46, pp. 18891–18896, 2012.

C. Picioreanu, M. van Loosdrecht, and J. Heijnen, Multidimensional Modeling of Biofilm Structure, Faculty of Applied Sciences, Delft University of Technology, 1999.

G. J. Kehrberger, Effect of recirculation on the performance of a trickling filter [Ph.D. thesis], Rice University, Houston, Tex, USA, 1969, http://hdl.handle.net/1911/14532.

R. C. Hoehn and A. D. Ray, “Effects of thickness on bacterial film,” Journal of the Water Pollution Control Federation, vol. 45, no. 11, pp. 2302–2320, 1973.

J. C. Merchuk and J. A. Asenjo, “The Monod equation and mass transfer,” Biotechnology and Bioengineering, vol. 45, no. 1, pp. 91–94, 1995.

O. Wanner and W. Gujer, “Competition in biofilms,” Water Science and Technology, vol. 17, no. 2-3, pp. 27–44, 1985.

B. E. Christensen and W. G. Characklis, “Physical and chemical properties of biofilms,” Biofilms, vol. 93, article 130, 1990.

C. S. Laspidou and B. E. Rittmann, “Modeling the development of biofilm density including active bacteria, inert biomass, and extracellular polymeric substances,” Water Research, vol. 38, no. 14-15, pp. 3349–3361, 2004.

H. T. Chang, B. E. Rittmann, D. Amar, R. Heim, O. Ehlinger, and Y. Lesty, “Biofilm detachment mechanisms in a liquid-fluidized bed,” Biotechnology and Bioengineering, vol. 38, no. 5, pp. 499–506, 1991.

W. K. Kwok, C. Picioreanu, S. L. Ong, M. C. M. Van Loosdrecht, W. J. Ng, and J. J. Heijnen, “Influence of biomass production and detachment forces on biofilm structures in a biofilm airlift suspension reactor,” Biotechnology and Bioengineering, vol. 58, no. 4, pp. 400–407, 1998.

H.-C. Flemming, T. R. Neu, and D. J. Wozniak, “The EPS matrix: the ‘house of biofilm cells’,” Journal of Bacteriology, vol. 189, no. 22, pp. 7945–7947, 2007.

H. W. Jannasch, “Growth of marine bacteria at limiting concentrations of organic carbon in seawater,” Limnology and Oceanography, vol. 12, no. 2, pp. 264–271, 1967.

B. Chen-Charpentier, “Numerical simulation of biofilm growth in porous media,” Journal of Computational and Applied Mathematics, vol. 103, no. 1, pp. 55–66, 1999.

M. S. Chae and H. Schraft, “Comparative evaluation of adhesion and biofilm formation of different Listeria monocytogenes strains,” International Journal of Food Microbiology, vol. 62, no. 1-2, pp. 103–111, 2000.

M. K. Borucki, J. D. Peppin, D. White, F. Loge, and D. R. Call, “Variation in biofilm formation among strains of Listeria monocytogenes,” Applied and Environmental Microbiology, vol. 69, no. 12, pp. 7336–7342, 2003.