An individual-based model for biofilm formation at liquid surfaces Maxime Ardré+,∗,o, Hervé Henry+, Carine Douarche∗, Mathis Plapp+ + Laboratoire PMC, Ecole Polytechnique CNRS, Palaiseau, France ∗ Laboratoire de Physique des Solides, Université Paris-Sud CNRS, Orsay, France o Laboratoire de Physique Statistique, ENS CNRS, Paris, France May 7, 2018 Abstract The bacterium Bacilus subtilis frequently forms biofilms at the interface between the culture medium and the air. We develop a mathematical model that couples a description of bacteria as individ- ual discrete objects to the standard advection-diffusion equations for the environment. The model takes into account two different bacte- rial phenotypes. In the motile state, bacteria swim and perform a run-and-tumble motion that is biased toward regions of high oxygen concentration (aerotaxis). In the matrix-producer state they excrete extracellular polymers, which allows them to connect to other bacteria and to form a biofilm. Bacteria are also advected by the fluid, and can trigger bioconvection. Numerical simulations of the model reproduce all the stages of biofilm formation observed in laboratory experiments. Finally, we study the influence of various model parameters on the dy- namics and morphology of biofilms. 1 Introduction arXiv:1512.04697v1 [q-bio.CB] 15 Dec 2015 Bacteria are unicellular prokaryotic microorganisms. They are ubiquitous and constitute a large part of the terrestrial biomass. Bacteria can live as individual cells during the planktonic phase. However, most of the time, they are part of self-organized communities of complex architecture adsorbed on interfaces: the biofilms. Besides the bacteria themselves, biofilms are mostly made of an extracellular matrix composed of macromolecules [1, 2, 3] that are produced by the bacteria and lead to cohesive interactions between them 1 [4, 5]. Most of the time biofilms appear in aqueous environments either on solid surfaces or at the water-air interfaces. Due to this multicellular organization, the bacteria have different properties than in the motile state. For example, bacteria trapped in biofilms can exhibit an increased resistance to antibiotics as well as to environmental stresses (desiccation, UV radiation, disinfecting agents, shear flow...) [6, 7]. Therefore, the association in biofilms is a crucial step both for survival and spreading of bacterial colonies [8]. Many environmental and genetic factors influence the development of biofilms [8]. Although biofilm formation is not understood in all details, a consistent picture of biofilm growth on solid substrates has been proposed: bacteria in the planktonic phase anchor preferentially on a stable surface, which initiates the nucleation of bacterial microcolonies [5, 8]. Bacteria con- stituting the microcolonies secrete an extracellular matrix in which they em- bed, and form a mature biofilm with a complex multiscale architecture [9, 10]. Later on, bacteria can also detach from the superficial layers, return to the planktonic state and spread to new parts of the surface[11]. Much less is known on biofilm formation at water-air interfaces. We have performed laboratory experiments on Bacillus subtilis (B. subtilis) and built a mathematical model that can reproduce the main experimental observations. The details of the experimental results will be published elsewhere; here, we briefly present a summary to motivate the construction of the mathematical model. Bacillus subtilis is a stricly aerobic bacterium which is able to form float- ing pellicles that can reach a thickness of several hundred microns [12] on the top of nutritive media. During a typical experiment, the nutritive medium (Luria broth supplemented with glycerol and MnO2) is initially inoculated with a small concentration of bacteria. Then, the cells divide and grow over several hours to yield a homogeneous suspension (see Fig. 1 (a)). Later (Fig. 1 (b)), bacteria start to accumulate close to the interface, and a rapid transition occurs (between (b) and (c)), which leads to the nucleation of bac- terial clusters on the interface (Fig. 1 (c)). At this stage, some macroscopic filaments, which can be seen inside the liquid, are advected by fluid motion. With time, the biofilm further develops into a mature floating pellicle that exhibits a typical wrinkled morphology (Fig. 1(d)). These experiments indicate several important phenomena that need be included in the model. First, B. subtilis is known to exhibit aerotaxis [13], that is, to migrate to areas in which the concentration of oxygen is high. In the experiments, the water-air interface acts as an oxygen source; therefore, aerotaxis is a possible explanation for the accumulation of bacteria close to the surface. Second, planktonic bacteria are slightly more dense than water. As a result, an accumulation of bacteria at the surface is unstable 2 Figure 1: Side view of B. subtilis pellicle developing on top of a liquid culture media at T=31◦C at : (a) 0 hours, (b) 15 hours (c) 19 hours 30 minutes, and (d) 67 hours after incubation. Images are acquired with a black and white CCD camera. The scale bar is 1 cm long. and should drive the development of bioconvection [14], which may hinder or help the biofilm formation. Finally, there is a well-defined change in bacterium phenotype from motile (planktonic phase) to matrix-producing (biofilm phase). This change is controlled by a genetic switch that is believed to be triggered by a quorum sensing mechanism [15, 10] (i.e. when the local density of bacteria exceeds a threshold value, the phenotype changes). Many different modelling strategies have already been proposed and used for biofilm formation on solid substrates. It is possible to average the con- tribution of individual bacteria and to write a full continuum model [16, 17]. It is also possible to consider the bacteria as individual objects in cellular automata models (CA) [18] or individual-based models (IdbM) [19]. The latter model can be hybridized with a continuum model to describe the con- tribution of the environment. Modeling is often accompanied by numerical simulations to study the processes which structure the biofilm [20, 21, 22]. Here, we formulate a detailed model for biofilm formation at liquid-air interfaces. The main point is that the model captures the transition between a “gas” of individual swimming bacteria and the biofilm pellicle (a soft solid). Since this transition involves a change in the connectivity between bacteria, a description using a full continuum model is difficult. We have chosen a hybrid approach, in which the bacteria are described as individual particles, whereas the local environment (oxygen concentration, fluid velocity) is described by continuous fields. In order to keep the model minimal, we only consider two bacterial phenotypes: in the motile state, bacteria are self-propelled swimmers that perform a standard “run-and-tumble” motion. The bacteria interact through a local repulsive potential to describe collisions between bacteria (hydrodynamic interactions are neglected). In the matrix-producer state, the bacteria stop moving actively, and produce macro-molecules that constitute the extracellular matrix. Due to the presence of this matrix, they are able to “connect” to other bacteria. In this case, the interaction between 3 bacteria is described by an attractive potential. The transition between the two phenotypes is triggered by a quorum-sensing mechanism. This model has been implemented in two dimensions, using the discrete- element method [23] to calculate the evolution of the bacteria, and the finite- difference method [24] for the continuous fields, and simulations involving up to 10000 bacteria were performed. This setting is sufficient to demonstrate that all the steps of biofilm formation that are observed in the experiments can be reproduced. The influence of the model parameters on these different stages was also studied, yielding suggestions for further model improvements and experiments. In the following, an overview of the model architecture is first given, followed by a detailed description of each ingredient (section 2). In section 3, the choice of the model parameters is discussed in detail. In Section 4, simulation results are presented, which demonstrate the ability of the model to properly describe biofilm formation. Finally, section 5 discusses the main conclusions and perspectives of this work. 2 Model description 2.1 Overview Our goal is to construct a minimal mathematical model that reproduces the different steps of biofilm formation, from motile bacteria (individuals swimming in the liquid) all the way to the mature biofilm (bacteria linked by extracellular matrix). In order to include the phenomena of aerotaxis and bioconvection, the local oxygen concentration c and the fluid velocity −→u are described as continuous fields which obey partial differential equations (PDE). Bacteria are represented as discrete objects, with interactions that de- pend on their internal state. Each bacterium is characterized by its position, its velocity, and internal variables that reproduce its behavior (aerotaxis, phenotype, cell cycle). One of these internal variables is the bacterial pheno- type. We take into account only two of them : motile and matrix producer. In the motile phenotype, the bacteria propel themselves with a constant ve- locity and change their direction with a frequency that is determined by the local concentration in oxygen (run-and-tumble motion). Each motile bac- terium increases its body size with time and divides into two cells with