Research Article  Open Access
Modeling of the Bacillus subtilis Bacterial Biofilm Growing on an Agar Substrate
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 diffusionlimited 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_{2}; 700 μM CaCl_{2}; 50 μM MnCl_{2}; 50 μM FeCl_{2}; 1 μM ZnCl_{2}; 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 matrixproducing 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_{600} = 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’ timelapse 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^{2}. The volume of the agar substrate is 0.5 cm^{3}.
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.


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 toHere, , 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
4. Results
4.1. Experimental Results
We make timelapse 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 BeerLambert 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 aswhere 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)
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 steadystate 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)
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 10dayold chemostatgrown P. aeruginosa biofilms are significantly more resistant to tobramycin and piperacillin than 2dayold 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.
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 steadystate 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.
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 timerateofchange 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)
(a)
(b)
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, matrixproducing, and sporulating types. We make timelapse 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). Matrixproducing 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 matrixproducing 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 particleparticle 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. View at: Publisher Site  Google Scholar
 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. View at: Google Scholar
 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. View at: Publisher Site  Google Scholar
 Q. Wang and T. Zhang, “Review of mathematical models for biofilms,” Solid State Communications, vol. 150, no. 2122, pp. 1009–1022, 2010. View at: Publisher Site  Google Scholar
 C. Picioreanu, M. C. M. Van Loosdrecht, and J. J. Heijnen, “Effect of diffusive and convective substrate transport on biofilm structure formation: a twodimensional modeling study,” Biotechnology and Bioengineering, vol. 69, no. 5, pp. 504–515, 2000. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 P. Cumsille, J. A. Asenjo, and C. Conca, “A novel model for biofilm growth and its resolution by using the hybrid immersed interfacelevel set method,” Computers and Mathematics with Applications, vol. 67, no. 1, pp. 34–51, 2014. View at: Publisher Site  Google Scholar
 T. Zhang, N. Cogan, and Q. Wang, “Phasefield models for biofilms II. 2D numerical simulations of biofilmflow interaction,” Communications in Computational Physics, vol. 4, no. 1, pp. 72–101, 2008. View at: Google Scholar
 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. 2122, pp. 5709–5721, 2012. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 E. Alpkvist, C. Picioreanu, M. C. M. van Loosdrecht, and A. Heyden, “Threedimensional biofilm model with individual cells and continuum EPS matrix,” Biotechnology and Bioengineering, vol. 94, no. 5, pp. 961–979, 2006. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 M. Henze, “Characterization of wastewater for modelling of activated sludge processes,” Water Science and Technology, vol. 25, no. 6, pp. 1–15, 1992. View at: Google Scholar
 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. View at: Google Scholar
 W. Gujer and M. Henze, “Activated sludge modelling and simulation,” Water Science & Technology, vol. 23, no. 4–6, pp. 1011–1023, 1991. View at: Google Scholar
 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. View at: Publisher Site  Google Scholar
 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 pulsedfield gradient nuclear magnetic resonance,” Biotechnology and Bioengineering, vol. 106, no. 6, pp. 928–937, 2010. View at: Publisher Site  Google Scholar
 L.S. Fan, R. LeyvaRamos, K. D. Wisecarver, and B. J. Zehner, “Diffusion of phenol through a biofilm grown on activated carbon particles in drafttube threephase fluidizedbed bioreactor,” Biotechnology and Bioengineering, vol. 35, no. 3, pp. 279–286, 1990. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 V. Gadani, P. Villon, J. Manem, and B. Rittmann, “A new method to solve a nonsteadystate multispecies biofilm model,” Bulletin of Mathematical Biology, vol. 55, no. 6, pp. 1039–1061, 1993. View at: Publisher Site  Google Scholar
 C. S. Laspidou and B. E. Rittmann, “Nonsteady state modeling of extracellular polymeric substances, soluble microbial products, and active and inert biomass,” Water Research, vol. 36, no. 8, pp. 1983–1992, 2002. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 K. Kovarova, A. J. Zehnder, and T. Egli, “Temperaturedependent growth kinetics of Escherichia coli ML 30 in glucoselimited continuous culture,” Journal of Bacteriology, vol. 178, pp. 4530–4539, 1996. View at: Google Scholar
 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. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 I. G. Duguid, E. Evans, M. R. W. Brown, and P. Gilbert, “Growthrateindependent killing by ciprofloxacin of biofilmderived Staphylococcus epidermidis; evidence for cellcycle dependency,” Journal of Antimicrobial Chemotherapy, vol. 30, no. 6, pp. 791–802, 1992. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 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. View at: Google Scholar
 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. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 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. View at: Google Scholar
 J. C. Merchuk and J. A. Asenjo, “The Monod equation and mass transfer,” Biotechnology and Bioengineering, vol. 45, no. 1, pp. 91–94, 1995. View at: Publisher Site  Google Scholar
 O. Wanner and W. Gujer, “Competition in biofilms,” Water Science and Technology, vol. 17, no. 23, pp. 27–44, 1985. View at: Google Scholar
 B. E. Christensen and W. G. Characklis, “Physical and chemical properties of biofilms,” Biofilms, vol. 93, article 130, 1990. View at: Google Scholar
 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. 1415, pp. 3349–3361, 2004. View at: Publisher Site  Google Scholar
 H. T. Chang, B. E. Rittmann, D. Amar, R. Heim, O. Ehlinger, and Y. Lesty, “Biofilm detachment mechanisms in a liquidfluidized bed,” Biotechnology and Bioengineering, vol. 38, no. 5, pp. 499–506, 1991. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 B. ChenCharpentier, “Numerical simulation of biofilm growth in porous media,” Journal of Computational and Applied Mathematics, vol. 103, no. 1, pp. 55–66, 1999. View at: Publisher Site  Google Scholar  MathSciNet
 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. 12, pp. 103–111, 2000. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
Copyright
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.